个人用dmap写了个小程序 生成给定几何模型的单元刚度矩阵和整体刚度矩阵
程序虽小,但是dmap编程从这可以进入 较有用处
SOL getmat
$
COMPILE getmat, LIST
$
SUBDMAP getmat
$ DECLARE PARAMETERS
TYPE PARM,,I,Y,GENSHU,GENSHUM,JDS,JDS1$
$
IFP1 /CASECC,PCDB,XYCDB,POSTCDB,FORCE/S,N,NOGOIFP1 $
XSORT FORCE,/BULK/S,N,NOGOXSRT $
IFP BULK/
GEOM1,EPT,MPT,EDT,DIT,DYNAMIC,
GEOM2,GEOM3,GEOM4,UNUSED1,UNUSED2,
MATPPOL,AXIC,PVT,DMI,DMINDX,DTI,DTINDX,
DEFUSET,EDOM,DEQATN,DEQIND,CONTACT,OINT,
UNUSED3/
S,N,NOGOIFP/S,N,RUNIFP3/S,N,RUNIFP4/
S,N,RUNIFP5/S,N,RUNIFP6/
S,N,RUNIFP7/S,N,RUNIFP8/S,N,RUNIFP9//
S,N,RUNMEPT $
PVT PVT,CASECC/ $ $$加入后可识别用户自定义参数param,此外需要type中设置y
GP1 GEOM1,GEOM2,,,,,,,/
GPL,EQEXIN,GPDT,CSTM,BGPDT,SIL,,,,/
S,N,LUSET/0/S,N,NOGPDT $
GP2 GEOM2,EQEXIN,EPT,,,BGPDT,CSTM,/ECT,/S,N,ACOUSTIC $
GP4 CASECC,GEOM4,EQEXIN,SIL,GPDT,BGPDT,CSTM,
,,,GEOM2,,ECT,,GEOM1/
RMG,YG1,USET0,GPSNTN/
LUSET/S,N,NOMSET/0/S,N,NOSSET/S,N,NOOSET/
S,N,NORSET/0/0/0/0/0////
//S,N,NEWGPSN// $
TA1 MPT,ECT,EPT,BGPDT,SIL,,CSTM,DIT,,,,/
EST,,GEI,GPECT,,VGFD,DITID,,NSMEST/
LUSET/-1/0/1/0/0 $
EMG EST,CSTM,MPT,DIT,,,,,DEQATN,DEQIND,BGPDT,,
ECT,EPT,,,,,,,,,,,,,,,,,,/
KELM,KDICT,MELM,MDICT,BELM,BDICT,,,,,/
S,N,NOKGG/S,N,NOMGG/S,N,NOBGG/S,N,NOK4GG//
//////////////////////S,N,MGE $
EMA GPECT,KDICT,KELM,BGPDT,SIL,CSTM,,/KJJLF, $
EMA GPECT,MDICT,MELM,BGPDT,SIL,CSTM,,/MJJLF, $
matprn KJJLF//
matprn KELM//
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$输出单元刚度矩阵和整体刚度矩阵
end
COMPILE SUBDMAP= PRTNSR, LIST, NOREF $
SUBDMAP PRTNSR ANYJJ,PARTU1,PARTU2/ANYJJ00,ANYJJ01,ANYJJ10,ANYJJ11/ $
$
PARTN ANYJJ,PARTU1,PARTU2/ANY00,ANY10,ANY01,ANY11/-1/ $
MERGE ,ANY00,,,,PARTU1,/ANYJJ00/-1/ $
MERGE ,,,,ANY11,PARTU1,/ANYJJ11/-1/ $
MERGE ,,ANY10,,,PARTU1,/ANYJJ10/-1/ $
MERGE ,,,ANY01,,PARTU1,/ANYJJ01/-1/ $
$
RETURN $
END $
cend
SUBTITLE = sample use
BEGIN BULK
PARAM AUTOSPC YES
PARAM PRTMAXIM YES
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
PARAM,GENSHU,4$
PARAM,JDS,134$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$
include'5.bdf' 我自己的模型
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
ENDDATA |