patran&nastran问题集锦 1/关于软件安装的问题去掉了 2、请问在PATRAN中输出图片能将黑色背景去掉 在display 下的color plat....下面调整。把上面得黑条,变成白得,点击apply就行了 方法二: 1. 用文字编辑器开启c:\MSC\patran2003r2\shareware\msc\unsupported \utilities\extra_files\bv_p3toolbar_ntgui.def 2. 选取所有的文字并复制。 3. 用文字编辑器开启 c:\MSC\patran2003r2\p3toolbar.def ,在最後面的位置贴上刚刚复制的文字, 最後存档离开 4. 复制c:\MSC\patran2003r2\shareware\msc\unsupported\utilities\ icons\*.bmp 到c:\MSC\patran2003r2\icons工具栏出现三个图标,背景颜色轻松改变。 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! 问题出在translation parameter里面的一个参数numerical,帮助文件里面说它用于比较两个数是否相等,其默认为1e-4.writing才是判断一个数是否近似为零,默认为1e-21.但实际上当你给出的强迫位移量小于numerical时,它就认为近似为0,在bdf文件中就给忽略掉了.你修改numerical为1e-5,你上面说的1e-5就可以算了。 5、自重怎么加到模型上去? 自重是在load/BCS里加的create->inertial load->element 在input data->load/BC set scale Factor [输入加速度的值一般取9.8] Trans Accel(a1 a2 a3)<0 -1 0> (力是沿Y轴向下 )后就ok了 tools下面有个mass properties是计算模型质量和惯量的,不知对你有没有帮助 /solu /output,mass,txt psolve,elform /output finish 将在当前目录下生成mass.txt文件里面会么都有。 6、共振时的应力 模态分析得到的位移是模态位移,没有输入载荷时的结果其值没有意义,同一模态下模态坐标位移之间的比值才有意义,因此,通过NORMALMODES分析的结果不表示真实的位移以及应力等。 需要计算FREQENCY RESPONSE或者 TRANSIENT RESPONSE才能获得真实应力。 7、体单元的节点没有转动自由度,只需考虑三个平移自由度。 8、我在用PATRAN做计算时,分析出的结果文件很大,往往提示磁盘空间不足,结果提示需要980M空间。 你试试不用 patran直接递交分析,进入nastran 用手工递交,可输入scr=yes 之后进行分析。因为直接递交会产生个dball文件 很大的 ,没什么大用,用手递交后,用scr=yes 就可以不用产生这个文件了, 另外所产生的dball文件可以delete它,它占的空间太大了 在patran下不进行分析,不用full run 用analysic deck 之后运行 nastran.ex桌面没有 在msc 文件夹下找到,用它开打你刚才生成的文件就行 在最底一行输入scr=yes 就行,之后和在patran下直接调用nastran的界面一样。 patran向nastran递交运算时的参数可以在nastxxxx.rcf文件中设定,这个文件可以用文本编辑软件修改,加一句scr=yes就行了,运算式通用的参数都可以在这里面定制。 9、做模态分析需要密度。 10、如何画椭圆 1) Create a circle, then use Transform/Curve/Scale to stretch it in one direction. The easiest way to create an ellipse is with the Utilities menu of Patran. Utilities--Geometry--Create curve by function Exemple : If you want to create an ellipse with a=110 and b=90 your equation will be: Y = SQRT(90**2*(1-(`X**2/110**2))) 2) 画圆,将圆旋转一个角度 具体旋转角度:设a为长轴长; b为短轴长 则满足关系 cos(angle)=(b/a) 旋转该角度后,将其投影到园形原来的所在平面即得所需椭圆 3) 利用局部坐标系 创建局部坐标系,选cylindrical 1,设置长短轴即可 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里的Pinned DOFs @ Node 1/2定义铰接。 4、Patran 如何把不小心Equivalence的node分开 用 Utility/FEM-Elements/Separate Elements 在equivalence时, 可以?选项切换为"List", 只选择特定某些节点作equivalence, 可以避免不小心的失锗。 5、Patran如何?Tri3单元转换为Tri6单元 在Patran -> Element -> Modify/Element/Edit , ?Type选项打勾, 在Shape中选Tri, New Shape 选Tri6, 最后选取想要改变的Tri3单元, 6、Patran 如何定义材料库 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 中出现刚才的材料名称。. 7、Patran 的完整信息输出: 在执行Patran的时候出现齐怪的错误信息时, 可以先把Patran关闭, 接着启动DOS窗口, 在DOS下直接输入patran -stdout ,执行 再重新启动Patran, 会多了一个信息窗口。 16、超单元 超单元的应用跟模态分析没有关系,它只是有限元中的一种分支方法,用于解决运算量过大的问题的,现在基本上已经不用了! 17、约束 我Patran/Nastran在做项目时,对结构应力进行分析,发现在结构被约束的部位往往产生比较大的应力,甚至是最大的应力值就在这些部位出现,好像这种结果并不太符合实际情况。请问约束该如何加,才能消除这种影响? 这种现象很正常的,只要你的约束反映了实际情况,约束部位的过大应力你可根据圣维南原理视而不见。 圣维南原理看看力学方面的书!如果一定要把约束区域分析准确的话,你的约束一定要尽可能地与真实情况相符,航空经常用的方法是加弹簧元,弹簧元的刚度是用的经验值或者试验值! 18、MSC.Flds 2004 计算飞行载荷-flightload and dynamics,主要用来计算气动弹性问题。 19、如何找到频谱分析的响应最大点? 在Tools--List---Create就可以将你想要的点的结果导入到一个组中,然后可以从结果文件中调出再在origin中编辑,或是直接显示。 20、自由模态 所谓自由模态,就是无约束模态 前6阶都是刚体模态 21、能否设置计算结束后不输出f06文件? 可以设置的,Solution Parameters里面Maximum Printed Lines = 就是f06输出的最大行数 22、一个利用DMAP提取质量、刚度矩阵的方法 对于一个结构静力的线性分析或非线性分析,总刚度矩阵是每次求解过程中所必须的。如果想再一次分析完成之后在结果中查看结构的总刚度矩阵和总质量矩阵,可以在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 23、MSC.Nastran可形成如下输出文件 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 即可 25、弹簧 弹簧的位置由两个节点定义,弹簧是只能承受指定自由度方向刚度的单元,所以在定义弹簧的时候,还必须指定弹簧的自由度方向(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。 模拟一个y方向可被压缩弹簧:在用SPRING属性的时候两个自由都选则UY,然后在边界条件里将你需要的节点全约束就可以了。STRING 是NASTRAN的使用的一个字符串。在这里指SPRING的方向。 零维弹簧指的是GROUND SPRING即地弹簧比如在结构的某个部位需要施加一个对地的弹性约束,就可以直接使用这种单元属性。当然用一维的同样可以模拟,用零维的比较方便而已。 定义弹簧元:property的1D单元中有一个项是弹簧,只要再输入这个弹簧单元的方向和刚度值即可。注意弹簧单元的性质只取决于作用方向和刚度系数这两个参数,而于单元的两个节点位置无关。 |
|
来自: 阿钢avckkrh3m1 > 《振动试验学习笔记》