求助如何在结果中输出单元刚度矩阵

maxiangsen

在nastran中如何实现，谢谢

FEM-Fan

我在好几年前干过，用DMAP编程。不过现在要重新寻找资料。

maxiangsen

帮帮忙,希望提供资料谢谢

yamette

2楼的高人，可否指点一下，我现在也需要导出结构的质量和刚度矩阵，但是在f06文件里只能看见正教规范化后的模态质量阵

重阳

在msc网站里有几个例子，可以打印刚阵。

FEM-Fan

在DMAP中适当的位置插入
MATPCH MGG,KGG,,,//$ 即可在.PCH文件里输出整体矩阵.或用OUTPUT2输出二进制文件。
至于maxiangsen要求的单元刚度矩阵, 我记得是用KDICT和KELM数据块来储存的,在DMAP中打印它们即可输出。
在情况控制卡中插入DIAG 14,(一般可在SOL***后面) 即可输出你所用的版本的NASTRAN的DMAP序列, 找到所需的数据块出现的位置, 在其后面一行插入上述的指令.

maxiangsen

好，我试试，感谢！！！！！！！

maxiangsen

这是Nastran上的一个例子，可我运行总是有错
USER FATAL MESSAGE 807 (XALTER)
     THERE IS A SYNTAX ERROR IN FIRST ARGUMENT OF ALTER STATEMENT.
  
SOL 101
TIME 20
DIAG 8
$
compile semg
alter ’kjjz.*stiffness’ $
matprn kjjz// $
$
compile sekr $
alter ’upartn.*kff’ $
matprn kff $
$
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
FORCE,1,2,,1000.,1.,0.,0.
SPC1,1,123456,1
ENDDATA

maxiangsen

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
这个终于可以了，我的是MSC.visualNastran for Windows，所以缺了以下这个是不行的
$ ***************************************************************************
$   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
$ ***************************************************************************