论坛 › MSC系列软件专区 › 怎么提取总刚度阵和总质量阵

北欧的冬天 发表于 2006-8-30 15:25

怎么提取总刚度阵和总质量阵

怎么提取总刚度阵和总质量阵

如果在patran里建了个模型(比如梁的模型),有什么办法在patran中提取到总刚度阵和总质量阵吗,或者用pcl语言能做到吗,谢谢

Von_Mises 发表于 2006-9-13 09:27

同问

北欧的冬天 发表于 2006-9-15 15:19

Von_Mises
你也碰到同样的问题吗能把你的qq给我吗,我们交流可以

superwolf 发表于 2006-9-15 22:00

加入代码：PARAM,EXTOUT,DMIGPCH,会在结果中得到一个.pch文件，里面包含了所需的总刚度矩阵和总质量矩阵。

北欧的冬天 发表于 2006-9-16 16:33

在什么地方加入代码呢?能详细点吗
我想得到总刚度矩阵和总质量矩阵后接着计算特征值和特征向量可以吗,

多情清秋 发表于 2006-9-16 18:48

质量阵直接可以在材料material菜单中查看。
刚度矩阵需要用到dmap输出

从傲雪论坛上抄来的dmap输出刚振
INIT MASTER(S)
ID C:\Mscvn4w2002\data\nh.MO,MSC.v
SOL SESTATICS
TIME 200
DIAG 8
$
compile semg
alter 'kjjz.*stiffness' $
matprn kjjz// $
$
compile sekr $
alter 'upartn.*kff' $
matprn kff $
$
CEND
TITLE = ma
SUBTITLE = ma
ECHO = NONE
DISPLACEMENT(PRINT) = ALL
SPC = 1
LOAD = 1
BEGIN BULK
$ ***************************************************************************
$ Written by : MSC.visualNastran for Windows
$ Version : 2002
$ Modeller : 8.10
$ Translator : MSC.Nastran
$ From Model : C:\Mscvn4w2002\data\nh.MOD
$ Date : Tue May 20 10:59:55 2003
$ ***************************************************************************
$
PARAM,POST,-1
PARAM,OGEOM,NO
PARAM,AUTOSPC,YES
PARAM,K6ROT,100.
PARAM,MAXRATIO,1.E+8
PARAM,GRDPNT,0
CORD2C 1 0 0. 0. 0. 0. 0. 1.+MSC.VC1
+MSC.VC1 1. 0. 1.
CORD2S 2 0 0. 0. 0. 0. 0. 1.+MSC.VC2
+MSC.VC2 1. 0. 1.
$ MSC.visualNastran for Windows Load Set 1 : NASTRAN 1
FORCE 1 2 0 1. 1000. 0. 0.
$ MSC.visualNastran for Windows Constraint Set 1 : NASTRAN SPC 1
SPC 1 1 123456 0.
$ MSC.visualNastran for Windows Property 1 : Untitled
PROD 1 1 0.2 0. 0. 0.
$ MSC.visualNastran for Windows Material 1 : Untitled
MAT1 1 1.E+73787879. 0.32 0. 0. 0.
GRID 1 0 0. 0. 0. 0
GRID 2 0 5. 0. 0. 0
GRID 3 0 10. 0. 0. 0
CROD 1 1 1 2
CROD 2 1 2 3
ENDDATA 015f757a

上面内容是转贴，本人不会用patran，是否好用希望大家验证一下

superwolf 发表于 2006-9-16 20:35

本帖最后由 wdhd 于 2016-5-11 10:22 编辑

　　原帖由 superwolf 于 2006-9-15 22:00 发表
　　加入代码：PARAM,EXTOUT,DMIGPCH ,会在结果中得到一个.pch文件，里面包含了所需的总刚度矩阵和总质量矩阵。
　　好久没搞了，都忘了，这里有个文件，希望能对你有帮助
　　一个利用DMAP提取质量、刚度矩阵的方法
　　对于一个结构静力的线性分析或非线性分析，总刚度矩阵是每次求解过程中所必须的。如果想再一次分析完成之后在结果中查看结构的总刚度矩阵和总质量矩阵，可以在nastran中加入下列代码：
　　PARAM,EXTOUT,DMIGPCH
　　在patran中直接输入，会在结果中得到一个.pch文件，里面包含了所需的总刚度矩阵和总质量矩阵。
　　若要在求解过程中输出总刚度矩阵或总质量矩阵，就需要用DMAP语言。下面给出一个简单的利用DMAP语言提取单元刚度矩阵，总质量矩阵和总质量矩阵的例子：
　　ID MSC cxh77 DIAG8TIME5DIAG8TIME5 BEGIN (开始)
　　SOL 100
　　MALTER 'MALTER:USERDMAP'
　　GP1 GEOM1,GEOM2,,,/GPL,EQEXIN,GPDT,CSTM,BGPDT,SIL,/S,N,LUSET/0/0 GP2GEOM2,EQEXIN,EPT,,/ECT,GP2GEOM2,EQEXIN,EPT,,/ECT,
　　PLTSET PCDB,BGPDT,ECT/PLTX,PLTPAR,GPSET,ELSET/S,N,NSLLS/S,N,JPLOT CONDP1S,JPLOTCONDP1S,JPLOT
　　LABEL P1S GP3GEOM3,bgpdt,GEOM2,,,,,/SLT,ETT/0/0/0GP3GEOM3,bgpdt,GEOM2,,,,,/SLT,ETT/0/0/0
　　TA1, ,ECT,EPT,BGPDT,SIL,ETT,CSTM,,,/EST,,GEI,GPECT,,/LUSET/-1/0/1/0/0 EMGEST,CSTM,MPT,DIT,GEOM2,,,,,,BGPDT,,,,/KELM,KDICT,MELM,MDICT,BELM,BDICT/S,N,NOKGGEMGEST,CSTM,MPT,DIT,GEOM2,,,,,,BGPDT,,,,/KELM,KDICT,MELM,MDICT,BELM,BDICT/S,N,NOKGG
　　KELM为所需要的单元刚度矩阵
　　EMA GPECT,KDICT,KELM,BGPDT,SIL,CSTM,,/KGG,/
　　组集单元刚度矩阵组集单元刚度矩阵
　　KGG为总刚度矩阵EMAGPECT,MDICT,MELM,BGPDT,SIL,CSTM,,/MGG,/−1/V,Y,WTMASSKGG为总刚度矩阵EMAGPECT,MDICT,MELM,BGPDT,SIL,CSTM,,/MGG,/−1/V,Y,WTMASS
　　MGG为总质量矩阵
　　matprn kelm// printelementstifinessmatricies(打印单元刚度矩阵)matprnkgg//printelementstifinessmatricies(打印单元刚度矩阵)matprnkgg// print global stifiness matrix (打印总刚度矩阵)
　　matprn mgg// printglobalmassmatrix(打印总质量矩阵)ENDALTERLINKUSERDMAPprintglobalmassmatrix(打印总质量矩阵)ENDALTERLINKUSERDMAP 下面为一段具体的结构模型
　　CEND
　　TITLE=
　　SUBTITLE=
　　SUBCASE 1
　　LOAD = 1
　　SPC = 1
　　DISP = ALL
　　STRESS = ALL
　　BEGIN BULK
　　$
　　GRID,1,,0.,0.,0.
　　GRID,2,,5.,0.,0.
　　GRID,3,,10.,0.,0.
　　CROD,1,1,1,2
　　CROD,2,1,2,3
　　PROD,1,1,.2
　　MAT1,1,1.+7,,.32,2.7
　　FORCE 1,2,,1000.,1.,0.,0.
　　SPC1,1,123456,1
　　ENDDATA
　　将之保存为matrix.bdf直接用nastran分析调用，在生成的.f06文件中可以查看输出矩阵：
　　一些片断：
　　0 MATRIX KELM (GINO NAME 101 ) IS A DB PREC 2 COLUMN X 78 ROW RECTANG MATRIX.
　　0COLUMN 1 ROWS 1 THRU 58 --------------------------------------------------
　　ROW
　　1) 4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00
　　ROW POSITIONS 11 THRU 50 NOT PRINTED - ALL ARE NULL.
　　51) 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 4.0000D+05
　　0COLUMN 2 ROWS 1 THRU 58 --------------------------------------------------
　　ROW
　　1) 4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00
　　ROW POSITIONS 11 THRU 50 NOT PRINTED - ALL ARE NULL.
　　51) 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 4.0000D+05
　　0THE NUMBER OF NON-ZERO TERMS IN THE DENSEST COLUMN = 3
　　0THE DENSITY OF THIS MATRIX IS 3.85 PERCENT.
　　1 DECEMBER 6, 2004 MSC.NASTRAN 9/23/03 PAGE 8
　　0
　　0 MATRIX KGG (GINO NAME 101 ) IS A DB PREC 18 COLUMN X 18 ROW SYMMETRC MATRIX.
　　0COLUMN 1 ROWS 1 THRU 7 --------------------------------------------------
　　ROW
　　1) 4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 -4.0000D+05
　　0COLUMNS 2 THRU 6 ARE NULL.
　　0COLUMN 7 ROWS 1 THRU 13 --------------------------------------------------
　　ROW
　　1) -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 8.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00
　　11) 0.0000D+00 0.0000D+00 -4.0000D+05
　　0COLUMNS 8 THRU 12 ARE NULL.
　　0COLUMN 13 ROWS 7 THRU 13 --------------------------------------------------
　　ROW
　　7) -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 4.0000D+05
　　0COLUMNS 14 THRU 18 ARE NULL.
　　0THE NUMBER OF NON-ZERO TERMS IN THE DENSEST COLUMN = 3
　　0THE DENSITY OF THIS MATRIX IS 2.16 PERCENT.
　　1 DECEMBER 6, 2004 MSC.NASTRAN 9/23/03 PAGE 9

ljy_0119 发表于 2006-9-16 21:24

本帖最后由 wdhd 于 2016-5-11 10:22 编辑

　　一个利用DMAP提取质量、刚度矩阵的方法
　　对于一个结构静力的线性分析或非线性分析，总刚度矩阵是每次求解过程中所必须的。如果想再一次分析完成之后在结果中查看结构的总刚度矩阵和总质量矩阵，可以在nastran中加入下列代码：
　　PARAM,EXTOUT,DMIGPCH
　　在patran中直接输入，会在结果中得到一个.pch文件，里面包含了所需的总刚度矩阵和总质量矩阵。
　　若要在求解过程中输出总刚度矩阵或总质量矩阵，就需要用DMAP语言。下面给出一个简单的利用DMAP语言提取单元刚度矩阵，总质量矩阵和总质量矩阵的例子：
　　ID MSC cxh77 DIAG8TIME5DIAG8TIME5 BEGIN (开始)
　　SOL 100
　　MALTER 'MALTER:USERDMAP'
　　GP1 GEOM1,GEOM2,,,/GPL,EQEXIN,GPDT,CSTM,BGPDT,SIL,/S,N,LUSET/0/0 GP2GEOM2,EQEXIN,EPT,,/ECT,GP2GEOM2,EQEXIN,EPT,,/ECT,
　　PLTSET PCDB,BGPDT,ECT/PLTX,PLTPAR,GPSET,ELSET/S,N,NSLLS/S,N,JPLOT CONDP1S,JPLOTCONDP1S,JPLOT
　　LABEL P1S GP3GEOM3,bgpdt,GEOM2,,,,,/SLT,ETT/0/0/0GP3GEOM3,bgpdt,GEOM2,,,,,/SLT,ETT/0/0/0
　　TA1, ,ECT,EPT,BGPDT,SIL,ETT,CSTM,,,/EST,,GEI,GPECT,,/LUSET/-1/0/1/0/0 EMGEST,CSTM,MPT,DIT,GEOM2,,,,,,BGPDT,,,,/KELM,KDICT,MELM,MDICT,BELM,BDICT/S,N,NOKGGEMGEST,CSTM,MPT,DIT,GEOM2,,,,,,BGPDT,,,,/KELM,KDICT,MELM,MDICT,BELM,BDICT/S,N,NOKGG
　　KELM为所需要的单元刚度矩阵
　　EMA GPECT,KDICT,KELM,BGPDT,SIL,CSTM,,/KGG,/
　　组集单元刚度矩阵组集单元刚度矩阵
　　KGG为总刚度矩阵EMAGPECT,MDICT,MELM,BGPDT,SIL,CSTM,,/MGG,/−1/V,Y,WTMASSKGG为总刚度矩阵EMAGPECT,MDICT,MELM,BGPDT,SIL,CSTM,,/MGG,/−1/V,Y,WTMASS
　　MGG为总质量矩阵
　　matprn kelm// printelementstifinessmatricies(打印单元刚度矩阵)matprnkgg//printelementstifinessmatricies(打印单元刚度矩阵)matprnkgg// print global stifiness matrix (打印总刚度矩阵)
　　matprn mgg// printglobalmassmatrix(打印总质量矩阵)ENDALTERLINKUSERDMAPprintglobalmassmatrix(打印总质量矩阵)ENDALTERLINKUSERDMAP 下面为一段具体的结构模型
　　CEND
　　TITLE=
　　SUBTITLE=
　　SUBCASE 1
　　LOAD = 1
　　SPC = 1
　　DISP = ALL
　　STRESS = ALL
　　BEGIN BULK
　　$
　　GRID,1,,0.,0.,0.
　　GRID,2,,5.,0.,0.
　　GRID,3,,10.,0.,0.
　　CROD,1,1,1,2
　　CROD,2,1,2,3
　　PROD,1,1,.2
　　MAT1,1,1.+7,,.32,2.7
　　FORCE 1,2,,1000.,1.,0.,0.
　　SPC1,1,123456,1
　　ENDDATA
　　将之保存为matrix.bdf直接用nastran分析调用，在生成的.f06文件中可以查看输出矩阵：
　　一些片断：
　　0 MATRIX KELM (GINO NAME 101 ) IS A DB PREC 2 COLUMN X 78 ROW RECTANG MATRIX.
　　0COLUMN 1 ROWS 1 THRU 58 --------------------------------------------------
　　ROW
　　1) 4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00
　　ROW POSITIONS 11 THRU 50 NOT PRINTED - ALL ARE NULL.
　　51) 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 4.0000D+05
　　0COLUMN 2 ROWS 1 THRU 58 --------------------------------------------------
　　ROW
　　1) 4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00
　　ROW POSITIONS 11 THRU 50 NOT PRINTED - ALL ARE NULL.
　　51) 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 4.0000D+05
　　0THE NUMBER OF NON-ZERO TERMS IN THE DENSEST COLUMN = 3
　　0THE DENSITY OF THIS MATRIX IS 3.85 PERCENT.
　　1 DECEMBER 6, 2004 MSC.NASTRAN 9/23/03 PAGE 8
　　0
　　0 MATRIX KGG (GINO NAME 101 ) IS A DB PREC 18 COLUMN X 18 ROW SYMMETRC MATRIX.
　　0COLUMN 1 ROWS 1 THRU 7 --------------------------------------------------
　　ROW
　　1) 4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 -4.0000D+05
　　0COLUMNS 2 THRU 6 ARE NULL.
　　0COLUMN 7 ROWS 1 THRU 13 --------------------------------------------------
　　ROW
　　1) -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 8.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00
　　11) 0.0000D+00 0.0000D+00 -4.0000D+05
　　0COLUMNS 8 THRU 12 ARE NULL.
　　0COLUMN 13 ROWS 7 THRU 13 --------------------------------------------------
　　ROW
　　7) -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 4.0000D+05
　　0COLUMNS 14 THRU 18 ARE NULL.
　　0THE NUMBER OF NON-ZERO TERMS IN THE DENSEST COLUMN = 3
　　0THE DENSITY OF THIS MATRIX IS 2.16 PERCENT.
　　1 DECEMBER 6, 2004 MSC.NASTRAN 9/23/03 PAGE 9
　　0

北欧的冬天 发表于 2006-9-17 11:28

谢谢各位的帮忙
请问有没有DMAP的资料,你们列出来的DMAP语言我看不明白

ljy_0119 发表于 2006-9-17 12:49

原帖由 北欧的冬天 于 2006-9-17 11:28 发表
谢谢各位的帮忙
请问有没有DMAP的资料,你们列出来的DMAP语言我看不明白

已经提供资料

北欧的冬天 发表于 2006-9-17 14:51

请问DMAP的程序是写在什么地方的,是和PCL一样写在记事本里,改个扩展名吗
ljy_0119 (aphorodite) (在线)
教研室主任
你有qq吗,我想有什么问题以后可以请教你

北欧的冬天 发表于 2006-9-17 14:56

原帖由 ljy_0119 于 2006-9-17 12:49 发表


已经提供资料

你给的是一些DMAP程序吧有没有介绍DMAP的基础资料,我还完全不会用呢
有没有解特征值问题的DMAP程序

superwolf 发表于 2006-9-17 21:43

http://www.mscsoftware.com/support/online_ex/Nastran/DMAP_ex.cfm       DMAP的基础资料下载

北欧的冬天 发表于 2006-9-17 22:50

你提供下载的这些就是楼上的朋友提供的,
我想要些讲基础的,语法什么的,有吗?还是就是这些

好象msc/nastran user'smanual

msc/nastran programmer's manual上有DMAP的介绍,有谁可以提供这些资料吗

[ 本帖最后由 北欧的冬天 于 2006-9-17 22:53 编辑 ]

ljy_0119 发表于 2006-9-17 23:08

原帖由 北欧的冬天 于 2006-9-17 22:50 发表
你提供下载的这些就是楼上的朋友提供的,
我想要些讲基础的,语法什么的,有吗?还是就是这些

好象msc/nastran user'smanual

msc/nastran programmer's manual上有DMAP的介绍,有谁可以提供这些资料吗


那试试这个，我在发一个，是dmap的guide，你可以参考，我已经发到下载中心了

[ 本帖最后由 ljy_0119 于 2006-9-17 23:10 编辑 ]

查看完整版本: 怎么提取总刚度阵和总质量阵