1、有限元静力学分析,王晓军 航空科学与工程学院固体力学研究所,第二章 ANSYS静力分析,主要内容 结构静力分析简介 结构线性静力分析基本步骤 建模、施加载荷并求解、计算结果的后处理与分析、结构线性静力分析实例 结构非线性静力分析实例 静力非线性计算的高级应用 屈曲分析的概念、结构线性屈曲分析一般步骤、线性屈曲分析实例、接触分析实例 思考题,2.1 结构静力分析简介,静力分析计算在固定载荷作用下结构的响应，它不考虑惯性和阻尼影响，如结构受随时间变化载荷作用的情况 静力分析可以计算那些固定不变的惯性载荷对结构的影响（如重力和离心力），以及那些可以近似为等价静力作用的随时间变化载荷（如通常在许多建筑
2、规范中所定义的等价静力风载）的作用 静力分析用于计算由那些不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位移、应力、应变和力。 固定不变的载荷和响应是一种假定，即假定载荷和结构相应随时间的变化非常缓慢。,2.1 结构静力分析简介,静力分析所施加的载荷类型有 外部施加的作用力和压力 稳态的惯性力 强迫位移 温度载荷 能流 静力分析既可以是线性的也可以是非线性的 非线性静力分析包括：大变形、塑性、蠕变、应力刚化、接触和超弹性等。 本节主要讨论线性静力分析，对非线性静力分析只作简单介绍,2.2 结构线性静力分析基本步骤,在线性静力分析中，主要有以下几个关键步骤 建模 施加载荷和边界条件并求解 结
3、果的后处理与分析,2.2 结构线性静力分析基本步骤,2.2.1 建模 为了建模，用户首先应指定作业名（Jobname） 、分析标题（Title）和单位（Unit），然后应用PREP7前处理程序定义单元类型、实常数、材料特性、模型的几何参数。 在线性静力分析中，单元类型同样可以使用线性或者是非线性的单元类型 GUI: Main MenuPreprocessorElement Type Add/Edit/Delete Command: ET,2.2 结构线性静力分析基本步骤,2.2.1 建模 选择的材料特性可以是线性或者是非线性，可以是各向同性或者各向异性材料，并且可以随温度变化或者与温度无关。
4、GUI: Main MenuPreprocessorMaterial Props-Constant-Isotropic/Orthotropic Command: MP,2.2 结构线性静力分析基本步骤,2.2.1 建模 另外在选择材料特性时，还需要注意以下几点： 材料的基本弹性性能参数必须定义：如杨氏模量EX、泊松比PRXY 如果施加惯性载荷（如重力），必须定义材料的密度DENS 创建几何实体模型并划分网格得到有限元模型 如果施加温度载荷时，必须定义材料的热膨胀系数ALPX 对应力敏感区域应划分较细的有限元网格,2.2 结构线性静力分析基本步骤,2.2.2 施加载荷并求解 进入ANSYS求解器
5、（Main MenuSolution） 定义分析类型并对求解进行控制 分析类型选New Analysis (GUI: SolutionNew Analysis)，进入求解控制对话框（GUI: SolutionSoln Contros） BasicAnalysis Options:选择Small Displacement Static Soln Option选项指定采用的求解器 实际上，求解控制对话框的绝大多数默认选项对于静力线性分析是合适的，用户只需要作很少的设置。,2.2 结构线性静力分析基本步骤,2.2.2 施加载荷并求解 施加载荷，设定载荷步 静力线性分析的载荷可以是DOF约束、集中载荷
6、、分布力、温度、重力和旋转惯性力 载荷步可以在求解控制对话框中设置，不设置载荷子步 保存数据文件 开始求解计算 退出求解器 Command: Finish GUI: Main MenuFinish,2.2 结构线性静力分析基本步骤,2.2.3 计算结果的后处理与分析 在结构的静力分析中，其计算结果将被写入结果文件Jobname.RST中，一般结果文件包含了一下数据： 基本数据：主要是关于节点的位移信息（UX, UY, UZ, ROTX, ROTY, ROTZ） 导出数据：包括节点和单元应力、节点和单元应变、单元集中力和节点支反力等,2.2 结构线性静力分析基本步骤,2.2.3 计算结果的后处理
7、与分析 后处理包括两个模块： 通用后处理其POST1 检查整个模型在指定时间步（或子步）下的计算结果 时间历程后处理器POST26 用于非线性分析中特定加载历史下的结果跟踪。 无论使用哪个后处理器，程序的数据库都必须包含求解时相同的计算模型，且计算结果文件Jobname.RST必须可用。,2.2 结构线性静力分析基本步骤,2.2.3 计算结果的后处理与分析 在结构线性静力分析中，后处理常用命令如下a. 查看结果 从数据文件中读取数据 GUI: Utility MenuFileResume from Command: Resume 读取选定的结果数据 GUI: Main MenuGeneral
8、PostprocRead ResultsBy Load Step Command: SET 使用SET命令可以用载荷步、子步或通过时间来选择数据读取，如果数据文件中没有制定时间点上的数据结果，程序会通过线性插值计算出该时间点上的结果,2.2 结构线性静力分析基本步骤,2.2.3 计算结果的后处理与分析 在结构线性静力分析中，后处理常用命令如下 b. 显示变形结果 GUI: Main MenuGeneral PostprocPlot ResultsDeformed Shape Command: PLDISP c. 列表出力和力矩的信息 列出节点的支反力和反力矩 GUI: Main MenuGen
9、eral PostprocList ResultsReaction Solu Command: PRRSOL,2.2 结构线性静力分析基本步骤,2.2.3 计算结果的后处理与分析 在结构线性静力分析中，后处理常用命令如下 列出节点单元的支反力和反力矩 Command: PRESOL, F(or M) 列出节点力和反力矩的和 GUI: Main MenuGeneral PostprocList ResultsEelement Solution GUI: Main MenuGeneral PostprocNodal CalcsTotal Force Sum Command: FSUM,2.2 结构
10、线性静力分析基本步骤,2.2.3 计算结果的后处理与分析 在结构线性静力分析中，后处理常用命令如下 d. 单元的结果处理 GUI: Main MenuGeneral PostprocElement TableDefine Table Command: ETABLE 对于线性单元体，例如梁单元、杆单元、管单元等可以通过此选项获得结果数据，如应力和应变等。,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 问题描述 具有圆孔的矩形板在拉伸状态下的应力分布。1.0 m2.0 m的矩形板，厚度为0.03 m，中心圆孔直径为0.25 m，弹性模量为207GPa，泊松比0.3，端部受拉
11、伸载荷600 N。,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 a. 建立几何模型 1) 定义任务标题：Utility MenuFileChange Title，输入“Tensile Loading of a Rectangular Plate with a Hole”，单击OK。 Command: /TITLE, Tensile Loading of a Rectangular Plate with a Hole 2)定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Delet，出现Element Ty
12、pes列表框。单击Add出现单元类型库对话框，左侧列表中选择Structural Shell，右侧列表中选择Elastic 4node63，单击OK确定。 Command: ET, 1, SOLID45,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 a. 建立几何模型 3)定义材料特性：Main MenuPreprocessorMaterial PropsMaterial Modes。在Define Material Model Behavior窗口中，双击StructuralLinearElastic Isotropic。在出现的对话框中，EX输入20
13、7E9，PRXY输入0.3，单击OK，Material Model Number 1出现在Material Models Defined窗口左侧列表中，选择菜单MaterialExit退出 Command: MPDATA, EX, 1,207e9 MPDATA, PRXY, 1,0.3,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 a. 建立几何模型 4)建立几何模型：Main MenuPreprocessorModelingCreateVolumesBlockBy Dimensions，出现Create Block by Dimensions对话框，X
14、1，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 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分
15、析步骤 a. 建立几何模型 4)建立几何模型：矩形板中心出现一个圆柱体。Main Menu PreprocessorModelingOperateBooleansSubtractVolume，出现Subtract Volumes选择框，先选择矩形板单击OK，点击subtract volume中的apply，然后选择圆柱体选择next，选择的模型颜色会发生变化，被切除部分为红色，单击OK ，单击Subtract Volumes 中的OK。最后得到的模型就是中心带有孔的矩形板，如图2.1所示。 Command: BLOCK, -1,1,-0.5,0.5,-0.015,0.015CYLIND,0.1
16、, ,-0.1,0.1,0,360,VSBV, 1, 2,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 a. 建立几何模型,图2.1 中心带孔的矩形板模型,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 b. 模型剖分 5)模型剖分：为了对应力集中区域进行较准确的捕捉，划分有限元网格之前，通常需要对几何模型进行适当的剖分，以利于网格的划分。选择Utility MenuWorkPlaneDisplay Working Plane，然后选择UtilityWorkPlaneOffset WP by Increments，在
17、Offset WP对话框的Degrees框中输入：0,-90,0然后点击OK确定。,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 b.模型剖分 5)模型剖分：然后选择Main MenuPreprocessor ModelingOperateBooleansDivideVolu by WrkPlane，选择Pick All，图形窗口中将显示模型被工作平面一分为二。1.工作平面和要分割的面垂直，所以要旋转工作平面（沿Z方向切绕X轴转；沿X轴方向切绕Y轴转；沿Y轴方向切绕Z轴转) 2.要移动工作平面和你要切的面的位置相交，并且交线不能是面的边线。类似地，通过移
18、动工作平面的位置，最后将几何模型剖分为如图2.2所示。最后，选择Main MenuPreprocessorModelingOperateBooleans GlueVolumes，在对话框中选择Pick All，将剖分开的各部分模型粘接在一起。,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 b.模型剖分,图2.2 几何模型的剖分示意图,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 c. 网格划分 6)有限元网格划分：选择Main MenuPreprocessorMeshingMesh Tool，出现如图2.3所示的M
19、esh Tool对话框。选择Size controlLinesset，将圆孔周边的线段和中线小正方型的线段都设定为10段，厚度方向的线段设定为6段，如图2.3所示。然后选择Mesh处下拉菜单为volume，shape设定为sweep，点击sweep，然后点击select all，然后点击OK确定，最后得到模型的有限元网格如图2.4所示。,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 c. 网格划分,图2.3 划分网格的参数设定,图2.4 带孔平板的网格划分,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 d. 施加载
20、荷及边界条件 7)载荷与边界条件：选择Main MenuSolutionDefine LoadsApplyDisplacementOn Areas，选择图形窗口最左端的面，约束这些面上所有节点UX方向的自由度。然后选择Main MenuSolutionDefine LoadsApplyPressureOn Ares，选择图形窗口最右端的所有面，输入压力为-2000，点击OK确定。最后得到模型上的边界条件和载荷如图2.5所示。,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 d. 施加载荷及边界条件,图2.5 划分网格的参数设定,2.2 结构线性静力分析基本
21、步骤,2.2.4 结构线性静力分析实例1 GUI分析步骤 e. 求解及后处理 8)求解：选择Main MenuSolutionSolveCurrent LS，点击OK开始计算。 9)结果分析：选择Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu，在Contour Nodal Solution Data对话框中选择Nodal SolutionStressvon Mises Stress，点击OK，图形窗口将出现如图2.6所示的von Mises应力云图分布。,2.2 结构线性静力分析基本步骤,2.2.4 结构线性静力分析实例1
22、 GUI分析步骤 e. 求解及后处理,图2.6 计算结果显示Mises应力云图,2.2 结构线性静力分析实例,2.2.5 结构线性静力分析实例 问题描述 如图所示，长度为L=2m的悬臂梁，自由端受集中力F=6000N作用，弹性模量E=200GPa，柏松比=0.25，梁的界面为工字型，并且已知梁的截面形式和具体尺寸如下：b=0.06m，h=0.08m，t1=t2=t3=0.01m，如图2.7所示。,图2.7 悬臂梁自由端受集中力的作用,2.2 结构线性静力分析实例,2.2.5 结构线性静力分析实例 求解步骤 1)建立有限元模型 定义分析标题：FileChange Title，输入“Bending
23、 of a I-Shaped Beam”，单击OK。 定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Delet，出现Element Types列表框。单击Add，出现单元类型库对话框，在左侧列表中选择Structural Beam，右侧列表中选择3D tapered 44，单击OK。,2.2 结构线性静力分析实例,2.2.5 结构线性静力分析实例 求解步骤 1)建立有限元模型 定义材料特性：Main MenuPreprocessorMaterial PropsMaterial Models，定义各向同性线弹性材料，弹性模量EX为200E9，P
24、RXY输入0.25，然后退出材料属性对话框。 定义截面形状：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 结构线性静力分析实例,2.2.5 结构线性静力分析实例 求解步骤 1)建立有限元模型 建立模型：Main MenuPreprocessorCreateKeypointsIn Active CS，依次建立坐标为(0,0,0) 和 (2
25、,0,0)的两个个关键点。执行Main Menu PreprocessorCreateLinesStraight Line，选择1、2关键点，生成一条直线段。 划分网格：Main MenuPreprocessorMesh Tool，将线段分为20段，然后在Mesh Tool对话框单击Mesh，选择以上的线段生成有限元单元。,2.2 结构线性静力分析实例,2.2.5 结构线性静力分析实例 求解步骤 2)施加载荷进行静力求解与分析 定义分析类型：Main MenuSolutionUnabridged MenuNew Analysis，选择Static，单击OK。 施加载荷：Main MenuSol
26、utionLoadsApplyDisplacementOn Nodes，将图形窗口中最左端的节点的所有自由度约束。执行Main MenuSolutionApplyForce/MomentOn Nodes，选择最右端节点，单击OK，在Lab选择FY，VALUE输入-6000，单击OK，在梁的自由端出现集中力标识。,2.2 结构线性静力分析实例,2.2.5 结构线性静力分析实例 求解步骤 2)施加载荷进行静力求解与分析 保存数据库文件：Save_DB 求解：Main MenuSolutionSolveCurrent LS，查看求解信息，关闭求解状态窗口，单击OK开始求解，求解完成后单击Close关
27、闭求解信息框。 3)查看分析结果 进入POST1，读入结果数据：Main MenuGeneral PostprocRead Results-LastSet。 观察结构变形：Main MenuGeneral PostprocPlot ResultsDeformed Shape，选择Def+undeformed，单击OK，结构变形前后的形状同时出现在图形窗口，如图2.8所示。,2.2 结构线性静力分析实例,2.2.5 结构线性静力分析实例 求解步骤 3)查看分析结果,图2.8 梁变形前后的形状改变图,2.2 结构线性静力分析实例,2.2.5 结构线性静力分析实例 求解步骤 3)查看分析结果 显示梁
28、单元空间形状：PlotCtrlsStyleSize and Shape，出现Size and Shape对话框。激活ESHAPE选项为ON，单击OK，图形窗口出现梁单元的空间三维效果。显示梁单元截面上的应力：Main MenuGeneral PostprocPlot Resultscontour plot Nodal Solution，在左侧Item列表中选择Stress，在右侧的Comp列表中选择SX，单击OK，梁单元上出现应力分布图，如图2.9所示。,2.2 结构线性静力分析实例,2.2.5 结构线性静力分析实例 求解步骤 3)查看分析结果,图2.9 梁上沿X方向应力分布图,2.2 结构线
29、性静力分析实例,2.2.5 结构线性静力分析实例 求解步骤 3)查看分析结果 列表显示反作用力：Main MenuGeneral PostprocList ResultsReaction Solu，在右侧列表中选All Items，单击OK，得到反作用力列表。 绘制弯矩图：首先定义单元表：Main MenuGeneral PostprocElement TableDefine Table，得到Element Table Data列表框。单击Add，在Lab框中输入MI，Items列表框选中By sequence num, Comp列表中选取SMISC，输入2，单击OK。,2.2 结构线性静力分
30、析实例,2.2.5 结构线性静力分析实例 求解步骤 3)查看分析结果 按照类似的方法将SMISC，15添加到单元表，并命名为MJ。执行Main MenuGeneral PostprocPlot ResultsLine Elem Res，LabI选择MI，LabJ选择MJ，单击OK，在图形窗口中得到弯矩图。 （6）图形输出：PlotCtrlsHardCopyTo File。,课后练习,一个口径为10mm的通用扳手，在它的末端承受100N的水平力和20N竖直向下的力。试确定该扳手在这两个载荷作用下的应力强度值。,2.4 静力非线性分析,2.4.1 屈曲分析的概念 屈曲分析是一种用于确定结构开始变得
31、不稳定时的临界载荷和屈曲模态形状的技术。 ANSYS在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Structural以及ANSYS/Professional中，提供两种结构屈曲载荷和屈曲模态的分析方法： 特征值（线性）屈曲 非线性屈曲分析,2.4 静力非线性分析,2.4.1 屈曲分析的概念 特征值（线性）屈曲分析 特征值屈曲分析用于预测一个理想弹性结构的理论屈曲强度。该方法相当于教科书中的弹性屈曲分析方法。例如，一个柱体结构的特征值屈曲分析的结果，将与经典欧拉解相当。但是，当初始缺陷和非线性使得很多实际结构都不是在其理论弹性屈曲强度处发生屈曲。因此，特
32、征值屈曲分析经常得出非常保守结果，通常不能用于实际的工程分析。 非线性屈曲分析。 非线性屈曲分析比线性屈曲分析更精确，故建议用于对实际结构的设计或计算。该方法用一种逐渐增加载荷的非线性静力分析技术来求得使结构开始变得不稳定时的临界载荷，如图所示。应用非线性技术，模型中就可以包括诸如初始缺陷、塑性、间隙、大变形响应等特征。此外，使用偏离控制加载，用户还可以跟踪结构的后屈曲行为。,2.4 静力非线性分析,2.4.2 结构线性屈曲分析的一般步骤 1. 建立有限元模型 定义工作文件名、标题和单元； 定义单元类型、单元实常数、材料性质。 2. 施加载荷并获得静力解 方法与一般静力分析相同，另外有几点需要
33、说明： 静力求解中必须激活预应力影响（PSTRES命令）； 由于屈曲分析计算出的特征值表示屈曲载荷系数，因此通常只需要施加单位载荷； 求解完成后使用FINISH命令退出求解器。,2.4 静力非线性分析,2.4.2 结构线性屈曲分析的一般步骤 3. 获得线性屈曲解 重新进入ANSYS求解器（GUI: Main MenuSolution）； 定义分析的类型对求解选项进行设置。 分析类型选择Eigen Buckling(特征值屈曲分析)（GUI: New Analysis） 屈曲分析选项（GUI: Analysis Options） Command: BUCOP, Method, NMODE, SH
34、IFT Method：指定特征值提取方法 NMODE：指定特征值提取的数目,2.4 静力非线性分析,2.4.2 结构线性屈曲分析的一般步骤 3. 获得线性屈曲解 定义载荷步选项 只需设置输出选项，控制结构输出。 保存数据库文件 开始求解 求解过程是输出特征值，若要得到屈曲模态形状还需要进行解的扩展。注：有些情况程序可能会计算出负的特征值，这表示结构在相反方向施加载荷也会发生屈曲。 退出求解器,2.4 静力非线性分析,2.4.2 结构线性屈曲分析的一般步骤 4. 扩展解 重新进入求解器 激活Expansion Pass（GUI: Main MenuSolutionExpansion Pass）
35、设置扩展模态选项（GUI：Main MenuSolutionExpand Modes） 进行输出控制 打印输出控制 数据库和结果文件输出控制 开始扩展处理 退出求解器（FINISH）,2.4 静力非线性分析,2.4.2 结构线性屈曲分析的一般步骤 5. 查看分析结果 进入后处理器POST1 列出屈曲载荷系数 显示屈曲变形 显示相对应力分布,2.4 静力非线性分析,2.4.3 线性屈曲分析实例 1)问题描述 一端固支，一端自由的柱型壳体，半径为200 mm，长度L=800 mm，弹性模量E=200 GPa，壁厚10 mm, 轴向载荷2000 KN。 2) 建模与参数设定 定义分析标题：Utili
36、ty MenuFileChange Title，输入“Buckling of a Shell”，单击OK。 定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Delete，出现Element Types列表框。单击Add出现单元类型库的对话框，在左侧列表中选择Structural Shell，右侧列表中选择SHELL63，单击OK。,2.4 静力非线性分析,2.4.3 线性屈曲分析实例 2) 建模与参数设定 定义实常数：Main MenuPreprocessorReal ConstantAdd/Edit/Delete，出现Real Constan
37、ts列表框。单击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输入20
38、0E9，单击OK，Material Model Number 1出现在Material Models Defined 窗口左侧列表中，执行MaterialExit退出。,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
39、关闭对话框。点击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确定，生成壳模型。,2.4 静力非线性分析,2.4.3 线性屈曲分析实例 2) 建模与参数
40、设定 划分有限元网格：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所示有限元网格。,2.4 静力非线性分
41、析,2.4.3 线性屈曲分析实例 2) 建模与参数设定,图2.10 圆壳的有限元网格,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。,2.4 静力非线性分析,2
42、.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对话框中，选择Direction of
43、force/mom为FZ，VALUE设定为-500，点击OK。,2.4 静力非线性分析,2.4.3 线性屈曲分析实例 3) 施加载荷进行静力分析 保存数据库文件，ANSYS ToolbarSAVE_DB。 求静力解：Main MenuSolutionSolve-Current LS，查看求解信息，单击OK开始求解，求解完成后单击Close.,2.4 静力非线性分析,2.4.3 线性屈曲分析实例 4)特征值屈曲分析 定义分析类型：Main MenuSolutionNew Analysis，忽略跳出的警告窗口，在New Analysis对话框中选取Eigen Buckling，单击OK。 屈曲分析
44、选项设置：Main MenuSolutionAnalysis Options, 出现Eigenvalue Buckling Options, 选定Block Lanczos提取模态，提取模态数输入4，单击Ok。,2.4 静力非线性分析,2.4.3 线性屈曲分析实例 4)特征值屈曲分析 模态扩展选项：Main MenuSolutionLoad Step OptsExpansionPass Single ExpandExpand Modes，在NMODE中输入4，单击OK。 开始特征值屈曲分析：Main MenuSolutionSolve-Current LS，单击Ok开始求解。 5) 查看求解结
45、果 读入结果数据：Main MenuGeneral PostprocRead ResultsFirst Set。 列出屈曲载荷系数：Main MenuGeneral PostprocList ResultsDetailed Summary。,2.4 静力非线性分析,2.4.3 线性屈曲分析实例 5) 查看求解结果 观察屈曲变形和应力云图：Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu，出现Contour Nodal Solution Data对话框。选择Nodal SolutionStressdisplacement ve
46、ctor sum，在Undisplaced shape key的下拉菜单中选择Deformed shape with undeformed edge，单击Ok，结构变形前后的形状同时出现在图形窗口，如图2.11所示。,2.4 静力非线性分析,2.4.3 线性屈曲分析实例 5) 查看求解结果,图2.11 一阶屈曲模态的von Mises应力分布图,2.4 静力非线性分析,2.4.4 接触分析实例 问题描述 经典的赫兹接触问题：衬底材料的弹性模量为60GPa，泊松比为0.3，尺寸如图2.12所示；刚性压头的半径为1 mm，压入载荷20KN。材料之间的摩擦系数为0.15。,图2.12 赫兹接触模型,
47、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对话框。,2.4 静力非线性分析,2.4.4 接触分析实例 接触分析过程 定义单元类型：选择Main MenuPreproces
48、sor 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为0.01，然后退出Real Constants对话框。,2.4 静力非线性分析,2.4.4 接触分析实例 接触分析
49、过程 定义材料模型：选择Main MenuPreprocessor Material PropsMaterial Models，定义材料模型1为各向同性材料，EX=60e9，PRXY=0.3。回到Define Material Model Behavior对话框，选择StructuralFriction Coefficient，设定MU=0.15。 建立几何模型：选择Main MenuPreprocessor ModelingCreateKeypointsIn Active CS，依次分别建立坐标为(0,0,0)、(Contact,0,0)和(Length,0,0)的三个关键点。,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,5，单击OK完成面的建立。,2.4 静力非线性分析,
