ANSYS实例分析内有结果

-.z.【ANSYS算例】3.4.2(1)基于图形界面的桁架桥梁结构分析(stepbystep)下面以一个简单桁架桥梁为例，以展示有限元分析的全过程。背景素材选自位于密执安的"OldNorthParkBridge"(1904-1988)，见图3-22。该桁架桥由型钢组成，顶梁及侧梁，桥身弦杆，底梁分别采用3种不同型号的型钢，结构参数见表3-6。桥长L=32m,桥高H=5.5m。桥身由8段桁架组成，每段长4m。该桥梁可以通行卡车，若这里仅考虑卡车位于桥梁中间位置，假设卡车的质量为4000kg，若取一半的模型，可以将卡车对桥梁的作用力简化为P1，P2和P3，其中P1=P3=5000N,P2=10000N，见图3-23。图3-22位于密执安的"OldNorthParkBridge"(1904-1988)图3-23桥梁的简化平面模型(取桥梁的一半)表3-6桥梁结构中各种构件的几何性能参数构件惯性矩m4横截面积m2顶梁及侧梁(Beam1)桥身弦梁(Beam2)底梁(Beam3)解答以下为基于ANSYS图形界面(GraphicUserInterface,GUI)的菜单操作流程。(1)进入ANSYS（设定工作目录和工作文件）程序→ANSYS→ANSYSInteractive→Workingdirectory（设置工作目录）→Initialjobname（设置工作文件名）:TrussBridge→Run→OK(2)设置计算类型ANSYSMainMenu：Preferences…→Structural→OK(3)定义单元类型ANSYSMainMenu：Preprocessor→ElementType→Add/Edit/Delete...→Add…→Beam:2delastic3→OK（返回到ElementTypes窗口）→Close(4)定义实常数以确定梁单元的截面参数ANSYSMainMenu:Preprocessor→RealConstants…→Add/Edit/Delete→Add…→selectType1Beam3→OK→inputRealConstantsSetNo.:1,AREA:2.19E-3，Izz:3.83e-6(1号实常数用于顶梁和侧梁)→Apply→inputRealConstantsSetNo.:2,AREA:1.185E-3，Izz:1.87E-6(2号实常数用于弦杆)→Apply→inputRealConstantsSetNo.:3,AREA:3.031E-3，Izz:8.47E-6(3号实常数用于底梁)→OK(backtoRealConstantswindow)→Close(theRealConstantswindow)(5)定义材料参数ANSYSMainMenu:Preprocessor→MaterialProps→MaterialModels→Structural→Linear→Elastic→Isotropic→inputE*:2.1e11,PR*Y:0.3(定义泊松比及弹性模量)→OK→Density(定义材料密度)→inputDENS:7800,→OK→Close（关闭材料定义窗口）(6)构造桁架桥模型生成桥体几何模型ANSYSMainMenu：Preprocessor→Modeling→Create→Keypoints→InActiveCS→NPTKeypointnumber：1，*，Y，ZLocationinactiveCS：0，0→Apply→同样输入其余15个特征点坐标（最左端为起始点，坐标分别为(4,0),(8,0),(12,0),(16,0),(20,0),(24,0),(28,0),(32,0),(4,5.5),(8,5.5),(12,5.5),(16.5.5),(20,5.5),(24,5.5),(28,5.5)）→Lines→Lines→StraightLine→依次分别连接特征点→OK网格划分ANSYSMainMenu:Preprocessor→Meshing→MeshAttributes→PickedLines→选择桥顶梁及侧梁→OK→selectREAL:1,TYPE:1→Apply→选择桥体弦杆→OK→selectREAL:2,TYPE:1→Apply→选择桥底梁→OK→selectREAL:3,TYPE:1→OK→ANSYSMainMenu：Preprocessor→Meshing→MeshTool→位于SizeControls下的Lines：Set→ElementSizeonPicked→Pickall→Apply→NDIV：1→OK→Mesh→Lines→Pickall→OK(划分网格)(7)模型加约束ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Displacement→OnNodes→选取桥身左端节点→OK→selectLab2:AllDOF(施加全部约束)→Apply→选取桥身右端节点→OK→selectLab2:UY(施加Y方向约束)→OK(8)施加载荷ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Force/Moment→OnKeypoints→选取底梁上卡车两侧关键点（*坐标为12及20）→OK→selectLab:FY，Value:-5000→Apply→选取底梁上卡车中部关键点（*坐标为16）→OK→selectLab:FY，Value:-10000→OK→ANSYSUtilityMenu：→Select→Everything(9)计算分析ANSYSMainMenu：Solution→Solve→CurrentLS→OK(10)结果显示ANSYSMainMenu：GeneralPostproc→PlotResults→Deformedshape→Defshapeonly→OK（返回到PlotResults）→ContourPlot→NodalSolu→DOFSolution,Y-ponentofDisplacement→OK（显示Y方向位移UY）(见图3-24(a))定义线性单元I节点的轴力ANSYSMainMenu→GeneralPostproc→ElementTable→DefineTable→Add→Lab:[bar_I],Bysequencenum:[SMISC,1]→OK→Close定义线性单元J节点的轴力ANSYSMainMenu→GeneralPostproc→ElementTable→DefineTable→Add→Lab:[bar_J],Bysequencenum:[SMISC,1]→OK→Close画出线性单元的受力图(见图3-24(b))ANSYSMainMenu→GeneralPostproc→PlotResults→ContourPlot→LineElemRes→LabI:[bar_I],LabJ:[bar_J],Fact:[1]→OK(11)退出系统ANSYSUtilityMenu：File→E*it→SaveEverything→OK(a)桥梁中部最大挠度值为0.003374m(b)桥梁中部轴力最大值为25380N图3.24桁架桥挠度UY以及单元轴力计算结果【ANSYS算例】3.4.2(2)基于命令流方式的桁架桥梁结构分析!%%%%%[ANSYS算例]3.4.2(2)%%%%%begin%%%%%%!注：命令流中的符号$，可将多行命令流写成一行/prep7!进入前处理/PLOPTS,DATE,0!设置不显示日期和时间!=====设置单元和材料ET,1,BEAM3!定义单元类型R,1,2.19E-3,3.83e-6,,,,, !定义1号实常数用于顶梁侧梁R,2,1.185E-3,1.87e-6,0,0,0,0, !定义2号实常数用于弦杆R,3,3.031E-3,8.47E-6,0,0,0,0, !定义3号实常数用于底梁MP,E*,1,2.1E11!定义材料弹性模量MP,PR*Y,1,0.30!定义材料泊松比MP,DENS,1,,7800!定义材料密度!定义几何关键点K,1,0,0,,$K,2,4,0,,$K,3,8,0,,$K,4,12,0,,$K,5,16,0,,$K,6,20,0,,$K,7,24,0,,$K,8,28,0,,$K,9,32,0,,$K,10,4,5.5,,$K,11,8,5.5,,$K,12,12,5.5,,$K,13,16,5.5,,$K,14,20,5.5,,$K,15,24,5.5,,$K,16,28,5.5,,!通过几何点生成桥底梁的线L,1,2$L,2,3$L,3,4$L,4,5$L,5,6$L,6,7$L,7,8$L,8,9!生成桥顶梁和侧梁的线L,9,16$L,15,16$L,14,15$L,13,14$L,12,13$L,11,12$L,10,11$L,1,10!生成桥身弦杆的线L,2,10$L,3,10$L,3,11$L,4,11$L,4,12$L,4,13$L,5,13$L,6,13$L,6,14$L,6,15$L,7,15$L,7,16$L,8,16!选择桥顶梁和侧梁指定单元属性LSEL,S,,,9,16,1,LATT,1,1,1,,,,!选择桥身弦杆指定单元属性LSEL,S,,,17,29,1,LATT,1,2,1,,,,!选择桥底梁指定单元属性LSEL,S,,,1,8,1,LATT,1,3,1,,,,!划分网格AllSEL,all!再恢复选择所有对象LESIZE,all,,,1,,,,,1!对所有对象进行单元划分前的分段设置LMESH,all!对所有几何线进行单元划分!=====在求解模块中，施加位移约束、外力，进行求解/soluNSEL,S,LOC,*,0!根据几何位置选择节点D,all,,,,,,ALL,,,,,!对所选择的节点施加位移约束AllSEL,all!再恢复选择所有对象NSEL,S,LOC,*,32!根据几何位置选择节点D,all,,,,,,,UY,,,,!对所选择的节点施加位移约束ALLSEL,all!再恢复选择所有对象!基于几何关键点施加载荷FK,4,FY,-5000$FK,6,FY,-5000$FK,5,FY,-10000/replot!重画图形Allsel,all!选择所有信息(包括所有节点、单元和载荷等)solve!求解!=====进入一般的后处理模块/post1!后处理PLNSOL,U,Y,0,1.0!显示Y方向位移PLNSOL,U,*,0,1.0!显示*方向位移!显示线单元轴力ETABLE,bar_I,SMISC,1ETABLE,bar_J,SMISC,1PLLS,BAR_I,BAR_J,0.5,1!画出轴力图finish!结束!%%%%%[ANSYS算例]3.4.2(2)%%%%%end%%%%%%【ANSYS算例】3.2.5(3)四杆桁架结构的有限元分析下面针对【典型例题】3.2.5(1)的问题，在ANSYS平台上，完成相应的力学分析。即如图3-8所示的结构，各杆的弹性模量和横截面积都为,，基于ANSYS平台，求解该结构的节点位移、单元应力以及支反力。图3-8四杆桁架结构解答对该问题进行有限元分析的过程如下。以下为基于ANSYS图形界面(graphicuserinterface，GUI)的菜单操作流程；注意：符号“→”表示针对菜单中选项的鼠标点击操作。关于ANSYS的操作方式见附录B。1．基于图形界面的交互式操作(stepbystep)(1)进入ANSYS(设定工作目录和工作文件)程序→Ansys→ANSYSInteractive→Workingdirectory(设置工作目录)→Initialjobname(设置工作文件名):planetruss→Run→OK(2)设置计算类型ANSYSMainMenu:Preferences…→Structural→OK(3)选择单元类型ANSYSMainMenu:Preprocessor→ElementType→Add/Edit/Delete…→Add…→Link：2Dspar1→OK(返回到ElementTypes窗口)→Close(4)定义材料参数ANSYSMainMenu:Preprocessor→MaterialProps→MaterialModels→Structural→Linear→Elastic→Isotropic:E*:2.95e11(弹性模量)，PR*Y:0(泊松比)→OK→鼠标点击该窗口右上角的“”来关闭该窗口(5)定义实常数以确定单元的截面积ANSYSMainMenu:Preprocessor→RealConstants…→Add/Edit/Delete→Add→Type1→OK→RealConstantSetNo:1(第1号实常数),AREA:1e-4(单元的截面积)→OK→Close(6)生成单元ANSYSMainMenu:Preprocessor→Modeling→Creat→Nodes→InActiveCS→Nodenumber1→*:0,Y:0,Z:0→Apply→Nodenumber2→*:0.4,Y:0,Z:0→Apply→Nodenumber3→*:0.4,Y:0.3,Z:0→Apply→Nodenumber4→*:0,Y:0.3,Z:0→OKANSYSMainMenu:Preprocessor→Modeling→Create→Elements→ElemAttributes（接受默认值）→Usernumbered→Thrunodes→OK→选择节点1，2→Apply→选择节点2，3→Apply→选择节点1，3→Apply→选择节点3，4→Apply→OK(7)模型施加约束和外载添加位移的约束，分别将节点1*和Y方向、节点2Y方向、节点4的*和Y方向位移约束。ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Displacement→OnNodes→用鼠标选择节点1→Apply→Lab2DOFs:U*，UY，VALUE：0→Apply→用鼠标选择节点2→Apply→Lab2DOFs:UY，VALUE：0→Apply→用鼠标选择节点4→Apply→Lab2DOFs:U*,UY，VALUE：0→OK加载集中力ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Force/moment→OnNodes→用鼠标选择结构节点2→Apply→F*，VALUE:20000→Apply→用鼠标选择结构节点3→Apply→FY，VALUE:－25000→OK(9)分析计算ANSYSMainMenu:Solution→Solve→CurrentLS→OK→ShouldTheSolvemandbeE*ecuted"Y→Close(Solutionisdone!)→关闭文字窗口(10)结果显示ANSYSMainMenu:GeneralPostproc→PlotResults→DeformedShape…→Def+Undeformed→OK(返回到PlotResults)→ContourPlot→NodalSolu…→DOFsolution→Displacementvectorsum(可以看到位移云图)ANSYSMainMenu:GeneralPostproc→ListResults→Nodalsolution→DOFsolution→Displacementvectorsum(弹出的文本文件显示各个节点的位移)ANSYSMainMenu:GeneralPostproc→ListResults→ReactionSolu→ALLitems→OK(弹出的文本文件显示各个节点反力)(11)退出系统ANSYSUtilityMenu:File→E*it…→SaveEverything→OK2．完整的命令流以下为命令流语句。注意：以“!”打头的文字为注释容，其后的文字和符号不起运行作用。关于命令流的调用方式见附录B。!%%%%%%%%[典型例题]3.2.5(1)%%%%begin%%%%%%/PREP7!进入前处理/PLOPTS,DATE,0!设置不显示日期和时间!=====设置单元、材料，生成节点及单元ET,1,LINK1!选择单元类型UIMP,1,E*,,,2.95e11,!给出材料的弹性模量R,1,1e-4,!给出实常数(横截面积)N,1,0,0,0,!生成1号节点,坐标(0,0,0)N,2,0.4,0,0,!生成2号节点,坐标(0.4,0,0)N,3,0.4,0.3,0,!生成3号节点,坐标(0.4,0.3,0)N,4,0,0.3,0,!生成4号节点,坐标(0,0.3,0)E,1,2!生成1号单元(连接1号节点和2号节点)E,2,3!生成2号单元(连接2号节点和3号节点)E,1,3!生成3号单元(连接1号节点和3号节点)E,4,3!生成4号单元(连接4号节点和3号节点)FINISH!前处理结束!=====在求解模块中，施加位移约束、外力，进行求解/SOLU!进入求解状态(在该状态可以施加约束及外力)D,1,ALL!将1号节点的位移全部固定D,2,UY，!将2号节点的Y方向位移固定D,4,ALL!将4号节点的位移全部固定F,2,F*,20000,!在2号节点处施加*方向的力(20000)F,3,FY,-25000,!在3号节点处施加Y方向的力(-25000)SOLVE!进行求解FINISH!结束求解状态!=====进入一般的后处理模块/POST1!进入后处理PLDISP,1!显示变形状况FINISH!结束后处理!%%%%%%%%[典型例题]3.2.5(1)%%%%end%%%%%%【ANSYS算例】3.3.7(3)三梁平面框架结构的有限元分析针对【典型例题】3.3.7(1)的模型，即如图3-19所示的框架结构，其顶端受均布力作用，用有限元方法分析该结构的位移。结构中各个截面的参数都为：，，，相应的有限元分析模型见图3-20。在ANSYS平台上，完成相应的力学分析。图3-19框架结构受一均布力作用（a）节点位移及单元编号（b）等效在节点上的外力图3-20单元划分、节点位移及节点上的外载解答对该问题进行有限元分析的过程如下。1．基于图形界面的交互式操作(stepbystep)(1)进入ANSYS(设定工作目录和工作文件)程序→Ansys→ANSYSInteractive→Workingdirectory(设置工作目录)→Initialjobname(设置工作文件名):beam3→Run→OK(2)设置计算类型ANSYSMainMenu:Preferences…→Structural→OK(3)选择单元类型ANSYSMainMenu:Preprocessor→ElementType→Add/Edit/Delete…→Add…→beam：2Delastic3→OK(返回到ElementTypes窗口)→Close(4)定义材料参数ANSYSMainMenu:Preprocessor→MaterialProps→MaterialModels→Structural→Linear→Elastic→Isotropic:E*:3e11(弹性模量)→OK→鼠标点击该窗口右上角的“”来关闭该窗口(5)定义实常数以确定平面问题的厚度ANSYSMainMenu:Preprocessor→RealConstants…→Add/Edit/Delete→Add→Type1Beam3→OK→RealConstantSetNo:1(第1号实常数),Cross-sectionalarea:6.8e-4(梁的横截面积)→OK→Close(6)生成几何模型生成节点ANSYSMainMenu:Preprocessor→Modeling→Creat→Nodes→InActiveCS→Nodenumber1→*:0,Y:0.96,Z:0→Apply→Nodenumber2→*:1.44,Y:0.96,Z:0→Apply→Nodenumber3→*:0,Y:0,Z:0→Apply→Nodenumber4→*:1.44,Y:0,Z:0→OK生成单元ANSYSMainMenu:Preprocessor→Modeling→Create→Element→AutoNumbered→ThruNodes→选择节点1，2(生成单元1)→apply→选择节点1，3(生成单元2)→apply→选择节点2，4(生成单元3)→OK(7)模型施加约束和外载左边加*方向的受力ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Force/Moment→OnNodes→选择节点1→apply→Directionofforce:F*→VALUE：3000→OK→上方施加Y方向的均布载荷ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Pressure→OnBeams→选取单元1(节点1和节点2之间)→apply→VALI：4167→VALJ：4167→OK左、右下角节点加约束ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Displacement→OnNodes→选取节点3和节点4→Apply→Lab:ALLDOF→OK(8)分析计算ANSYSMainMenu:Solution→Solve→CurrentLS→OK→ShouldtheSolvemandbeE*ecuted"Y→Close(Solutionisdone!)→关闭文字窗口(9)结果显示ANSYSMainMenu:GeneralPostproc→PlotResults→DeformedShape…→Def+Undeformed→OK(返回到PlotResults)(10)退出系统ANSYSUtilityMenu:File→E*it…→SaveEverything→OK(11)计算结果的验证与MATLAB支反力计算结果一致。2．完全的命令流!%%%%%%%%%%[典型例题]3.3.7(3)%%%begin%%%%%/PREP7!进入前处理ET,1,beam3!选择单元类型R,1,6.5e-7,6.8e-4!给出实常数(横截面积、惯性矩)MP,E*,1,3e11!给出材料的弹性模量N,1,0,0.96,0!生成4个节点,坐标(0,0.96,0)，以下类似N,2,1.44,0.96,0N,3,0,0,0N,4,1.44,0,0E,1,2!生成单元(连接1号节点和2号节点)，以下类似E,1,3E,2,4D,3,ALL!将3号节点的位移全部固定D,4,ALL!将4号节点的位移全部固定F,1,F*,3000!在1号节点处施加*方向的力(3000)SFBEAM,1,1,PRESS,4167!施加均布压力FINISH!结束前处理状态/SOLU!进入求解模块SOLVE!求解FINISH!结束求解状态/POST1!进入后处理PLDISP,1!显示变形状况FINISH!结束后处理!%%%%%%%%%%[典型例题]3.3.7(3)%%%end%%%%%实验一衍架的结构静力分析一、问题描述图1所示为由9个杆件组成的衍架结构，两端分别在1，4点用铰链支承，3点受到一个方向向下的力Fy,衍架的尺寸已在图中标出，单位：m。试计算各杆件的受力。图1衍架结构简图图1衍架结构简图 弹性模量（也称扬式模量）E=206GPa；泊松比μ=0.3； 作用力Fy=-1000N；杆件的横截面积A=0.125m2.显然，该问题属于典型的衍架静力分析问题，通过理论求解方法（如节点法或截面法）也可以很容易求出个杆件的受力，但这里为什么要用ANSYS软件对其分析呢？三、结果演示通过使用ANSYS8.0软件对该衍架结构进行静力分析，其分析结果与理论计算结果如表1所示。表1ANSYS分析结果与理论计算结果的比较杆件序号123456789ANSYS分析结果/N333.33333.33666.67-471.400-666.67471.40666.67-947.81理论计算结果/N333.333333.333666.667-471.4050-666.667471.405666.667-942.809误差/%0.30.30.30.500.30.50.30.1比较结果表明，使用ANSYS分析的结果与理论计算结果的误差不超过0.5%，因此，利用ANSYS软件分析来替代理论计算是完全可行的。四、实训步骤（一）ANSYS8.0的启动与设置1． 启动。点击：开始>所有程序>ANSYS8.0>ANSYS，即可进入ANSYS图形用户主界面。如图2所示。其中，几个常用的部分有应用菜单，命令输入栏，主菜单，图形显示区和显示调整工具栏，分别如图2所示。图3Preference参数设置对话框2． 功能设置。电击主菜单中的“Preference”菜单，弹出“参数设置”对话框，选中“Structural”复选框，点击“OK”按钮，关闭对话框，如图3所示。本步骤的目的是为了仅使用该软件的结构分析功能，以简化主菜单中各级子菜单的结构。图3Preference参数设置对话框图4单元类型库对话框图5单元类型对话框3． 系统单位设置。由于ANSYS软件系统默认的单位为英制，因此，在分析之前，应将其设置成国际公制单位。在命令输入栏中键入“/UNITS，SI”图4单元类型库对话框图5单元类型对话框（二）单元类型，几何特性及材料特性定义图6单元类型对话框图7实常数对话框图8材料特性对话框1．定义单元类型。电击主菜单中的“Preference>ElementType>Add/Edit/Delete”，弹出对话框，点击对话框中的“Add图6单元类型对话框图7实常数对话框图8材料特性对话框的“Link”和“2Dspar1”选项，点击“OK”，关闭图4对话框，返回至上一级对话框，此时，对话框中出现刚才选中的单元类型：LINK1，如图5所示。点击“Close”，关闭图5所示对话框。注：LINK1属于二维平面杆单元，即我们常说的二力杆，只承受拉压，不考虑弯矩。2．定义几何特性。在ANSYS中主要是实常数的定义：点击主菜单中的“Preprocessor>RealContants>Add/Edit/Delete”,弹出对话框，点击“Add…”按钮，第二（1）步定义的LINK1单元出现于该对话框中，点击“OK”，弹出下一级对话框，如图6所示。在AREA一栏杆件的截面积0.125，点击“OK”，回到上一级对话框，如图7所示。点击“Close”，关闭图7所示对话框。3．定义材料特性。点击主菜单中的“Preprocessor>MaterialProps>MaterialModels”,弹出对话框，如图8所示,逐级双击右框中“Structural,Linear,Elastic,Isotropic”前图标，弹出下一级对话框，在弹性模量文本框中输入：206E9，在泊松比文本框中输入：0.3，如图9所示，点击“OK”返回上一级对话框，并点击“关闭”按钮，关闭图8所示对话框。图9材料特性参数对话框（三）衍架分析模型的建立图9材料特性参数对话框图11生成节点显示图10节点生成参数输入对话框1． 生成节点。图10所示衍架中共有6个节点，其坐标根据已知条件容易求出如下：1（0，0，0），2（1，0，0），3（2，0，0），4（3，0，0），5（1，1，0），6（2，1，0）。点击主菜单中的“Preprocessor>Modeling>Create>Nodes>InActiveCS”,弹出对话框.在“Nodenumber”一栏中输入节点号1，在“*YZLocation”一栏中输入节点1的坐标（0，0，0），如图10所示，点击“Apply”图11生成节点显示图10节点生成参数输入对话框按钮，在生成1节点的同时弹出与图10一样的对话框，同理将2-6点的坐标输入，以生成其余5个节点。此时，在显示窗口上显示所生成的6个节点的位置，如图11所示。图12节点选择对话框2．生成单元格。点击主菜单中的“Preprocessor>Modeling>Create>Elements>AutoNumbered>ThruNodes”,弹出“节点选择”对话框，如图12所示。依次点选节点1、2，点击Apply按钮，既可生成①单元。同理，分别点击2、3；3、4；1、5；2、5；5、6；3、5；3、6；4、6可生成其余8个单元。生成后的单元如图13所示。图12节点选择对话框（四）施加载荷图13生成单元显示1．施加位移约束。点击主菜单中的“Preprocessor>Loads>Apply>Structural>Displacement>OnNodes”,弹出与图12所示类似的“节点选择”对话框，点选1节后，然后点击“Apply”按钮，弹出对话框如图14所示，选择右上列表框中的“AllDOF”,并点击“Apply”按钮，弹出对话框如图14所示，选择右上列表框中的UY，并点击“OK”按钮，即可完成对节点4沿y方向的位移约束。图13生成单元显示图14节点1的位移约束2．施加集中力载荷。点击主菜单中的“Preprocessor>Loads>DefineLoads>Apply>Structural>Force/Moment>OnNodes”,弹出对话框如图15所示，在“Directionofforce/mom”一项中选择：“FY”，在“Force/Momentvalue”一项中输入：-1000（注：负号表示力的方向与Y的正向相反），然后点击“OK”按钮关闭对话框，这样，就在节点3处给桁架结构施加了一个竖直向下的集中载荷。图14节点1的位移约束图15施加载荷说明：根据图1所示有限元分析的基本过程，到此为止，有限元分析的前置处理部分已经结束。但在使用ANSYS软件进行分析的过程中，施加载荷这一步骤往往既可以在前置处理中完成（如本实训一样），也可以在求解器中完成（如点击主菜单中的“Solution>DefineLoads>Apply>Stuctural…”,实现过程完全一样）。图15施加载荷（五）开始求解图16求解对话框点击主菜单中的“Solution>Solve>CurrentLS”,弹出对话框（图16），点击“OK”按钮，开始进行分析求解。分析完成后，又弹出一信息窗口（图17）提示用户已完成求解，点击“Close”按钮关闭对话框即可。至于在求解时产生的STATUSmand窗口，点击“File>Close”关闭即可。图16求解对话框图17求解完成说明：到此为止，有限元分析的求解器计算部分已经结束。图17求解完成（六）分析结果显示图18显示变形图设置图2-3用户主界面1．显示变形图。点击主菜单中的“GeneralPostproc>PlotResults>DeformedShape”,弹出对话框如图18所示。选中“Def+undeformed”选项，并点击“OK”按钮，即可显示本实训桁架结构变形前后的结果，如图19所示。图18显示变形图设置图2-3用户主界面图20变形动画参数设置图19用户主界面2显示变形动画。点击应用菜单（UtilityMenu）中的PlotCtrls>Animate>DeformedShape…,弹出对话框如图20所示。选中Def+undeformed”选项，并在“Timedelay”文本框中输入：0.1,然后点击“OK”按钮，即可显示本实训桁架结构的变性动画。由于集中力FY作用在3节点上，因此，3节点产生的位移最大。图21是动画片、显示桁架受力变形的过程，右边窗口是动画显示的控制窗口，可以暂停，也可以拖动显示进度条。图20变形动画参数设置图19用户主界面图22显示支反力参数设置对话框图21动画仿真控制对话框3．列举支反力计算结果。点击主菜单中的“GeneralPostproc>ListResults>ReactionSolu”，弹出对话框如图22所示。接受缺省设置，点击“OK”按钮关闭对话框，并弹出一列表窗口，显示了两铰链点（1、4节点）所受的支反力情况，如图23所示。图22显示支反力参数设置对话框图21动画仿真控制对话框4.列举各杆件的轴向力计算结果。点击主菜单中的“GeneralPostproc>ListResult>ElementSolution”，弹出对话框如图24所示，在中间列表框中移动滚动条至最后，选择“BySequencenum”选项，右上列表框中选择“SMISC”选项，右下文本框中输入“SMISC，1”，点击“OK”按钮关闭对话框，并弹出一列表窗口，显示了9个杆单元所受的轴向力，如图25所示，此外，还给出了最大、最小力及其发生位置。说明：到此为止，有限元分析的后置处理部分就可以结束了。实际上，ANSYS软件的后置处理功能非常强大，除了能显示上述结果外，还可以显示其他许多结果，比如，各节点产图25显示轴向力列表图23显示支反力列表生的位移（GeneralPostproc>ListResults>NodalSolution）和各节点所受载荷（GeneralPostproc>ListResults>NodalLoads）等等，有兴趣的读者可自行尝试。图25显示轴向力列表图23显示支反力列表（七）ANSYS软件的保存与退出1．保存。点击应用菜单中的“File>SaveasJobname.db,ANSYS”自动将结果保存，缺省的文件名为：下次打开时，可直接点击“file.db”2．退出。点击应用菜单中的“File>E*it…”，弹出保存对话框，选中“SaveEverything”，点击“OK”按扭，即可退出ANSYS。图24单元求解结果参数选择对话框图24单元求解结果参数选择对话框实验二三维实体结构的分析介绍三维实体结构的有限元分析。一、问题描述图25所示为一工字钢梁，两端均为固定端，其截面尺寸为，。试建立该工字钢梁的三维实体模型，并在考虑重力的情况下对其进行结构静力分析。其他已知参数如下：图25工字钢结构示意图弹性模量（也称式模量）E=206GPa；泊松比；图25工字钢结构示意图材料密度；重力加速度；作用力Fy作用于梁的上表面沿长度方向中线处，为分布力，其大小Fy=-5000N三、结果演示图27单元类型对话框图26单元类型库对话框使用ASSYS8。0软件对该工字钢梁进行结构静力分析，显示其节点位移云图。图27单元类型对话框图26单元类型库对话框四、实训步骤（一）ASSYS8.0的启动与设置与实训1第一步骤完全相同，请参考。（二）单元类型、几何特性及材料特性定义图28材料特性参数对话框1定义单元类型。点击主菜单中的“Preprocessor>ElementType>Add/Edit/Delete”，弹出对话框，点击对话框中的“Add…”按钮，又弹出一对话框（图26），选中该对话框中的“Solid”和“Brick8node45”选项，点击“OK”，关闭图26对话框，返回至上一级对话框，此时，对话框中出现刚才选中的单元类型：Solid45，如图27所示。点击“Close”，关闭图27所示对话框。注：Solid45单元用于建立三维实体结构的有限元分析模型，该单元由8个节点组成，每个节点具有*、Y、Z方向的三个移动自由度。图28材料特性参数对话框图28材料特性对话框2．定义材料特性。点击主菜单中的“Preprocessor>MaterialProps>MaterialModels”，弹出窗口如图28所示，逐级双击右框中“Structural\Linear\Elastic\Isotropic”前图标，弹出下一级对话框，在“弹性模量”（E*）文本框中输入：2.06e11，在“泊松比”（PR*Y）文本框中输入：0.3，如图29所示，点击“OK”按钮，回到上一级对话框，然后，双击右框中的“Density”选项，在弹出对话框的“DENS”一栏中输入材料密度：7800，点击“OK”按钮关闭对话框。最后，点击图2-31所示窗口右上角“关闭”该窗口。图28材料特性对话框（三）工字钢三维实体模型的建立图30节点生成参数输入对话框1．生成关键点。图30所示的工字钢梁的横截面由12个关键点连线而成，其各点坐标分别为：1（-0.08,0,0）、2（0.08,0,0）、3（0.08,0.02,0）、4（0.015,0.02,0）、5（0.015,0.18,0）、6（0.08,0.18,0）、7（0.08,0.20,0）、8（-0.08,0.20,0）、9（-0.08,0.18,0）、10（-0.015,0.18,0）、11（-0.015,0.02,0）、12（-0.08,0.02,0）。点击主菜单中的“Preprocessor>Modeling>Create>图30节点生成参数输入对话框Keypoints>InActiveCS”，弹出对话框。在“Keypointnumber”一栏中输入关键点号1，在“*YZLocation”一栏中输入关键点1的坐标（-0.08，0，0），如图30所示，点击“Apply”按钮，同理将2~12点的坐标输入，此时，在显示窗口上显示所生成的12个关键点的位置。图31生成直线显示图32面选择对话框2．生成直线。点击主菜单中的“Preprocessor>Modeling>Create>Lines>Lines>StraightLine”，弹出关键点选择对话框，依次点选关键点1、2，点击“Apply”图31生成直线显示图32面选择对话框图33平面拉伸成体的参数设置3．生成平面。点击主菜单中的“Preprocessor>Modeling>Create>Areas>Arbitrary>ByLines”，弹出“直线选择”对话框，依次点选1~12直线，点击“OK”按钮关闭对话框，即可生成工字钢的横截面。图33平面拉伸成体的参数设置4．生成三维实体。点击主菜单中的“Preprocessor>Modeling>>Operate>E*trude>Areas>AlongNormal”，弹出平面选择对话框32，点选上一步骤生成的平面，点击“OK”按钮。之后弹出另一对话框33，在“DIST”一栏中输入：1（工字钢梁的长度），其他保留缺省设置，点击OK按钮关闭对话框，即可生成工字钢梁的三维实体模型。如图34所示。图34工字钢梁三维实体模型（四）网络划分图34工字钢梁三维实体模型图36划分网格后的工字钢梁模型1．设定单元大小。点击主菜单中的“Preprocessor>Meshing>MeshTool”,弹出对话框，在“SizeControl”标签中的Global一栏点击Set按钮，弹出“网格尺寸设置”对话框，在SIZE一栏中输入：0.02，其他保留缺省设置，点击OK按钮关闭对话框。图36划分网格后的工字钢梁模型2．接着上一步，在图35的划分网格的对话框中，选中单选框“He*”和“Sweep”，其他保留缺省设置，然后点击“Sweep”按钮，弹出体选择对话框，点选34中的工字钢梁实体，并点击OK按钮，即可完成对整个实体结构的网格划分，其结果如36所示。（五）施加载荷1．施加位移约束。点击主菜单的“Preprecessor>Loads>DefineLoads>Apply>Structuaral>Displacemengt>OnAreas”，弹出面选择对话框，点击该工字梁的左端面，点击“OK”按钮，弹出对话框如图示，选择右上列表框中的“AllDOF”，并点击“OK”按钮，即可完成对左端面的位移约束，相当于梁的固定端。同理，对工字钢的右端面进行固定端约束。2．施加分布力（Fy）载荷。（1）选择施力节点。点击应用菜单中的“Select>Entities...”，弹出对话框如图37所示，在第一个列表框中选择“Nodes”选项，第二个列表框中选择“ByLocation”选项，选中“Zcoordinates”单选框，并在“Min,Ma*”参数的文本框中输入：0.5（表示选择工字钢梁沿的中间横截面上的所有节点），其他参数保留缺省设置，点击“Apply”按钮完成选择。点击“Plot”按钮，在显示窗口上显示出工字钢梁中间横截面上的所有节点。然后，在图37所示对话框中选中“Ycoordinates”单选框，在“Min,Ma*”参数文本框中输入：0.2（表示工字钢梁的上表面），选中“Reselect”（表示在现有活动节点——即上述选择的中间横截面中，再选择y坐标等于0.2的节点为活动节点）单选框，其他参数保留缺省设置（参见图37），然后依次点击“Apply”和“Plot”按钮，即可在显示窗口上显示出工字钢梁上表面沿长度方向中线处的一组节点，这组节点即为施力节点。图35网格划分对话框（2）施加载荷。点击主菜单中的“Preprocessor>Loads>DefineLoads>Apply>Structural>Force/Moment>OnLoads”，弹出“节点选择”对话框，点击“PickAll”按钮，即可选中（1）中所选择的这组需要施力的节点，之后弹出另一个对话框，在该对话框中的“Directionofforce/mom”一项中选择：“FY,在Force/momentvalue”一项中输入：-5000（注：负号表示力的方向与Y的方向相反），其他保留缺省装置，如图38所示，然后点击“OK”按钮关闭对话框，这样，通过在该组节点上施加与Y向相反的作用力，就可以模拟该实训中所要求的分布力Fy=-5000N。图35网格划分对话框（3）恢复选择所有节点。在求解之前必须选择所有已创建的对象为活动对象（如点、线、面、体、单元等），否则求解会出错。因此，点击应用菜单中的“Select>Everything”，即可完成该项工作。图37节点选择需要注意的是，此时显示窗口仅显示施力节点及作用力的方向箭头。若要显示整个工字钢梁的网络模型，可点击应用菜单中的“Plot>Elements”即可。图37节点选择3．施加重力载荷。点击主菜单中的“Preprocessor>Loads>DefineLoads>Apply>Structural>Inertia>Gravity”,在弹出对话框的“ACELY”一栏中输入：9.8（表示沿Y方向的重力加速度为9.8m/s，系统会自动利用密度等参数进行分析计算），其他保留缺省设置，点击“OK”关闭对话框。到此为止，有限元分析的前置处理部分已经结束。（六）求解点击主菜单中的“Solution>Solve>CurrentLS”，在弹出对话框中点击“OK”按钮，开始进行分析求解。分析完成后，又弹出一信息窗口提示用户已完成求解，点击“Close”按钮关闭对话框即可。至于在求解时产生的STATUSmand窗口，点击“File>Close”关闭即可。说明：到此为止，有限元分析的求解部分已经结束。（七）分析结果浏览1．绘制节点位移云图。点击主菜单中的“GeneralPostproc>PlotResults>ContourPlot>NodalSolu”，弹出对话框，选中右上列表框“Translation”栏中的“UY”选项，其他保留缺省设置。点击“OK”按钮，即可显示本实训工字钢梁各节点在重力和Fy作用下的位移云图，如图40所示。同理，通过在图39所示对话框中选择不同的选项，也可以绘制各节点的应力以及沿其他方向的云图。图38施加载荷时的参数设置图39节点位移显示设置2．列举各节点的位移解。点击主菜单中的GeneralPostproc>PlotResults>ContourPlot>NodalSolu”，弹出对话框如图39所示，全部保留缺省设置，点击“OK”按钮关闭对话框，并弹出一列表窗口，显示了该工字钢梁各节点的位移情况，显然，由于受力方向为y方向，因此，从窗口数据看出，各节点沿y的位移最大。图38施加载荷时的参数设置图39节点位移显示设置3．显示变形动画。点击应用菜单(Utility>Menu)中的“PlotCtrls>Animate>DeformedResults...)，在弹出的对话框中的“Timedelay”文本框中输入：0.1，并选中右列表框中的“UY”选项，其他保留缺省设置，点击“OK”按钮关闭对话框，即可显示本实训工字钢梁的变形动画。由于分布力Fy作用于梁中间，可以看出Fy对梁的局部作用过程。图40节点位移云图说明：到此为止，有限元分析的后置部分已经结束。与实验一一样，本实验还可以显示其他许多结果。例如，各单元沿各方向产生的应力、应变等容，有兴趣的读者可自行尝试。图40节点位移云图（八）ANSYS软件的保存与退出实验三变截面高速轴的最优化设计图41轴的简图一、 问题描述图41轴的简图 在图41所示的变截面轴上，安装有一个质量为Q的轮子。轴的长度为3,各段长度已知，要求解决的问题是确定在满足动力稳定性条件下，轴的质量最小时的直径和。 该问题的分析过程如下：选择设计变量。此最优化问题有两个设计变量，即要求确定的轴直径和。建立目标函数。以轴的质量作为目标函数，选用材料的密度为，则目标函数如下： 3.确定约束条件。当轴的旋转速度达到其临界转速时的角速度或横向振动（弯曲振动）的固有频率时，轴便处于共振状态。在大多数情况下，需要进行这种动力稳定性设计的轴，其质量总小于轮子的质量Q，为了简化计算，在确定时时常忽略，而简化为单自由度的振动问题。因此，轴的横向振动的固有频率为：式中g为重力加速度；u为轴的中间截面处的静挠度。 按照图41给出的条件，根据材料力学可求得： 式中E为材料的弹性模量。 为保证轴在工作时的动力稳定性，应使 式中K为大于1的安全系数。 将u的表达式及代入上式，经过整理后得到动力稳定性所要求的等式约束条件为： 另外，根据结构尺寸要求的条件为： 因此，装有一个重轮的变截面高速旋转直轴，当以轴的质量为目标函数时，其最优化设计的数学模型为：图43LINGO程序主界面二、实训目的图43LINGO程序主界面 理解“最优化设计”的概念和方法；了解“最优化设计”的建模方法；了解“最优化设计”问题的解决过程；了解LINGO的编程、方法和求解过程。图42求解结果三、结果演示图42求解结果图42为此问题的求解结果。四、实验步骤 选择材料密度，材料弹性模量，应用LINGO程序对问题进行最优化求解，具体操作过程如下：运行LINGO程序。LINGO程序主界面如图43所示。2.点击下拉菜单“File”，选择“New”或者单击工具栏中的按钮，新建一个优化实例模型LINGOModel(*.lg4).3.在LINGOModel窗口建立模型，输入以下程序语句：（注意：每一句后面要以分号“；”作为语句的结束。冒号“：”和分号“；”要在英文状态下输入。）model:L=0.4;Q=2;!轮子的质量p=7.8*10^3;!材料密度;Mg=(3.14*L*p/4)*(2*d1^2+d2^2);!轴的质量;min=Mg;!目标函数;!直径约束条件如下;图44求解状态窗口d11=0.01;图44求解状态窗口d12=0.4;d1>=d11;d1<=d12;!直径约束条件如下;d21=0.01;d22=0.4;d2>=d21;d2<=d22;E=200*10^9;!料性弹材模量;g=9.8;!重力加速度;w=3000;!轴的转速;k=1.3;!安全系数;(3.14*E*g)/(10.67*Q*l^3*w^2*k^2)-((1/d1^4)+(2.38/d2^4))=0;!约束条件;END4.模型建立完成后要对模型进行求解：单击下拉菜单“LINGO”选择“Solve”，或者单击工具栏中的模型求解命令按钮。LINGO程序会弹出“LINGOSolveStatus”窗口和“SolutionReport”窗口。在“LINGOSolveStatus”窗口中列出了求解器的有关信息，如图44所示，在“SolutionReport”窗口中列出了完整的模型求解过程，如图42所示。实验四压杆的最优化设计一、问题描述 图45所示的空心压杆两端受轴向外载荷P。轴的径为，外径为，支承间距尺寸为。试确定压杆的结构尺寸、和，以保证在压杆不产生屈服并且不破坏压杆稳定性条件下，压杆的体积和重量最小。 该问题的分析过程如下： 压杆为细长直杆，承受轴向压力，会因轴向压力达到临界值时突然弯曲而失去稳定性。设计压杆，除应使其压力不超过材料的弹性极限外，还必须使其承受的轴向压力小于压杆的临界载荷。 压杆在机械装置中应用的例子较多，例如在液压机构中当活塞的行程足够大时，会导致活塞杆足够长，这种细长的活塞杆便是压杆。 根据欧拉压杆公式，对与两端均为铰支的压杆，其临界载荷为图45压杆机构简图图45压杆机构简图式中E为压杆材料的弹性模量；J为压杆横截面的最小惯性矩，EJ为抗弯刚度；L为压杆长度。 将欧拉公式推广到端部不同约束的压杆，则上式变为式中为长度折算系数，其值将随压杆两端约束形式的不同而异。当两端铰支时取；一端固定，另一端自由时取；一端固定，另一端铰支时取；两端均固定时，取。 由欧拉公式可知，与成正比。合理的设计压杆截面形状，使其材料尽量远离形心分布，就能使J增大而提高压更的抗弯刚度EJ，增大临界载荷。所以在相同截面面积的条件下，管状压杆比实心压杆有更大的临界载荷。 以管状压杆的径，外径，长度作为设计变量，以其体积或重量作为目标函数，以压杆不产生屈服和不破坏轴向稳定性以及其尺寸约束条件，则管状压杆最优化设计的数学模型为式中分别为大于1的安全系数；P为设计给定的外载荷。二、实训目的 理解“最优化设计”的概念和方法；了解“最优化设计”的建模方法；了解“最优化设计”问题的解决过程；了解LINGO的编程、方法和求解过程。三、实验步骤 1.运行LINGO程序。LINGO程序主界面如图3-3所示。2.点击下拉菜单，选择或者单击工具栏中的按钮，新建一个优化实例模型LINGOModel(*.lg4)。 3.在“LINGOModel”窗口建立模型，输入以下程序语句：Model:min=3.14*(d2^2-d1^2)*L/4;!目标函数； !径约束条件如下;d1min=0.02; d1ma*=0.03; d1>=dimin; d1<=d1ma*;!外径约束条件如下; d2min=0.02; d2ma*=0.03; d2>=d2min; d2<=d2ma*;!长度约束条件如下; Lmin=0.1; Lma*=1; L>=Lmin; L<=Lma*;!定义安全系数; a=1.2; b=1.3;!外加载荷20kN; P=20*10^3; E=200*10^9; !弹性模量;C=350*10^6; !材料抗压极限; (4*a*P)/(3.14*(d2^2-d1^2))<=C; !不产生屈服的约束条件;u=0.5; !长度折算系数; J=3.14*(d2^4-d1^4)/64; !刚度系数;图46求解状态窗口 P<=(3.14^2*E*J)/(b*(u*L)^2);图46求解状态窗口END4.模型建立完成后要对模型进行求解：单击下拉菜单“LINGO”选择“Sovle”，或者单击工具栏中的模型求解命令按钮。LINGO程序会弹出“LINGOSolveStatus”窗口和“SolutionReport”窗口。在“SolutionReport”窗口中会列出完整的模型求解过程，在“LINGOSolveStatus”窗口中列出了求解器的有关信息，如图46所示。实验五悬臂梁的可靠度分析一、问题描述如图47所示，长度L=2m的悬臂梁，自由端受集中力F=50KN作用。弹性模量E=3.0e7kPa，泊松比v=0.17。悬臂梁截面：面积二、实验步骤1)建立有限元模型图47自由端受集中力F=50KN作用(1)定义分析标题:UtilityMenu->File->ChangeTitle,输入”ProbabilisticDesignofabeam”图47自由端受集中力F=50KN作用，单击OK。(2)定义参数初始值：UtilityMenu->Parameters->ScalarParameters,打开数值参数对话框。在Selection编辑栏中输入F＝50，按Enter键：输入E＝3．0E7，按Enter键；单击C1ose关闭数值参数对话框。(3)定义单元类型：MainMenu->Preprocessor->ElementType->Add/Edit/Delete,出现ElementTypes列表框。单击Add出现单元类型库对话框，左侧列表中选择StructuralBeam，右侧列表中选择2DElastic3，单击OK。(4)定义实常数：MainMenu->Preprocessor->RealConstants->Add/Edit/Delete,出现RealConstants列表框。单击Add出现E1ementTypeforRealConstants对话框，单击OK，出现RealConstantSetNumber1,ForBeam3对话框。在AREA中输入0.04，IZZ中输入0.0003，HEIGHT中输入0.2，单击OK。单击Close关闭RealConstants列表框.(5)定义材料特性：MainMenu->Preprocessor->MaterialProps->MaterialModels,在DefineMaterialModelBehavior窗口中，双击Structural->Linear->Elastic->Isotropic。在出现的对话框中，在E*中输入3.0E7，在PR*Y中输入0.17，单击OK,MaterialModell出现在MaterialModelsDefined窗口左侧列表中，选择菜单Material->E*it退出。(6)建立模型：MainMenu->Preprocessor->Create->Keypoints->InActiveCS,出现CreateKeypointsinActiveCoordinateSystem对话框。在LocationinActiveCS中输入0，0，0，单击Apply,输入2，0，0，单击OK,图形窗口出现2个关键点。执行MainMenu->Preprocessor->Create->Lines->StraightLine，出现CreateStraightLine选择框。在图形窗口中依次选择l、2关键点，单击OK生成一直线段。(7)定义线的属性：MainMenu->Preprocessor->Define->AllLines，出现Lineattributes对话框。将PickOrientationKeypoint激活为YES，单击OK，出现Lineattributes选择框，在选择框单击OK。(8)划分网格：Menu->Preprocessor->->MeshTool,出现MeshTool对话框。单击Lines右侧的Set，出现ElementSizeonPicked选择框，单击PickA11，出现ElementSizeonPickedlines对话框。在NDIV中输入20，单击OK，在图形窗口中每段弧线以20段均分线段表示。在MeshTool对话框中单击Mesh，出现MeshLines选择框，单击PickAll生成单元。(9)定义分析类型：Menu->Solution->UnabridgedMenu->NewAnalysis．出现NewAnalysis对话框，选样Static,单击OK。(10)施加约束与荷载：Menu->Solution->Loads->Apply->Displacement->OnNodes，出现ApplyUROTonNodes选择框。在图形窗口中选择左端节点，在选择框中单击Ok出现ApplyUROTonNodes对话框。选中A11DOF单击OK。约束标识出现在图形窗口。执行Menu->Solution->Loads->Apply->Force/Moment->OnNodes，出现ApplyF/MonNodes选择框，选择右端节点，单击OK，出现ApplyF/MonNodes对话框。在Lab中选样FY，在Value中输入“F”，单击OK，集中力标识出现在梁的自由端。(11)保存数据库文件：UtilityMenu->File->SaveasJobname.db。(12)求解：Menu->Solution->Loads->Solve->CurrentLS，查看求解信息，关闭求解状态窗口，单击OK开始求解。求解完成后单击Close关闭求解信息框。(13)读出数据：MainMenu->GeneralPostproc->ElementTable->DefineTable，得到ElementTableData对话框.单击Add，得到DefineAdditionElementaryTableItem对话框.在UserLabel编辑框中输入SMIN_I，在Item列表框中选择BySequenceNum，征p列表框中选择NMISC,在p列表框下方的编辑框中输入“NMISC，2”，单击Apply；再次在UserLabel编辑框中输入SMIN_J，在Item列表框中选择BySequenceNum，征p列表框中选择NMISC,在p列表框下方的编辑框中输入“NMISC，4”，单击OK关闭对话框，在ElementTableData对话框中单击Close关闭对话框。(14)单元表应力数据排序：MainMenu->GeneralPostproc->->ListResult->SortElems,得到SortElement对话框。在ORDER中选择Descendingorder，在KABS中激活Yes选项，Item，p中选择SMIN—I，单击OK，对单元表数据SMIN—I按绝对值降序排列。提取应力数据：在输入对话框中输入*GET,SMINI，SORT，，MA*，提取单元I端绝对值最大应力数据。执行类似的操作：MainMenu->GeneralPostproc->->ListResult->SortElems,得到SortElement对话框。在ORDER中选择Descendingorder，在KABS中激活Yes选项，Item，p中选择SMIN—J，单击OK，对单元表数据SMIN—J按绝对值降序排列。提取应力数据：在输入对话框中输入*GET,SMINJ，SORT，，MA*，提取单元J端绝对值最大应力数据。(15)定义应力极值：UtilityMenu->Parameters->ScalarParameters，打开数值参数对话框。在Selection编辑栏中输入SMIN＝ABS(SMINI)-ABS(SMINJ)，按Enter键，单击Close关闭数值参数对话框。(16)生成优化分析文件：UtilityMenu->File->WriteDBLogFile，打开WriteDBLog对话框，在WriteDBLogto编