个人用dmap写了个小程序 生成给定几何模型的单元刚度矩阵和整体刚度矩阵
程序虽小,但是dmap编程从这可以进入 较有用处
SOL getmat
COMPILEgetmat,LIST
SUBDMAP getmat
DECLAREPARAMETERSTYPEPARM,,I,Y,GENSHU,GENSHUM,JDS,JDS1
IFP1/CASECC,PCDB,XYCDB,POSTCDB,FORCE/S,N,NOGOIFP1
XSORT FORCE,/BULK/S,N,NOGOXSRT IFPBULK/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中设置yGP1GEOM1,GEOM2,,,,,,,/GPL,EQEXIN,GPDT,CSTM,BGPDT,SIL,,,,/S,N,LUSET/0/S,N,NOGPDT
GP2 GEOM2,EQEXIN,EPT,,,BGPDT,CSTM,/ECT,/S,N,ACOUSTIC GP4CASECC,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 EMGEST,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, EMAGPECT,MDICT,MELM,BGPDT,SIL,CSTM,,/MJJLF,
matprn KJJLF//
matprn KELM//
输出单元刚度矩阵和整体刚度矩阵endCOMPILESUBDMAP=PRTNSR,LIST,NOREF$SUBDMAPPRTNSRANYJJ,PARTU1,PARTU2/ANYJJ00,ANYJJ01,ANYJJ10,ANYJJ11/$$PARTNANYJJ,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$cendSUBTITLE=sampleuseBEGINBULKPARAMAUTOSPCYESPARAMPRTMAXIMYESPARAM,GENSHU,4
PARAM,JDS,134$include′5.bdf′我自己的模型$
ENDDATA |