ANSYS静力分析.ppt

《ANSYS静力分析.ppt》由会员分享，可在线阅读，更多相关《ANSYS静力分析.ppt（76页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

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

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

4、P7前处理程序定义单元类型、实常数、材料特性、模型的几何参数。 n在线性静力分析中，单元类型同样可以使用线性或者是非线性的单元类型pGUI: Main MenuPreprocessorElement Type Add/Edit/DeletepCommand: ET2.2 结构线性静力分析基本步骤p2.2.1 建模 n选择的材料特性可以是线性或者是非线性，可以是各向同性或者各向异性材料，并且可以随温度变化或者与温度无关。pGUI: Main MenuPreprocessorMaterial Props-Constant-Isotropic/OrthotropicpCommand: MP2.2 结

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

6、对话框（GUI: SolutionSoln Contros）pBasicAnalysis Options:选择Small Displacement StaticpSoln Option选项指定采用的求解器p实际上，求解控制对话框的绝大多数默认选项对于静力线性分析是合适的，用户只需要作很少的设置。2.2 结构线性静力分析基本步骤p2.2.2 施加载荷并求解 n施加载荷，设定载荷步p静力线性分析的载荷可以是DOF约束、集中载荷、分布力、温度、重力和旋转惯性力p载荷步可以在求解控制对话框中设置，不设置载荷子步n保存数据文件n开始求解计算n退出求解器pCommand: FinishpGUI: Main

7、 MenuFinish2.2 结构线性静力分析基本步骤p2.2.3 计算结果的后处理与分析n在结构的静力分析中，其计算结果将被写入结果文件Jobname.RST中，一般结果文件包含了一下数据：p基本数据：主要是关于节点的位移信息（UX, UY, UZ, ROTX, ROTY, ROTZ）p导出数据：包括节点和单元应力、节点和单元应变、单元集中力和节点支反力等2.2 结构线性静力分析基本步骤p2.2.3 计算结果的后处理与分析n后处理包括两个模块：p通用后处理其POST1 检查整个模型在指定时间步（或子步）下的计算结果p时间历程后处理器POST26 用于非线性分析中特定加载历史下的结果跟踪。n无

8、论使用哪个后处理器，程序的数据库都必须包含求解时相同的计算模型，且计算结果文件Jobname.RST必须可用。 2.2 结构线性静力分析基本步骤p2.2.3 计算结果的后处理与分析n在结构线性静力分析中，后处理常用命令如下p a. 查看结果 从数据文件中读取数据 GUI: Utility MenuFileResume fromCommand: Resume 读取选定的结果数据GUI: Main MenuGeneral PostprocRead ResultsBy Load StepCommand: SET使用使用SETSET命令可以用载荷步、子步或通过时间来选择数据读命令可以用载荷步、子步或通

9、过时间来选择数据读取，如果数据文件中没有制定时间点上的数据结果，程序会取，如果数据文件中没有制定时间点上的数据结果，程序会通过线性插值计算出该时间点上的结果通过线性插值计算出该时间点上的结果2.2 结构线性静力分析基本步骤p2.2.3 计算结果的后处理与分析n在结构线性静力分析中，后处理常用命令如下pb. 显示变形结果 GUI: Main MenuGeneral PostprocPlot ResultsDeformed Shape Command: PLDISPpc. 列表出力和力矩的信息 列出节点的支反力和反力矩GUI: Main MenuGeneral PostprocList Resul

10、tsReaction SoluCommand: PRRSOL2.2 结构线性静力分析基本步骤p2.2.3 计算结果的后处理与分析n在结构线性静力分析中，后处理常用命令如下p列出节点单元的支反力和反力矩 Command: PRESOL, F(or M)p列出节点力和反力矩的和 GUI: Main MenuGeneral PostprocList ResultsEelement Solution GUI: Main MenuGeneral PostprocNodal CalcsTotal Force Sum Command: FSUM2.2 结构线性静力分析基本步骤p2.2.3 计算结果的后处理与

11、分析n在结构线性静力分析中，后处理常用命令如下nd. 单元的结果处理pGUI: Main MenuGeneral PostprocElement TableDefine TablepCommand: ETABLEp对于线性单元体，例如梁单元、杆单元、管单元等可以通过此对于线性单元体，例如梁单元、杆单元、管单元等可以通过此选项获得结果数据，如应力和应变等。选项获得结果数据，如应力和应变等。 2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1n问题描述p具有圆孔的矩形板在拉伸状态下的应力分布。1.0 m2.0 m的矩形板，厚度为0.03 m，中心圆孔直径为0.25 m，弹性模量为

12、207GPa，泊松比0.3，端部受拉伸载荷600 N。2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤pa. 建立几何模型1) 定义任务标题：Utility MenuFileChange Title，输入“Tensile Loading of a Rectangular Plate with a Hole”，单击OK。Command: /TITLE, Tensile Loading of a Rectangular Plate with a Hole2)定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Dele

13、t，出现Element Types列表框。单击Add出现单元类型库对话框，左侧列表中选择Structural Shell，右侧列表中选择Elastic 4node63，单击OK确定。Command: ET, 1, SOLID452.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤pa. 建立几何模型3)定义材料特性：Main MenuPreprocessorMaterial PropsMaterial Modes。在Define Material Model Behavior窗口中，双击StructuralLinearElastic Isotropic。在出现的

14、对话框中，EX输入207E9，PRXY输入0.3，单击OK，Material Model Number 1出现在Material Models Defined窗口左侧列表中，选择菜单MaterialExit退出 Command: MPDATA, EX, 1,207e9 MPDATA, PRXY, 1,0.3 2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤pa. 建立几何模型4)建立几何模型：Main MenuPreprocessorModelingCreateVolumesBlockBy Dimensions，出现Create Block by Dimen

15、sions对话框，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，弹出Create Cylinder by Dimensions对话框，RAD1输入0.1，Z1,Z2,输入-0.1,0.1，THETA1输入0，THETA2输入360，单击OK退出。2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分

16、析实例1nGUI分析步骤na. 建立几何模型p4)建立几何模型：矩形板中心出现一个圆柱体。Main Menu PreprocessorModelingOperateBooleansSubtractVolume，出现Subtract Volumes选择框，先选择矩形板单击OK，点击subtract volume中的apply，然后选择圆柱体选择next，选择的模型颜色会发生变化，被切除部分为红色，单击OK ，单击Subtract Volumes 中的OK。最后得到的模型就是中心带有孔的矩形板，如图2.1所示。 pCommand: BLOCK, -1,1,-0.5,0.5,-0.015,0.015

17、 CYLIND,0.1, ,-0.1,0.1,0,360, VSBV, 1, 2 2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤pa. 建立几何模型图图2.1 2.1 中心带孔的矩形板模型中心带孔的矩形板模型 2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤nb. 模型剖分p5)模型剖分：为了对应力集中区域进行较准确的捕捉，划分有限元网格之前，通常需要对几何模型进行适当的剖分，以利于网格的划分。选择Utility MenuWorkPlaneDisplay Working Plane，然后选择UtilityWorkPlan

18、eOffset WP by Increments，在Offset WP对话框的Degrees框中输入：0,-90,0然后点击OK确定。2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤nb.模型剖分p5)模型剖分：然后选择Main MenuPreprocessor ModelingOperateBooleansDivideVolu by WrkPlane，选择Pick All，图形窗口中将显示模型被工作平面一分为二。1.工作平面和要分割的面垂直，所以要旋转工作平面（沿Z方向切绕X轴转；沿X轴方向切绕Y轴转；沿Y轴方向切绕Z轴转)p2.要移动工作平面和你要切的面

19、的位置相交，并且交线不能是面的边线。类似地，通过移动工作平面的位置，最后将几何模型剖分为如图2.2所示。最后，选择Main MenuPreprocessorModelingOperateBooleans GlueVolumes，在对话框中选择Pick All，将剖分开的各部分模型粘接在一起。2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤pb.模型剖分图图2.2 2.2 几何模型的剖分示意图几何模型的剖分示意图 2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤nc. 网格划分p6)有限元网格划分：选择Main MenuPr

20、eprocessorMeshingMesh Tool，出现如图2.3所示的Mesh Tool对话框。选择Size controlLinesset，将圆孔周边的线段和中线小正方型的线段都设定为10段，厚度方向的线段设定为6段，如图2.3所示。然后选择Mesh处下拉菜单为volume，shape设定为sweep，点击sweep，然后点击select all，然后点击OK确定，最后得到模型的有限元网格如图2.4所示。 2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤nc. 网格划分图图2.3 2.3 划分网格的参数设定划分网格的参数设定 图图2.4 2.4 带孔平

21、板的网格划分带孔平板的网格划分 2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤nd. 施加载荷及边界条件p7)载荷与边界条件：选择Main MenuSolutionDefine LoadsApplyDisplacementOn Areas，选择图形窗口最左端的面，约束这些面上所有节点UX方向的自由度。然后选择Main MenuSolutionDefine LoadsApplyPressureOn Ares，选择图形窗口最右端的所有面，输入压力为-2000，点击OK确定。最后得到模型上的边界条件和载荷如图2.5所示。 2.2 结构线性静力分析基本步骤p2.2

22、.4 结构线性静力分析实例1nGUI分析步骤nd. 施加载荷及边界条件图图2.5 2.5 划分网格的参数设定划分网格的参数设定 2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤ne. 求解及后处理p8)求解：选择Main MenuSolutionSolveCurrent LS，点击OK开始计算。p9)结果分析：选择Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu，在Contour Nodal Solution Data对话框中选择Nodal SolutionStressvon Mises

23、Stress，点击OK，图形窗口将出现如图2.6所示的von Mises应力云图分布。2.2 结构线性静力分析基本步骤p2.2.4 结构线性静力分析实例1nGUI分析步骤ne. 求解及后处理图图2.6 2.6 计算结果显示计算结果显示MisesMises应力云图应力云图 2.2 结构线性静力分析实例p2.2.5 结构线性静力分析实例n问题描述n如图所示，长度为L=2m的悬臂梁，自由端受集中力F=6000N作用，弹性模量E=200GPa，柏松比=0.25，梁的界面为工字型，并且已知梁的截面形式和具体尺寸如下：b=0.06m，h=0.08m，t1=t2=t3=0.01m，如图2.7所示。图图2.7

24、 2.7 悬臂梁自由端受集中力的作用悬臂梁自由端受集中力的作用 2.2 结构线性静力分析实例p2.2.5 结构线性静力分析实例n求解步骤n1)建立有限元模型 p定义分析标题：FileChange Title，输入“Bending of a I-Shaped Beam”，单击OK。p定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Delet，出现Element Types列表框。单击Add，出现单元类型库对话框，在左侧列表中选择Structural Beam，右侧列表中选择3D tapered 44，单击OK。2.2 结构线性静力分析实例p2.2

25、.5 结构线性静力分析实例n求解步骤n1)建立有限元模型 p定义材料特性：Main MenuPreprocessorMaterial PropsMaterial Models，定义各向同性线弹性材料，弹性模量EX为200E9，PRXY输入0.25，然后退出材料属性对话框。 p定义截面形状：Main MenuPreprocessorSectionsbeamCommon Sections，出现Beam Tool对话框，在Sub-Type下拉菜单中选择“工”字型梁截面，输入截面参数：W1=0.06，W2=0.06，W3=0.08，t1=0.01，t2=0.01，t3=0.01，单击OK。2.2 结构

26、线性静力分析实例p2.2.5 结构线性静力分析实例n求解步骤n1)建立有限元模型 p建立模型：Main MenuPreprocessorCreateKeypointsIn Active CS，依次建立坐标为(0,0,0) 和 (2,0,0)的两个个关键点。执行Main Menu PreprocessorCreateLinesStraight Line，选择1、2关键点，生成一条直线段。p划分网格：Main MenuPreprocessorMesh Tool，将线段分为20段，然后在Mesh Tool对话框单击Mesh，选择以上的线段生成有限元单元。2.2 结构线性静力分析实例p2.2.5 结构

27、线性静力分析实例n求解步骤n2)施加载荷进行静力求解与分析 p定义分析类型：Main MenuSolutionUnabridged MenuNew Analysis，选择Static，单击OK。p施加载荷：Main MenuSolutionLoadsApplyDisplacementOn Nodes，将图形窗口中最左端的节点的所有自由度约束。执行Main MenuSolutionApplyForce/MomentOn Nodes，选择最右端节点，单击OK，在Lab选择FY，VALUE输入-6000，单击OK，在梁的自由端出现集中力标识。2.2 结构线性静力分析实例p2.2.5 结构线性静力分析

28、实例n求解步骤n2)施加载荷进行静力求解与分析 p保存数据库文件：Save_DBp求解：Main MenuSolutionSolveCurrent LS，查看求解信息，关闭求解状态窗口，单击OK开始求解，求解完成后单击Close关闭求解信息框。 n3)查看分析结果 p进入POST1，读入结果数据：Main MenuGeneral PostprocRead Results-LastSet。p观察结构变形：Main MenuGeneral PostprocPlot ResultsDeformed Shape，选择Def+undeformed，单击OK，结构变形前后的形状同时出现在图形窗口，如图2.

29、8所示。2.2 结构线性静力分析实例p2.2.5 结构线性静力分析实例n求解步骤n3)查看分析结果 图图2.8 2.8 梁变形前后的形状改变图梁变形前后的形状改变图2.2 结构线性静力分析实例p2.2.5 结构线性静力分析实例n求解步骤n3)查看分析结果 p显示梁单元空间形状：PlotCtrlsStyleSize and Shape，出现Size and Shape对话框。激活ESHAPE选项为ON，单击OK，图形窗口出现梁单元的空间三维效果。显示梁单元截面上的应力：Main MenuGeneral PostprocPlot Resultscontour plot Nodal Solution

30、，在左侧Item列表中选择Stress，在右侧的Comp列表中选择SX，单击OK，梁单元上出现应力分布图，如图2.9所示。2.2 结构线性静力分析实例p2.2.5 结构线性静力分析实例n求解步骤n3)查看分析结果 图图2.9 2.9 梁上沿梁上沿X X方向应力分布图方向应力分布图2.2 结构线性静力分析实例p2.2.5 结构线性静力分析实例n求解步骤n3)查看分析结果 p列表显示反作用力：Main MenuGeneral PostprocList ResultsReaction Solu，在右侧列表中选All Items，单击OK，得到反作用力列表。p绘制弯矩图：首先定义单元表：Main Me

31、nuGeneral PostprocElement TableDefine Table，得到Element Table Data列表框。单击Add，在Lab框中输入MI，Items列表框选中By sequence num, Comp列表中选取SMISC，输入2，单击OK。2.2 结构线性静力分析实例p2.2.5 结构线性静力分析实例n求解步骤n3)查看分析结果 p按照类似的方法将SMISC，15添加到单元表，并命名为MJ。执行Main MenuGeneral PostprocPlot ResultsLine Elem Res，LabI选择MI，LabJ选择MJ，单击OK，在图形窗口中得到弯矩图

32、。p（6）图形输出：PlotCtrlsHardCopyTo File。 课后练习p一个口径为10mm的通用扳手，在它的末端承受100N的水平力和20N竖直向下的力。试确定该扳手在这两个载荷作用下的应力强度值。2.4 静力非线性分析 p2.4.1 屈曲分析的概念n屈曲分析是一种用于确定结构开始变得不稳定时的临界载荷和屈曲模态形状的技术。nANSYS在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Structural以及ANSYS/Professional中，提供两种结构屈曲载荷和屈曲模态的分析方法：p特征值（线性）屈曲p非线性屈曲分析2.4 静力非线性分析p

33、2.4.1 屈曲分析的概念p特征值（线性）屈曲分析特征值屈曲分析用于预测一个理想弹性结构的理论屈曲强度。该方法相当于教科书中的弹性屈曲分析方法。例如，一个柱体结构的特征值屈曲分析的结果，将与经典欧拉解相当。但是，当初始缺陷和非线性使得很多实际结构都不是在其理论弹性屈曲强度处发生屈曲。因此，特征值屈曲分析经常得出非常保守结果，通常不能用于实际的工程分析。p非线性屈曲分析。非线性屈曲分析比线性屈曲分析更精确，故建议用于对实际结构的设计或计算。该方法用一种逐渐增加载荷的非线性静力分析技术来求得使结构开始变得不稳定时的临界载荷，如图所示。应用非线性技术，模型中就可以包括诸如初始缺陷、塑性、间隙、大变形

34、响应等特征。此外，使用偏离控制加载，用户还可以跟踪结构的后屈曲行为。2.4 静力非线性分析p2.4.2 结构线性屈曲分析的一般步骤n1. 建立有限元模型p定义工作文件名、标题和单元；p定义单元类型、单元实常数、材料性质。n2. 施加载荷并获得静力解p方法与一般静力分析相同，另外有几点需要说明：静力求解中必须激活预应力影响（PSTRES命令）；由于屈曲分析计算出的特征值表示屈曲载荷系数，因此通常只需要施加单位载荷；求解完成后使用FINISH命令退出求解器。 2.4 静力非线性分析p2.4.2 结构线性屈曲分析的一般步骤n3. 获得线性屈曲解p重新进入ANSYS求解器（GUI: Main Menu

35、Solution）；p定义分析的类型对求解选项进行设置。分析类型选择Eigen Buckling(特征值屈曲分析)（GUI: New Analysis）屈曲分析选项（GUI: Analysis Options）Command: BUCOP, Method, NMODE, SHIFTMethod：指定特征值提取方法NMODE：指定特征值提取的数目 2.4 静力非线性分析p2.4.2 结构线性屈曲分析的一般步骤n3. 获得线性屈曲解p定义载荷步选项只需设置输出选项，控制结构输出。p保存数据库文件p开始求解求解过程是输出特征值，若要得到屈曲模态形状还需要进行解的扩展。 注：有些情况程序可能会计算出负

36、的特征值，这表示结构在注：有些情况程序可能会计算出负的特征值，这表示结构在相反方向施加载荷也会发生屈曲。相反方向施加载荷也会发生屈曲。p退出求解器2.4 静力非线性分析p2.4.2 结构线性屈曲分析的一般步骤n4. 扩展解p重新进入求解器p激活Expansion Pass（GUI: Main MenuSolutionExpansion Pass）p设置扩展模态选项（GUI：Main MenuSolutionExpand Modes）p进行输出控制打印输出控制数据库和结果文件输出控制开始扩展处理退出求解器（FINISH） 2.4 静力非线性分析p2.4.2 结构线性屈曲分析的一般步骤n5. 查看

37、分析结果p进入后处理器POST1p列出屈曲载荷系数p显示屈曲变形p显示相对应力分布2.4 静力非线性分析p2.4.3 线性屈曲分析实例 n1)问题描述p一端固支，一端自由的柱型壳体，半径为200 mm，长度L=800 mm，弹性模量E=200 GPa，壁厚10 mm, 轴向载荷2000 KN。 n2) 建模与参数设定p定义分析标题：Utility MenuFileChange Title，输入“Buckling of a Shell”，单击OK。p定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Delete，出现Element Types列表框

38、。单击Add出现单元类型库的对话框，在左侧列表中选择Structural Shell，右侧列表中选择SHELL63，单击OK。 2.4 静力非线性分析p2.4.3 线性屈曲分析实例n2) 建模与参数设定p定义实常数：Main MenuPreprocessorReal ConstantAdd/Edit/Delete，出现Real Constants列表框。单击Add，出现Element Type for Real Constants 对话框。单击OK，出现Real Constants for Set number 1 Shell63对话框。在Shell Thickness中输入10，单击Ok。单

39、击Close关闭Real Constants列表框。 p定义材料特性：Main MenuPreprocessorMaterial PropsMaterial Models, 在Define Material Model Behavior窗口中，双击StrucutralLinearElastic Isotropic。在出现的对话框中，EX输入200E9，单击OK，Material Model Number 1出现在Material Models Defined 窗口左侧列表中，执行MaterialExit退出。2.4 静力非线性分析p2.4.3 线性屈曲分析实例n2) 建模与参数设定p通过线扫略

40、生成壳模型：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关闭对话框。点击Main MenuProprocessorModeling CreateLines LinesStraight Line，分别选取点3和点4，点击OK关闭对话框。p点击Main MenuProprocessorModelingOperate

41、ExtrudeLinesAbout Axis，拾取line 1，点击Ok确定，然后分别拾取Keypoint 1和Keypoint 2，点击OK确定，出现Sweep lines about axis对话框在ARC输入360， NSEG文本框输入4，点击OK确定，生成壳模型。2.4 静力非线性分析p2.4.3 线性屈曲分析实例n2) 建模与参数设定p划分有限元网格：Main MenuPreprocessorMeshingMesh Tool，跳出MeshTool对话框，点击Size Controls: Lines后的Set按钮，然后选取壳沿轴向的任意一条线段，点击OK确定，出现Elements Si

42、zes on Picked Lines对话框，然后在NDIV文本框输入30。采用统一的办法将任意一端圆周上的线段均分为10段。选择Main MenuPreprocessorMeshing MeshTool，在MeshTool对话框中shape选中Quad和map，然后点击sweep按钮，生成如图2.10所示有限元网格。 2.4 静力非线性分析p2.4.3 线性屈曲分析实例n2) 建模与参数设定图图2.10 2.10 圆壳的有限元网格圆壳的有限元网格2.4 静力非线性分析p2.4.3 线性屈曲分析实例n3) 施加载荷进行静力分析p定义分析类型：Main MenuSolutionAbridged

43、MenuAnalysis TypeNew Analysis，出现New Analysis对话框。选择Static，单击OK。执行Main MenuSolution Analysis Options，出现Static or Steady-State Analysis对话框。在SSTIFPSTRES中选定Prestress ON，单击OK。 2.4 静力非线性分析p2.4.3 线性屈曲分析实例n3) 施加载荷进行静力分析p施加载荷：Main MenuSolutionDefine LoadsApplyStrucutral DisplacementOn Lines，出现Apply U, ROT on

44、Lines对话框，在图形窗口中选择Z=0的所有圆周上的四条弧线，单击OK，在对话框中选择ALL DOF，单击OK。执行Main MenuSolutionDefine LoadsApplyStructural Force/MomentOn Nodes，选择Z=800的所有节点（共40个），在Apply F/M on Nodes对话框中，选择Direction of force/mom为FZ，VALUE设定为-500，点击OK。 2.4 静力非线性分析p2.4.3 线性屈曲分析实例n3) 施加载荷进行静力分析p保存数据库文件，ANSYS ToolbarSAVE_DB。p求静力解：Main Menu

45、SolutionSolve-Current LS，查看求解信息，单击OK开始求解，求解完成后单击Close.2.4 静力非线性分析p2.4.3 线性屈曲分析实例n4)特征值屈曲分析p定义分析类型：Main MenuSolutionNew Analysis，忽略跳出的警告窗口，在New Analysis对话框中选取Eigen Buckling，单击OK。p屈曲分析选项设置：Main MenuSolutionAnalysis Options, 出现Eigenvalue Buckling Options, 选定Block Lanczos提取模态，提取模态数输入4，单击Ok。 2.4 静力非线性分析p

46、2.4.3 线性屈曲分析实例n4)特征值屈曲分析p模态扩展选项：Main MenuSolutionLoad Step OptsExpansionPass Single ExpandExpand Modes，在NMODE中输入4，单击OK。p开始特征值屈曲分析：Main MenuSolutionSolve-Current LS，单击Ok开始求解。 n5) 查看求解结果p读入结果数据：Main MenuGeneral PostprocRead ResultsFirst Set。p列出屈曲载荷系数：Main MenuGeneral PostprocList ResultsDetailed Summa

47、ry。2.4 静力非线性分析p2.4.3 线性屈曲分析实例n5) 查看求解结果p观察屈曲变形和应力云图：Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu，出现Contour Nodal Solution Data对话框。选择Nodal SolutionStressdisplacement vector sum，在Undisplaced shape key的下拉菜单中选择Deformed shape with undeformed edge，单击Ok，结构变形前后的形状同时出现在图形窗口，如图2.11所示。 2.4 静力非线性

48、分析p2.4.3 线性屈曲分析实例n5) 查看求解结果图图2.11 一阶屈曲模态的一阶屈曲模态的von Mises应力分布图应力分布图 2.4 静力非线性分析p2.4.4 接触分析实例n问题描述p经典的赫兹接触问题：衬底材料的弹性模量为60GPa，泊松比为0.3，尺寸如图2.12所示；刚性压头的半径为1 mm，压入载荷20KN。材料之间的摩擦系数为0.15。 图图2.12 赫兹接触模型赫兹接触模型2.4 静力非线性分析p2.4.4 接触分析实例n接触分析过程p定义任务名称：启动ANSYS，选择材料FileChange Title.，输入“Contact Analysis”，单击OK退出。p设定

49、参数变量：选择主菜单下ParametersScalar Parameters.，弹出Scalar Parameters对话框。在Section下分别定义Contact=1.5、Radius=1、Thick=2和Length=12，然后退出Scalar Parameters对话框。 2.4 静力非线性分析p2.4.4 接触分析实例n接触分析过程p定义单元类型：选择Main MenuPreprocessor Element TypeAdd/Edit/ Delete，分别添加Solid下的4node 182单元、Contact下的2D target 169和2 nd surf171单元，点击OK退出

50、该对话框。p定义实常数：选择Main MenuPreprocessor Real Constants，选择Type 2 TARGE169，不修改任何数据，单击OK退出。单击Add.，选择Type 3 CONTA171，设定参数FKN为1，FTOLN为0.1，PMAX为0.1，PMIN为0.01，然后退出Real Constants对话框。 2.4 静力非线性分析p2.4.4 接触分析实例n接触分析过程p定义材料模型：选择Main MenuPreprocessor Material PropsMaterial Models，定义材料模型1为各向同性材料，EX=60e9，PRXY=0.3。回到De