OptiStruct は予荷重または予張力を考慮した固有値の解析ができます。弦楽器の弦を張れば張るほど高い音が出るあれです。
次のモデルを使って説明します。
やり方は簡単で、固有値解析サブケースにて STATSUB(PRELOAD)=予荷重として使いたい線形静解析サブケース ID を打ち込むだけです。
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ 線形制解析サブケース
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
SUBCASE 1
LABEL static
ANALYSIS STATICS
SPC = 1
LOAD = 2
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ 固有値解析サブケース
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
SUBCASE 2
LABEL eigrl
ANALYSIS MODES
SPC = 1
METHOD(STRUCTURE) = 1
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
STATSUB(PRELOAD) = 1
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
実際に計算すると固有値は 2 となります。Eigenvalue を見てください。
Generalized Generalized
Subcase Mode Frequency Eigenvalue Stiffness Mass
2 1 2.250791E-01 2.000000E+00 2.000000E+00 1.000000E+00
こちらが正しいか確認してみましょう。こちらの資料がよくまとまっていてありがたいです。P1 に K_L (もともとの剛性) と K_G (予荷重による剛性) の行列式があります。
http://civilyarou.web.fc2.com/WANtaroHP_html5_win/f90_FEM/f90_FDA/TeX_FDA_frame.pdf
今回のモデルは X 自由度以外を拘束しているので、行列式の一番左上だけ見れば OK です。すべての変数を 1 にしてあるので K_L=1, K_G=1 です。全体剛性はK= K_L + K_G = 2 です。したがって、固有値は K/M=2 となります。
ということで正しく計算できたことがわかりました。
今回は固有値解析としましたが、周波数応答解析と過渡応答解析も同様に利用可能です。お試しください。
この入力モデルのダウンロード: model01.fem