Workbench有限元静力学分析演示幻灯片

有限元静力学分析,王晓军航空科学与工程学院固体力学研究所,第二章ANSYS静力分析,主要内容结构静力分析简介结构线性静力分析基本步骤建模、施加载荷并求解、计算结果的后处理与分析、结构线性静力分析实例结构非线性静力分析实例静力非线性计算的高级应用屈曲分析的概念、结构线性屈曲分析一般步骤、线性屈曲分析实例、接触分析实例思考题,2.1结构静力分析简介,静力分析计算在固定载荷作用下结构的响应，它不考虑惯性和阻尼影响，如结构受随时间变化载荷作用的情况静力分析可以计算那些固定不变的惯性载荷对结构的影响（如重力和离心力），以及那些可以近似为等价静力作用的随时间变化载荷（如通常在许多建筑规范中所定义的等价静力风载）的作用静力分析用于计算由那些不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位移、应力、应变和力。固定不变的载荷和响应是一种假定，即假定载荷和结构相应随时间的变化非常缓慢。,2.1结构静力分析简介,静力分析所施加的载荷类型有外部施加的作用力和压力稳态的惯性力强迫位移温度载荷能流静力分析既可以是线性的也可以是非线性的非线性静力分析包括：大变形、塑性、蠕变、应力刚化、接触和超弹性等。本节主要讨论线性静力分析，对非线性静力分析只作简单介绍,2.2结构线性静力分析基本步骤,在线性静力分析中，主要有以下几个关键步骤建模施加载荷和边界条件并求解结果的后处理与分析,2.2结构线性静力分析基本步骤,2.2.1建模为了建模，用户首先应指定作业名（Jobname）、分析标题（Title）和单位（Unit），然后应用PREP7前处理程序定义单元类型、实常数、材料特性、模型的几何参数。在线性静力分析中，单元类型同样可以使用线性或者是非线性的单元类型GUI:MainMenuPreprocessorElementTypeAdd/Edit/DeleteCommand:ET,2.2结构线性静力分析基本步骤,2.2.1建模选择的材料特性可以是线性或者是非线性，可以是各向同性或者各向异性材料，并且可以随温度变化或者与温度无关。GUI:MainMenuPreprocessorMaterialProps-Constant-Isotropic/OrthotropicCommand:MP,2.2结构线性静力分析基本步骤,2.2.1建模另外在选择材料特性时，还需要注意以下几点：材料的基本弹性性能参数必须定义：如杨氏模量EX、泊松比PRXY如果施加惯性载荷（如重力），必须定义材料的密度DENS创建几何实体模型并划分网格得到有限元模型如果施加温度载荷时，必须定义材料的热膨胀系数ALPX对应力敏感区域应划分较细的有限元网格,2.2结构线性静力分析基本步骤,2.2.2施加载荷并求解进入ANSYS求解器（MainMenuSolution）定义分析类型并对求解进行控制分析类型选NewAnalysis(GUI:SolutionNewAnalysis)，进入求解控制对话框（GUI:SolutionSolnContros）BasicAnalysisOptions:选择SmallDisplacementStaticSolnOption选项指定采用的求解器实际上，求解控制对话框的绝大多数默认选项对于静力线性分析是合适的，用户只需要作很少的设置。,2.2结构线性静力分析基本步骤,2.2.2施加载荷并求解施加载荷，设定载荷步静力线性分析的载荷可以是DOF约束、集中载荷、分布力、温度、重力和旋转惯性力载荷步可以在求解控制对话框中设置，不设置载荷子步保存数据文件开始求解计算退出求解器Command:FinishGUI:MainMenuFinish,2.2结构线性静力分析基本步骤,2.2.3计算结果的后处理与分析在结构的静力分析中，其计算结果将被写入结果文件Jobname.RST中，一般结果文件包含了一下数据：基本数据：主要是关于节点的位移信息（UX,UY,UZ,ROTX,ROTY,ROTZ）导出数据：包括节点和单元应力、节点和单元应变、单元集中力和节点支反力等,2.2结构线性静力分析基本步骤,2.2.3计算结果的后处理与分析后处理包括两个模块：通用后处理其POST1检查整个模型在指定时间步（或子步）下的计算结果时间历程后处理器POST26用于非线性分析中特定加载历史下的结果跟踪。无论使用哪个后处理器，程序的数据库都必须包含求解时相同的计算模型，且计算结果文件Jobname.RST必须可用。,2.2结构线性静力分析基本步骤,2.2.3计算结果的后处理与分析在结构线性静力分析中，后处理常用命令如下a.查看结果从数据文件中读取数据GUI:UtilityMenuFileResumefromCommand:Resume读取选定的结果数据GUI:MainMenuGeneralPostprocReadResultsByLoadStepCommand:SET使用SET命令可以用载荷步、子步或通过时间来选择数据读取，如果数据文件中没有制定时间点上的数据结果，程序会通过线性插值计算出该时间点上的结果,2.2结构线性静力分析基本步骤,2.2.3计算结果的后处理与分析在结构线性静力分析中，后处理常用命令如下b.显示变形结果GUI:MainMenuGeneralPostprocPlotResultsDeformedShapeCommand:PLDISPc.列表出力和力矩的信息列出节点的支反力和反力矩GUI:MainMenuGeneralPostprocListResultsReactionSoluCommand:PRRSOL,2.2结构线性静力分析基本步骤,2.2.3计算结果的后处理与分析在结构线性静力分析中，后处理常用命令如下列出节点单元的支反力和反力矩Command:PRESOL,F(orM)列出节点力和反力矩的和GUI:MainMenuGeneralPostprocListResultsEelementSolutionGUI:MainMenuGeneralPostprocNodalCalcsTotalForceSumCommand:FSUM,2.2结构线性静力分析基本步骤,2.2.3计算结果的后处理与分析在结构线性静力分析中，后处理常用命令如下d.单元的结果处理GUI:MainMenuGeneralPostprocElementTableDefineTableCommand:ETABLE对于线性单元体，例如梁单元、杆单元、管单元等可以通过此选项获得结果数据，如应力和应变等。,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1问题描述具有圆孔的矩形板在拉伸状态下的应力分布。1.0m2.0m的矩形板，厚度为0.03m，中心圆孔直径为0.25m，弹性模量为207GPa，泊松比0.3，端部受拉伸载荷600N。,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤a.建立几何模型1)定义任务标题：UtilityMenuFileChangeTitle，输入“TensileLoadingofaRectangularPlatewithaHole”，单击OK。Command:/TITLE,TensileLoadingofaRectangularPlatewithaHole2)定义单元类型：MainMenuPreprocessorElementTypeAdd/Edit/Delet，出现ElementTypes列表框。单击Add出现单元类型库对话框，左侧列表中选择StructuralShell，右侧列表中选择Elastic4node63，单击OK确定。Command:ET,1,SOLID45,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤a.建立几何模型3)定义材料特性：MainMenuPreprocessorMaterialPropsMaterialModes。在DefineMaterialModelBehavior窗口中，双击StructuralLinearElasticIsotropic。在出现的对话框中，EX输入207E9，PRXY输入0.3，单击OK，MaterialModelNumber1出现在MaterialModelsDefined窗口左侧列表中，选择菜单MaterialExit退出Command:MPDATA,EX,1,207e9MPDATA,PRXY,1,0.3,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤a.建立几何模型4)建立几何模型：MainMenuPreprocessorModelingCreateVolumesBlockByDimensions，出现CreateBlockbyDimensions对话框，X1，X2X-coordinates分别填入-1,1，Y1，Y2Y-coordinates分别填入-0.5,0.5，Z1，Z2Z-coordinates分别填入-0.015,0.015，单击OK确定。然后MainMenuPreprocessorModelingCreateVolumesCylinderByDimensions，弹出CreateCylinderbyDimensions对话框，RAD1输入0.1，Z1,Z2,输入-0.1,0.1，THETA1输入0，THETA2输入360，单击OK退出。,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤a.建立几何模型4)建立几何模型：矩形板中心出现一个圆柱体。MainMenuPreprocessorModelingOperateBooleansSubtractVolume，出现SubtractVolumes选择框，先选择矩形板单击OK，点击subtractvolume中的apply，然后选择圆柱体选择next，选择的模型颜色会发生变化，被切除部分为红色，单击OK，单击SubtractVolumes中的OK。最后得到的模型就是中心带有孔的矩形板，如图2.1所示。Command:BLOCK,-1,1,-0.5,0.5,-0.015,0.015CYLIND,0.1,-0.1,0.1,0,360,VSBV,1,2,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤a.建立几何模型,图2.1中心带孔的矩形板模型,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤b.模型剖分5)模型剖分：为了对应力集中区域进行较准确的捕捉，划分有限元网格之前，通常需要对几何模型进行适当的剖分，以利于网格的划分。选择UtilityMenuWorkPlaneDisplayWorkingPlane，然后选择UtilityWorkPlaneOffsetWPbyIncrements，在OffsetWP对话框的Degrees框中输入：0,-90,0然后点击OK确定。,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤b.模型剖分5)模型剖分：然后选择MainMenuPreprocessorModelingOperateBooleansDivideVolubyWrkPlane，选择PickAll，图形窗口中将显示模型被工作平面一分为二。1.工作平面和要分割的面垂直，所以要旋转工作平面（沿Z方向切绕X轴转；沿X轴方向切绕Y轴转；沿Y轴方向切绕Z轴转)2.要移动工作平面和你要切的面的位置相交，并且交线不能是面的边线。类似地，通过移动工作平面的位置，最后将几何模型剖分为如图2.2所示。最后，选择MainMenuPreprocessorModelingOperateBooleansGlueVolumes，在对话框中选择PickAll，将剖分开的各部分模型粘接在一起。,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤b.模型剖分,图2.2几何模型的剖分示意图,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤c.网格划分6)有限元网格划分：选择MainMenuPreprocessorMeshingMeshTool，出现如图2.3所示的MeshTool对话框。选择SizecontrolLinesset，将圆孔周边的线段和中线小正方型的线段都设定为10段，厚度方向的线段设定为6段，如图2.3所示。然后选择Mesh处下拉菜单为volume，shape设定为sweep，点击sweep，然后点击selectall，然后点击OK确定，最后得到模型的有限元网格如图2.4所示。,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤c.网格划分,图2.3划分网格的参数设定,图2.4带孔平板的网格划分,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤d.施加载荷及边界条件7)载荷与边界条件：选择MainMenuSolutionDefineLoadsApplyDisplacementOnAreas，选择图形窗口最左端的面，约束这些面上所有节点UX方向的自由度。然后选择MainMenuSolutionDefineLoadsApplyPressureOnAres，选择图形窗口最右端的所有面，输入压力为-2000，点击OK确定。最后得到模型上的边界条件和载荷如图2.5所示。,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤d.施加载荷及边界条件,图2.5划分网格的参数设定,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤e.求解及后处理8)求解：选择MainMenuSolutionSolveCurrentLS，点击OK开始计算。9)结果分析：选择MainMenuGeneralPostprocPlotResultsContourPlotNodalSolu，在ContourNodalSolutionData对话框中选择NodalSolutionStressvonMisesStress，点击OK，图形窗口将出现如图2.6所示的vonMises应力云图分布。,2.2结构线性静力分析基本步骤,2.2.4结构线性静力分析实例1GUI分析步骤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)建立有限元模型定义分析标题：FileChangeTitle，输入“BendingofaI-ShapedBeam”，单击OK。定义单元类型：MainMenuPreprocessorElementTypeAdd/Edit/Delet，出现ElementTypes列表框。单击Add，出现单元类型库对话框，在左侧列表中选择StructuralBeam，右侧列表中选择3Dtapered44，单击OK。,2.2结构线性静力分析实例,2.2.5结构线性静力分析实例求解步骤1)建立有限元模型定义材料特性：MainMenuPreprocessorMaterialPropsMaterialModels，定义各向同性线弹性材料，弹性模量EX为200E9，PRXY输入0.25，然后退出材料属性对话框。定义截面形状：MainMenuPreprocessorSectionsbeamCommonSections，出现BeamTool对话框，在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)建立有限元模型建立模型：MainMenuPreprocessorCreateKeypointsInActiveCS，依次建立坐标为(0,0,0)和(2,0,0)的两个个关键点。执行MainMenuPreprocessorCreateLinesStraightLine，选择1、2关键点，生成一条直线段。划分网格：MainMenuPreprocessorMeshTool，将线段分为20段，然后在MeshTool对话框单击Mesh，选择以上的线段生成有限元单元。,2.2结构线性静力分析实例,2.2.5结构线性静力分析实例求解步骤2)施加载荷进行静力求解与分析定义分析类型：MainMenuSolutionUnabridgedMenuNewAnalysis，选择Static，单击OK。施加载荷：MainMenuSolutionLoadsApplyDisplacementOnNodes，将图形窗口中最左端的节点的所有自由度约束。执行MainMenuSolutionApplyForce/MomentOnNodes，选择最右端节点，单击OK，在Lab选择FY，VALUE输入-6000，单击OK，在梁的自由端出现集中力标识。,2.2结构线性静力分析实例,2.2.5结构线性静力分析实例求解步骤2)施加载荷进行静力求解与分析保存数据库文件：Save_DB求解：MainMenuSolutionSolveCurrentLS，查看求解信息，关闭求解状态窗口，单击OK开始求解，求解完成后单击Close关闭求解信息框。3)查看分析结果进入POST1，读入结果数据：MainMenuGeneralPostprocReadResults-LastSet。观察结构变形：MainMenuGeneralPostprocPlotResultsDeformedShape，选择Def+undeformed，单击OK，结构变形前后的形状同时出现在图形窗口，如图2.8所示。,2.2结构线性静力分析实例,2.2.5结构线性静力分析实例求解步骤3)查看分析结果,图2.8梁变形前后的形状改变图,2.2结构线性静力分析实例,2.2.5结构线性静力分析实例求解步骤3)查看分析结果显示梁单元空间形状：PlotCtrlsStyleSizeandShape，出现SizeandShape对话框。激活ESHAPE选项为ON，单击OK，图形窗口出现梁单元的空间三维效果。显示梁单元截面上的应力：MainMenuGeneralPostprocPlotResultscontourplotNodalSolution，在左侧Item列表中选择Stress，在右侧的Comp列表中选择SX，单击OK，梁单元上出现应力分布图，如图2.9所示。,2.2结构线性静力分析实例,2.2.5结构线性静力分析实例求解步骤3)查看分析结果,图2.9梁上沿X方向应力分布图,2.2结构线性静力分析实例,2.2.5结构线性静力分析实例求解步骤3)查看分析结果列表显示反作用力：MainMenuGeneralPostprocListResultsReactionSolu，在右侧列表中选AllItems，单击OK，得到反作用力列表。绘制弯矩图：首先定义单元表：MainMenuGeneralPostprocElementTableDefineTable，得到ElementTableData列表框。单击Add，在Lab框中输入MI，Items列表框选中Bysequencenum,Comp列表中选取SMISC，输入2，单击OK。,2.2结构线性静力分析实例,2.2.5结构线性静力分析实例求解步骤3)查看分析结果按照类似的方法将SMISC，15添加到单元表，并命名为MJ。执行MainMenuGeneralPostprocPlotResultsLineElemRes，LabI选择MI，LabJ选择MJ，单击OK，在图形窗口中得到弯矩图。（6）图形输出：PlotCtrlsHardCopyToFile。,课后练习,一个口径为10mm的通用扳手，在它的末端承受100N的水平力和20N竖直向下的力。试确定该扳手在这两个载荷作用下的应力强度值。,2.4静力非线性分析,2.4.1屈曲分析的概念屈曲分析是一种用于确定结构开始变得不稳定时的临界载荷和屈曲模态形状的技术。ANSYS在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Structural以及ANSYS/Professional中，提供两种结构屈曲载荷和屈曲模态的分析方法：特征值（线性）屈曲非线性屈曲分析,2.4静力非线性分析,2.4.1屈曲分析的概念特征值（线性）屈曲分析特征值屈曲分析用于预测一个理想弹性结构的理论屈曲强度。该方法相当于教科书中的弹性屈曲分析方法。例如，一个柱体结构的特征值屈曲分析的结果，将与经典欧拉解相当。但是，当初始缺陷和非线性使得很多实际结构都不是在其理论弹性屈曲强度处发生屈曲。因此，特征值屈曲分析经常得出非常保守结果，通常不能用于实际的工程分析。非线性屈曲分析。非线性屈曲分析比线性屈曲分析更精确，故建议用于对实际结构的设计或计算。该方法用一种逐渐增加载荷的非线性静力分析技术来求得使结构开始变得不稳定时的临界载荷，如图所示。应用非线性技术，模型中就可以包括诸如初始缺陷、塑性、间隙、大变形响应等特征。此外，使用偏离控制加载，用户还可以跟踪结构的后屈曲行为。,2.4静力非线性分析,2.4.2结构线性屈曲分析的一般步骤1.建立有限元模型定义工作文件名、标题和单元；定义单元类型、单元实常数、材料性质。2.施加载荷并获得静力解方法与一般静力分析相同，另外有几点需要说明：静力求解中必须激活预应力影响（PSTRES命令）；由于屈曲分析计算出的特征值表示屈曲载荷系数，因此通常只需要施加单位载荷；求解完成后使用FINISH命令退出求解器。,2.4静力非线性分析,2.4.2结构线性屈曲分析的一般步骤3.获得线性屈曲解重新进入ANSYS求解器（GUI:MainMenuSolution）；定义分析的类型对求解选项进行设置。分析类型选择EigenBuckling(特征值屈曲分析)（GUI:NewAnalysis）屈曲分析选项（GUI:AnalysisOptions）Command:BUCOP,Method,NMODE,SHIFTMethod：指定特征值提取方法NMODE：指定特征值提取的数目,2.4静力非线性分析,2.4.2结构线性屈曲分析的一般步骤3.获得线性屈曲解定义载荷步选项只需设置输出选项，控制结构输出。保存数据库文件开始求解求解过程是输出特征值，若要得到屈曲模态形状还需要进行解的扩展。注：有些情况程序可能会计算出负的特征值，这表示结构在相反方向施加载荷也会发生屈曲。退出求解器,2.4静力非线性分析,2.4.2结构线性屈曲分析的一般步骤4.扩展解重新进入求解器激活ExpansionPass（GUI:MainMenuSolutionExpansionPass）设置扩展模态选项（GUI：MainMenuSolutionExpandModes）进行输出控制打印输出控制数据库和结果文件输出控制开始扩展处理退出求解器（FINISH）,2.4静力非线性分析,2.4.2结构线性屈曲分析的一般步骤5.查看分析结果进入后处理器POST1列出屈曲载荷系数显示屈曲变形显示相对应力分布,2.4静力非线性分析,2.4.3线性屈曲分析实例1)问题描述一端固支，一端自由的柱型壳体，半径为200mm，长度L=800mm，弹性模量E=200GPa，壁厚10mm,轴向载荷2000KN。2)建模与参数设定定义分析标题：UtilityMenuFileChangeTitle，输入“BucklingofaShell”，单击OK。定义单元类型：MainMenuPreprocessorElementTypeAdd/Edit/Delete，出现ElementTypes列表框。单击Add出现单元类型库的对话框，在左侧列表中选择StructuralShell，右侧列表中选择SHELL63，单击OK。,2.4静力非线性分析,2.4.3线性屈曲分析实例2)建模与参数设定定义实常数：MainMenuPreprocessorRealConstantAdd/Edit/Delete，出现RealConstants列表框。单击Add，出现ElementTypeforRealConstants对话框。单击OK，出现RealConstantsforSetnumber1Shell63对话框。在ShellThickness中输入10，单击Ok。单击Close关闭RealConstants列表框。定义材料特性：MainMenuPreprocessorMaterialPropsMaterialModels,在DefineMaterialModelBehavior窗口中，双击StrucutralLinearElasticIsotropic。在出现的对话框中，EX输入200E9，单击OK，MaterialModelNumber1出现在MaterialModelsDefined窗口左侧列表中，执行MaterialExit退出。,2.4静力非线性分析,2.4.3线性屈曲分析实例2)建模与参数设定通过线扫略生成壳模型：MainMenuProprocessorModelingCreateKeypointsInActiveCS，出现CreateKeypointsinActiveCoordinateSystem对话框，分别建立四个关键点(0,0,0)，(0,0,800)，(200,0,0)，(200,0,800)，点击OK关闭对话框。点击MainMenuProprocessorModelingCreateLinesLinesStraightLine，分别选取点3和点4，点击OK关闭对话框。点击MainMenuProprocessorModelingOperateExtrudeLinesAboutAxis，拾取line1，点击Ok确定，然后分别拾取Keypoint1和Keypoint2，点击OK确定，出现Sweeplinesaboutaxis对话框在ARC输入360，NSEG文本框输入4，点击OK确定，生成壳模型。,2.4静力非线性分析,2.4.3线性屈曲分析实例2)建模与参数设定划分有限元网格：MainMenuPreprocessorMeshingMeshTool，跳出MeshTool对话框，点击SizeControls:Lines后的Set按钮，然后选取壳沿轴向的任意一条线段，点击OK确定，出现ElementsSizesonPickedLines对话框，然后在NDIV文本框输入30。采用统一的办法将任意一端圆周上的线段均分为10段。选择MainMenuPreprocessorMeshingMeshTool，在MeshTool对话框中shape选中Quad和map，然后点击sweep按钮，生成如图2.10所示有限元网格。,2.4静力非线性分析,2.4.3线性屈曲分析实例2)建模与参数设定,图2.10圆壳的有限元网格,2.4静力非线性分析,2.4.3线性屈曲分析实例3)施加载荷进行静力分析定义分析类型：MainMenuSolutionAbridgedMenuAnalysisTypeNewAnalysis，出现NewAnalysis对话框。选择Static，单击OK。执行MainMenuSolutionAnalysisOptions，出现StaticorSteady-StateAnalysis对话框。在SSTIFPSTRES中选定PrestressON，单击OK。,2.4静力非线性分析,2.4.3线性屈曲分析实例3)施加载荷进行静力分析施加载荷：MainMenuSolutionDefineLoadsApplyStrucutralDisplacementOnLines，出现ApplyU,ROTonLines对话框，在图形窗口中选择Z=0的所有圆周上的四条弧线，单击OK，在对话框中选择ALLDOF，单击OK。执行MainMenuSolutionDefineLoadsApplyStructuralForce/MomentOnNodes，选择Z=800的所有节点（共40个），在ApplyF/MonNodes对话框中，选择Directionofforce/mom为FZ，VALUE设定为-500，点击OK。,2.4静力非线性分析,2.4.3线性屈曲分析实例3)施加载荷进行静力分析保存数据库文件，ANSYSToolbarSAVE_DB。求静力解：MainMenuSolutionSolve-CurrentLS，查看求解信息，单击OK开始求解，求解完成后单击Close.,2.4静力非线性分析,2.4.3线性屈曲分析实例4)特征值屈曲分析定义分析类型：MainMenuSolutionNewAnalysis，忽略跳出的警告窗口，在NewAnalysis对话框中选取EigenBuckling，单击OK。屈曲分析选项设置：MainMenuSolutionAnalysisOptions,出现EigenvalueBucklingOptions,选定BlockLanczos提取模态，提取模态数输入4，单击Ok。,2.4静力非线性分析,2.4.3线性屈曲分析实例4)特征值屈曲分析模态扩展选项：MainMenuSolutionLoadStepOptsExpansionPassSingleExpandExpandModes，在NMODE中输入4，单击OK。开始特征值屈曲分析：MainMenuSolutionSolve-CurrentLS，单击Ok开始求解。5)查看求解结果读入结果数据：MainMenuGeneralPostprocReadResultsFirstSet。列出屈曲载荷系数：MainMenuGeneralPostprocListResultsDetailedSummary。,2.4静力非线性分析,2.4.3线性屈曲分析实例5)查看求解结果观察屈曲变形和应力云图：MainMenuGeneralPostprocPlotResultsContourPlotNodalSolu，出现ContourNodalSolutionData对话框。选择NodalSolutionStressdisplacementvectorsum，在Undisplacedshapekey的下拉菜单中选择Deformedshapewithundeformededge，单击Ok，结构变形前后的形状同时出现在图形窗口，如图2.11所示。,2.4静力非线性分析,2.4.3线性屈曲分析实例5)查看求解结果,图2.11一阶屈曲模态的vonMises应力分布图,2.4静力非线性分析,2.4.4接触分析实例问题描述经典的赫兹接触问题：衬底材料的弹性模量为60GPa，泊松比为0.3，尺寸如图2.12所示；刚性压头的半径为1mm，压入载荷20KN。材料之间的摩擦系数为0.15。,图2.12赫兹接触模型,2.4静力非线性分析,2.4.4接触分析实例接触分析过程定义任务名称：启动ANSYS，选择材料FileChangeTitle.，输入“ContactAnalysis”，单击OK退出。设定参数变量：选择主菜单下ParametersScalarParameters.，弹出ScalarParameters对话框。在Section下分别定义Contact=1.5、Radius=1、Thick=2和Length=12，然后退出ScalarParameters对话框。,2.4静力非线性分析,2.4.4接触分析实例接触分析过程定义单元类型：选择MainMenuPreprocessorElementTypeAdd/Edit/Delete，分别添加Solid下的4node182单元、Contact下的2Dtarget169和2ndsurf171单元，点击OK退出该对话框。定义实常数：选择MainMenuPreprocessorRealConstants，选择Type2TARGE169，不修改任何数据，单击OK退出。单击Add.，选择Type3CONTA171，设定参数FKN为1，FTOLN为0.1，PMAX为0.1，PMIN为0.01，然后退出RealConstants对话框。,2.4静力非线性分析,2.4.4接触分析实例接触分析过程定义材料模型：选择MainMenuPreprocessorMaterialPropsMaterialModels，定义材料模型1为各向同性材料，EX=60e9，PRXY=0.3。回到DefineMaterialModelBehavior对话框，选择StructuralFrictionCoefficient，设定MU=0.15。建立几何模型：选择MainMenuPreprocessorModelingCreateKeypointsInActiveCS，依次分别建立坐标为(0,0,0)、(Contact,0,0)和(Length,0,0)的三个关键点。,2.4静力非线性分析,2.4.4接触分析实例接触分析过程建立几何模型：选择MainMenuPreprocessorModelingCopyKeypoints，弹出CopyKeypoints对话框，选定pickall。设定DY=thick，单击OK退出。选择MainMenuPreprocessorModelingCreateAreasArbitraryThroughKeypoints，依次单击1,2,5,4关键点，点击Apply，然后单击2,3,6,5，单击OK完成面的建立。,2.4静力非线性分析,2.4.4接触分析实例接触分析过程建立几何模型：选择UtilityMenuWorkPlaneOffsetWPtoXYZLocations+，在GlobalCartesian栏输入0，thick+radius，单击OK。选择UtilityMenuWorkPlaneLocalCoordianteSystemsCreateLocalCSAtWPOrigin.，弹出CreateLocalCSatWPOrigin对话框，在KCS栏中输入Cylindrical1，单击OK退出。选择MainMenuPreprocessorModelingCreateKeypointsInActiveCS，依次建立编号为20-23的四个关键点，其坐标分别为(radius,180,0)、(radius,-90,0)、(radius,0,0)和(0,0,0)。选择MainMenuPreprocessorModelingCreateLinesInActiveCoord，分别将关键点20,21以及21,22连成线。,2.4静力非线性分析,2.4.4接触分析实例接触分析过程划分有限元网格：选择MainMenuPreprocessorMeshingDefaultAttribs,选择1PLANE1