ANSYS载荷施加及后处理

ANSYS稳态热分析的基本过程ANSYS热分析可分为三个步骤：•前处理：建模、材料和网格•分析求解：施加载荷计算•后处理：査看结果1、建模①、确定jobname> title、unit；②、进入PREP7前处理，定义单元类型，设定单元选项；③、定义单元实常数；④、定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定的，也可以随温度变化；⑤、创建儿何模型并划分网格，请参阅《ANSYS Modeling and Meshing Guide》。

2、施加载荷计算①、定义分析类型•如果进行新的热分析：Command: ANTYPE, STATIC, NEWGUI: Main menu>Solution>-Analysis Type->New Analysis>Steady-state•如果继续上一次分析，比如增加边界条件等：Command: ANTYPE, STATIC, RESTGUI: Main menu>Solution>Analysis Type->Restart②、施加载荷可以直接在实体模型或单元模型上施加五种载荷（边界条件）：a、恒定的温度通常作为自由度约束施加于温度已知的边界上。

Command Family: DGUI： Main Menu>Solution>-Loads-Apply>-Thermal-Temperatureb、热流率热流率作为节点集中载荷，主要用于线单元模型中（通常线单元模型不能施加对流或热流密度载荷），如果输入的值为正，代表热流流入节点，即单元获取热量。

如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。

注意：如果在实体单元的某一节点上施加热流率，则此节点周用的单元要密一些，在两种导热系数差别很大的两个单元的公共节点上施加热流率时, 尤其要注意:。

此外，尽可能使用热生成或热流密度边界条件，这样结果会更精确些。

ANSYS稳态热分析的基本过程ANSYS热分析可分为三个步骤：•前处理：建模、材料和网格•分析求解：施加载荷计算•后处理：查看结果1、建模①、确定jobname、title、unit；②、进入PREP7前处理，定义单元类型，设定单元选项；③、定义单元实常数；④、定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定的，也可以随温度变化；⑤、创建几何模型并划分网格，请参阅《ANSYS Modeling and Meshing Guide》。

2、施加载荷计算①、定义分析类型●如果进行新的热分析：Command: ANTYPE, STATIC, NEWGUI: Main menu>Solution>-Analysis Type->New Analysis>Steady-state●如果继续上一次分析，比如增加边界条件等：Command: ANTYPE, STATIC, RESTGUI: Main menu>Solution>Analysis Type->Restart②、施加载荷可以直接在实体模型或单元模型上施加五种载荷(边界条件) ：a、恒定的温度通常作为自由度约束施加于温度已知的边界上。

Command Family: DGUI：Main Menu>Solution>-Loads-Apply>-Thermal-Temperatureb、热流率热流率作为节点集中载荷，主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷)，如果输入的值为正，代表热流流入节点，即单元获取热量。

如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。

注意：如果在实体单元的某一节点上施加热流率，则此节点周围的单元要密一些，在两种导热系数差别很大的两个单元的公共节点上施加热流率时，尤其要注意。

此外，尽可能使用热生成或热流密度边界条件，这样结果会更精确些。

ANSYS后处理1.ANSYS后处理时如何按灰度输出云图？1）你可以到utilitymenu-plotctrls-style-colors-window colors试试2）直接utilitymenu-plotctrls－redirect plots2 将云图输出为JPG菜单->PlotCtrls->Redirect Plots->To JPEG Files3.怎么在计算结果实体云图中切面?命令流/cplane/type图形界面操作<1.设置工作面为切面<2.PlotCtrls-->Style-->Hidden line Options将[/TYPE]选项选为section将[/CPLANE]选项选为working plane4.非线性计算过程中收敛曲线实时显示solution>load step opts>output ctrls>grph solu track>on5.运用命令流进行计算时,一个良好的习惯是:使用SELECT COMMEND后.........其后再加上ALLSEL.........6.应力图中左侧的文字中，SMX与SMN分别代表最大值和最小值如你plnsolv,s,eqv则 SMX与SMN分别代表最大值等效应力和最小值等效应力如你要看的是plnsolv,u则SMX与SMN分别代表位移最大值和位移最小值不要被S迷惑mx(max)mn(min)7.在非线性分析中，如何根据ansys的跟踪显示来判断收敛？在ansys output windows 有 force convergence value值和 criterion 值当前者小于后者时，就完成一次收敛你自己可以查看两条线的意思分别是:FL2：不平衡力的2范数 FCRIT：不平衡力的收敛容差，如果前者大于后者说明没有收敛，要继续计算，当然如果你以弯矩M为收敛准则那么就对应 M L2 和 M CRIT希望你现在能明白8.两个单元建成公共节点，就成了刚性连接，不是接触问题了。

在任意面施加任意方向任意变化的压力在某些特殊的应用场合，可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷，当然，这在一定程度上可以通过 ANSYS表面效应单元实现。

如果利用 ANSYS的参数化设计语言，也可以非常完美地实现此功能，下面通过一个小例子描述此方法。

!!!在执行如下加载命令之前 ,请务必用选择命令 asel 将需要加载的几何面选择出来!!!finish/prep7et,500,shell63press=100e6amesh, allesla, snsla,s,1!如果载荷的反向是一个特殊坐标系的方向 ,可在此建立局部坐标系 ,并将 !所有节点坐标系旋转到局部坐标系下 .*get,enmax,elem,,num,maxdofsel,s,fx,fy,fzfcum,add !!! 将力的施加方式设置为 " 累加 ",而不是缺省的 " 替代 "*do,i,1,enmax*if,esel,eq,1,then*get,ae,elem,i,area! 此命令用单元真实面积，如用投影面积，请用下几条命令!*get,ae,elem,i,aproj,x ! 此命令用单元 X 投影面积，如用真实面积，请用上一条命令!*get,ae,elem,i,aproj,y ! 此命令用单元 Y投影面积!*get,ae,elem,i,aproj,z ! 此命令用单元 Z 投影面积xe=centrx !单元中心 X 坐标 (用于求解压力值 )ye=centry !单元中心 Y 坐标 (用于求解压力值 )ze=centrz !单元中心 Z 坐标 (用于求解压力值 )!下面输入压力随坐标变化的公式 ,本例的压力随 X 和 Y 坐标线性变化.p_e=(xe-10)*press+(ye-5)*pressf_tot=p_e*aeesel,s,elem,,insle,s,corner*get,nn,node,,countf_n=f_tot/nn*do,j,1,nnf,nelem(i,j),fx,f_n ! 压力的作用方向为X 方向!f,nelem(i,j),fy,f_n ! 压力的作用方向为 Y 方向!f,nelem(i,j),fz,f_n ! 压力的作用方向为 Z 方向*enddo*endifesla,s*enddoaclear,allfcum,repl !!! 将力的施加方式还原为缺省的"替代 "dofsel,allallsel说明：本信息在任意面施加任意方向任意变化的压力在某些特殊的应用场合，可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷，当然，这在一定程度上可以通过 ANSYS表面效应单元实现。

ANSYS结构分析基础篇一、总体介绍进行有限元分析的基本流程：1.分析前的思考1)采用哪种分析静态,模态,动态...2)模型是零件还是装配件零件可以form a part形成装配件,有时为了划分六面体网格采用零件,但零件间需定义bond接触3)单元类型选择线单元,面单元还是实体单元4)是否可以简化模型如镜像对称,轴对称2.预处理1)建立模型2)定义材料3)划分网格4)施加载荷及边界条件3.求解4.后处理1)查看结果位移,应力,应变,支反力2)根据标准规范评估结构的可靠性3)优化结构设计高阶篇：一、结构的离散化将结构或弹性体人为地划分成由有限个单元,并通过有限个节点相互连接的离散系统;这一步要解决以下几个方面的问题:1、选择一个适当的参考系,既要考虑到工程设计习惯,又要照顾到建立模型的方便;2、根据结构的特点,选择不同类型的单元;对复合结构可能同时用到多种类型的单元,此时还需要考虑不同类型单元的连接处理等问题;3、根据计算分析的精度、周期及费用等方面的要求,合理确定单元的尺寸和阶次;4、根据工程需要,确定分析类型和计算工况;要考虑参数区间及确定最危险工况等问题;5、根据结构的实际支撑情况及受载状态,确定各工况的边界约束和有效计算载荷;二、选择位移插值函数1、位移插值函数的要求在有限元法中通常选择多项式函数作为单元位移插值函数,并利用节点处的位移连续性条件,将位移插值函数整理成以下形函数矩阵与单元节点位移向量的乘积形式;位移插值函数需要满足相容协调条件,采用多项式形式的位移插值函数,这一条件始终可以满足;但近年来有人提出了一些新的位移插值函数,如：三角函数、样条函数及双曲函数等,此时需要检查是否满足相容条件;2、位移插值函数的收敛性完备性要求：1 位移插值函数必须包含常应变状态;2位移插值函数必须包含刚体位移;3、复杂单元形函数的构造对于高阶复杂单元,利用节点处的位移连续性条件求解形函数,实际上是不可行的;因此在实际应用中更多的情况下是利用形函数的性质来构造形函数;形函数的性质：1相关节点处的值为 1,不相关节点处的值为 0;2形函数之和恒等于 1;这里我们称为的相关节点, 为的相关节点,其它点均为不相关节点;三、单元分析目的：计算单元弹性应变能和外力虚功;使用最小势能原理,需要计算结构势能,由弹性应变能和外力虚功两部分构成;结构已经被离散,弹性应变能可以由单元弹性应变能叠加得到,外力虚功中的体力、面力都是分布在单元上的,也可以采用叠加计算;2、计算单元外力功从前面推导可以看出：单元弹性应变能可计算的部分只有单元刚度矩阵,单元外力虚功可计算的部分只有单元等效体力载荷向量和等效面力载荷向量;在实际分析时并不需要进行上述推导,只需要将假定的位移插值函数代入本节推导得出的单元刚度矩阵、等效体力载荷向量和等效面力载荷向量的计算公式即可;所以我们说有限元分析的第三步是计算单元刚度矩阵、等效体力载荷向量和等效面力载荷向量;几点说明：1单元刚度矩阵具有正定性、奇异性和对称性三各重要特性;所谓正定性指所有对角线元素都是正数,其物理意义是位移方向与载荷方向一致；奇异性是说单元刚度矩阵不满秩是奇异矩阵,其物理意义是单元含有刚体位移；对称性是说单元刚度矩阵是对称矩阵,程序设计时可以充分利用;2按照本节公式计算的单元等效体力载荷向量和等效面力载荷向量称为一致载荷向量;实际分析时有时也采用静力学原理计算单元等效体力载荷向量和等效面力载荷向量,实际应用表明在大多数情况下,这样做可以简化计算,同时又基本上不影响分析结果;二、预处理总述1、实体分析可是3D或2D,3D分析采用的高阶单元SOLID186或SOLID187划分的四面体TET 或六面体HEX单元,2D分析采用的高阶单元的三角形TRI或四边形QUA单元,2D分析时需要在创建项目时在GEOMETRY的分析类型项选择2D,实体分析得每个节点结构上只有3个自由度,如下图所示2、面体分析几何上是2D但离散元是3D,通常面体厚度给予赋值;面体网格划分采用壳单元,具有6个自由度;3、线体几何上是1D,离散元是3D,截面形状可通过line body进行设置,线体网格划分采用梁单元,具有6个自由度;4、同个part下的所有body共享相交边界,网格划分时共用交界上的节点,不需要设置接触;5、NameSelection的使用技巧,在model模块下,可点击右键insert NameSelection,一般Nameselection的选择方法可用几何选取,直接在模型上鼠标点选;另一种实用的选取方法为Worksheet,可以添加多种条件进行筛选,模型划分网格后,可以精确到对每个单元的选取;三、网格划分1、relevance选项控制网格的精度,值在-100到100间,越小越粗糙,越大越精密;relevance center 控制relevance中间点的精度,element size控制整个模型的最大单元尺寸;2、网格的高级尺寸控制a)接近度和曲度结合控制b)曲度c)邻近度d)固定尺寸曲度对于一些含曲线特征的几何体,可以控制其划分网格的精密度邻近度可以控制某个区域两个邻近的几何特征间的网格划分密度2、网格的高级选项形状检测：标准力学-线性分析、模态和热分析进阶力学-大变形分析、材料的非线性分析3、局部网格划分控制Method 选择Automatic 首先若能SWEEP则选用sweep划分HEX网格,否则选用patchconforming划分TET网格;四面体TET网格划分有两种方法：patchconforming和patch independent;对于不能通过sweep得到六面体的几何体可以选用Hex dominant或者Multizone划分方法4、尺寸控制Sizing可以通过element size单元最大尺寸、Number of divisions每个边的单元数量、Sphere of influence控制影响区,可设置影响半径来调节网格划分尺寸;Contact sizing可设置接触面的尺寸;5、其他设置element refinement可设置选择几何体的网格密度加密倍数；mapped face meshing 可设置映射面生成结构化网格；可通过side、corner、end点的定义来设置映射策略;inflation control设置膨胀层,主要用于流体分析的边界层划分;pinch 可以移出一些不必要的小的几何特征,划分网格时可以去掉一些小的凸起部分;划分网格前有个小圆台采用pinch划分网格后没有凸台Master选择蓝色线,Slave选择红色线,tolerance的值要比凸台的高度大;6、虚拟拓扑的应用虚拟拓扑有助于优化几何模型,可以合并面,分割面或边来提高网格划分质量;虚拟拓扑可以自动控制虚拟拓扑合并面虚拟拓扑分割边虚拟拓扑设置：behavior可以设置拓扑搜寻深度;7、子模型的应用当原几何模型较大,网格数量有限,为了对模型局部进行更精确的计算分析,可以采用子模型;子模型的一般创建方法：先对整体模型项目A进行分析计算,然后copy原项目得到项目B,对项目B中几何进行切割细化网格,将项目A的solution栏拖到Setup栏,最好在B项目求解设置下的submodeling 插入边界条件,子模型的切割边界应远离高应力区;四、静力学分析线弹性静力分析假设：a)各向同性线弹性材料b)小变形理论c)无时间、无阻尼效应1、point mass,质量点可以通过坐标或选择几何面、线、点加载在几何体上,质量点只受Acceleration,Standard earth gravity,Rotational velocity影响;2、求解设置可设置求解步数,定义每步的终止时间,静力分析中的time只是一个跟踪量求解器选择：自动,直接求解Direct,迭代求解Iterative弱弹簧的使用：为了满足静止约束,程序可自动添加弱弹簧,可以在结果中查看弱弹簧的反力,应该是一个很小的值,并不影响结构的应力分析;惯性释放：当物体受力不平衡产生加速度时,利用惯性释放可以产生一个惯性力进行静力分析,惯性释放只能用于线性结构分析;惯性释放下的应力：静力平衡下的应力3、施加载荷加速度、角速度、压力、力,静水压力模拟水压轴承力Bearing Load,施加在整个圆柱面上;remote force定义力的作用点螺栓预紧力Bolt Pretension施加在圆柱面上,可以定义预紧力或伸长量;Thermal condition,计算热应力,需要设置reference temperature4、施加约束Fix support 约束点、线、面的所有自由度;Displacement 位移约束Elastic Support 无摩擦的弹性支持面Frictionless Support,约束面的法向运动,作用在平面上等同于对称边界条件作用在圆柱面上约束径向运动cylindrical support 只作用在圆柱面上,可以设轴向,径向,切向三个自由度compression only support 基于罚函数方法对目标面建立一个刚性接触面simply supported 作用于点或边,面体或线体,约束所有平动除了转动自由度Fixed rotation 约束转动,放开平动nodal load and support 必须通过name selection 来选取nodetools-Solve process settings可以设置求解用的计算机CPU数五、接触基本设置接触是一种高度非线性特征,接触一般通过接触对描述,包括接触面contact和目标面target,程序一共有5种接触方式,其主要特征如下：Bonded 和 no separation 都是线性接触,bonded使两个接触面固定在一起,无间隙不能相对滑动而no separation 允许有较小的滑动,其他接触都是非线性;contact 接触行为behavior分为对称和非对称两种行为;接触面的处理interface treatment：adjust to touch程序自动取消两个接触面的间隙; add offset 可以设置偏移量,正值使两个接触面靠近可以模拟过盈配合,负值使两个接触面远离;Pinball region 可以设置判断接触区域的大小,当两个面都进入pinball region时程序则判定为发生接触;mesh connections建立网格连接connection worksheet表格查看连接信息joint 定义约束副,共有九种约束形式来约束body-body 或者body-ground;定义joint时需要定义reference和mobile regions,几何窗口左边显示的自由度,其中灰色的是被约束的,彩色的是自由的joint configure可以定义约束的初始状态Set定义初设状态,revert恢复原始状态;对于旋转面或圆柱面的约束类型,可以定义扭转刚度和扭转阻尼;大多数joints都可以通过stops来定义他的运动区域spring and beam：spring可以通过弹簧来连接body,可以定义初始值和弹簧刚度,beam可以定义材料和圆形截面半径;六、remote边界条件1、Remote boundary conditions provide a means to apply a condition whose center of action is not located where the condition is scoped , “remotely”.Remote 边界条件包括 point mass,springs,joints,remote displacement,remote force and moment loads;所有的remote边界条件都是采用MPC约束方程进行计算,几何行为可以设置为rigid,deformable and coupled,remote计算更耗时;设置remote边界一般先定义remote point,可以直接选择几何特征或给定坐标定义,也可以在定义remote边界条件时通过右键“promote remote point”定义;2、behavior controlrigid,deformable and coupled3、pinball control 可以通过pinball大小来定义约束方程的数量4、constraint equations 可以多个remote point间的相互约束关系;七、MultiStep的设置应用1、对于多步分析中的每一步,软件都作为一个独立的分析过程,载荷约束都可以单独设置;对于某些载荷或约束可以通过右键激活或抑制该步当查看计算结果选择两个载荷步之间的时间节点时,如0与1步的,则程序通过线性插值的方式得到的计算结果;2、Solution Combination结果组合Solution Combination可以通过不同的计算环境共享几何网格进行组合Solution Combination也可以通过同一计算环境的不同载荷步进行组合八、模态分析自由振动其中K-刚度矩阵和M-质量矩阵是常量,忽略阻尼C和外力F,应用线弹性材料和小变形理论,结构可以是约束的或非约束的,φ为模态坐标是个相对量;1.结构载荷和热载荷步,非线性接触不适用于模态分析,但可以施加约束或预应力;2.可以定义求解阶数和频率范围;3.由于并没有外部激励,模态变形只是一个相对量,并且是一个质量归一化的量;4.拉伸预应力将会增大自然频率,而压缩预应力将会降低自然频率;九、稳态热分析1.不考虑瞬态影响,K和Q可以是常量也可以是温度的函数,可以施加固定温度的边界条件;壳单元不考虑厚度方向的温度变化,线单元不考虑截面上的温度变化;接触中热传递：如果接触是bonded或no separation,热传递将会发生在pinball区域内的表面热接触通过以下公式进行传热：TCC默认被设为一个较大的数值用来模拟完美传热,同样可以人为设置较低的数值来模拟热阻;2.边界条件heat flow 热流量j/s,可应用于点、线、面heat flux 热通量j/m2/s,只能应用于面2D时可用于线internal heat generation 热源j/m3/s 只能用于实体perfectly insulated 绝热,默认应用于所有未设置边界条件的地方temperature 恒定温度,应用于点、线、面、实体convection 对流只能应用于面,其中h-对流传热系数,Tam-环境温度,用户可以自己设置;radiation 热辐射其中σ-玻尔兹曼常数,程序自动给定；ε-发射率,用户输入;F-form factor角系数,当correlation设为To ambient-F=1,即所有的辐射能都与周围环境进行交换当correlation设为 surface to surface ,辐射能只参与面面之间的交换,这时你可以设置Enclosure每个辐射面应该设置相同的enclosure number和Enclosure type可设为open 或perfect,如果计算报错可将其设为open;十、结果处理1.编辑legendPlane可以通过鼠标左键拖曳生产剖切面,也可以通过局部坐标系的XY平面生产剖切面 Tool 可以通过Geometry selection查看选择几何特征的计算结果,也可以先定义一个局部坐标系,再通过coordinate system查看具体某点局部坐标系的原点的计算结果;chart and Table可以对多个计算结果进行图表分析,Alert可以设置报警值,如强度极限;Geometry可以添加path和surface,path可以通过局部坐标系,边,点来定义,surface可以通过局部坐标系定义;查看edge的结果可以通过鼠标右键Convert to path result转换成基于path的计算结果,把X轴设为S即可绘制关于位置的图表;另外利用path结果可以得到应力线性化用于应力评判;error可以通过高的能量差异区来鉴别几何网格的合理性;可以通过Convergence来判断网格是否足够8.应力奇点,结构分析时由于几何模型、载荷施加等因素常常会导致应力奇点,影响计算结果的准确性,我们通过审查收敛结果来避免应力奇点;如果应力奇异区并不是我们感兴趣的区域,我们可以只对感兴趣区域的计算结果定义收敛控制,如下图所示;ANSYS结构动态分析篇一、简介动态分析包括以下模块：模态分析,谐响应分析,随机振动分析,响应谱分析及瞬态分析;动态分析中结构的惯性、阻尼都扮演着重要角色;自由振动：结构的自然频率和振型激励振动：曲柄轴和其他的旋转机械地震冲击载荷：地震工况,爆炸随机振动：火箭发射,道路交通时间载荷：汽车碰撞,汽锤、水锤等以上每种情况都可以选择相应的动态模块进行分析;1、模态分析模态分析是用来确定结构的振动特性,如自然频率和振型,通常也是进行其他动态分析的先决条件;如汽车的固有频率应发动机频率,叶片在预应力下的振动特性;2、谐响应分析谐响应分析常用来分析结构在持续的简谐载荷下的响应,如转动机械的响应;3、响应谱分析响应谱分析通常用来分析建筑结构在地震工况下的响应;4、随机振动分析宇宙空间站、航天飞机等一般都要进行随机振动分析,以便能承受一段时间内不同频率下的随机载荷;5、瞬态分析动态分析各模块的特点如下：基本方程如下：其中只有瞬态分析允许非线性,包括几何非线性、接触非线性、材料非线性;二、阻尼概述阻尼定义：阻尼是导致振动不断减弱甚至停止的一种能量耗散机制;阻尼一般与材料性质,运动速度,振动频率有关;阻尼分为以下类型：粘性阻尼-缓冲器、减震器材料/固体/滞后阻尼-内摩擦库伦或干摩擦阻尼-滑动摩擦数值阻尼-人工阻尼1、瞬态分析和阻尼模态分析中结构阻尼矩阵C的完整表达式如下：α和β阻尼用来确定瑞利阻尼对于大多数结构来说,α阻尼可以忽略,这时因此对于给定的β,低频率阻尼小,高频率阻尼大;而对于给定的α,低频率阻尼大,高频率阻尼小;α和β阻尼可以通过定义材料时输入：也可以通过全局阻尼输入：2、在谐响应分析中的材料/固体/滞后阻尼全函数的谐响应分析和模态叠加法分析中的结构阻尼矩阵C的完整表达式为：同样,α,β,g可以通过定义材料输入也可以通过求解设置输入：3、模态叠加法分析模态叠加法中的阻尼控制在谐响应分析、瞬态分析、响应谱分析及随机振动分析中都支持以下表达式：4、数值阻尼数值阻尼并不是真实的阻尼,是人工抑制由高频结构产生的数值噪声;默认值为用来过滤掉虚假的高频模态;使用较小的值来过滤掉对最终结果影响较小的非物理响应;注意：数值阻尼只适用于瞬态分析;三、模态分析应用模态分析用来分析结构的振动特性自然频率和振型,是大多数动态分析得基础;假设和限制：结构是线性的M和K是常量.线性无阻尼的自由振动方程：假设{u}为简谐运动,则有因此求解行列式的特征值和特征向量;注意,{φ}为振型反应结构振幅的比例关系,可对质量矩阵进行正则化2、参与因子与有效质量参与因子：,其中{D}是笛卡尔坐标系中各个坐标轴单位位移响应;测量各个模态在各个方向运动的总质量,较大的值意味着该模态在该方向容易被激励;有效质量：理论上,各个方向的有效质量的总和应该等于结构的总质量,但取决于模态展开的数量;3、模态展开方法接触：由于模态分析时线性分析,只允许Bonded和No separation,其他接触程序视为无接触;4、阻尼模态分析特征值是复杂的,特征值的虚数部分表示自然频率,而实数部分衡量系统的稳定性,正值不稳定,负值稳定;模态展开方法：四、谐响应分析应用输入条件：简谐变化的载荷力,压力和位移,多个载荷应具有同样得频率,力和位移可以是同相或异相;假设和限制：结构具有固定的或与频率相关的刚度,阻尼,质量,不允许非线性；所有的载荷位移按相同频率做简谐变化;当施加的载荷的频率接近结构的自然频率时,发生共振；增加阻尼降低响应的振幅；阻尼较小的变化都会导致共振区响应的大幅变化；谐响应的运动方程如下：求解方法有两种：1、全函数法,直接求解矩阵方程;该方法求解准确,但速度慢于MSUP且耗资源,支持几乎所有的载荷和边界条件,其中加速度、轴承载荷、力矩相角只能为0;2、模态叠加法MUSP,对方程进行坐标变换{u}={φ}{y},将{M}和{K}变换成对角矩阵进行解耦,再求解n个解耦的方程{y},其中{C}必须是是对称矩阵,此方法需先进行模态分析;模态叠加法是一种近似求解,准确度取决了模态的展开阶数,一般比FULL法快;基本设置：cluster results-include residual vector-在模态叠加分析中,当施加的载荷激励高阶模态时,动态响应将会很粗糙;因此采用residual vector方法,除了采用模态的特征向量,还利用附件的模态转换向量来计算高阶频率;五、响应谱分析响应谱分析主要用来替代时程分析来确定结构对时间变化载荷的响应：如地震载荷,风载,海浪载荷,活塞载荷,火箭发动机振动等;对于多自由度长时程的分析往往通过响应谱分析来近似快速的求解最大响应;1、响应谱响应谱一般是单自由度系统在给定时程内的最大响应,该响应可以是位移,速度和加速度;多个不同频率相同阻尼的单自由度振荡器K,C,M就可以绘制响应谱,其中阻尼已经包含在响应谱中,也可以给定其他的阻尼绘制相应的响应谱;位移,速度,加速度响应谱之间是可以相互转换的,转换公式如下：2、分析类型响应谱分析分为单点响应谱SPRS分析和多点响应谱分析MPRS.SPRS-已知激励方向和频率的响应谱作用在所有的支撑点上,通常用来分析建筑结构的地震载荷;参与因子γ是对给定自然频率结构响应的量度,表征每个模态对特定方向的响应贡献多少;对于每个特征频率ω,谱值S都可以通过对数插值从响应谱中得到,但超过响应谱频率不会进行插值,而是取最近点的谱值;模态系数A,定义为放大系数来乘以特征向量来给出每个模态的实际位移,计算公式如下;响应R,计算公式如下如果系统有多个模态,那么应该对各个模态下的响应R进行叠加组合响应谱分析计算最大的位移和应力响应,它不能准确计算实际响应,因此有以下3种叠加方法SRSS,CQC和ROSE;SRSS:以下情况,SRSS法不再适用：1)考虑近间距自然频率的模态2)考虑部分或全刚度响应的模态3)包含未展开的高阶频率4、如果各阶模态频率有足够的间距,可以使用SRSS法叠加;评判各阶模态是否是近距频率,对于不同的阻尼比有不同的评判准则;对于阻尼比ζ≤2%,如果fi<fj,且fj≤,则是近距频率；对于阻尼比ζ＞2%,如果fi<fj,且fj≤1+5ζfi,则是近距频率;对于近距频率模态,可选用CQC或ROSE进行叠加,其中纠正系数0≤ε≤1,ε=0,不纠正；ε=1,全纠正；0<ε<1,部分纠正;CQC和ROSE计算公式中ε是基于模态的频率和阻尼计算得到;CQC计算公式如下ROSE计算公式如下5、响应谱中有两个特征频率fsp峰值频率和fzpa0周期加速度区域低频区<fsp,不考虑模态纠正除非有近距频率,可用SRSS,CQC或ROSE;中频区在fsp和fzpa之间,由周期区向刚性区转变,模态包含周期部分和刚性部分,通常用系数α将响应分为周期部分和刚性部分;α=0,周期；α=1,刚性；0<α<1,部分周期部分刚性;高频区>fzpa,刚性区,模态需要完全纠正;计算α有两种方法：Lindley-Yow和Gupta;Lindley-Yow法：α=αSa, α=ZPA/Sa,ZPA-0周期的加速度,Sa第i阶频率的加速度;当Sa<ZPA,α=0；Sa=ZPA, α=1；Sa>ZPA,随着Sa的减小α增大;Gupta法：α=αf,Lindley-Yow法中刚性响应影响所有的模态其对应的频率响应Sa>ZPA,但不应该用于其模态频率f<fsp；Gupta法中刚性响应影响所有的模态只有其频率f>f1=fsp,因此Gupta法适用大部分情况,应优先选用;6、刚性响应计算首先如前面描述的单独进行各个模态的响应计算,当打开刚性响应影响Rigid Response Effect时,这些模态响应R就不再是进行直接组合,而是分为周期Rp和刚性部分Rr;刚性响应系数α可选择Gupta或Lindley-Yow法计算;周期部分和刚性部分响应计算如下：然后分别进行组合叠加,对于周期部分响应Rp可用SRSS,CQC或ROSE方法进行叠加,如果含有近距频率模态时需要纠正不能使用SRSS法;刚性部分响应Rr进行代数和叠加即可最后将周期部分响应和刚性部分响应进行组合得到总的响应Rt7、缺省质量响应进行模态分析时,我们不可能展开所有模态来考虑结构100%的质量,因此我们关心的模态中所有质量占总质量的百分比即为有效质量比率,但展开的最高模态频率因远大于响应谱的fzpa,才能得到较为准确的分析结果;有时需要展开的模态阶数太多,我们可以通过模态分析计算缺省的质量将其进行额外的响应分析Missing Mass Response,这样就不必展开的模态频率要远大于fzpa;当f>fzpa,加速度响应是刚性的,因此可以进行静态的加速度分析;1)首先可以计算频率大于fzpa总的惯性力FT2)计算各个模态的惯性力3)计算各模态惯性力的合力。

ANSYSWorkbench使⽤中99%的时候都会⽤到的操作本⽂源⽂来⾃公众号CAD初学者结合个⼈经验，介绍⼀些ANSYS Workbench使⽤过程中的⼀些实⽤操作，主要包括：印记⾯建⽴、局部⽹格信息读取、求解设置（载荷步、并⾏计算、求解过程信息查看）以及结果后处理（节点结果输出、Surface、Path等）。

1 模型中的印迹⾯经典版的ANSYS中，可以直接施加载荷在节点上从⽽实现某个局部范围上的载荷施加，但在ANSYSWorkbench中就不怎么⽅⾯。

Workbench中有⼀个功能可以实现在局部区域施加载荷，即创建Imprint face（印记⾯功能）。

该功能须在Geometry中进⾏编辑，随后在Mechanical中将载荷局部施加在所创建的印记⾯上。

对于外部导⼊的模型，geometry编辑时，先对操作对象进⾏解冻（Unfreeze），若为geometry所建模型则⽆需此操作。

根据需求，在所需平⾯内绘制载荷施加形状，这⾥为圆。

在modeling中对该草图进⾏拉伸，在拉伸选项中选择Imprint Face并generate。

完成印记⾯的添加如下：2 ⽹格2.1 ⽹格质量检查在Mesh→Statistics→Mesh metric中，可选择不同项对单元⽹格质量进⾏综合评估。

常⽤的包括单元质量（ElementQuality）、单元长宽⽐（AspectRatio）、雅克⽐（JacobianRatio）以及最⼤⾓度（MaximumCorner Angle）等。

通过合理的⽹格划分⽅法，综合考虑这⼏项单元质量指标，有助于计算过程的顺利进⾏（尤其是遇到⾮线性求解）。

⽹格质量：单元长宽⽐：单元雅克⽐：单元最⼤顶⾓：2.2 局部⽹格信息输出对于局部区域的⽹格信息，可通过建⽴Named selection导出信息。

右键选择Named Selection选项，选择Export，导出txt⽂件，即可得到该区域的⽹格及其节点信息，包括单元编号、单元类型、节点编号。

ansys施加⾯荷载1、如果是线性变化的，可以采⽤⽔压⽅式定义；1、如果可以⽤函数表⽰，则可以⽤函数来定义，2、也可以使⽤表⾯效应单元来定义；在ANSYS中如果要在⼀个⾯上施加沿某个⽅向变化的⾯荷载，需要有两步来完成：这⾥以⼀个在圆筒内表⾯加内⽔压⼒的例⼦进⾏说明。

第⼀步，设置⾯荷载变化规律。

如果⾯荷载沿Z向变化，后⾯指定⾯荷载从Z=100开始变化，并按斜率为-9800进⾏变化，可⽤如下语句sfgrad,pres,,z,100,-9800 !也就是准备在⾼100⽶的圆柱加内⽔压⼒吧第⼆步，施加⾯荷载。

在指定的⾯上施加按第⼀步设置的⾯荷载变化规律的⾯荷载。

SFA,P51X,1,PRES,0这个语句相当于在指定⾯上施加法向荷载(选圆筒体内表⾯)，在Z=100时荷载值为0，随Z坐标变化荷载值以变化率-9800进⾏变化，这样在Z=0时荷载值为-9800*100每次⽤sfgrad进⾏设置后仅对随后的sfa命令有效，直倒下次再⽤sfgrad进⾏设置。

在⾯上施加荷载后，对模型剖分后可以执⾏以下命令来查看加的⾯荷载是否正确/PSF,PRES,NORM,2,0,1 以箭头⽅式显⽰⾯荷载sftran 将⾯荷载转化到有限元模型上本⽂摘⾃《ANSYS⼯程分析进阶实例》---王呼佳、陈洪军主编，在此对本书作者表⽰感谢！⼀般可以通过两种⽅法施加⾯荷载，⼀是在表⾯上覆盖⼀层表⾯效应单元SURF153或SURF154；⼆是通过apdl语⾔编程施加。

基本思路如下：⼈为将⾯上压⼒荷载换算成集中⼒并施加到节点上。

施加集中⼒时，将合⼒分解为X,Y,Z⽅向的分⼒。

（1）选中所要施加压⼒的表⾯，在⾯上⽣成⼀层shell63单元。

（2）对⽣成的shell63单元，使⽤循环语句逐步进⾏以下操作。

（3）得到每个单元的⾯积及单元中⼼的X,Y,Z坐标值。

（4）将坐标值代⼊压⼒随坐标变化的函数式，得到单元中⼼点处的压⼒值，并乘以⾯积得到单元所受的合⼒。

ANSYS为最著名的FEA软件，工程设计人员可用它完成以下工作：建立模型或转换结构，产品，零件或系统等的CAD模型。

施加工作载荷或其它设计执行条件。

研究物理反应，例如应力水平，温度分布，或电磁场在开发过程的早期优化设计以降底生产成本。

由于恶劣或不可能的环境（如生化应用）而用其他方式不可行时，可用它进行模拟试验。

一个典型的ANSYS分析有三个步骤：1，建立模型。

2，施加载荷和得到求解。

3，后处理查看结果。

1，建立模型：1) 指定一个工作名和分析标题。

（如不指明系统默认工作名为file）工作名:运行intractive在initial jobname后把file改为你的工作名。

或者进入后用以下方法：命令输入: /FILNAMEGUI: Utility Menu>File>Change Jobname分析标题:命令: /TITLE (副标题/STITLE)GUI: Utility Menu>File>Change Title指定单位系统：/UNITS2)定义单元类型ET，1，......Main Menu>Preprocessor>Create>Elements>Elem Attributes3)定义单元常量REALMain Menu>Preprocessor>Create>Elements>Elem Attributes4)定义材料性质命令: MATGUI: Main Menu>Preprocessor>-Attributes->Define>Default Attribs线性材料命令: MPGUI: Main Menu>Preprocessor>Material Props>非线性材料命令: TBPLOT, TBLISTGUI:Main Menu>Preprocessor>Material Props>Data T ables>GraphMain Menu>Preprocessor>Material Props>Data T ables>List5)生成几何模型2，施加载荷和求解1）定义分析类型和选项命令: ANTYPEMain Menu>Preprocessor>Loads>New Analysis 或Main Menu> Preprocessor>Loads>Restart2）加载3）设置载荷步选项Main Menu>Solution>-Load Step Options4）求解问题命令: SOLVEGUI:Main Menu>Solution>Current LSMain Menu>Solution>solution_method多步求解命令: LSSOLVEGUI: Main Menu>Solution>From LS Files3. 后处理命令: /POST1Main Menu>General Postproc命令:/POST26Main Menu>TimeHist Postpro。

ansys中施加梯形载荷总结前段时间做的一个项目中，大量使用了梯形荷载，尤其是在柱坐标系下定义渐变荷载，查阅了一些资料，现将所学心得贴出，希望对您能有益处。

(希望斑竹加分，呵呵)1、使用格式sfgrad实验室slkcn sldir slzer斜坡lab：有效的表面荷载标签，如pres、conv、hflux等slkcn：斜率坐标系统的参考编号，默认为0。

sldir：斜率的方向。

Slzer：坡度基准值为零的坐标位置。

slope：每单位长度或每单位角度的载荷值。

然后可以使用命令SF、SFE、SFL和SFA再次应用曲面载荷，每个节点上的载荷值为：cvalue=value+（slope*（coord slzer））2。

如果取消先前定义的渐变，则可以在没有指定值的情况下定义sfgrad。

3.在笛卡尔坐标系中使用：sfgrad,pres,0,y,0,-25!斜率为-25nsel，!选择压力施加的节点sf，全部，压力，500！y=0时为500，y=10时为250，y=20时为04、在柱坐标系下应遵循的规则(在柱坐标系下施加渐变荷载必须遵守这两条规则)(1)slzer以度表示，slope以荷载/度表示。

（2）设置cscir，使要加载的曲面不会通过坐标系的奇点。

(3)选择slzer，使之与cscir设置一致。

如果奇异点在180度(cscir，kcn，0，默认)，slzer应在-180-180之间。

如果奇异点在0度处(cscir，kcn，1)，slzer应在0度-360度之间。

5.在柱坐标系中使用的示例。

因为做这个比较多而且相对在笛卡儿坐标系下复杂些，因此说的较多些作用在外部的楔形压力被施加到位于局部圆柱坐标系11中的半圆形壳体上，压力位置从-90位置的400逐渐变化到90度位置的580。

缺省情况下，奇异点位于柱坐标系中的180度，因此壳的坐标范围从-90-90度，施加命令流如下：当地人，11岁！定义局部圆柱坐标系sfgrad，pres，11，y，-90，1!指定压力作用于-90度，斜率为1个单位/度sf，全体，主席，400！在-90度时为400，在0度时为490，在90度时为580。

Ansys中轴孔处正弦、随位置变化载荷的施加方法边界条件：力为在下半圆弧施加大小为50sinθ，随圆弧变化压力载荷。

这里主要谈一下上图圆孔下半园面上随位置变化的正弦或余弦载荷的添加方法。

1：Parameters----functions----define/edit2：弹出界面如下：我们这里是一段函数，所以function type选择默认single equation，这里我们的载荷是随角度变化的，所以CSYS：0那里该选择1，（全局柱坐标系）。

柱坐标系的y为角度，所以在TIME 那里选择Y输入公式。

最后结果如下图所示。

注意，ansys里的角度是度，而进行计算的时候要转换成弧度。

因为全局柱坐标系正好在孔中心，故这里不用新建局部坐标系，如果不在孔中心，就必须先在孔中心建立一个局部坐标系。

另外0度从x那里开始的，所以要减去一个PI。

上面只是定义了一个函数，然后保存为sin.func。

3，定义table选择刚才定义的sin.func。

Table parameter name那里随便输入aaa作为名字，确保下面的坐标系为1。

4在下半面施加压力。

如下图选择existing table。

其他的不用管。

然后在这里选择刚才定义的AAA就行了。

5，计算结束后，想查看的话采用如下命令。

Plotctrls---symbols--Symbols对话框surface load symbols选择pressures就行了。

如果想改变显示效果，show pres and convect as选择arrows，就是箭头形式了。

6,最后效果如下图结果：位移云图如下：Von misses stress云图如下：在90度角的地方，应力集中严重，上图中红色圈中局部放大图如下修改云图标签最大值为300MPa，如下图：二：左端完全固定，单元为solid186，尺寸为5mm，施加压力和上面的一样为50sinθ。

修改标签最大应力为300MPa，von misses应力结果如下：。

ansys接触设置、施加载荷、后处理有限元
原理
有限元分析是现代工程设计中必不可少的工具，而在有限元分析软件中，ANSYS更是被广泛使用。

在进行ANSYS分析时，接触设置、施加载荷和后处理是三个非常重要的环节，下面我们来进行一些介绍：
1、接触设置：
接触是指两个物体表面之间的接触，包括摩擦和各种变形。

在ANSYS中，接触问题的处理是非常复杂的。

一般通过指定接触面和接触类型等参数来进行接触分析。

同时，还可以通过设置接触初始距离、接触材料和接触力等参数来精确模拟接触情况，从而得出精确的分析结果。

2、施加载荷：
在进行有限元分析时，施加载荷是必不可少的。

在ANSYS中，可以通过多种方式来施加载荷，比如输入位移、力、压力、温度等等。

同时，在施加载荷时需要考虑不同类型的加载方式，比如静载、动载和瞬态载荷等等。

为了得到准确的分析结果，需要根据实际情况选择合适的加载方式进行分析。

3、后处理：
有限元分析并不是最终的目标，最终的目标是得出准确的分析结果。

然而，对于大多数工程师来说，如何准确分析和解释计算结果是非常困难的。

因此，ANSYS提供了强大的后处理功能，可以对分析结果进行可视化处理和分析。

例如，绘制应力云图、变形云图和动态图等来展现分析结果，并进一步分析数据以得到更准确的结果。

总之，ANSYS是一种非常强大的有限元分析软件，然而，要想得到准确的结果，需要进行准确的接触设置、施加载荷和后处理等细节工作。

因此，在进行ANSYS分析时，需要根据实际情况选择合适的模型和参数，并进行合理的分析处理，以得到更加精确的结果。

2.1 载荷概述有限元分析的主要目的是检查结构或构件对一定载荷条件的响应。

因此，在分析中指定合适的载荷条件是关键的一步。

在ANSYS程序中，可以用各种方式对模型加载，而且借助于载荷步选项，可以控制在求解中载荷如何使用。

2.2 什么是载荷在ANSYS术语中，载荷（loads）包括边界条件和外部或内部作用力函数，如图2-1所示。

不同学科中的载荷实例为：结构分析：位移，力，压力，温度（热应变），重力热分析：温度，热流速率，对流，内部热生成，无限表面磁场分析：磁势，磁通量，磁场段，源流密度，无限表面电场分析：电势（电压），电流，电荷，电荷密度，无限表面流体分析：速度，压力图2-1 “载荷”包括边界条件以及其它类型的载荷载荷分为六类：DOF约束，力（集中载荷），表面载荷，体积载荷、惯性力及耦合场载荷。

·DOF constraint（DOF约束）将用一已知值给定某个自由度。

例如，在结构分析中约束被指定为位移和对称边界条件；在热力分析中指定为温度和热通量平行的边界条件。

·Force（力）为施加于模型节点的集中载荷。

例如，在结构分析中被指定为力和力矩；在热力分析中为热流速率；在磁场分析中为电流段。

·Surface load（表面载荷）为施加于某个表面上的分布载荷。

例如，在结构分析中为压力；在热力分析中为对流和热通量。

·Body load（体积载荷）为体积的或场载荷。

例如，在结构分析中为温度和fluences；在热力分析中为热生成速率；在磁场分析中为流密度。

·Inertia loads（惯性载荷）由物体惯性引起的载荷，如重力加速度，角速度和角加速度。

主要在结构分析中使用。

·Coupled-field loads（耦合场载荷）为以上载荷的一种特殊情况，从一种分析得到的结果用作为另一分析的载荷。

例如，可施加磁场分析中计算出的磁力作为结构分析中的力载荷。

其它与载荷有关的术语的定义在下文中出现。

ANSYS/LS-DYNA软件压力载荷加载总结隐式分析：1、进入Solution项，找到Pressure项进行施加压力载荷，选取作用面（Areas）上。

2、设置Plotctrls/Symbal/定义载荷显示形式，如Pressure（压力）以Arrow（箭头表示方向），这样设置好，可以帮助判断施加载荷位置及方向是否正确。

3、一旦施加载荷后，压力的方向将以箭头形式出现。

显式分析：1、只有Element形式加载压力载荷。

2、定义压力载荷载体Element组元。

3、定义时间-压力载荷曲线，没有正负之分。

4、设置Plotctrls/Symbal/定义载荷显示形式，如Pressure（压力）以Arrow（箭头表示方向）。

可以帮助判断施加载荷位置及方向是否正确。

一旦施加载荷后，压力的方向将以箭头形式出现。

5、一般要以映射网格划分有限元模型，正常状态下单元有6个面，并且分配了面号。

施加压力载荷：Specify Load/PRES，面号需要通过Arrow方向来以手工方式逐个试出来（因为每个每次模型的方向和顺序是不一样的），选时间和载荷。

6、加载完毕。

隐式-显式分析步骤：1. 求解分析的隐式部分，从而得到预载求解分析的隐式部分(预加载荷)2. 改变现在的文件名进行显式求解部分Utility Menu > File > Change Jobname > Jobname2 （由原来的Jobname1改为Jobname2）…为了防止显式求解的结果覆盖隐式求解的结果3. 将隐式单元改为相应的显式单元Preprocessor > Element Type > Switch Elem Type …> Implic to Explic如果使用了非对应的单元，使用ETCHG, ITE 命令不能自动地将它们转变，而是用EMODIF命令手动将它们转变。

Preprocessor > Move/Modify > -Elements- Modify Attrib >Select elements to be modified > Elem Type –TYPE (STLOC field) > TYPE参考号与显式单元关联若LINK160, BEAM161, and LINK167 单元都需要第三个节点（方向点）, 所以如果相应的隐式单元只定义了端点，那么必须增加第三个节点。