/MAT/LAW19 で DETERMINANT OF MATERIAL MATRIX IS LESS THAN 0. が出る訳
altair_fukuoka
Altair Employee
Altair Employeeさっそくですが、この材料モデルはスターターで
/MAT/LAW19/1 law19 # Rho 1e-9 # E11 E22 NU12 10000 150000 .3 # G12 G23 G31 6000 6000 3000 # RE ZeroStress Fsclae_Por sens_ID
このエラーを出します。
ERROR ID : 307 ** INPUT ERROR IN MATERIAL LAW DESCRIPTION : -- MATERIAL ID: 1 -- MATERIAL TITLE: law19 DETERMINANT OF MATERIAL MATRIX IS LESS THAN 0. CHECK CONSTITUTIVE LAW
なぜでしょうか？
実はこちらの材料マトリクス C-1 の行列式 (Determination) が負になるからです。
では実際に計算してみます。
とんでもないポアソン比とともに、行列式が負になっています。行列式が負というのは、大雑把に言えば、膨らませる応力を掛けたら縮んでしまった、ということになります。計算できなくもないのかもしれませんが、RADIOSS はこれをエラーとしています。
実はこのデータ NU21=0.3, NU12=0.02 に対して、 NU12=0.3 と言うように入力を間違えてしまったのです。では間違いに気が付いたので、NU12=0.02 に直して計算し直してみます。
これで RADIOSS のエラーも出なくなります。
Compose スクリプトも張り付けておきます。
clc; clear; % 入力値 E11=10000.; E22=150000.; NU12=0.02; G12=6000.; G23=6000.; G31=3000.; % 自動計算 NU21=NU12*E22/E11 CREV=[ 1./E11, -NU21/E22, 0., 0., 0.; -NU12/E11, 1./E22, 0., 0., 0.; 0., 0., 1./G12, 0., 0.; 0., 0., 0., 1./G23, 0.; 0., 0., 0., 0., 1./G31]; DET=det(CREV)
おそらく一番多い間違いが NU12 と NU21 の取り間違いだと思います。
アルテアジャパン公式製品ページ
https://www.altairjp.co.jp/radioss/
https://www.altairjp.co.jp/compose/
2