|
本帖最后由 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
|
评分
-
1
查看全部评分
-
|