第2讲 有限元软件功能概览.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,*,第二讲,有限元（,FEA,）软件介绍,1,、软件功能,2,、前处理、求解器和后处理,3,、有限元建模,4,、计算结果的评价,MSC.PATRAN&.NASTRAN,大型通用,FEA,分析系统,MSC.PATRAN&.NASTRAN,是美国,MSC,公司生产的大型通用,FEA,分析系统,为工业标准的,FEA,原代码程序及国际招标中工程分析和校验的工具,是船舶结构分析的主流软件,1,、有限元软件功能,结构分析,热分析,电磁分析,流体分析,(,CFD),耦合场分析,-,多物理场,4/8/2026,3,结构分析,概览,结构分析的类型,:,静力分析,-,用于静态载荷,.,可以考虑结构的线性及非线性行为，例如,:,大变形、大应变、应力刚化、接触、塑性、超弹及蠕变等,.,模态分析,-,计算线性结构的自振频率及振形,.,谱分析,是模态分析的扩展，用于计算由于随机振动引起的结构应力和应变,(,也叫作,响应谱或,PSD,).,结构分析用于确定结构的变形、应变、应力及反作用力等,.,4/8/2026,4,结构分析,概览,(,续,),谐响应分析,-,确定线性结构对随时间按正弦曲线变化的载荷的响应,.,瞬态动力学分析,-,确定结构对随时间任意变化的载荷的响应,.,可以考虑与静力分析相同的结构非线性行为,.,特征屈曲分析,-,用于计算线性屈曲载荷并确定屈曲模态形状,.(,结合瞬态动力学分析可以实现非线性屈曲分析,.),专项分析,:,断裂分析,复合材料分析，疲劳分析,Courtesy:Sikorsky Aircraft,4/8/2026,5,结构分析,概览,(,续,),它的显式方程求解冲击、碰撞、快速成型等问题，是目前求解这类问题最有效的方法,.,用于模拟非常大的变形，惯性力占支配地位，并考虑所有的非线性行为,.,4/8/2026,6,热分析概览,热分析之后往往进行结构分析,计算由于热膨胀或收缩不均匀引起的应力,.,功能,:,相变,(,熔化及凝固,),内热源,(,例如电阻发热等,),三种热传递方式,(,热传导、热对流、热辐射,),热分析计算物体的稳态或瞬态温度分布,，,以及热量的获取或损失、热梯度、热通量等,.,4/8/2026,7,流体分析,概览,流体分析,用于确定流体的流动及热行为,.,流体分析分以下几类,:,CFD,提供强大的计算流体动力学分析功能，包括不可压缩或可压缩流体、层流及湍流，以及多组份流等,.,声学分析,-,考虑流体介质与周围固体的相互作用,进行声波传递或水下结构的动力学分析等,.,容器内流体,分析,-,考虑容器内的非流动流体的影响,.,可以确定由于晃动引起的静水压力,.,流体动力学耦合,分析,-在考虑流体约束质量的动力响应基础上，,在结构动力学分析中使用流体耦合单元,.,4/8/2026,8,耦合场分析概览,耦合场分析,考虑两个或多个物理场之间的相互作用。如果两个物理场之间相互影响，单独求解一个物理场是不可能得到正确结果的，因此你需要一个能够将两个物理场组合到一起求解的分析软件。,耦合场分析的典型情况有：,热,应力分析,流体,结构相互作用,两根热膨胀系数不同的棒焊接在一起，图示为加热后的变形,.,4/8/2026,9,2,、前处理、求解器和后处理,前处理,有限元的建模工具,功能强大的有限元网格生成功能,自动与半自动网格剖分,网格优化功能,网格加密功能,属性定义,加载,约束处理,载荷分类,自由度,DOF,-,定义节点的自由度（,DOF,）值,(,结构分析,_,位移热分析,_,温度、电磁分析,_,磁势等,),集中载荷,-,点载荷,(,结构分析,_,力、热分析,_,热导率、电磁分析,_ magnetic current segments,),面载荷,-,作用在表面的分布载荷,(,结构分析,_,压力、热分析,_,热对流、电磁分析,_magnetic Maxwell surfaces,等,),体积载荷,-,作用在体积或场域内,(,热分析,_,体积膨胀、内生成热、电磁分析,_,magnetic current density,等),惯性载荷,-,结构质量或惯性引起的载荷,(,重力、角速度等,),11,加载,可在实体模型或,FEA,模型,(,节点和单元,),上加载,.,在关键点处约束,实体模型,沿线均布的压力,在关键点加集中力,在节点处约束,FEA,模型,沿单元边界均布的压力,在节点加集中力,12,求解器,有限元的核心,求解结果保存在数据库中并输出到结果文件,结果文件,结果数据,数据库,求解器,结果,输入数据,13,后处理,结果的绘图和列表,静力分析结果后处理的步骤主要包括,:,1.,绘变形图,2.,变形动画,3.,支反力列表,4.,应力等值线图,5.,网格密度检查,Guidelines,后处理,绘变形图,绘出结构在静力作用下的变形结果,:,后处理,变形动画,以动画方式模拟结构在静力作用下的变形过程,:,Utility Menu:PlotCtrls Animate Deformed Shape.,后处理,应力等值线,应力等值线方法可清晰描述一种结果在整个模型中的变化，可以快速确定模型中的“危险区域”。,显示应力等值线,:,Main Menu:General Postprocessor Plot Results -Contour Plot-Nodal Solution.,后处理,应力等值线动画,结果动画,:,Utility Menu:PlotCtrls Animate Deformed Results,后处理,应力平均,FEA,的计算结果包括通过计算直接得到的初始量和导出量。,任一节点处的,DOF,结果,(UX,、,UY,、,TEMP,等,),是初始量。它们只是在每个节点计算出来的初始值。其它量，如应力应变，是由,DOF,结果通过单元计算导出而得到的,。,因此，在给定节点处，可能存在不同的应力值。这是由以与此节点相连的不同单元计算而产生的。,“,节点结果”（,nodal solution,）,画出,的是在节点处导出量的平均值，而,“,单元结果,”,（,element solution,）,画出非平均量。,在同一个节点处，相临单元计算出的结果不同。,后处理,应力平均,(,续,),在弹性模量不同的材料交界处，应力分量会不连续。,(,PowerGraphics,自动考虑到这一点并对此界面不进行平均处理。,),在多数情况下，画出平均应力图，但有时要画出：,非平均单元应力显示不连续的应力,平均的节点应力显示连续的应力,20,后处理,应力平均,(,续,),在不同厚度的壳单元的交界处，大多数应力会不连续。,(,PowerGraphics,自动考虑到这一点并对此界面不进行平均处理。,),非平均单元应力显示不连续的应力,平均的节点应力显示掩盖了应力的不连续性,21,后处理,应力平均,(,续,),在壳单元构成的尖角或连接处，某些应力分量不连续,。,非平均单元应力显示不连续的应力,平均的节点应力显示掩盖了应力的不连续性,22,3,、有限元建模,制订分析方案是很重要的。一般考虑下列问题：,分析领域,分析目标,线性,/,非线性问题,静力,/,动力问题,分析细节的考虑,几何模型对称性,奇异,单元类型,网格密度,单位制,材料特性,载荷,求解器,制订的分析方案好坏直接影响分析的,精度,和,成本,（人耗工时，计算机资源等），但通常情况下精度和成本是相互,冲突，,特别是分析较大规模和具有切割边界的模型时更为明显。一个糟糕的分析方案可能导致分析资源紧张和分析方式受得限制。,June 3,1996,23,确定合适的分析学科领域,实体运动，承受压力，或实体间存在接触,施加热、高温或存在温度变化,恒定的磁场或磁场,电流（直流或交流）,气（液）体的运动，或受限制的气体,/,液体,以上各种情况的耦合,结构,热,磁,流体,电,耦合场,准则,June 3,1996,24,有限元建模,-,分析目的,分析目的直接决定分析近似模型的确定。分析目的，就是这样一个问题的答案：“,利用,FEA,我想研究结构哪些方面的情况？,”,结构分析,:,要想得到极高精度的应力结果，必须保证影响精度的任何结构部位有理想 的单元网格，不对几何形状进行细节上的简化。应力收敛应当得到保证，而任何位置所作的任何简化都可能引起明显误差。,在忽略细节的情况下，用相对较粗糙的单元网格计算转角和法向应力。,复杂的模型要求具有较好的均匀单元网格，并允许忽略细节因素。,准则,June 3,1996,25,有限元建模,-,分析目的,(,续,),模态分析,:,简单模态振型和频率可以忽略细节因素而使用相对较粗糙的单元网格进行分析计算。,热分析：,温度分布梯度变化不大时可以忽略细节，划分均匀且相对稀疏的单元网格。,当温度场梯度较大时，在梯度较大的方向划分细密的单元网格。梯度越大，单元划分就越细密。,利用一个能同时模拟两个物理场的模型求解温度和热耗散应力，但热和应力模型都是相对独立的。,June 3,1996,26,在建立分析模型之前必须制订好建模方案：,必须考虑那些细节问题？,对称,/,反对称,/,轴对称？,模型中存在应力奇异？,选用那种类型的单元？,线单元,壳单元,X,Y,平面单元,平面应力或应变单元,轴对称单元,谐单元,实体单元,专用单元,线性单元,/,高阶单元,/P,单元,四边形单元,/,三角形单元，块单元,/,四面体单元,June 3,1996,27,有限元建模,-,细节处理,对于分析不重要的细节不应当包含在分析模型中。当从,CAD,系统传一个模型到,ANSYS,程序中时往往可以作大量的简化处理。,然而，诸如倒角或孔等细节可以是最大应力出现的位置，这些细节对于你的分析目的是十分重要的。,带倒角,不带倒角,准则,June 3,1996,28,有限元建模,-,对称性模型,对称,当物理系统的形状、材料和载荷具有对称性时，就可以只对实际结构中具有,代表性,的部分或截面进行建模分析，再将结果映射到整个模型上，就能获得,相同精度的结果,。,物理系统对称分析要求具有以下对称性条件几何结构对称,材料特性对称,具有零位移约束,存在非零位移约束,定义,June 3,1996,29,有限元建模,-,对称性模型,(,续,),对称类型,轴对称,即绕某一轴线存在对称性，这类结构如：电灯泡，直管，圆锥体，圆盘和圆屋顶。对称面就是旋转形成结构的横截面，它可以在任何位置。大多数轴对称分析求解必须假定非零约束（边界），集中力、压力和体截荷均具有轴对称。然而，如果截荷不存在轴对称性，并且是线性分析，可以将截荷分成简谐成分，进行独立求解（然后进行叠加）,。,定义,June 3,1996,30,有限元建模,-,对称性模型,(,续,),对称类型,(,续,),旋转对称即结构由绕轴分布的几个重复部分组成，诸如涡轮叶片这类物体。大多数旋转对称分析求解要求非零位移约束（边界），集中力、压力和体载荷应具有对称性。然而，如果载荷不对称分布，并且如果是线性分析，它们可以利用周期对称求解。,定义,June 3,1996,31,有限元建模,-,对称性模型,(,续,),对称类型,(,续,),平面,或,镜面对称,即结构的一半与另一半成镜面映射关系，对称位置（镜面）称为对称平面。大多数平面对称分析求解要求非零位移约束（边界），集中力、压力和体力应当对称。但是，如果这些载荷不对称，并且是线性分析，它们可以分成,对称,或,反对称,问题进行独立求解,。,定义,该图显示了镜面对称和旋转对称,June 3,199632,有限元建模,-,对称性模型,(,续,),对称类型,(,续,),重复,或,平移对称,即结构是由沿一直线分布的重复部分组成，诸如带有均匀分布冷却节的长管等结构。该对称要求非零位移约束，集中力、压力和体载荷应具有对称性。,一个结构可能由多个对称平面，这样就可以利用对称性建立一个很小的等效分析模型。,定义,图示模型具有镜面对称,(2X),和,重复对称,June 3,1996,33,有限元建模,-,对称性模型,(,续,),在实际当中，可以利用对称模型进行分析能获得,更好的分析结果,，因为可以建立更精确、综合考虑各细节的模型。,准则,June 3,1996,34,求解时模型是否准备就绪,?,在求解初始化前，应进行分析数据检查，包括下面内容,:,统一的单位,单元类型和选项,材料性质参数,考虑惯性时应输入材料密度,热应力分析时应输入材料的热膨胀系数,实常数,(,单元特性,),单元实常数和材料类型的设置,实体模型的质量特性,模型中不应存在的缝隙,壳单元的法向,节点坐标系,集中、体积载荷,面力方向,温度场的分布和范围,热膨胀分析的参考温度,(,与,ALPX,材料特性协调,?),35,进行求解,(,续,),没有获得结果的原因是什么,?,往往是求解输入的模型不完整或存在错误，典型原因有,:,约束不够,!(,通常出现的问题,),。,材料性质参数有负值,如密度或瞬态热分析时的比热值。,未约束铰接结构,如两个水平运动的梁单元在竖直方向没有约束。,屈曲,-,当应力刚化效应为负（压）时，在载荷作用下整个结构刚度弱化。如果刚度减小到零或更小时，求解存在奇异性，因为整个结构已发生屈曲。,4,、分析结果评价,验证分析的结果，在任何有限元分析中无疑是最为重要的步骤,.,在开始任何分析以前，应该至少对分析的结果有粗略的估计（来自经验、试验、标准考题等）,.,如果结果与预期的不一样，应该研究差别的原因,.,37,识别无效的结果,应该知道所分析的对象的一些基本的行为,:,重力方向总是竖直向下的,离心力总是沿径向向外的,物体受热一般要膨胀,没有一种材料能抵抗,1,000,MPa,的应力,轴对称的物体几乎没有为零的 环向应力,弯曲载荷造成的应力使一侧受压，另一侧受拉,如果只有一个载荷施加在结构上，检验结果比较容易,.,如果有多个载荷，可单独施加一个或几个载荷分别检验，然后施加所有载荷检验分析结果,.,Objective,识别分析中无效的结果,.,38,变形,/,温度,/,应力,检验求解的自由度及应力,.,确认施加在模型上的载荷环境是合理的,.,确认模型的运动行为与预期的相符,-,无刚体平动、无刚体转动、无裂缝等,.,确认位移和应力的分布与期望的相符，或者利用物理学或数学可以解释,.,Objective,检验求解的变形、温度及应力,.,39,反作用力或节点力,模型所有的反作用力应该与施加的点力、压力和惯性力平衡,.,在所有约束节点的竖直方向的反作用力,.,在所有约束节点水平方向的反作用力必须与水平方向的载荷平衡,.,所有约束节点的反作用力矩必须与施加的载荷平衡,.,注意包含在约束方程中自由度的反力，不包括由这个约束方程传递的力,.,必须与施加的竖直方向的载荷平衡,Objective,检验整个或部分模型的反作用力或节点力,.,40,检查网格精度,由于网格密度影响分析结果的精度，因此有必要验证网格的精度是否足够。,有三种方法进行网格精度检查：,1,.,观察（,Visual inspection,）,2.,误差估计,3,.,将网格加密一倍，重新求解并比较两者结果。注意,:,有些情况下这种做法不适用。,误差估计,误差估计是依据沿单元内边界的应力或热流的不连续性，是平均与未平均节点应力间的差值,.,Elem 1,Elem 2,s,avg,=1100,s,=1200,s,=1000,s,avg,=1200,s,=1300,s,=1100,(,节点的,s,s,是积分点的外插）,Objective,定义误差估计的基础,.,42,验证足够的网格密度,(,续,),能量百分比误差,能量百分比误差是对所选择的单元的位移、应力、温度或热流密度的粗略估计,.,它可以用于比较承受相似载荷的相似结构的相似模型,.,这个值的通常应该在,10,%,以下,.,如果不选择其他单元，而只选择在节点上施加点载荷或应力集中处的单元，误差值有时会达到,50%,或以上,.,在绘制变形图时，百分比误差在图形右边的文本中以“,SEPC”,表示,.,也可以在通用后处理的,List Results,PerCent,Error,菜单中列出,.,SEPC 2%,43,验证足够的网格密度,(,续,),应力偏差,要检验某个位置的网格离散应力误差，可以列出或绘制应力偏差,.,某一个单元的应力偏差是此单元上全部节点的六个应力分量值与此节点的平均应力值之差的最大值,.,所关心位置上的应力,偏差,值,450,psi,(30,000,psi,应力的,1.5%),44,验证足够的网格密度,(,续,),能量误差,每个单元的另一种误差值是能量误差,.,它与单元上节点应力差值有关的,用于计算选择的单元的能量百分比误差,.,能量误差的单位是能量的单位,.,能量误差在通用后处理的菜单中得到,.,45,验证足够的网格密度,(,续,),应力上下限,应力上下限可以帮助确定由于网格离散误差对模型的应力最大值的影响,.,显示或列出的应力上下限包括,:,估计的上限,-,SMXB,估计的下限,-,SMNB,应力上下限限并不是估计实际的最高或最小应力,.,对于有些情况，,SMXB,过于保守,.,而有些情况比实际的要小,.,应力上下限定义了一个”确信范围,.”,如果没有其他的确凿的验证，就不能认为实际的最大应力低于,SMXB.,X stress SMAX 32,750 psi,SMXB 33,200 psi,(difference 450 psi 1.5%),46,使用,p,方法验证足够的网格精度,如果使用,p-,方法,进行结构分析，可以依靠,p,单元自动调整单元多项式阶数，达到收敛到设定的精度,.,对这种方法的相信程度，与使用经验有关,.,规定,0.1%,局部应力误差，使用,p,方法计算的最大,X,方向应力约为,34,700,psi,(,比普通,h,方法高出大约,5%,),47,第三讲,船舶结构的有限元分析,1,、克令吊机座结构分析,2,、锚机机座结构分析,3,、舱长超过,30m,的船舶结构分析,4,、船长大于,150m,的船舶结构分析,