PATRAN的一些精华小技巧.docx
《PATRAN的一些精华小技巧.docx》由会员分享,可在线阅读,更多相关《PATRAN的一些精华小技巧.docx(26页珍藏版)》请在冰豆网上搜索。
![PATRAN的一些精华小技巧.docx](https://file1.bdocx.com/fileroot1/2022-12/1/327ba5a0-3551-4ac7-9af3-1f845b36982f/327ba5a0-3551-4ac7-9af3-1f845b36982f1.gif)
PATRAN的一些精华小技巧
PATRAN的一些精华小技巧
1、在Patran里如何Move一组Points的位置,而不改变这组Points的ID编号?
Group/Transform/Translate的功能,
这样不但编号不会变,连property跟边界条件都会保留
2、Patran如何执行多次Undo?
所有Patran的操作步骤,都记录在最新的一个patran.ses.xx中,
如果需要多次undo,可以刪除最后不需要的步骤指令行,
再利用File->Session->Play的方式,执行改过的patran.ses.xx,
这样可以无限制的undo。
3、Patran中如何定义杆件之间的铰接?
用的是rod元素,不需要定义铰接,因为rod元间本身就是以铰接形式连接。
如果用bar或是beam,需在properties里的PinnedDOFs@Node1/2定义铰接。
4、Patran如何把不小心Equivalence的node分开
用Utility/FEM-Elements/SeparateElements
在equivalence时,可以將选项切换为"List",
只选择特定某些节点作equivalence,可以避免不小心的失誤。
5、Patran如何將Tri3单元转换为Tri6单元
在Patran->Element->Modify/Element/Edit,將Type选项打勾,
在Shape中选Tri,NewShape选Tri6,最后选取想要改变的Tri3单元,
6、Patran如何定义材料库
Patran除了可以直接读取MSC.Mvision的材料材料库外,
还可利用执行SessionFile的方式,直接使用以前已经定义过的材料。
编辑patran.ses.xx,将定义材料的PCL指令剪下,
粘贴到另外一个文件中(如mat.ses)。
之后便可直接由Patran的File/Session/Play来加入此一材料的定义。
也可以直接加入Patran菜单的中:
把刚刚定义的sessionfile复制到
C:
\MSC\patran2001r3\shareware\msc\unsupported\utilities\data_files\bv_material_data,
之后就会在Utilities/Material/MaterialSessionFileLibrary中出现刚才的材料名称。
.
7、Patran的完整信息输出:
在执行Patran的时候出现齐怪的错误信息时,
可以先把Patran关闭,接着启动DOS窗口,
在DOS下直接输入patran-stdout,执行
再重新启动Patran,会多了一个信息窗口。
3、计算完毕后,只想显示应力超过某个值的单元,而其它单元不想显示,如何设置?
tool-list-creat,方法选attribute,设f>你要求的应力,apply以后选中在list里面的即为你要求的.再用plot/erase不显示你不要的单元。
4、一个四边形平板,一端的边上节点6个自由度全约束住,另外一端上几个节点上施加z方向强迫位移<,,1E-5>,没有别的条件。
计算完以后看F06文件,看看那些节点的位移是否加上了!
用的是loads/BCs中的creat-diaplacement,我很奇怪的是:
我试了几次这个强迫位移值,如0.1,0.01,0.001,0.0001,f06文件中显示正确,节点位移值确实就是输入值!
但是这个值在变小时如1e-5,1e-6,F06文件显示结果为0感觉好像是nastran的识别问题,把10的-5次方一下的数默认为0!
问题出在translationparameter里面的一个参数numerical,帮助文件里面说它用于比较两个数是否相等,其默认为1e-4.writing才是判断一个数是否近似为零,默认为1e-21.但实际上当你给出的强迫位移量小于numerical时,它就认为近似为0,在bdf文件中就给忽略掉了.你修改numerical为1e-5,你上面说的1e-5就可以算了。
5、自重怎么加到模型上去?
自重是在load/BCS里加的create->inertialload->element在inputdata->load/BCsetscaleFactor[输入加速度的值一般取9.8]TransAccel(a1a2a3)<0-10>(力是沿Y轴向下)后就ok了
tools下面有个massproperties是计算模型质量和惯量的,不知对你有没有帮助
/solu
/output,mass,txt
psolve,elform
/output
finish
将在当前目录下生成mass.txt文件里面会么都有。
6、共振时的应力
模态分析得到的位移是模态位移,没有输入载荷时的结果其值没有意义,同一模态下模态坐标位移之间的比值才有意义,因此,通过NORMALMODES分析的结果不表示真实的位移以及应力等。
需要计算FREQENCYRESPONSE或者TRANSIENTRESPONSE才能获得真实应力。
7、体单元的节点没有转动自由度,只需考虑三个平移自由度。
8、我在用PATRAN做计算时,分析出的结果文件很大,往往提示磁盘空间不足,结果提示需要980M空间。
你试试不用patran直接递交分析,进入nastran用手工递交,可输入scr=yes之后进行分析。
因为直接递交会产生个dball文件很大的,没什么大用,用手递交后,用scr=yes就可以不用产生这个文件了,另外所产生的dball文件可以delete它,它占的空间太大了
在patran下不进行分析,不用fullrun用analysicdeck之后运行nastran.ex桌面没有在msc文件夹下找到,用它开打你刚才生成的文件就行在最底一行输入scr=yes就行,之后和在patran下直接调用nastran的界面一样。
patran向nastran递交运算时的参数可以在nastxxxx.rcf文件中设定,这个文件可以用文本编辑软件修改,加一句scr=yes就行了,运算式通用的参数都可以在这里面定制。
9、做模态分析需要密度。
10、如何画椭圆
1) Createacircle,thenuseTransform/Curve/Scaletostretchitinonedirection.
TheeasiestwaytocreateanellipseiswiththeUtilitiesmenuofPatran.Utilities--Geometry--Createcurvebyfunction
Exemple:
Ifyouwanttocreateanellipsewitha=110andb=90yourequationwillbe:
Y=SQRT(90**2*(1-(`X**2/110**2)))
2) 画圆,将圆旋转一个角度
具体旋转角度:
设a为长轴长;b为短轴长
则满足关系cos(angle)=(b/a)
旋转该角度后,将其投影到园形原来的所在平面即得所需椭圆
3) 利用局部坐标系
创建局部坐标系,选cylindrical1,设置长短轴即可
4) 用curve?
一?
?
,然後在move/scale?
不同比例,?
?
curve就可以.
5) 据说通过creat-curve-conic可以画椭圆
11、MPC
MPC是一个很好的工具
用来拟和多点自由度的
比如
要对一个截面上的多个点施加位移或者力
特别是对整个截面施加弯矩
就可以使用MPC
把多个点的位移
力
弯矩施加到一个点上就可以了
特别是对弯矩
你知道,对实体模型的横截面施加一个弯矩
如果没有MPC的话
会有多么地困难
12、模态分析的时候是不能够分析应力的,分析应力要单独用别的求解类型再分析。
14、注意转速的单位是n转/s,不是弧度
15、一些小技巧
1、在Patran里如何Move一组Points的位置,而不改变这组Points的ID编号?
Group/Transform/Translate的功能,
这样不但编号不会变,连property跟边界条件都会保留
2、Patran如何执行多次Undo?
所有Patran的操作步骤,都记录在最新的一个patran.ses.xx中,
如果需要多次undo,可以刪除最后不需要的步骤指令行,
再利用File->Session->Play的方式,执行改过的patran.ses.xx,
这样可以无限制的undo。
3、Patran中如何定义杆件之间的铰接?
用的是rod元素,不需要定义铰接,因为rod元间本身就是以铰接形式连接。
如果用bar或是beam,需在properties里的PinnedDOFs@Node1/2定义铰接。
4、Patran如何把不小心Equivalence的node分开
用Utility/FEM-Elements/SeparateElements
在equivalence时,可以將选项切换为"List",
只选择特定某些节点作equivalence,可以避免不小心的失誤。
5、Patran如何將Tri3单元转换为Tri6单元
在Patran->Element->Modify/Element/Edit,將Type选项打勾,
在Shape中选Tri,NewShape选Tri6,最后选取想要改变的Tri3单元,
6、Patran如何定义材料库
Patran除了可以直接读取MSC.Mvision的材料材料库外,
还可利用执行SessionFile的方式,直接使用以前已经定义过的材料。
编辑patran.ses.xx,将定义材料的PCL指令剪下,
粘贴到另外一个文件中(如mat.ses)。
之后便可直接由Patran的File/Session/Play来加入此一材料的定义。
也可以直接加入Patran菜单的中:
把刚刚定义的sessionfile复制到
C:
\MSC\patran2001r3\shareware\msc\unsupported\utilities\data_files\bv_material_data,
之后就会在Utilities/Material/MaterialSessionFileLibrary中出现刚才的材料名称。
.
7、Patran的完整信息输出:
在执行Patran的时候出现齐怪的错误信息时,
可以先把Patran关闭,接着启动DOS窗口,
在DOS下直接输入patran-stdout,执行
再重新启动Patran,会多了一个信息窗口。
16、超单元
超单元的应用跟模态分析没有关系,它只是有限元中的一种分支方法,用于解决运算量过大的问题的,现在基本上已经不用了!
17、约束
我Patran/Nastran在做项目时,对结构应力进行分析,发现在结构被约束的部位往往产生比较大的应力,甚至是最大的应力值就在这些部位出现,好像这种结果并不太符合实际情况。
请问约束该如何加,才能消除这种影响?
这种现象很正常的,只要你的约束反映了实际情况,约束部位的过大应力你可根据圣维南原理视而不见。
圣维南原理看看力学方面的书!
如果一定要把约束区域分析准确的话,你的约束一定要尽可能地与真实情况相符,航空经常用的方法是加弹簧元,弹簧元的刚度是用的经验值或者试验值!
18、MSC.Flds2004
计算飞行载荷-flightloadanddynamics,主要用来计算气动弹性问题。
19、如何找到频谱分析的响应最大点?
在Tools--List---Create就可以将你想要的点的结果导入到一个组中,然后可以从结果文件中调出再在origin中编辑,或是直接显示。
20、自由模态
所谓自由模态,就是无约束模态
前6阶都是刚体模态
21、能否设置计算结束后不输出f06文件?
可以设置的,SolutionParameters里面MaximumPrintedLines=
就是f06输出的最大行数
22、一个利用DMAP提取质量、刚度矩阵的方法
对于一个结构静力的线性分析或非线性分析,总刚度矩阵是每次求解过程中所必须的。
如果想再一次分析完成之后在结果中查看结构的总刚度矩阵和总质量矩阵,可以在nastran中加入下列代码:
PARAM,EXTOUT,DMIGPCH
在patran中直接输入,会在结果中得到一个.pch文件,里面包含了所需的总刚度矩阵和总质量矩阵。
若要在求解过程中输出总刚度矩阵或总质量矩阵,就需要用DMAP语言。
下面给出一个简单的利用DMAP语言提取单元刚度矩阵,总质量矩阵和总质量矩阵的例子:
IDMSCcxh77$
DIAG8
TIME5
$BEGIN(开始)
SOL100
MALTER'MALTER:
USERDMAP'
GP1GEOM1,GEOM2,,,/GPL,EQEXIN,GPDT,CSTM,BGPDT,SIL,/S,N,LUSET/0/0$
GP2GEOM2,EQEXIN,EPT,,/ECT,$
PLTSETPCDB,BGPDT,ECT/PLTX,PLTPAR,GPSET,ELSET/S,N,NSLLS/S,N,JPLOT$
CONDP1S,JPLOT$
LABELP1S$
GP3GEOM3,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,NOKGG$
$
$KELM为所需要的单元刚度矩阵
EMAGPECT,KDICT,KELM,BGPDT,SIL,CSTM,,/KGG,/$
$
$组集单元刚度矩阵
$
$KGG为总刚度矩阵
EMAGPECT,MDICT,MELM,BGPDT,SIL,CSTM,,/MGG,/-1/V,Y,WTMASS$
$
$MGG为总质量矩阵
matprnkelm//$printelementstifinessmatricies(打印单元刚度矩阵)
matprnkgg//$printglobalstifinessmatrix(打印总刚度矩阵)
matprnmgg//$printglobalmassmatrix(打印总质量矩阵)
ENDALTER
LINKUSERDMAP
$下面为一段具体的结构模型
CEND
TITLE=
SUBTITLE=
SUBCASE1
LOAD=1
SPC=1
DISP=ALL
STRESS=ALL
BEGINBULK
$
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
FORCE1,2,,1000.,1.,0.,0.
SPC1,1,123456,1
ENDDATA
将之保存为matrix.bdf直接用nastran分析调用,在生成的.f06文件中可以查看输出矩阵:
一些片断:
0MATRIXKELM(GINONAME101)ISADBPREC2COLUMNX78ROWRECTANGMATRIX.
0COLUMN1ROWS1THRU58--------------------------------------------------
ROW
1)4.0000D+050.0000D+000.0000D+000.0000D+000.0000D+000.0000D+00-4.0000D+050.0000D+000.0000D+000.0000D+00
ROWPOSITIONS11THRU50NOTPRINTED-ALLARENULL.
51)0.0000D+000.0000D+000.0000D+000.0000D+000.0000D+000.0000D+000.0000D+004.0000D+05
0COLUMN2ROWS1THRU58--------------------------------------------------
ROW
1)4.0000D+050.0000D+000.0000D+000.0000D+000.0000D+000.0000D+00-4.0000D+050.0000D+000.0000D+000.0000D+00
ROWPOSITIONS11THRU50NOTPRINTED-ALLARENULL.
51)0.0000D+000.0000D+000.0000D+000.0000D+000.0000D+000.0000D+000.0000D+004.0000D+05
0THENUMBEROFNON-ZEROTERMSINTHEDENSESTCOLUMN=3
0THEDENSITYOFTHISMATRIXIS3.85PERCENT.
1DECEMBER6,2004MSC.NASTRAN9/23/03PAGE8
0
0MATRIXKGG(GINONAME101)ISADBPREC18COLUMNX18ROWSYMMETRCMATRIX.
0COLUMN1ROWS1THRU7--------------------------------------------------
ROW
1)4.0000D+050.0000D+000.0000D+000.0000D+000.0000D+000.0000D+00-4.0000D+05
0COLUMNS2THRU6ARENULL.
0COLUMN7ROWS1THRU13--------------------------------------------------
ROW
1)-4.0000D+050.0000D+000.0000D+000.0000D+000.0000D+000.0000D+008.0000D+050.0000D+000.0000D+000.0000D+00
11)0.0000D+000.0000D+00-4.0000D+05
0COLUMNS8THRU12ARENULL.
0COLUMN13ROWS7THRU13--------------------------------------------------
ROW
7)-4.0000D+050.0000D+000.0000D+000.0000D+000.0000D+000.0000D+004.0000D+05
0COLUMNS14THRU18ARENULL.
0THENUMBEROFNON-ZEROTERMSINTHEDENSESTCOLUMN=3
0THEDENSITYOFTHISMATRIXIS2.16PERCENT.
1DECEMBER6,2004MSC.NASTRAN9/23/03PAGE9
0
23、MSC.Nastran可形成如下输出文件
Job_name.DBALL包含数据库运行的永久数据
Job_name.MASTER数据库运行的总词典
Job_name.F04包含数据库文件信息和模块执行摘要
Job_name.F06包含MSC.Nastran的分析结果,为文本文件。
Job_name.LOG包含系统信息和系统错误信息
Job_name.OP2MSC.Nastran的分析结果,为二进制文件
Job_name.XDBMSC.Nastran的分析结果,为二进制文件
24、GRIDPOINTIDDEGREEOFFREEDOMMATRIX/FACTORDIAGONALRATIOMATRIXDIAGONAL
102T1-2.78284E+143.28730E+01
102T2-9.16132E+133.28730E+01
113R3-8.35181E+082.72727E-03
USERFATALMESSAGE9050(SEKRRS)
RUNTERMINATEDDUETOEXCESSIVEPIVOTRATIOSINMATRIXKLL
USERACTION:
CONSTRAINMECHANISMSWITHSPCIORSUPORTIENTRIESORSPECIFYPARAM,BAILOUT,-1TO
在bdf文件中,找到solutionparameter字段,加入“PARAM,BAILOUT,-1”再提交计算!
不过这种方法用于调试局部模型可以,调试总体模型是不能用的。
我们工作的规定文件中明确说了,计算总体模型不能够用bailout卡,也不能打开automaticconstraint选项!
你的加载板面的1,2,6自由度没有刚度(即加载板没有约束刚体位移)
你可以任意选择一点加
spc1,1,126,**
同时在BEGINBULK卡下一行加卡片
PARAM,K6ROT,100.0
即可
25、弹簧
弹簧的位置由两个节点定义,弹簧是只能承受指定自由度方向刚度的单元,所以在定义弹簧的时候,还必须指定弹簧的自由度方向(UX,UY等)比如弹簧只能拉压是在x方向,则定义DOF为UX。
上述所说的自由度依赖于你定义节点时用的坐标系,也就是节点若定义在总体坐标下,则DOF也是总体系下的方向。
上述所说的自由度依赖于你定义节点时用的坐标系,也就是节点若定义在总体坐标下,则DOF也是总体系下的方向。
这个节点的坐标系是可以改的,可以用element->modify->node->edit之下就有改变坐标系的命令。
对于一维的弹簧单元,两端头的自由度选择,我认为是为它选的可