|
楼主 |
发表于 2005-8-24 08:54
|
显示全部楼层
回复:(ecass)[分享]Patran技巧及其常见问题-Patran...
本帖最后由 wdhd 于 2016-3-23 11:01 编辑
31
Q:Patran 如何定义材料库
A:Patran除了可以直接读取MSC.Mvision的材料材料库外,还可利用执行Session File的方式,直接使用以前已经定义过的材料。 编辑patran.ses.xx,将定义材料的PCL指令剪下,
粘贴到另外一个文件中(如mat.ses)。 之后便可直接由Patran的File/ Session/ Play来加入此一材料的定义。
也可以直接加入Patran菜单的中:把刚刚定义的session file 复制到 C:\MSC\patran2001r3\shareware\msc\unsupported\utilities\data_files\bv_material_data , 之后就会在Utilities/ Material/ Material Session File Library 中出现刚才的材料名称。.
32
Q:Patran 的完整信息输出:
A:在执行Patran的时候出现齐怪的错误信息时, 可以先把Patran关闭, 接着启动DOS窗口,
在DOS下直接输入patran -stdout ,执行再重新启动Patran, 会多了一个信息窗口。
33
Q:超单元
A:超单元的应用跟模态分析没有关系,它只是有限元中的一种分支方法,用于解决运算量过大的问题的,现在基本上已经不用了!
34
Q:约束
A:我Patran/Nastran在做项目时,对结构应力进行分析,发现在结构被约束的部位往往产生比较大的应力,甚至是最大的应力值就在这些部位出现,好像这种结果并不太符合实际情况。请问约束该如何加,才能消除这种影响?
这种现象很正常的,只要你的约束反映了实际情况,约束部位的过大应力你可根据圣维南原理视而不见。
圣维南原理看看力学方面的书!如果一定要把约束区域分析准确的话,你的约束一定要尽可能地与真实情况相符,航空经常用的方法是加弹簧元,弹簧元的刚度是用的经验值或者试验值!
35
Q:MSC.Flds 2004
A:计算飞行载荷-flightload and dynamics,主要用来计算气动弹性问题。
36
Q:如何找到频谱分析的响应最大点?
A:在Tools--List---Create就可以将你想要的点的结果导入到一个组中,然后可以从结果文件中调出再在origin中编辑,或是直接显示。科研中国SciEi.com整理.
37
Q:自由模态
A:所谓自由模态,就是无约束模态,前6阶都是刚体模态,
38
Q:能否设置计算结束后不输出f06文件?
A:可以设置的,Solution Parameters里面Maximum Printed Lines =就是f06输出的最大行数
38
Q:一个利用DMAP提取质量、刚度矩阵的方法
A:对于一个结构静力的线性分析或非线性分析,总刚度矩阵是每次求解过程中所必须的。如果想再一次分析完成之后在结果中查看结构的总刚度矩阵和总质量矩阵,可以在nastran中加入下列代码:PARAM,EXTOUT,DMIGPCH
在patran中直接输入,会在结果中得到一个.pch文件,里面包含了所需的总刚度矩阵和总质量矩阵。
若要在求解过程中输出总刚度矩阵或总质量矩阵,就需要用DMAP语言。下面给出一个简单的利用DMAP语言提取单元刚度矩阵,总质量矩阵和总质量矩阵的例子:
ID MSC cxh77 $
DIAG 8
TIME 5
$ BEGIN (开始)
SOL 100
MALTER 'MALTER:USERDMAP'
GP1 GEOM1,GEOM2,,,/GPL,EQEXIN,GPDT,CSTM,BGPDT,SIL,/S,N,LUSET/0/0 $
GP2 GEOM2,EQEXIN,EPT,,/ECT, $
PLTSET PCDB,BGPDT,ECT/PLTX,PLTPAR,GPSET,ELSET/S,N,NSLLS/S,N,JPLOT $
COND P1S,JPLOT $
LABEL P1S $
GP3 GEOM3,bgpdt,GEOM2,,,,,/SLT,ETT/0/0/0 $
TA1, ,ECT,EPT,BGPDT,SIL,ETT,CSTM,,,/EST,,GEI,GPECT,,/LUSET/-1/0/1/0/0 $
EMG EST,CSTM,MPT,DIT,GEOM2,,,,,,BGPDT,,,,/
KELM,KDICT,MELM,MDICT,BELM,BDICT/
S,N,NOKGG $
$
$ KELM为所需要的单元刚度矩阵
EMA GPECT,KDICT,KELM,BGPDT,SIL,CSTM,,/KGG,/ $
$
$ 组集单元刚度矩阵
$
$ KGG为总刚度矩阵
EMA GPECT,MDICT,MELM,BGPDT,SIL,CSTM,,/MGG,/-1/V,Y,WTMASS $
$
$ MGG为总质量矩阵
matprn kelm// $ print element stifiness matricies (打印单元刚度矩阵)
matprn kgg// $ print global stifiness matrix (打印总刚度矩阵)
matprn mgg// $ print global mass matrix (打印总质量矩阵)
ENDALTER
LINK USERDMAP
$ 下面为一段具体的结构模型
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
39
Q:MSC.Nastran可形成如下输出文件
A:
Job_name.DBALL 包含数据库运行的永久数据
Job_name.MASTER 数据库运行的总词典
Job_name.F04 包含数据库文件信息和模块执行摘要
Job_name.F06 包含MSC.Nastran的分析结果,为文本文件。
Job_name.LOG 包含系统信息和系统错误信息
Job_name.OP2 MSC.Nastran的分析结果,为二进制文件
Job_name.XDB MSC.Nastran的分析结果,为二进制文件
24、GRID POINT ID DEGREE OF FREEDOM MATRIX/FACTOR DIAGONAL RATIO MATRIX DIAGONAL
102 T1 -2.78284E+14 3.28730E+01
102 T2 -9.16132E+13 3.28730E+01
113 R3 -8.35181E+08 2.72727E-03
USER FATAL MESSAGE 9050 (SEKRRS)
RUN TERMINATED DUE TO EXCESSIVE PIVOT RATIOS IN MATRIX KLL
USER ACTION: CONSTRAIN MECHANISMS WITH SPCI OR SUPORTI ENTRIES OR SPECIFY PARAM,BAILOUT,-1 TO
在bdf文件中,找到solution parameter字段,加入“PARAM,BAILOUT,-1”再提交计算!
不过这种方法用于调试局部模型可以,调试总体模型是不能用的。我们工作的规定文件中明确说了,计算总体模型不能够用bailout卡,也不能打开automatic constraint选项!
你的加载板面的1,2,6自由度没有刚度(即加载板没有约束刚体位移)
你可以任意选择一点加
spc1,1,126,**
同时在BEGIN BULK卡下一行加卡片
PARAM,K6ROT,100.0
即可
40
Q:弹簧
A:弹簧的位置由两个节点定义,弹簧是只能承受指定自由度方向刚度的单元,所以在定义弹簧的时候,还必须指定弹簧的自由度方向(UX,UY等)比如弹簧只能拉压是在x方向,则定义DOF为UX。上述所说的自由度依赖于你定义节点时用的坐标系,也就是节点若定义在总体坐标下,则DOF也是总体系下的方向。
上述所说的自由度依赖于你定义节点时用的坐标系,也就是节点若定义在总体坐标下,则DOF也是总体系下的方向。这个节点的坐标系是可以改的,可以用element->modify->node->edit之下就有改变坐标系的命令。
对于一维的弹簧单元,两端头的自由度选择,我认为是为它选的可以释放的(即可自由的)的自由度。可能描述的不是很清楚,你去IDEAS9.0的帮助文件里面看看,关于弹簧单元的解释比较详细。
在定义弹簧单元的property时,有两个框框dof at node 1:UX,dof at node 2:UX:节点一的自由度UX(约束X方向的平动),节点2的自由度UX约束;除了x方向约束(受刚度影响)外,其他自由
CURVE是有方向性的,,每条CURVE生成的时候方向是固定的,这个可以从GEOMETRY里面进行确认。以CURVE生成网格的时候,沿着方向依次是NODE1和NODE2。如果直接生成一维单元先选的是NODE1后选的是NODE2。科研中国SciEi.com整理.
模拟一个y方向可被压缩弹簧:在用SPRING属性的时候两个自由都选则UY,然后在边界条件里将你需要的节点全约束就可以了。STRING 是NASTRAN的使用的一个字符串。在这里指SPRING的方向。
零维弹簧指的是GROUND SPRING即地弹簧比如在结构的某个部位需要施加一个对地的弹性约束,就可以直接使用这种单元属性。当然用一维的同样可以模拟,用零维的比较方便而已。
定义弹簧元:property的1D单元中有一个项是弹簧,只要再输入这个弹簧单元的方向和刚度值即可。注意弹簧单元的性质只取决于作用方向和刚度系数这两个参数,而于单元的两个节点位置无关。
|
|