MSC.patran,nastran中输出总体刚度阵和质量阵的方法

1、从 NASTRAN 中导出质量阵和刚度阵1.将总体刚度阵和总体质量阵输出到.f06 文件中COMPILE SEMG $ALTER 'KJJZ.*STIFFNESS' $DIAGON (30) $MATPRN KJJZ / $DIAGOFF (30) $ALTER 'MJJX,.*MASS' $ DIAGON (30) $ MATPRN MJJX / $DIAGOFF (30) $把以上的代码插入到 bdf 文件中 3 个$的位置，如下图 SOL 101TIME 5diag 8$ (INSERT DMAP HERE)$CEND然后就可以在f06文件中找到刚度阵kj

2、jz和质量阵mjjx2.将刚度、质量阵输出到外部文件中output41)用0UTPUT2、0UTPUT4将刚度阵和质量阵分别输出到 outout2 (二进制码) (数据)文件中：assign output2=output2.dat',unit=12,formatted,deleteassign output4='output4.dat',unit=13,formatted,deleteS0L 103$ Direct Text Input for Executive ControlDIAG 8C0MPILE SEMG $ALTER 'KJJZ.*STIFFNESS

3、' $ALTER 'MJJX,.*MASS' $0UTPUT2 MJJX,KJJZ, , /12 $0UTPUT4 MJJX,KJJZ, , /13/0 $CEND在 output4 文件中，数据基本格式为：列 行一个节点具有6个自由度.定义22个节点，故去醉总大小为132叫321321326 2KJJZ lPEW.g1 17 第一列中，笫1行后的7牛数据15.3666666t57E+13 O.ODOODOOOOE+OO O.OOODOOOOOE+OO D OODOOODOOE+OO C.OOOODOOODE+OOD 0 0 DOO D DO 0E+0. 8 6(5(5

4、66(5 S7E + 132 211第二列中，笫2行后的11个数据2.479324333E+D6 O.DOOOOOODOE+M O.OODOOOOOOE+CO D.ODOOODOOOE+OD 1.85943325OE+D7D Q00000000E+ 00-2.475)324333E+O6 0 OODOOOOOOE+OO O.OOOOOOOODE+QOi O.DOOOOOODOE+OO1.859493250E+072) MATPCH KJJZ,MJJX, , / $将矩阵中的非零数据输入到* pch文件中，具体操作如下所示：SOL 103$ Direct Text In put for Exec

5、utive Con trolDIAG 8$COMPILE SEMG $ALTER 'KJJZ.*STIFFNESS' $ALTER 'MJJX,.*MASS' $MATPCH KJJZ, MJJX, , / $CEND其格式如下DM1KJJZ0S10矩阵大小132132笫一列的第1行+第7行数据DMI*KIIZ11 6.86666(5 (53E+13*7-(5r86(56fi663E+13第一列的第2行+第6行+第呂行+第辽行数据DMPKJJZ222.47932425E+0C*6 1.8594P320E+078 2479卫425E十D6*12 L85P43320

6、E+07格式相关参考资料有矩阵MAT，利用DAMP语句：MATPCH MAT/ $得到如下格式1.00.06.00.00.00.00.00.07.00.00.00.0MAT-204.00.00.00.00.0005.00,00.000903.00,08.00.00,00.0DMIMAT0Z1256DMIMATi11000000E 003z.oiknoeii w5扌 ooouftjt onDMI*MAT231 CMJOWlK Od5.000000E00DMI'316.00000JL 007.0000001( 005BOOOUOOlt DOMATg49 DoomnE Od0UTPUT2、0

7、UTPUT4 和 MATPCH 可任选3.用NASTRAN自带语句直接得到*pch文件在NASTRAN中加入下列代码:PARAM,EXTOUT, DMIGPCH 得到.PCH 文件格式如下:DMIG缩减矩阵KAAXDMIG*+KAAXT L-DMKTKAAXDMIG”«KAAX1DMIG*KAAX12对鮒辭0 6对应矩阵排列_ _ . ?_ .EC666efil57D+13L3 2.479324JJ 30+015L3-1.8594332500+075 1.8594932 5 0LH0 321-6 866666667D+131 1,刃刃站333D+W相应列位置1个节点的第个自由度相应行

8、位晝-笫1亍节点的拿1个自由度值为6.3666 6156(5704-13相应列应置T1个节点的第3个自由度相应行位置第I个节点的第了个自由由值为1.A79324333D+G6相应列位養Tl个节点的第于,自由度相应行位置一笫1个节点的第3个自由度值为-1朗姒血0D+O7 相应行位置韩个节点的第5个自由度值为1.8594?3250>08 相应列位宣2个节点的衛个自由度相应行位置-笫1个节点的第1个自由度值为-6 866666667D+13 相应行位置第2个节点的第1个自由度值为1.373333333D+14第1个节点的石个自由度212.47P3JJ4333D+0d第2伞节点的疔个自由度123

9、_ _A _L 1-I 科P493彷UD+CfT1.3594?32500+083456NASTRAN 中关于 DMIG 的介绍（ DMIGdirect matrixinput at points定义与网格点相关的直接输入的矩阵，矩阵由一列或多列定义，列输入中无 0 元素。表头格式： DMIG ， NAME ， “0” , IFO, TIN, TOUT, POLAR, , NCOL, NAME:矩阵名IFO：IFO=6 指定的矩阵是由 K2GC， M2GC，B2GC 情形控制语句指定的1= 正方形； 9 或 2= 矩形； 6= 对称形TIN :矩阵输入形式1= 实数，单精度（每一元素用一个 fi

10、eld ）2= 实数，双精度（每一元素用一个 field ）3= 复数，单精度（每一元素用两个 field ）4= 复数，双精度（每一元素用两个 field ）TOUT:生成的矩阵类型0= 由精确的系统单元设定（默认）1= 实数，单精度2= 实数，双精度3= 复数，单精度4= 复数，双精度POLAR Ai,Bi的形式（实根为0或black，为实根虚根形式；实根大于 0为振幅相位形式） NCOL矩形矩阵的列数，仅在IF0=9时有效数据格式： DMIG, NAME, GJ, CJ, , G1, C1, A1, B1,G2 , C2,A2,B2, -etc-,GJ：Grid, scalar or extra point identification number for column index （列中网格点、标量点的辨识号）CJ：Component number for grid point GJ. （0 < Integer < 6; blank or zero if GJ is a scalar or extra poin