声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 9104|回复: 28

[Patran专区] 怎么提取总刚度阵和总质量阵

[复制链接]
发表于 2006-8-30 15:25 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

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

使用道具 举报

发表于 2006-9-13 09:27 | 显示全部楼层
同问
 楼主| 发表于 2006-9-15 15:19 | 显示全部楼层
Von_Mises
你也碰到同样的问题吗能把你的qq给我吗,我们交流可以
发表于 2006-9-15 22:00 | 显示全部楼层
加入代码:PARAM,EXTOUT,DMIGPCH  ,会在结果中得到一个.pch文件,里面包含了所需的总刚度矩阵和总质量矩阵。

点评

赞成: 3.0
赞成: 3
  发表于 2014-3-28 11:07

评分

1

查看全部评分

 楼主| 发表于 2006-9-16 16:33 | 显示全部楼层
在什么地方加入代码呢?能详细点吗
我想得到总刚度矩阵和总质量矩阵后接着计算特征值和特征向量可以吗,
发表于 2006-9-16 18:48 | 显示全部楼层
质量阵直接可以在材料material菜单中查看。
刚度矩阵需要用到dmap输出

从傲雪论坛上抄来的dmap输出刚振
  1. INIT MASTER(S)  
  2. ID C:\Mscvn4w2002\data\nh.MO,MSC.v  
  3. SOL SESTATICS  
  4. TIME 200  
  5. DIAG 8  
  6. $  
  7. compile semg  
  8. alter 'kjjz.*stiffness' $  
  9. matprn kjjz// $  
  10. $  
  11. compile sekr $  
  12. alter 'upartn.*kff' $  
  13. matprn kff $  
  14. $  
  15. CEND  
  16.   TITLE = ma  
  17.   SUBTITLE = ma  
  18.   ECHO = NONE  
  19.   DISPLACEMENT(PRINT) = ALL  
  20.   SPC = 1  
  21.   LOAD = 1  
  22. BEGIN BULK  
  23. $ ***************************************************************************  
  24. $ Written by : MSC.visualNastran for Windows  
  25. $ Version : 2002  
  26. $ Modeller : 8.10  
  27. $ Translator : MSC.Nastran  
  28. $ From Model : C:\Mscvn4w2002\data\nh.MOD  
  29. $ Date : Tue May 20 10:59:55 2003  
  30. $ ***************************************************************************  
  31. $  
  32. PARAM,POST,-1  
  33. PARAM,OGEOM,NO  
  34. PARAM,AUTOSPC,YES  
  35. PARAM,K6ROT,100.  
  36. PARAM,MAXRATIO,1.E+8  
  37. PARAM,GRDPNT,0  
  38. CORD2C 1 0 0. 0. 0. 0. 0. 1.+MSC.VC1  
  39. +MSC.VC1 1. 0. 1.  
  40. CORD2S 2 0 0. 0. 0. 0. 0. 1.+MSC.VC2  
  41. +MSC.VC2 1. 0. 1.  
  42. $ MSC.visualNastran for Windows Load Set 1 : NASTRAN 1  
  43. FORCE 1 2 0 1. 1000. 0. 0.  
  44. $ MSC.visualNastran for Windows Constraint Set 1 : NASTRAN SPC 1  
  45. SPC 1 1 123456 0.  
  46. $ MSC.visualNastran for Windows Property 1 : Untitled  
  47. PROD 1 1 0.2 0. 0. 0.  
  48. $ MSC.visualNastran for Windows Material 1 : Untitled  
  49. MAT1 1 1.E+73787879. 0.32 0. 0. 0.  
  50. GRID 1 0 0. 0. 0. 0  
  51. GRID 2 0 5. 0. 0. 0  
  52. GRID 3 0 10. 0. 0. 0  
  53. CROD 1 1 1 2  
  54. CROD 2 1 2 3  
  55. ENDDATA 015f757a
复制代码


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

点评

赞成: 2.0
赞成: 2
  发表于 2014-3-26 16:56
发表于 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

点评

赞成: 3.0
赞成: 3
  发表于 2014-3-26 16:55

评分

1

查看全部评分

发表于 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

点评

赞成: 5.0
赞成: 5
  发表于 2014-3-26 16:56

评分

1

查看全部评分

 楼主| 发表于 2006-9-17 11:28 | 显示全部楼层
谢谢各位的帮忙
请问有没有DMAP的资料,你们列出来的DMAP语言我看不明白
发表于 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程序
发表于 2006-9-17 21:43 | 显示全部楼层
 楼主| 发表于 2006-9-17 22:50 | 显示全部楼层
你提供下载的这些就是楼上的朋友提供的,
我想要些讲基础的,语法什么的,有吗?还是就是这些

好象msc/nastran user'smanual

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

[ 本帖最后由 北欧的冬天 于 2006-9-17 22:53 编辑 ]
发表于 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 编辑 ]
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-10 15:53 , Processed in 0.068815 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表