patrannastran问题集锦Word文档格式.docx

patrannastran问题集锦Word文档格式.docx

《patrannastran问题集锦Word文档格式.docx》由会员分享，可在线阅读，更多相关《patrannastran问题集锦Word文档格式.docx（28页珍藏版）》请在冰点文库上搜索。

我试了几次这个强迫位移值，如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文件，里面包含了所需的总刚度矩阵和总质量矩阵。

若要在求解过程中输出总刚度矩阵或总质量矩阵，就需要用ＤＭＡＰ语言。

下面给出一个简单的利用ＤＭＡＰ语言提取单元刚度矩阵，总质量矩阵和总质量矩阵的例子：

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

将之保存为ｍａｔｒｉｘ．ｂｄｆ直接用ｎａｓｔｒａｎ分析调用，在生成的．ｆ０６文件中可以查看输出矩阵：

一些片断：

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--------------------------------------------------

0THENUMBEROFNON-ZEROTERMSINTHEDENSESTCOLUMN=3

0THEDENSITYOFTHISMATRIXIS3.85PERCENT.

1DECEMBER6,2004MSC.NASTRAN9/23/03PAGE8

0

0MATRIXKGG（GINONAME101）ISADBPREC18COLUMNX18ROWSYMMETRCMATRIX.

0COLUMN1ROWS1THRU7--------------------------------------------------

1）4.0000D+050.0000D+000.0000D+000.0000D+000.0000D+000.0000D+00-4.0000D+05

0COLUMNS2THRU6ARENULL.

0COLUMN7ROWS1THRU13--------------------------------------------------

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--------------------------------------------------

7）-4.0000D+050.0000D+000.0000D+000.0000D+000.0000D+000.0000D+004.0000D+05

0COLUMNS14THRU18ARENULL.

0THEDENSITYOFTHISMATRIXIS2.16PERCENT.

1DECEMBER6,2004MSC.NASTRAN9/23/03PAGE9

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也是总体系下的方向。

这个节点的坐标系是可以改的，可以用element-&

modify-&

node-&

edit之下就有改变坐标系的命令。

对于一维的弹簧单元，两端头的自由度选择，我认为是为它选的可以释放的（即可自由的）的自由度。

可能描述的不是很清楚，你去IDEAS9.0的帮助文件里面看看，关于弹簧单元的解释比较详细。

在定义弹簧单元的property时，有两个框框dofatnode1:

UX,dofatnode2:

UX：

节点一的自由度UX（约束X方向的平动），节点2的自由度UX约束；

除了x方向约束（受刚度影响）外，其他自由

CURVE是有方向性的，，每条CURVE生成的时候方向是固定的，这个可以从GEOMETRY里面进行确认。

以CURVE生成网格的时候，沿着方向依次是NODE1和NODE2。

如果直接生成一维单元先选的是NODE1后选的是NODE2。

模拟一个y方向可被压缩弹簧：

在用SPRING属性的时候两个自由都选则UY，然后在边界条件里将你需要的节点全约束就可以了。

STRING是NASTRAN的使用的一个字符串。

在这里指SPRING的方向。

零维弹簧指的是GROUNDSPRING即地弹簧比如在结构的某个部位需要施加一个对地的弹性约束，就可以直接使用这种单元属性。

当然用一维的同样可以模拟，用零维的比较方便而已。

定义弹簧元：

property的1D单元中有一个项是弹簧，只要再输入这个弹簧单元的方向和刚度值即可。

注意弹簧单元的性质只取决于作用方向和刚度系数这两个参数，而于单元的两个节点位置无关。

Q：

对于有限元或MSC.NASTRAN的初学者，MSC.NASTRAN有哪些入门手册？

A：

MSC.NASTRANGettingStarte