ANSYS静力分析PPT教学课件

1、第二章 ANSYS静力分析 主要内容 结构静力分析简介 结构线性静力分析基本步骤 建模、施加载荷并求解、计算结果的后处理与分析、结构线性静力分析实例 结构非线性静力分析实例 静力非线性计算的高级应用 屈曲分析的概念、结构线性屈曲分析一般步骤、线性屈曲分析实例、接触分析实例 思考题第1页/共75页2.1 结构静力分析简介 静力分析计算在固定载荷作用下结构的响应，它不静力分析计算在固定载荷作用下结构的响应，它不考虑惯性和阻尼影响考虑惯性和阻尼影响，如结构受随，如结构受随时间变化载荷作用的情况时间变化载荷作用的情况 静力分析可以计算那些静力分析可以计算那些固定不变的惯性载荷固定不变的惯性载荷对结构的

2、影响（如重力和离心力），以及对结构的影响（如重力和离心力），以及那些可以近似为那些可以近似为等价静力作用等价静力作用的随时间变化载荷（如通常在许多建筑规范中所定义的的随时间变化载荷（如通常在许多建筑规范中所定义的等价静力风载）的作用等价静力风载）的作用 静力分析用于计算由那些不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位静力分析用于计算由那些不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位移、应力、应变和力。移、应力、应变和力。 固定不变的载荷和响应是一种假定，即假定载荷和结构相应随时间的变化非常缓慢。固定不变的载荷和响应是一种假定，即假定载荷和结构相应随时间的变化非常缓慢。第2页/

3、共75页2.1 结构静力分析简介 静力分析所施加的载荷类型有 外部施加的作用力和压力 稳态的惯性力 强迫位移 温度载荷 能流 静力分析既可以是线性的也可以是非线性的 非线性静力分析包括：大变形、塑性、蠕变、应力刚化、接触和超弹性等。 本节主要讨论线性静力分析，对非线性静力分析只作简单介绍第3页/共75页2.2 结构线性静力分析基本步骤 在线性静力分析中，主要有以下几个关键步骤 建模 施加载荷和边界条件并求解 结果的后处理与分析第4页/共75页2.2 结构线性静力分析基本步骤 2.2.1 建模 为了建模，用户首先应指定作业名（Jobname） 、分析标题（Title）和单位（Unit），然后应用

4、PREP7前处理程序定义单元类型、实常数、材料特性、模型的几何参数。 在线性静力分析中，单元类型同样可以使用线性或者是非线性的单元类型 GUI: Main MenuPreprocessorElement Type Add/Edit/Delete Command: ET第5页/共75页2.2 结构线性静力分析基本步骤 2.2.1 建模 选择的材料特性可以是线性或者是非线性，可以是各向同性或者各向异性材料，并且可以随温度变化或者与温度无关。 GUI: Main MenuPreprocessorMaterial Props-Constant-Isotropic/Orthotropic Command

5、: MP第6页/共75页2.2 结构线性静力分析基本步骤 2.2.1 建模 另外在选择材料特性时，还需要注意以下几点： 材料的基本弹性性能参数必须定义：如杨氏模量EX、泊松比PRXY 如果施加惯性载荷（如重力），必须定义材料的密度DENS 创建几何实体模型并划分网格得到有限元模型 如果施加温度载荷时，必须定义材料的热膨胀系数ALPX 对应力敏感区域应划分较细的有限元网格 第7页/共75页2.2 结构线性静力分析基本步骤 2.2.2 施加载荷并求解 进入ANSYS求解器（Main MenuSolution） 定义分析类型并对求解进行控制 分析类型选New Analysis (GUI: Solut

6、ionNew Analysis)，进入求解控制对话框（GUI: SolutionSoln Contros） BasicAnalysis Options:选择Small Displacement Static Soln Option选项指定采用的求解器 实际上，求解控制对话框的绝大多数默认选项对于静力线性分析是合适的，用户只需要作很少的设置。第8页/共75页2.2 结构线性静力分析基本步骤 2.2.2 施加载荷并求解 施加载荷，设定载荷步 静力线性分析的载荷可以是DOF约束、集中载荷、分布力、温度、重力和旋转惯性力 载荷步可以在求解控制对话框中设置，不设置载荷子步 保存数据文件 开始求解计算 退

7、出求解器 Command: Finish GUI: Main MenuFinish第9页/共75页2.2 结构线性静力分析基本步骤 2.2.3 计算结果的后处理与分析 在结构的静力分析中，其计算结果将被写入结果文件Jobname.RST中，一般结果文件包含了一下数据： 基本数据：主要是关于节点的位移信息（UX, UY, UZ, ROTX, ROTY, ROTZ） 导出数据：包括节点和单元应力、节点和单元应变、单元集中力和节点支反力等第10页/共75页2.2 结构线性静力分析基本步骤 2.2.3 计算结果的后处理与分析 后处理包括两个模块： 通用后处理其POST1 检查整个模型在指定时间步（或子

8、步）下的计算结果 时间历程后处理器POST26 用于非线性分析中特定加载历史下的结果跟踪。 无论使用哪个后处理器，程序的数据库都必须包含求解时相同的计算模型，且计算结果文件Jobname.RST必须可用。 第11页/共75页2.2 结构线性静力分析基本步骤 2.2.3 计算结果的后处理与分析 在结构线性静力分析中，后处理常用命令如下 a. 查看结果 从数据文件中读取数据 GUI: Utility MenuFileResume from Command: Resume 读取选定的结果数据 GUI: Main MenuGeneral PostprocRead ResultsBy Load Step

9、 Command: SET 使用使用SETSET命令可以用载荷步、子步或通过时间来选择数据读取，如果数据文件命令可以用载荷步、子步或通过时间来选择数据读取，如果数据文件中没有制定时间点上的数据结果，程序会通过线性插值计算出该时间点上的结中没有制定时间点上的数据结果，程序会通过线性插值计算出该时间点上的结果果第12页/共75页2.2 结构线性静力分析基本步骤 2.2.3 计算结果的后处理与分析 在结构线性静力分析中，后处理常用命令如下 b. 显示变形结果 GUI: Main MenuGeneral PostprocPlot ResultsDeformed Shape Command: PLDIS

10、P c. 列表出力和力矩的信息 列出节点的支反力和反力矩 GUI: Main MenuGeneral PostprocList ResultsReaction Solu Command: PRRSOL第13页/共75页2.2 结构线性静力分析基本步骤 2.2.3 计算结果的后处理与分析 在结构线性静力分析中，后处理常用命令如下 列出节点单元的支反力和反力矩 Command: PRESOL, F(or M) 列出节点力和反力矩的和 GUI: Main MenuGeneral PostprocList ResultsEelement Solution GUI: Main MenuGeneral P

11、ostprocNodal CalcsTotal Force Sum Command: FSUM第14页/共75页2.2 结构线性静力分析基本步骤 2.2.3 计算结果的后处理与分析 在结构线性静力分析中，后处理常用命令如下 d. 单元的结果处理 GUI: Main MenuGeneral PostprocElement TableDefine Table Command: ETABLE 对于线性单元体，例如梁单元、杆单元、管单元等可以通过此选项获得结果数据，如应力和应变对于线性单元体，例如梁单元、杆单元、管单元等可以通过此选项获得结果数据，如应力和应变等。等。 第15页/共75页2.2 结构线

12、性静力分析基本步骤 2.2.4 结构线性静力分析实例1 问题描述 具有圆孔的矩形板在拉伸状态下的应力分布。1.0 m2.0 m的矩形板，厚度为0.03 m，中心圆孔直径为0.25 m，弹性模量为207GPa，泊松比0.3，端部受拉伸载荷600 N。第16页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤a. 建立几何模型1) 定义任务标题：Utility MenuFileChange Title，输入“Tensile Loading of a Rectangular Plate with a Hole”，单击OK。Command: /TITLE, Te

13、nsile Loading of a Rectangular Plate with a Hole2)定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Delet，出现Element Types列表框。单击Add出现单元类型库对话框，左侧列表中选择Structural Shell，右侧列表中选择Elastic 4node63，单击OK确定。Command: ET, 1, SOLID45第17页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤a. 建立几何模型3)定义材料特性：Main MenuPreproc

14、essorMaterial PropsMaterial Modes。在Define Material Model Behavior窗口中，双击StructuralLinearElastic Isotropic。在出现的对话框中，EX输入207E9，PRXY输入0.3，单击OK，Material Model Number 1出现在Material Models Defined窗口左侧列表中，选择菜单MaterialExit退出 Command: MPDATA, EX, 1,207e9 MPDATA, PRXY, 1,0.3 第18页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静

15、力分析实例1GUI分析步骤a. 建立几何模型4)建立几何模型：Main MenuPreprocessorModelingCreateVolumesBlockBy Dimensions，出现Create Block by Dimensions对话框，X1，X2 X-coordinates分别填入-1,1，Y1，Y2 Y-coordinates分别填入-0.5,0.5，Z1，Z2 Z-coordinates分别填入-0.015,0.015，单击OK确定。然后Main Menu PreprocessorModelingCreateVolumesCylinderBy Dimensions，弹出Crea

16、te Cylinder by Dimensions对话框，RAD1输入0.1，Z1,Z2,输入-0.1,0.1，THETA1输入0，THETA2输入360，单击OK退出。第19页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤a. 建立几何模型4)建立几何模型：矩形板中心出现一个圆柱体。Main Menu PreprocessorModelingOperateBooleansSubtractVolume，出现Subtract Volumes选择框，先选择矩形板单击OK，点击subtract volume中的apply，然后选择圆柱体选择next，选择的

17、模型颜色会发生变化，被切除部分为红色，单击OK ，单击Subtract Volumes 中的OK。最后得到的模型就是中心带有孔的矩形板，如图2.1所示。 Command: BLOCK, -1,1,-0.5,0.5,-0.015,0.015 CYLIND,0.1, ,-0.1,0.1,0,360, VSBV, 1, 2 第20页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤a. 建立几何模型图2.1 中心带孔的矩形板模型 第21页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤b. 模型剖分5)模型剖分：为了对

18、应力集中区域进行较准确的捕捉，划分有限元网格之前，通常需要对几何模型进行适当的剖分，以利于网格的划分。选择Utility MenuWorkPlaneDisplay Working Plane，然后选择UtilityWorkPlaneOffset WP by Increments，在Offset WP对话框的Degrees框中输入：0,-90,0然后点击OK确定。第22页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤b.模型剖分5)模型剖分：然后选择Main MenuPreprocessor ModelingOperateBooleansDivideV

19、olu by WrkPlane，选择Pick All，图形窗口中将显示模型被工作平面一分为二。1.工作平面和要分割的面垂直，所以要旋转工作平面（沿Z方向切绕X轴转；沿X轴方向切绕Y轴转；沿Y轴方向切绕Z轴转)2.要移动工作平面和你要切的面的位置相交，并且交线不能是面的边线。类似地，通过移动工作平面的位置，最后将几何模型剖分为如图2.2所示。最后，选择Main MenuPreprocessorModelingOperateBooleans GlueVolumes，在对话框中选择Pick All，将剖分开的各部分模型粘接在一起。第23页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性

20、静力分析实例1GUI分析步骤b.模型剖分图2.2 几何模型的剖分示意图 第24页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤c. 网格划分6)有限元网格划分：选择Main MenuPreprocessorMeshingMesh Tool，出现如图2.3所示的Mesh Tool对话框。选择Size controlLinesset，将圆孔周边的线段和中线小正方型的线段都设定为10段，厚度方向的线段设定为6段，如图2.3所示。然后选择Mesh处下拉菜单为volume，shape设定为sweep，点击sweep，然后点击select all，然后点击OK确

21、定，最后得到模型的有限元网格如图2.4所示。 第25页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤c. 网格划分图2.3 划分网格的参数设定 图2.4 带孔平板的网格划分 第26页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤d. 施加载荷及边界条件7)载荷与边界条件：选择Main MenuSolutionDefine LoadsApplyDisplacementOn Areas，选择图形窗口最左端的面，约束这些面上所有节点UX方向的自由度。然后选择Main MenuSolutionDefine Load

22、sApplyPressureOn Ares，选择图形窗口最右端的所有面，输入压力为-2000，点击OK确定。最后得到模型上的边界条件和载荷如图2.5所示。 第27页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤d. 施加载荷及边界条件图2.5 划分网格的参数设定 第28页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤e. 求解及后处理8)求解：选择Main MenuSolutionSolveCurrent LS，点击OK开始计算。9)结果分析：选择Main MenuGeneral PostprocPlot

23、ResultsContour PlotNodal Solu，在Contour Nodal Solution Data对话框中选择Nodal SolutionStressvon Mises Stress，点击OK，图形窗口将出现如图2.6所示的von Mises应力云图分布。第29页/共75页2.2 结构线性静力分析基本步骤2.2.4 结构线性静力分析实例1GUI分析步骤e. 求解及后处理图2.6 计算结果显示Mises应力云图 第30页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例问题描述如图所示，长度为L=2m的悬臂梁，自由端受集中力F=6000N作用，弹性模量E=2

24、00GPa，柏松比=0.25，梁的界面为工字型，并且已知梁的截面形式和具体尺寸如下：b=0.06m，h=0.08m，t1=t2=t3=0.01m，如图2.7所示。图2.7 悬臂梁自由端受集中力的作用 第31页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤1)建立有限元模型 定义分析标题：FileChange Title，输入“Bending of a I-Shaped Beam”，单击OK。定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Delet，出现Element Types列表框。单击Add，出现单元类型库

25、对话框，在左侧列表中选择Structural Beam，右侧列表中选择3D tapered 44，单击OK。第32页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤1)建立有限元模型 定义材料特性：Main MenuPreprocessorMaterial PropsMaterial Models，定义各向同性线弹性材料，弹性模量EX为200E9，PRXY输入0.25，然后退出材料属性对话框。 定义截面形状：Main MenuPreprocessorSectionsbeamCommon Sections，出现Beam Tool对话框，在Sub-Type下拉菜单中选

26、择“工”字型梁截面，输入截面参数：W1=0.06，W2=0.06，W3=0.08，t1=0.01，t2=0.01，t3=0.01，单击OK。第33页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤1)建立有限元模型 建立模型：Main MenuPreprocessorCreateKeypointsIn Active CS，依次建立坐标为(0,0,0) 和 (2,0,0)的两个个关键点。执行Main Menu PreprocessorCreateLinesStraight Line，选择1、2关键点，生成一条直线段。划分网格：Main MenuPreprocesso

27、rMesh Tool，将线段分为20段，然后在Mesh Tool对话框单击Mesh，选择以上的线段生成有限元单元。第34页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤2)施加载荷进行静力求解与分析 定义分析类型：Main MenuSolutionUnabridged MenuNew Analysis，选择Static，单击OK。施加载荷：Main MenuSolutionLoadsApplyDisplacementOn Nodes，将图形窗口中最左端的节点的所有自由度约束。执行Main MenuSolutionApplyForce/MomentOn Nodes

28、，选择最右端节点，单击OK，在Lab选择FY，VALUE输入-6000，单击OK，在梁的自由端出现集中力标识。第35页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤2)施加载荷进行静力求解与分析 保存数据库文件：Save_DB求解：Main MenuSolutionSolveCurrent LS，查看求解信息，关闭求解状态窗口，单击OK开始求解，求解完成后单击Close关闭求解信息框。 3)查看分析结果 进入POST1，读入结果数据：Main MenuGeneral PostprocRead Results-LastSet。观察结构变形：Main MenuGen

29、eral PostprocPlot ResultsDeformed Shape，选择Def+undeformed，单击OK，结构变形前后的形状同时出现在图形窗口，如图2.8所示。第36页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤3)查看分析结果 图2.8 梁变形前后的形状改变图第37页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤3)查看分析结果 显示梁单元空间形状：PlotCtrlsStyleSize and Shape，出现Size and Shape对话框。激活ESHAPE选项为ON，单击OK，图形窗口出现梁单元的空间三

30、维效果。显示梁单元截面上的应力：Main MenuGeneral PostprocPlot Resultscontour plot Nodal Solution，在左侧Item列表中选择Stress，在右侧的Comp列表中选择SX，单击OK，梁单元上出现应力分布图，如图2.9所示。第38页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤3)查看分析结果 图2.9 梁上沿X方向应力分布图第39页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤3)查看分析结果 列表显示反作用力：Main MenuGeneral PostprocList

31、ResultsReaction Solu，在右侧列表中选All Items，单击OK，得到反作用力列表。绘制弯矩图：首先定义单元表：Main MenuGeneral PostprocElement TableDefine Table，得到Element Table Data列表框。单击Add，在Lab框中输入MI，Items列表框选中By sequence num, Comp列表中选取SMISC，输入2，单击OK。第40页/共75页2.2 结构线性静力分析实例2.2.5 结构线性静力分析实例求解步骤3)查看分析结果 按照类似的方法将SMISC，15添加到单元表，并命名为MJ。执行Main Me

32、nuGeneral PostprocPlot ResultsLine Elem Res，LabI选择MI，LabJ选择MJ，单击OK，在图形窗口中得到弯矩图。（6）图形输出：PlotCtrlsHardCopyTo File。 第41页/共75页课后练习 一个口径为10mm的通用扳手，在它的末端承受100N的水平力和20N竖直向下的力。试确定该扳手在这两个载荷作用下的应力强度值。第42页/共75页2.4 静力非线性分析 2.4.1 屈曲分析的概念屈曲分析是一种用于确定结构开始变得不稳定时的临界载荷和屈曲模态形状的技术。ANSYS在ANSYS/Multiphysics、ANSYS/Mechanic

33、al、ANSYS/Structural以及ANSYS/Professional中，提供两种结构屈曲载荷和屈曲模态的分析方法：特征值（线性）屈曲非线性屈曲分析第43页/共75页2.4 静力非线性分析2.4.1 屈曲分析的概念特征值（线性）屈曲分析特征值屈曲分析用于预测一个理想弹性结构的理论屈曲强度。该方法相当于教科书中的弹性屈曲分析方法。例如，一个柱体结构的特征值屈曲分析的结果，将与经典欧拉解相当。但是，当初始缺陷和非线性使得很多实际结构都不是在其理论弹性屈曲强度处发生屈曲。因此，特征值屈曲分析经常得出非常保守结果，通常不能用于实际的工程分析。非线性屈曲分析。非线性屈曲分析比线性屈曲分析更精确，

34、故建议用于对实际结构的设计或计算。该方法用一种逐渐增加载荷的非线性静力分析技术来求得使结构开始变得不稳定时的临界载荷，如图所示。应用非线性技术，模型中就可以包括诸如初始缺陷、塑性、间隙、大变形响应等特征。此外，使用偏离控制加载，用户还可以跟踪结构的后屈曲行为。第44页/共75页2.4 静力非线性分析2.4.2 结构线性屈曲分析的一般步骤1. 建立有限元模型定义工作文件名、标题和单元；定义单元类型、单元实常数、材料性质。2. 施加载荷并获得静力解方法与一般静力分析相同，另外有几点需要说明：静力求解中必须激活预应力影响（PSTRES命令）；由于屈曲分析计算出的特征值表示屈曲载荷系数，因此通常只需要

35、施加单位载荷；求解完成后使用FINISH命令退出求解器。 第45页/共75页2.4 静力非线性分析2.4.2 结构线性屈曲分析的一般步骤3. 获得线性屈曲解重新进入ANSYS求解器（GUI: Main MenuSolution）；定义分析的类型对求解选项进行设置。分析类型选择Eigen Buckling(特征值屈曲分析)（GUI: New Analysis）屈曲分析选项（GUI: Analysis Options）Command: BUCOP, Method, NMODE, SHIFTMethod：指定特征值提取方法NMODE：指定特征值提取的数目 第46页/共75页2.4 静力非线性分析2.

36、4.2 结构线性屈曲分析的一般步骤3. 获得线性屈曲解定义载荷步选项只需设置输出选项，控制结构输出。保存数据库文件开始求解求解过程是输出特征值，若要得到屈曲模态形状还需要进行解的扩展。 注：有些情况程序可能会计算出负的特征值，这表示结构在相反方向施加载荷也会发生屈曲。注：有些情况程序可能会计算出负的特征值，这表示结构在相反方向施加载荷也会发生屈曲。退出求解器第47页/共75页2.4 静力非线性分析2.4.2 结构线性屈曲分析的一般步骤4. 扩展解重新进入求解器激活Expansion Pass（GUI: Main MenuSolutionExpansion Pass）设置扩展模态选项（GUI：M

37、ain MenuSolutionExpand Modes）进行输出控制打印输出控制数据库和结果文件输出控制开始扩展处理退出求解器（FINISH） 第48页/共75页2.4 静力非线性分析2.4.2 结构线性屈曲分析的一般步骤5. 查看分析结果进入后处理器POST1列出屈曲载荷系数显示屈曲变形显示相对应力分布第49页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例 1)问题描述一端固支，一端自由的柱型壳体，半径为200 mm，长度L=800 mm，弹性模量E=200 GPa，壁厚10 mm, 轴向载荷2000 KN。 2) 建模与参数设定定义分析标题：Utility MenuFile

38、Change Title，输入“Buckling of a Shell”，单击OK。定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Delete，出现Element Types列表框。单击Add出现单元类型库的对话框，在左侧列表中选择Structural Shell，右侧列表中选择SHELL63，单击OK。 第50页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例2) 建模与参数设定定义实常数：Main MenuPreprocessorReal ConstantAdd/Edit/Delete，出现Real Constants列表框。

39、单击Add，出现Element Type for Real Constants 对话框。单击OK，出现Real Constants for Set number 1 Shell63对话框。在Shell Thickness中输入10，单击Ok。单击Close关闭Real Constants列表框。 定义材料特性：Main MenuPreprocessorMaterial PropsMaterial Models, 在Define Material Model Behavior窗口中，双击StrucutralLinearElastic Isotropic。在出现的对话框中，EX输入200E9，单击

40、OK，Material Model Number 1出现在Material Models Defined 窗口左侧列表中，执行MaterialExit退出。第51页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例2) 建模与参数设定通过线扫略生成壳模型：Main MenuProprocessorModeling Create KeypointsIn Active CS，出现Create Keypoints in Active Coordinate System对话框，分别建立四个关键点(0,0,0)，(0, 0, 800)，(200, 0,0)，(200, 0,800)，点击OK关

41、闭对话框。点击Main MenuProprocessorModeling CreateLines LinesStraight Line，分别选取点3和点4，点击OK关闭对话框。点击Main MenuProprocessorModelingOperateExtrudeLinesAbout Axis，拾取line 1，点击Ok确定，然后分别拾取Keypoint 1和Keypoint 2，点击OK确定，出现Sweep lines about axis对话框在ARC输入360， NSEG文本框输入4，点击OK确定，生成壳模型。第52页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例2) 建

42、模与参数设定划分有限元网格：Main MenuPreprocessorMeshingMesh Tool，跳出MeshTool对话框，点击Size Controls: Lines后的Set按钮，然后选取壳沿轴向的任意一条线段，点击OK确定，出现Elements Sizes on Picked Lines对话框，然后在NDIV文本框输入30。采用统一的办法将任意一端圆周上的线段均分为10段。选择Main MenuPreprocessorMeshing MeshTool，在MeshTool对话框中shape选中Quad和map，然后点击sweep按钮，生成如图2.10所示有限元网格。 第53页/共7

43、5页2.4 静力非线性分析2.4.3 线性屈曲分析实例2) 建模与参数设定图2.10 圆壳的有限元网格第54页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例3) 施加载荷进行静力分析定义分析类型：Main MenuSolutionAbridged MenuAnalysis TypeNew Analysis，出现New Analysis对话框。选择Static，单击OK。执行Main MenuSolution Analysis Options，出现Static or Steady-State Analysis对话框。在SSTIFPSTRES中选定Prestress ON，单击OK。

44、 第55页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例3) 施加载荷进行静力分析施加载荷：Main MenuSolutionDefine LoadsApplyStrucutral DisplacementOn Lines，出现Apply U, ROT on Lines对话框，在图形窗口中选择Z=0的所有圆周上的四条弧线，单击OK，在对话框中选择ALL DOF，单击OK。执行Main MenuSolutionDefine LoadsApplyStructural Force/MomentOn Nodes，选择Z=800的所有节点（共40个），在Apply F/M on Nodes

45、对话框中，选择Direction of force/mom为FZ，VALUE设定为-500，点击OK。 第56页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例3) 施加载荷进行静力分析保存数据库文件，ANSYS ToolbarSAVE_DB。求静力解：Main MenuSolutionSolve-Current LS，查看求解信息，单击OK开始求解，求解完成后单击Close.第57页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例4)特征值屈曲分析定义分析类型：Main MenuSolutionNew Analysis，忽略跳出的警告窗口，在New Analysis对

46、话框中选取Eigen Buckling，单击OK。屈曲分析选项设置：Main MenuSolutionAnalysis Options, 出现Eigenvalue Buckling Options, 选定Block Lanczos提取模态，提取模态数输入4，单击Ok。 第58页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例4)特征值屈曲分析模态扩展选项：Main MenuSolutionLoad Step OptsExpansionPass Single ExpandExpand Modes，在NMODE中输入4，单击OK。开始特征值屈曲分析：Main MenuSolutionS

47、olve-Current LS，单击Ok开始求解。 5) 查看求解结果读入结果数据：Main MenuGeneral PostprocRead ResultsFirst Set。列出屈曲载荷系数：Main MenuGeneral PostprocList ResultsDetailed Summary。第59页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例5) 查看求解结果观察屈曲变形和应力云图：Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu，出现Contour Nodal Solution Data对话框。

48、选择Nodal SolutionStressdisplacement vector sum，在Undisplaced shape key的下拉菜单中选择Deformed shape with undeformed edge，单击Ok，结构变形前后的形状同时出现在图形窗口，如图2.11所示。 第60页/共75页2.4 静力非线性分析2.4.3 线性屈曲分析实例5) 查看求解结果图2.11 一阶屈曲模态的von Mises应力分布图 第61页/共75页2.4 静力非线性分析2.4.4 接触分析实例问题描述经典的赫兹接触问题：衬底材料的弹性模量为60GPa，泊松比为0.3，尺寸如图2.12所示；刚性

49、压头的半径为1 mm，压入载荷20KN。材料之间的摩擦系数为0.15。 图2.12 赫兹接触模型第62页/共75页2.4 静力非线性分析2.4.4 接触分析实例接触分析过程定义任务名称：启动ANSYS，选择材料FileChange Title.，输入“Contact Analysis”，单击OK退出。设定参数变量：选择主菜单下ParametersScalar Parameters.，弹出Scalar Parameters对话框。在Section下分别定义Contact=1.5、Radius=1、Thick=2和Length=12，然后退出Scalar Parameters对话框。 第63页/共

50、75页2.4 静力非线性分析2.4.4 接触分析实例接触分析过程定义单元类型：选择Main MenuPreprocessor Element TypeAdd/Edit/ Delete，分别添加Solid下的4node 182单元、Contact下的2D target 169和2 nd surf171单元，点击OK退出该对话框。定义实常数：选择Main MenuPreprocessor Real Constants，选择Type 2 TARGE169，不修改任何数据，单击OK退出。单击Add.，选择Type 3 CONTA171，设定参数FKN为1，FTOLN为0.1，PMAX为0.1，PMIN

51、为0.01，然后退出Real Constants对话框。 第64页/共75页2.4 静力非线性分析2.4.4 接触分析实例接触分析过程定义材料模型：选择Main MenuPreprocessor Material PropsMaterial Models，定义材料模型1为各向同性材料，EX=60e9，PRXY=0.3。回到Define Material Model Behavior对话框，选择StructuralFriction Coefficient，设定MU=0.15。建立几何模型：选择Main MenuPreprocessor ModelingCreateKeypointsIn Acti

52、ve CS，依次分别建立坐标为(0,0,0)、(Contact,0,0)和(Length,0,0)的三个关键点。第65页/共75页2.4 静力非线性分析2.4.4 接触分析实例接触分析过程建立几何模型：选择Main MenuPreprocessor ModelingCopyKeypoints，弹出Copy Keypoints对话框，选定pick all 。设定DY=thick，单击OK退出。选择Main MenuPreprocessorModeling CreateAreas ArbitraryThrough Keypoints，依次单击1,2,5,4关键点，点击Apply，然后单击2,3,6

53、,5，单击OK完成面的建立。第66页/共75页2.4 静力非线性分析2.4.4 接触分析实例接触分析过程建立几何模型：选择Utility MenuWorkPlane Offset WP toXYZ Locations +，在Global Cartesian栏输入0，thick+radius，单击OK。选择Utility MenuWorkPlaneLocal Coordiante SystemsCreate Local CSAt WP Origin.，弹出Create Local CS at WP Origin对话框，在KCS栏中输入Cylindrical 1，单击OK退出。选择Main MenuPreprocessorModelingCreate KeypointsIn Active CS，依次建立编号为20-23的四个关键点，其坐标分别为(radius,180,0)、(radius,-90,0)、(radius,0,0)和(0,0,0)。选择Main MenuPreprocessorModelingCreateLinesIn Active Coord，分别将关键点20,21以及21,22连成线。第67页/共75页2.4 静力非线性分析2.4.4 接触分析实例接触分析过程划分有限元网格：选择Main Menu PreprocessorMeshingDefault Attri