さっそくですが、この材料モデルはスターターで
/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/