ANSYS结构分析指南.docx

1、ANSYS结构分析指南ANSYS结构分析指南（上）第一章 结构分析概述1.1 结构分析定义结构分析是有限元分析方法最常用的一个应用领域。结构这个术语是一个广义的概念，它包括土木工程结构如桥梁和建筑物，汽车结构如车身骨架，海洋结构如船舶结构，航空结构如飞机机身，还包括机械零部件如活塞、传动轴等。1.2 结构分析的类型在 ANSYS 产品家族中有七种结构分析的类型。结构分析中计算得出的基本未知量(节点自由度)是位移。其他的一些未知量，如应变、应力和反力可通过节点位移导出。包含结构分析功能的ANSYS产品有：ANSYS/Multiphysics,ANSYS/Mechanical, ANSYS/Str

2、uctural和ANSYS/Professional。下面简单列出了这七种类型的结构分析：静力分析-用于求解静力载荷作用下结构的位移和应力等。包括线性和非线性分析。非线性分析涉及塑性、应力刚化、大变形、大应变、超弹性、接触面和蠕变等。模态分析-用于计算结构的固有频率和模态。提供了不同的模态提取方法。谐波分析-用于确定结构在随时间正弦变化的载荷作用下的响应。瞬态动力分析-用于计算结构在随时间任意变化的载荷作用下的响应，并且可计及上述静力分析中提到的所有的非线性特性。谱分析-是模态分析的扩展，用于计算由于响应谱或 PSD 输入(随机振动)引起的应力和应变。屈曲分析-用于计算曲屈载荷和确定曲屈模态。

3、ANSYS 可进行线性(特征值)屈曲和非线性曲屈分析。显式动力分析-ANSYS/LS-DYNA 可用于计算高度非线性动力学问题和复杂的接触问题。此外，除前面提到的七种分析类型外，还可以进行如下的特殊分析：断裂力学复合材料疲劳分析p-Method梁分析1.3 结构分析所使用的单元从简单的杆单元和梁单元，一直到较为复杂的层合壳单元和大应变实体单元，绝大多数的 ANSYS 单元类型都可用于结构分析。注意-显式动力分析只能采用显式动力单元(LINK160、BEAM161、PLANE162、SHELL163、SOLID164、COMBI165、MASS166、LINK167)。表1-1 结构单元类型分类

4、单元名说明杆LINK1,LINK8,LINK180LINK10梁BEAM3,BEAM4BEAM54,BEAM44BEAM23,BEAM24BEAM188,BEAM189管PIPE16,PIPE17,PIPE18PIPE59PIPE20,PIPE602-D实体PLANE42,PLANE82,PLANE182,PLANE183PLANE2HYPER84,HYPER56,HYPER74VISCO88VISCO106,VISCO108PLANE83,PLANE25PLANE145,PLANE146p-单元,63-D实体SOLID45,SOLID95,SOLID185,SOLID186SOLID92,S

5、OLID187SOLID46,SOLID191SOLID64,SOLID65HYPER86,HYPER58,HYPER158VISCO89VISCO107SOLID147,SOLID148p-单元,6壳SHELL93,SHELL63,SHELL41,SHELL43,SHELL181SHELL51,SHELL61SHELL91,SHELL99SHELL28SHELL150p-单元,6接触CONTAC48,CONTAC49,CONTA171,CONTA172,CONTA173,CONTA174CONTAC12,CONTAC52CONTAC26TARGE169,TARGE170耦合场FLUID29,

6、FLUID30,FLUID129,FLUID130,INFIN110,INFIN111PLANE13,SOLID5,SOLID98PLANE13,SOLID5,SOLID98PLANE13,SOLID5,SOLID62,SOLID98FLUID38,FLUID79,FLUID80,FLUID81FLUID116特殊COMBIN14,COMBIN40,COMBIN39MASS21COMBIN37SURF153,SURF154COMBIN7LINK11MATRIX27,MATRIX50显式动力分析LINK160BEAM161PLANE162SHELL163SOLID164COMBI165MASS1

7、66LINK1671.4 材料模式界面对于本书论述的分析，如果采用GUI交互式操作，用户可以通过直观的“材料模式交互界面”来定义材料特性。这种方法采用树状结构的材料分类，使用户在分析中选择合适的材料模式变得更加简单。具体方法见ANSYS Basic Analysis Guide1.2.4.4。对于显式动力分析(ANSYS/LS-DYNA)，材料定义见ANSYS /LS-DYNA Users Guide7.1。1.5 求解方法在 ANSYS 产品中，求解结构问题有两种方法：h-方法和p-方法。h-方法可用于任何类型的结构分析，而p-方法只能用于线性结构静力分析。根据所求的问题，h-方法通常需要比

8、p-方法更密的网格。p-方法在应用较粗糙的网格时，提供了求得适当精度的一种很好的途径。本书主要讨论h-方法，而6 则专门研究 p-方法。第二章 结构线性静力分析2.1 静力分析的定义静力分析计算在固定不变载荷作用下结构的响应，它不考虑惯性和阻尼影响-如结构受随时间变化载荷作用的情况。可是，静力分析可以计算那些固定不变的惯性载荷对结构的影响(如重力和离心力)，以及那些可以近似为等价静力作用的随时间变化载荷(如通常在许多建筑规范中所定义的等价静力风载和地震载荷)的作用。静力分析用于计算由那些不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位移、应力、应变和力。固定不变的载荷和响应是一种假定，即假

9、定载荷和结构响应随时间的变化非常缓慢。静力分析所施加的载荷包括：外部施加的作用力和压力稳态的惯性力(如重力和离心力)强迫位移温度载荷(对于温度应变)能流(对于核能膨胀)关于载荷，还可参见2.3.4。2.2 线性静力分析与非线性静力分析静力分析既可以是线性的也可以是非线性的。非线性静力分析包括所有类型的非线性：大变形、塑性、蠕变、应力刚化、接触(间隙)单元、超弹性单元等。本章主要讨论线性静力分析。对非线性静力分析只作简单介绍，其详细论述见ANSYS Structural Analysis Guide8。2.3 静力分析的求解步骤2.3.1建模首先用户应指定作业名和分析标题，然后通过 PREP7

10、前处理程序定义单元类型、实常数、材料特性、模型的几何元素。这些步骤是大多数分析类型共同的，并已在ANSYS Basic Analysis Guide1.2 论述。有关建模的进一步论述，见ANSYS Modeling and Meshing Guide。2.3.1.1注意事项在进行静力分析时，要注意如下内容：1、可以采用线性或非线性结构单元。2、材料特性可以是线性或非线性，各向同性或正交各向异性，常数或与温度相关的：必须按某种形式定义刚度(如弹性模量 EX，超弹性系数等)。对于惯性载荷(如重力等)，必须定义质量计算所需的数据，如密度DENS。对于温度载荷，必须定义热膨胀系数 ALPX。3、对于网

11、格密度，要注意：应力或应变急剧变化的区域(通常是用户感兴趣的区域)，需要比应力或应变近乎常数的区域较密的网格：在考虑非线性的影响时，要用足够的网格来得到非线性效应。如塑性分析需要相当的积分点密度，因而在高塑性变形梯度区需要较密的网格。2.3.2设置求解控制设置求解控制包括定义分析类型、设置一般分析选项、指定载荷步选项等。当进行结构静力分析时，可以通过“求解控制对话框”来设置这些选项。该对话框对于大多数结构静力分析都已设置有合适的缺省值，用户只需作很少的设置就可以了。我们推荐采用该对话框进行设置。如用户不喜欢采用求解控制对话框，则可应用 ANSYS 的标准求解命令集和相应的菜单(Main Men

12、uSolutionUnabridged Menuoption)来设置求解控制选项。关于求解控制对话框的总体情况，见ANSYS Basic Analysis Guide3.11。2.3.2.1进入求解控制对话框用户可通过选择(Main MenuSolution -Analysis Type- Soln Control)进入求解控制对话框。下面诸小节简要论述该对话框中各标签的选项。关于如何设置这些选项，可在按该标签的 Help 按钮进入帮助系统，得到详细介绍。2.3.2.2Basic标签在求解控制对话框中共有五个标签，这些标签按从基本到高级的顺序排列。根据这种排列方式，可使求解设置较为平顺。在进入

13、求解控制对话框时，缺省激活的是Basic标签。Basic标签中的设置，提供了分析中所需的最少数据。一旦在Basic标签中的设置满足以后，就不需要设置其他标签中的选项，除非因为要进行高级控制而修改其他缺省设置。按 OK 按钮以后，设置存储到ANSYS数据库，并关闭对话框。用户可以在Basic标签中设置的选项如 表2-1 所示。有关详细说明见该标签的Help帮助系统。表2-1选项详细信息指定分析类型ANTYPE,NLGEOMANSYS Basic Analysis Guide1.2.6.1ANSYS Basic Analysis Guide8ANSYS Basic Analysis Guide3.

14、16控制时间设置，包括载荷步结束的时间TIME,自动时间步AUTOTS,在一个载荷步中的子步数NSUBST或DELTIMANSYS Basic Analysis Guide2.4ANSYS Basic Analysis Guide2.7.1设置写到数据库中的结果数据OUTRESANSYS Basic Analysis Guide2.7.4在静力分析中需要特别注意的选项主要有：在设置 ANTYPE 和 NLGEOM 时，如进行一个新的分析并忽略大变形效应(如大挠度、大转角、大应变)时，请选择“Small Displacement Static”项。如预期有大挠度(如弯曲的长细杆)或大应变(如金属

15、成形问题)，则选择“Large Displacement Static”。如想重启动一个失败的非线性分析，或者用户已进行了完整的静力分析，而想指定其他载荷，则选择“Restart Current Analysis”项。在设置 TIME 时，记住这个载荷步选项指定该载荷步结束的时间，缺省值为1。对于后续的载荷步，缺省为1加上前一个载荷步指定的时间。虽然在静力分析(除蠕变、粘塑性或其他率相关材料行为外)中，时间没有物理意义，但可以用于追踪时间步和子步，见ANSYS Basic Analysis Guide2。在设置OUTRES时，请记住：缺省时只有1,000个结果集记录到结果文件(Jobname.

16、RST)中，如果超过这一数目(基于用户的 OUTRES 设置)，程序将出错停机。可以通过 /CONFIG,NRES 命令来增大这一限值，见ANSYS Basic Analysis Guide20。2.3.2.3Transient标签Transient 标签设置瞬态分析控制，只有在 Basic 标签中选择了瞬态分析时才能激活这一标签，如果在 Basic 标签中选择了静态分析，则这一标签不能设置。所以在这里暂不讨论。2.3.2.4Soln Options标签Soln Options 标签用于设置 表2-2 所列的选项。详细说明可从Help按钮进入帮助系统而得到。表2-2选项参见ANSYS Basi

17、c Analysis Guide指定方程求解器EQSLV3.23.10对于多重启动指定参数RESCONTROL3.16.2在静力分析中设置 EQSLV 时，选择下列求解器之一：程序选择求解器(ANSYS 将根据问题的领域自动选择一个求解器)；稀疏矩阵求解器(对线性和非线性、静力和完全瞬态分析，为缺省项)；PCG求解器(对于大模型/高波前，巨形结构推荐使用)；AMG的求解器(其应用与PCG求解器相同，但提供并行算法；在用于多处理器环境时，转向更快)；DDS求解器，通过网络在多处理器系统中提供并行算法；迭代求解器(自动选择；只适用于线性静力/完全瞬态结构分析，或稳态温度分析；推荐)；波前直接求解器

18、。注意-AMG 和 DDS 求解器，是ANSYS并行算法的一部分，需要单独购买。见ANSYS Advanced Analysis Techniques Guide9。2.3.2.5Nonlinear标签Nonlinear 标签用于设置 表2-3 所列的选项。详细内容可通过Help按钮进入帮助系统。表2-3选项参见ANSYS Basic Analysis Guide激活线性搜索LNSRCH8.5.2.8.58.10.2.3激活DOF解的预测PRED8.5.2.8.4指定每个子步的最大迭代次数NEQIT8.5.2.8.3指明是否包括蠕变计算RATE8.3.1.58.5.3.2.1设置收敛准则CNV

19、TOL8.5.2.8.2控制二分CUTCONTROL8.5.2.8.62.3.2.6Advanced NL标签Advanced NL 标签用于设置 表2-4 所列的选项。详细内容见该标签中的Help帮助系统。表2-4选项参见ANSYS Basic Analysis Guide指定分析结束准则NCNV8.5.2.8.3激活和终止弧长法控制ARCLEN,ARCTRM8.10.2.42.3.3设置其他求解选项本节讨论求解的其他选项的设置。由于很少用到这些选项，并且一般都采用其缺省设置，因此这些选项没有出现在求解控制对话框中。本节中许多选项是非线性选项，将在ANSYS Structural Analy

20、sis Guide8进一步讨论。2.3.3.1 应力刚度效应一些单元，如 18X 族单元，不论 SSTIF 如何，都包括了应力刚度效应。为了确定一个单元是否包括应力刚度，请见ANSYS Element Reference Manual说明。在缺省时，如果 NLGEOM 为 ON 的话，应力刚度效应为 ON。在下面的这些特殊情况下，用户可能会关闭应力刚度效应：应力刚度仅与非线性分析相关。如果进行线性分析NLGEOM,OFF，则可以关闭应力刚度。在分析之前，用户知道结构不会因屈曲(分叉或跳跃屈曲)而破坏。通常，包括应力刚度效应时，可以加速非线性分析收敛。请记住上面所述的各点，用户可能对一些看起来收

21、敛困难的特殊问题，选择关闭应力刚度效应，如局部破坏。命令：SSTIFGUI：Main MenuSolutionUnabridged MenuAnalysis Options2.3.3.2Newton-Raphson选项这一选项只能用于非线性分析中，它说明在求解时如何修正切线刚度矩阵。用户可以选下列选项之一：程序选择；完全；修正；初始刚度；完全并且非对称矩阵。命令：NROPT2.3.3.3预应力效应计算通过这一选项在同一模型中执行预应力分析，如预应力模态的分析。缺省值为 OFF。注意-应力刚度效应和预应力效应计算都控制应力刚度矩阵的生成，因此在一个分析中不能同时采用。如二者都指定，则最后选项将覆

22、盖前者。命令：PSTRESSGUI：Main MenuSolutionUnabridged MenuAnalysis Options2.3.3.4质量矩阵公式通过该选项在结构中施加惯性载荷(如重力或旋转载荷)。可以指定下列选项之一：缺省(与单元类型有关)；集中质量近似。注意-对于静力分析，用户所用的质量矩阵并不明显影响求解精度(假设网格密度足够)。然而，如果想在同一模型上作预应力动力分析，选择质量矩阵公式就很重要；参见动力分析的有关章节。命令：LUMPMGUI：Main MenuSolutionUnabridged MenuAnalysis Options2.3.3.5参考温度这个载荷步选项适

23、用于温度应变计算。可用MP,REFT命令来设置材料相关的参考温度。命令：TREFGUI：Main MenuSolution-Load Step Opts-OtherReference Temp2.3.3.6模态数这个载荷步选项用于轴对称简谐单元。命令：MODEGUI：Main MenuSolution-Load Step Opts-OtherFor Harmonic Ele2.3.3.7蠕变准则这个非线性载荷步选项为自动时间步指定蠕变准则。命令：CRPLIMGUI：Main MenuSolutionUnabridged Menu-Load Step Opts-NonlinearCreep Cr

24、iterion2.3.3.8输出选项这个载荷步选项用于指定在输出文件(Jobname.out)中包括任意结果数据。命令：OUTPRGUI：Main MenuSolutionUnabridged Menu-Load Step Opts-Output CtrlsSolu Printout 警告-应用多个 OUTPR 命令时，有时可能会有一些冲突，见ANSYS Basic Analysis Guide2.7.4。2.3.3.9外插应用这个载荷步选项，可以通过把单元积分点结果拷贝到节点上，而不是通过外插(存在材料非线性时，这是缺省设置)。命令：ERESXGUI：Main MenuSolutionUna

25、bridged Menu-Load Step Opts-Output CtrlsIntegration Pt2.3.4施加载荷设置了求解选项以后，就可以对模型施加载荷了。2.3.4.1载荷类型下面列出的所有载荷类型，都可应用于静力分析中。2.3.4.1.1位移(UX,UY,UZ,ROTX,ROTY,ROTZ)这些自由度约束常施加到模型边界上，用以定义刚性支承点。它们也可以用于指定对称边界条件以及已知运动的点。由标号指定的方向是按照节点座标系定义的。2.3.4.1.2力(FX,FY,FZ) 和力矩(MX,MY,MZ)这些集中力通常在模型的外边界上指定。其方向是按节点座标系定义的。2.3.4.1.

26、3压力 (PRES)这是表面载荷，通常作用于模型的外部。正压力为指向单元面（起到压缩的效果）。2.3.4.1.4温度(TEMP)温度用于研究热膨胀或热收缩(即温度应力)。如果要计算热应变的话，必须定义热膨胀系数。用户可以从热分析LDREAD中读入温度，或者直接指定温度(通过 BF 族命令)。2.3.4.1.5流(FLUE)用于研究膨胀(由于中子流或其他原因而引起的材料膨胀)或蠕变的效应。只在输入膨胀或蠕变方程时才能使用。2.3.4.1.6重力、旋转等整个结构的惯性载荷。如果要计算惯性效应，必须定义密度(或某种形式的质量)。2.3.4.2在模型上施加载荷除了与模型无关的惯性载荷以外，用户可以在模

27、型的几何实体(关键点、线、面)或在有限元模型(节点和单元)上定义载荷。用户还可以通过 TABLE 类型的数组参数(见2.3.4.2.1)施加边界条件或作为函数的边界条件(见2.6.15)。表2-5 汇总了静力分析可以使用的载荷。在一个分析中，可以施加、删除、操作或列表载荷。表2-5载荷类型分类定义这些荷在的命令和菜单位移(UX,UY,UZ,ROTX,ROTY,ROTZ)约束ANSYSBasicAnalysisGuide2.6.3力、力矩(FX,FY,FZ,MX,MY,MZ)力ANSYSBasicAnalysisGuide2.6.6压力(PRES)面载荷ANSYSBasicAnalysisGui

28、de2.6.7温度(TEMP)、流(FLUE)体载荷ANSYSBasicAnalysisGuide2.6.8重力、旋转等惯性载荷ANSYSBasicAnalysisGuide2.6.92.3.4.2.1应用TABLE类数组参数施加载荷用户可以通过 TABLE 类数组参数施加载荷。对于应用表格边界条件，参见ANSYS Basic Analysis Guide2.6.14。在结构分析中，有效的初变量有时间(TIME)、温度(TEMP)和位置(X,Y,Z)。定义随时间(TIME)变化的一维表时，时间必须为升序排列。用户可以通过命令或交互式定义表格数组参数。具体操作参见ANSYS APDL Progr

29、ammers Guide。2.3.4.3计算惯性解除用户可以通过静力分析来执行惯性解除计算，即计算与施加载荷反向平衡的加速度。用户可以把惯性解除想象成一个等价自由体分析。要在 SOLVE 命令之前应用这一命令作为惯性载荷命令的一部分。模型应当满足下面的要求：模型不应当包括轴对称单元、子结构、或非线性。不推荐使用混合2D和3D单元的模型。对于梁单元(BEAM23、BEAM24、BEAM44和BEAM54)以及分层单元(SHELL91、SHELL99、SOLID46和SOLID191)，忽略偏置和楔形效应。也忽略层状单元的不对称分层效应。把楔形变截面单元分解成数个单元将得出更精确的结果。必须提供质

30、量计算所需的数据，如密度。提供所需的最少位移约束，即保证不发生刚体运动即可。对于2D单元需要三个约束(根据单元类型，可能更少)，对于3D单元只需要6个约束(根据单元类型，可能更少)。附加的约束，如对称边界条件也是允许的，但必须对所有约束检查0反力，以确保在惯性解除分析中不出现过约束。应当指定对于惯性解除计算合适的载荷。命令：IRLE,1GUI：Main MenuSolution-Load Step Opts-OtherInertia Relief2.3.4.3.1惯性解除的输出通过 IRLIST 命令来打印惯性解除计算的输出。该输出包括平衡施加载荷所需要的平移和转动加速度，而且可用于其他程序来进行运动学研究。质量和惯性矩列表汇总是精确解