ANSYS经典应用实例(结构分析详解).ppt

5.1 5.1 结结构静力学分析中用到的构静力学分析中用到的单单元元 ANSYS软件中结构静力分析用来分析由于稳态外载荷引起的系统或部件的位移、应变、应力和力。稳态外载荷包括稳定的惯性力（如重力、旋转件所受的离心力）和能够等效为静载荷的随时间变化的载荷。这种分析类型有很广泛的应用，如确定确定结构应力集中程度，预测结构最大应力等。类别形状和特性单元类型1.杆普通LINK1,LINK8双线性LINK102.梁普通BEAM3,BEAM4截面渐变BEAM54,BEAM44塑性BEAM23,BEAM24考虑剪切变形BEAM188,BEAM189类别形状和特性单元类型3.管普通普通PIPE16，PIPE17，PIPE18浸入浸入PIPE59塑性塑性PIPE20，PIPE604.二维实体三角形三角形PLANE2四边形四边形PLANE42，PLANE82，PLANE182超弹性单元超弹性单元HYPER84，HYPER56，HYPER74粘弹性粘弹性VISCO88大应变大应变VISO106，VISO108谐单元谐单元PLANE83，PLANE25P单元单元PLANE145，PLANE146类别形状和特性单元类型5.三维实体块块SOLID45，SOLID95，SOLID73，SOLID185四面体四面体SOLID92，SOLID72层层SOLID46各向异性各向异性SOLID64，SOLID65超弹性单元超弹性单元HYPER86，HYPER58，HYPER158粘弹性粘弹性VISO89大应变大应变VISO107P单元单元SOLID147，SOLID1486.壳 四边形四边形SHELL93，SHELL63，SHELL41，SHELL43，SHELL181轴对称轴对称SHELL51，SHELL61层层SHELL91，SHELL99剪切板剪切板SHELL28P单元单元SHELL150类别 形状和特性单元类型7.接触单元面面TARGET169，TARGET170，SURF171，SURF172，SURF173点面SURF174点点CONTAC48，CONTAC49刚性表面CONTAC12，CONTAC52，CONTAC268.专业单元 弹簧COMBIN14，COMBIN39，COMBIN40质量MASS21控制单元COMBIN37表面效应单元SURF19，SURF22，SURF153，SURF154铰COMBIN7线性激发器LINK11矩阵MATRIX27，MATRIX50类别 形状和特性单元类型9.耦合场声学声学TARGET169，TARGET170，SURF171，SURF172，SURF173压电压电热应力热应力SURF174磁结构磁结构CONTAC48，CONTAC49流体结构流体结构CONTAC12，CONTAC52，CONTAC265.2 5.2 结结构静力学分析的构静力学分析的类类型型静力分析线性静力分析非线性静力分析非线性静力分析允许有大变形、蠕变、应力刚化、接触单元、超弹性单元等。线性静力分析根据结构特征和所受外载荷的形式通常归结为如下几个问题：平面问题、轴对称问题、周期对称问题、三维问题等。5.2.1 5.2.1 平面平面问题问题平面问题平面应变平面应力特点：物体沿某坐标轴（如Z轴）的尺寸远大于其它两个坐标轴的尺寸；垂直于Z轴各截面的形状和尺寸均相同；所有外力与Z轴垂直且不随Z坐标变化；物体的约束条件不随Z坐标而变化。在这种情况下，可以认为物体沿Z轴方向各截面没有z向位移，而沿x和y向的位移对各截面均相同（与z坐标无关），各截面内将产生平面应变。任何机械零部件，一般说来都是空间结构。但是在某些条件下，它们可以简化为平面问题来处理。平面问题包括平面应变问题和平面应力问题两类。特点：物体沿某坐标轴（如Z轴）的尺寸远小于其它两个坐标轴的尺寸；外力沿周边作用且与XY平面平行，且体积力也垂直于z轴；由于物体在z方向厚度很小，故外载的表面力和体积力都可看成是沿z向不变化的。约束条件在XY平面内。在这种情况下，可以认为物体沿Z轴方向无应力，所有应力都发生在XY平面内。5.2.2 5.2.2 轴对轴对称称结结构构问题问题与周期与周期对对称称结结构构问题问题轴对称结构问题周期对称结构问题特点：（1）结构为回转体（截面绕它的回转中心轴旋转而形成）；（2）载荷关于轴心线对称。ANSYS软件中提供了专门的分析方法对这类问题进行求解，与普通方法相比可以节约大量的人力和计算机资源，大大提高求解问题的效率。特点：如果结构绕其轴旋转一个角度，结构（包括材料常数）与旋转前完全相同，则将这种结构称为周期对称结构（循环对称结构）。符合这一条件的最小旋转角5.3 杆系结构分析实例问题描述：图5-1所示为一平面桁架，长度L=0.1m，各杆横截面面积均为A=1 10-4m2 ，力P=2000N，计算各杆的轴向力Fa、轴向应力a。图5-1 平面桁架 实实例例1 1 平面桁架分析平面桁架分析杆轴向力Fa/N轴向应力a/MPa10001001000100-1414.2-141.400-1414.2-141.4思考:根据静力平衡条件，很容易计算出轴向力Fa、轴向应力a，如表5-1所示。表5-1 各杆的轴向力和轴向应力求解步骤：求解步骤：创建单元类型选择Structusral Beam类的2D elastic 3。定义实常数定义AREA=1E-4。定义材料特性输入Ex=2e11，PRXY=0.3。创建节点GUI：Main MenuPreprocessorModelingCreateNodesIn Active cs在NODE文本框中依次输入节点1、2、3、4的X、Y、Z坐标节点1：0，0，0；节点2：0.1，0，0；节点3：0.2，0，0；节点4：0.1，0.1，0；显示节点号、单元号GUI：Utility MenuPlotCtrlsNumbering。在弹出的“Plot Numbering Controls”对话框，将节点号和单元号打开。创建单元GUI：Main MenuPreprocessorModelingCreateEleMentsAuto NumberedThru Nodes。弹出拾取窗口，拾取节点1和2，单击拾取窗口的“Apply”按钮，于是在节点1和2之间创建了一个单元。重复以上过程，在节点2和3、1和4、3和4间分别创建单元，建模图形如左图5-2所示。图5-2 桁架建模图 施加约束GUI：Main Menu SolutionDefine LoadesApplyStructuralDisplacementOn Node。弹出拾取窗口，拾取节点1，单击ok按钮，在弹出的对话框列表中选择“UX、UY”，单击Apply按钮。再次弹出拾取窗口，拾取节点3，单击 ok按钮，在弹出的对话框列表中选择“UY”，单击ok按钮，完成对模型的约束施加。施加载荷GUI：Main Menu SolutionDefine LoadesApplyStructuralForce/MomentOn Nodes。弹出拾取窗口，拾取节点4，单击ok按钮，在对话框中选择Lab为“FY”，在VALUE文本框中输入-2000，单击ok按钮。求解GUI：Main Menu SolutionSolvCurrent LS。单击“Solution Current Load Step”对话框ok按钮。出现“solution is done！”提示时，求解结束，即可查看求解结果。结果显示：定义单元列表GUI：Main MenuGeneral PostrocElement TableDefine Table。弹出“Element Table Data”对话框，单击Add按钮，在Lab文本框中输入FA，在“Item,Comp”两个列表中分别选择“By sequn-ence num”、“SMISC”,在右侧列表下方文本框输入SMISC，1，单击Apply按钮，于是定义了单元表FA，该单元列表保存了各单元的轴向力；在Lab文本框中输入SA，在“Item,Comp”两个列表中分别选择“By sequnence num”、“Ls”,在右侧列表下方文本框输入LS，1，单击ok按钮，于是定义了单元表SA，该单元列表保存了各单元的轴向应力。列表单元表数据GUI：Main MenuGeneral PostrocElement TableListElem Table。在弹出的对话框的列表中选择FA、SA，单击ok按钮，即显示出求解结果。与表5-1对照，二者完全一致。5.4 5.4 梁结构分析实例梁结构分析实例 实实例例11简简支梁支梁载载荷工况荷工况组组合合实实例例图5-3 简支梁问题描述及解析解：图5-3（a）所示为一圆截面简支梁，跨度L=1m，圆截面直径D=30mm，作用在梁上的集中力为P=1000N,作用点距支座A距a=0.2m，已知弹性模量E=2e11.由材料力学知识可得：梁截面的惯性距为最大挠度思路分析:当进行线性分析时，简支梁的应力、应变和变形等于如图5-3（b）、（c）所示两个简支梁的结果叠加。如图5-3（b)所示的简支梁结构和载荷均对称于梁的中点O，故应力、应变和变形也对称于梁的中点O，进行有限元分析时，可简化为如图5-4(a)所示的模型。如图5-3(c）所示的简支梁结构、载荷反对称，故应力、应变和变形也反对称于梁的中点O，因此可简化为如图5-3（b）所示的模型。对如图5-4（b）、（c）所示的模型进行有限元分析，将结果分别进行相加和相减，即可分别得到如图？所示的简支梁中点左、右两半部分的结果。图5-3 简支梁图5-4 简支梁的简化模型 创建单元类型选择Structusral Beam类的2D elastic 3 单元。定义实常数定义AREA=7.069e-4、IZZ=3.976e-8、HEIGHT=0.03。定义材料特性输入Ex=2e11（弹性模量），PRXY=0.3（泊松比）。创建关键点GUI：MainMenuPreprocessorModelingCreatKeypointsIn Active cs。在弹出对话框的NPT文本框中输入1，在“X、Y、Z”文本框中分别输入0，0，0.单击Apply按钮；在NPT文本框中输入2，“X、Y、Z”文本框中分别输入0.5，0，0.单击ok按钮.求解步骤：显示关键点号GUI：Utility MenuPlotCtrlsNumbering。在弹出的对话框中，将关键点号打开，单击ok按钮。创建直线GUI：MainMenuPreprocessorModelingCreateLines LinesStraight Line弹出拾取窗口，拾取关键点1和2，单击ok按钮。创建硬点GUI:MainMenuPreprocessorModelingCreateKeypointsHard PT On line Hard PT by ratio。弹出拾取窗口，拾取直线，单击ok按钮，在文本框中输入0.4，单击ok按钮。（为了施加集中载荷，作用点处创建了一个硬点，该硬点属于直线）划分单元GUI：Main MenuPreprocessorMeshingMeshTool。弹出MeshTool对话框，单击“Size Controls”区域中的“Line”后的set按钮，弹出拾取窗口，拾取直线，单击ok按钮，在NDIV文本框中输入50，单击Mesh按钮，弹出拾取窗口，拾取直线，单击ok按钮。显示点、线、单元GUI：UtilityPlotMulti-Plots。施加载荷施加第一个载荷步 施加第一个载荷步的位移载荷GUI:MainMenuSolutionDefineLoadsApplyStructuralDisplacementOn Keypoint。对关键点1加载UX、UY方向的载荷，关键点2加载UX、ROTZ方向的载荷。（注意：约束关键点2的转动自由度，相当于在该处时加一个弯矩）施加第一个载荷步的力载荷GUI:MainMenuSolutionDefineLoadsApplyStructuralForce/MomentOn Keyopints。弹出拾取窗口，拾取关键点3，单击ok按钮，选择Lab为FY,在VALUE文本框中输入-500，单击ok按钮。写第一个载荷步文件GUI：MainMenuGeneral PostprocLoad CaseRead Load Case。在LCNO文本框中输入1，单击ok按钮。删除位移载荷GUI：MainMenuSolutionDefine LoadsDeleteStructuralDisplacementOn Keypoint。弹出拾取窗口，拾取关键点2，单击ok按钮，弹出一个对话框，然后单击ok按钮。施加第二个载荷步 施加第二个载荷步的位移载荷GUI:MainMenuSolutionDefineLoadsApplyStructuralDisplacementOn Keypoint。对关键点2加载UX、UY方向的载荷。写第二个载荷步文件GUI：MainMenuSolutionLoad Step OptsWrite LS File.在LSNUM文本框中输入2，单击ok按钮。求解GUI：Main Menu SolutionSolveFrom LS Files。在LSMIN文本框中输入1,在LSMAX文本框中输入2,在LSINC文本框中输入1，单击ok按钮，出现“solution is done！”提示时时，求解结束，即可查看求解结果。结果显示：创建载荷工况GUI：MainMenuGeneral PostprocLoad CaseCreate Load Case。单击ok按钮，在LCNO文本框中输入1，在LSTEP文本框中输入1，单击Apply按钮，再次单击ok按钮，再在LCNO文本框中输入2，在LSTEP文本框中输入2，单击ok按钮。于是，创建了两个载荷工况。其中，载荷工况1对应着载荷步1，载荷工况2对应着载荷步2。将载荷工况1读入内存GUI：MainMenuGeneral PostprocLoad CaseRead Load Case。在LCNO文本框中输入1，单击ok按钮。对载荷工况进行加运算GUI：MainMenuGeneral PostprocLoad CaseAdd。在LCASE1文本框中输入2，单击ok按钮。查看结果，显示变形GUI：MainMenuGeneral PostprocPlot ResultsDeformed Shape。在弹出的对话框中，选中Def Shape only，单击ok按钮。即得到简支梁左半部分的结果。写载荷工况组合到结果文件GUI：MainMenuGeneral PostprocWrite Results。在弹出的对话框的LSTEP文本框中输入9999，单击ok按钮。查看简支梁有半部分的结果重复以上-的过程，在步骤中，拾取菜单MainMenuGeneral PostprocLoad CaseSubtract，对载荷工况进行减运算，可得到简支梁右半部分的结果。（注意：与简支梁的实际情况呈对称关系）。5.5 平面问题分析实例问题描述：图5-2所示为一厚壁圆筒，其内径r1=50mm,外径r2=100mm,作用在内孔上的压力p=10MPa,无轴向压 力，轴向长度视为无穷。要求计算后壁圆筒的径向应力r和切向应力t沿半径r方向的分布。图5-2厚壁圆筒问题 实实例例11厚壁厚壁圆圆筒筒问题问题思路:根据材料力学的知识，r、t沿r方向的分布的解析解为:根据对称性，可取圆筒的四分之一并施加垂直于对称面的约束进行分析。(5-1)求解步骤：求解步骤：创建单元类型选择Structusral Solid类的8node 183 单元。单击options按钮，选择K3为Plane strain(平面应变）。定义材料特性输入Ex=2e11（弹性模量），PRXY=0.3（泊松比）。创建实体模型GUI：Main MenuPreprocessorModelingCreateAreas Circle By Dimensions。在RAD1、RAD2、THETA2文本框中分别输入0.1、0.05和90，单击ok按钮。划分单元GUI：Main MenuPreprocessorMeshingMeshTool。弹出MeshTool对话框，单击“Size Controls”区域中的“Line”后的set按钮，弹出拾取窗口，拾取面的任意直线边，单击ok按钮，在NDIV文本框中输入6，单击Apply按钮，再次弹出拾取窗口，拾取面的任意圆弧边，单击ok按钮，在NDIV文本框中输入8，单击OK按钮。在Mesh区域，选择单元形状Quad（四边形），选择划分单元的方法为Mapped（映射）。单击Mesh按钮，弹出拾取窗口，单击ok按钮。施加约束GUI：Main MenuSolutionDefine LoadsApplyStructural DisplacementOn Lines。弹出拾取窗口，拾取水平直线边，单击ok按钮，在列表中选择UY，单击Apply按钮，再次弹出拾取窗口，拾取面的垂直直线边，单击ok按钮，在列表中选择UX，单击ok按钮。施加载荷GUI：MainMenuSolutionDefineLoadsApplyStructuralPressureOn Lines。弹出拾取窗口，拾取面的内侧圆弧边（较短的一条圆弧），单击ok按钮，在VALUE文本框中输入10e6，单击ok按钮。求解GUI：Main Menu SolutionSolveCurrent LS。单击“Solution Current Load Step”对话框ok按钮。出现“solution is done！”提示时，求解结束，即可查看求解结果。结果显示：显示节点GUI：Utility MenuPlotNodes。定义路径GUI：MainMenuGeneralPostprocPathOperationsDefine PathBy Nodes。弹出拾取窗口，由内向外依次拾取与X轴平行的水平直线边上的所有节点，单击ok按钮，在Name文本框中输入p1，单击ok按钮。将数据映射到路径上GUI：Main MenuGeneral PostprocPath OperationsMap onto Path。在Lab文本框中输入SR，在“Item,Comp”两个列表中分别选择“Stress”、“X-direction SX”,单击Apply按钮，在Lab文本框中输入ST，在“Item,Comp”两个列表中分别选择“Stress”、“X-direction SX”,单击ok按钮。（注：该路径上各节点X、Y放上的应力即径向应力r和切向应力t）作路径图GUI：Main MenuGeneral PostprocPath OperationsPlot Path ItemOn Graph。在弹出的对话框中的列表中选择SR、ST，单击ok按钮。即可显示径向应力r和切向应力t 随半径的分布情况。5.6 5.6 空间问题实例分析空间问题实例分析 实例实例11 各种坐标系下的扭转问题各种坐标系下的扭转问题设等直圆轴的圆截面直径D=50mm，长度L=120mm，作用在圆轴两端上的转矩Mn=1.5103N M,求圆柱的最大剪切应力和转角。问题描述及解析解：由材料力学知识可得：圆截面对圆心的极惯性矩为圆截面的抗扭截面模量为圆截面上任意一点的剪应力与该点半径成正比，在圆截面的边缘上有最大值等直圆轴距离为0.09m的两截面间的相对转角求解步骤：创建单元类型选择Structusral Solid类的Quad 4node 42和Brick 8node 45单元。定义材料特性输入Ex=2e11（弹性模量），PRXY=0.3(泊松比）。创建矩形面GUI：MainMenuPreprocessorModelingCreatAreasBy Dimensions。在X1、X2文本框中输入0，0.025；在Y1、Y2文本框中输入0，0.12.划分单元GUI：Main MenuPreprocessorMeshingMeshTool弹出MeshTool对话框，单击“Size Controls”区域中的“Line”后的set按钮，弹出拾取窗口，拾取矩形面的任一短边，单击ok按钮，在NDIV文本框中输入4，单击Apply按钮，再次弹出拾取窗口，拾取矩形面的任一长边，再次弹出单元尺寸对话框，在NDIV文本框中输入8，单击ok按钮。在Mesh区域，选择单元形状Quad（四边形），选择划分单元的方法为Mapped，单击Mesh按钮，弹出拾取窗口，拾取面，单击ok按钮。设定挤出选项GUI：MainMenuPreprocessorModelingOperateExtrudeElem Ext opts。在VAL1文本框中输入3（挤出段数），选定ACLEAR为Yes（清除矩形面上的单元），单击ok按钮。由面旋转挤出体GUI：MainMenuPreprocessorModelingOperateExtrudeAreasAbout Axis弹出拾取窗口，拾取矩形面，单击ok按钮；再次弹出拾取窗口，拾取矩形面在Y轴上两个关键点，单击ok按钮；在随后出现的对话框中的ARC文本框中输入360，单击ok按钮。显示单元GUI：Utility MenuPlotElements。改变视点GUI：Utility MenuPlotCtrlsPan zoom Roate。在弹出的对话框中，依次单击ISO、Fit按钮。旋转工作平面GUI：Utility MenuWorkPlaneOffset WP By Increment。在XY、YZ、ZX Angles文本框中输入0，90，单击ok按钮。创建局部坐标系GUI：Utility MenuLocal Coordinate SystemCreate LocalCSAt WP Origin。在KCN文本框中输入11，选择KCS为Cylindrical 1，单击ok按钮。即创建了一个代号为11、类型为圆柱坐标系的局部坐标系，并激活其成为当前坐标系。选中圆柱面上的所有节点GUI：Utility MenuSlectEntity。在弹出的对话框的下拉列表框、文本框、单选按钮中依次选择或输入Nodes、By Location、X Coordinates、0.025、From Full，单击Apply按钮。旋转节点坐标系到当前坐标系GUI：MainMenuPreprocessorModelingMove/ModifyRoate Node CSTo Active CS。弹出拾取窗口，单击Pick All按钮。施加约束GUI：MainMenuSolutionDefineLoadsApplyStructuralDisplacementOn Nodes。弹出拾取窗口，选择Pick All按钮。在弹出的对话框中的Lab2列表框中选择UX，单ok按钮。选中圆柱面上最上端的所有节点GUI：Utility MenuSlectEntity。在弹出的对话框的下拉列表框、文本框、单选按钮中依次选择或输入Nodes、By Location、Z Coordinates、-0.12、Reselect，单击Apply按钮。施加载荷GUI：MainMenuSolutionDefineLoadsApplyStructuralForce/MomentOn Nodes。弹出拾取窗口，单击Pick All按钮。在弹出的对话框中的Lab下拉列表框中选择FY，在VALUE文本框中输入5000，单击ok按钮。这样在结构上一共施加了12个大小为5000N的集中力，它们对圆心的力矩和为1500N.m 选择所有GUI：Utility MenuSlectEverything。显示体GUI：Utility MenuPlotVolumes。施加约束GUI：MainMenuSolutionDefineLoadsApplyStructuralDisplacementOn Areas。弹出拾取窗口，拾取圆柱体下侧底面（由四部分组成），单击ok按钮。在弹出的对话框中的Lab2列表框中选择All DOF，单ok按钮。求解GUI：Main Menu SolutionSolveCurrent LS。单击“Solution Current Load Step”对话框ok按钮。出现“solution is done！”提示时，求解结束，即可查看求解结果。结果显示：显示变形GUI：Main MenuGeneral PostprocPlot ResultsDeformed Shape。在弹出的对话框中，选中Def+undeformed（变形+未变形的单元边界），单击ok按钮。改变结果坐标系为局部坐标系GUI：Main MenuGeneral PostprocOptions for Outp。在RSYS下拉列表框中选择Local system，在Local system Reference no文本框中输入11，单击ok按钮。选择z=-0.09m的所有节点GUI：Utility MenuSlectEntity。在弹出的对话框的下拉列表框、文本框、单选按钮中依次选择或输入Nodes、By Location、Z Coordinates、-0.09、From Full，单击Apply按钮。再在下拉列表框、文本框、单选按钮中依次选择或输入Nodes、By Location、Y Coordinates、0、Reselect，单击Apply按钮。列表显示节点位移GUI：Main MenuGeneral PostprocList ResultsNodal Solution。在左侧列表中选择DOF Solution，在右侧列表中选择UY，单击ok按钮。选择单元GUI：Utility MenuSlectEntity。在弹出的对话框的下拉列表框、文本框、单选按钮中依次选择或输入Nodes、By Location、Z Coordinates、-0.1，0、From Full，单击Apply按钮。再在下拉列表框、文本框、单选按钮中依次选择或输入Elements、Attachedto、Nodes all、Reselect，单击Apply按钮。这样做的目的是，在下一步显示应力时，不包含集中力作用点附近的单元，以得到更好的计算结果。查看结果，用等高线显示剪应力GUI：Main MenuGeneral PostprocPlot ResultsContour PlotElements Solu。在Item,Comp两个列表中分别选择Stress、“YZ-shear SYZ”,即可查看结果。如图所示轴承座结构，尺寸单位为mm。已知沉孔径向内柱面承受外推压力0.0069MPa，轴承孔柱面下部分承受向下压力0.0345MPa，计算轴承的最大应力及其位置。实例2:轴承座分析分析：约束方式：四个通孔为轴承座安装固定孔，必须给予对称约束；除此之外，还必须限制轴承座垂直方向的移动，避免该方向的刚体位移。1.几何建模（1）建立底座：VolumesBlockBy Dimension 轮子二维结构如图所示（尺寸单位为英寸），现要分析该轮仅承受绕其中心轴旋转角速度的作用下，轮的受力及其变形情况。已知角速度：w=525rad/s,材料属性：弹性模量E30106psi 泊松比为0.3 密度为0.000731bf-2/in4实例3轮子分析思路分析：这是一个空间问题的静力学分析，由于结构的复杂性，不能采用由2D网格拖拉生成3D网格的方式，宜采用自由网格和映射网格相结合的网格生成方式，根据该轮的对称性，在分析式只要分析其中的一部分即可。1）先建立三维几何模型；2）根据三维模型的特点，对三维模型进行分割，使其中一部分模型能够采用映射网格方式划分，其余部分则可以用自由网格划分。因此必须选用20节点六面体单元（Solid95)，它可以进行映射网格划分，同时也可以变成20节点的四面体网格进行自由网格划分。因此能够实现不同网格形状的连接。3）对划分网格的连接处进行网格转换。即将20节点的四面体自由网格转换为10节点四面体网格，以减六面体单元（Solid95)，它可以进行映射网格划分，同时也可以变成20节点的四面体网格进行自由网格划分。因此能够实现不同网格形状的连接。4）对划分网格的连接处进行网格转换。即将20节点的四面体自由网格转换为10节点四面体网格，以减少解题规模。5）施加对称约束和角加速度。6）求解分析。7）显示结果。GUI操作方式（1）定义工作名”Wheel_Anal“和工作标题”The Stress calculating of Wheel by angular velocity”（2）定义单元类型和材料属性选择Structural Solid类型的Brick 20node 95两种单元；材料属性：在EX框输入30e6，PRXY框输入0.3，DENS框输入0.00073（3）建立2D模型（如下页所示）Main MenuPreprocessorModelingCreateAreasRectangleBy Dimension 生成左边矩形面：X11.0,X21.5，Y10,Y2=5.0 生成右边矩形面：X13.25,X23.75，Y10.5,Y2=3.75 生成中间矩形面：X11.0,X23.75，Y11.75,Y2=2.5面叠分布尔操作Overlap，Pick All线倒角，生成由倒角线围成的面生成弧线的中心关键点（共两个）：X=3.5,Y=0.475X=3.5,Y=3.505生成圆弧线(如图所示）ArcBy End Kps&RadR=0.35生成由圆弧线围成的面（共两处）面相加（Pick All）线相加以减少线的数量（可不作）：L17L12L13；L22L10L20压缩编号操作，并重新显示线保存结果数据：存为Wheel_Anal_2D（4）通过拖拉生成3D模型生成轴线的关键点：(0,0)-(0,5)2D拖拉成3DPreprocessorModelingOperateExtrudeAbout Axis在Arc length in degree框中输入22.5度，输入NSEG1，即生成的实体由一块体积组成关闭线号显示改变视图方向：ISO结果如下页所示。（5）生成一个圆柱孔激活工作平面旋转工作平面：WorkplaneOffset by increments，在”XY，YZ，ZX Angles”下面的输入栏中输入”0，90“,工作平面的”WZ”将指向上方；生成一个圆柱体：CreateCylinderSolid Cylinder，出现一个对话框，输入WP X=2.375WP Y=0，Radius0.45，Depth2.6(注：2.5)体相减得到的3D几何模型如图所示另存盘为：Save as“Wheel_Anal_3D”（6）生成网格平移工作平面到关键点9上。用工作平面切分体：OperateDivideVolu by WrkPlane，拾取“Pick All”。显示体平移工作平面到关键点14上。用工作平面切分体：OperateDivideVolu by WrkPlane，拾取下部分体（编号为V4）。关闭工作平面结果如上图所示设置单元尺寸：用Mesh Tool工具，单击Gloabal上的Set，输入Size0.25采用映射网格生成单元：在MeshTool工具条下的“Shape”下选择”Hex”六面体网格类型和“Mapped”映射生成方式，单击Mesh，拾取V1，V2，V3，V5。生成的网格单元如下图所示。采用自由网格划分单元：在MeshTool工具条上，选择Shape下的“Tet”四面体网格类型和“Free”自由网格生成方式，单击Mesh，拾取V6，则完成自由网格划分。注：20节点六面体网格到10节点四面体网格之间的连接，必须在会合面上通过五面体（20节点的退化形式）过渡，即20节点六面体20节点五面体20节点四面体，退化是通过节点重合实现的。转变的单元类型：MeshingModify MeshChange Tets，出现一个对话框，在“Change from”后面的框中选择”95 to 92”。即20节点95退化四面体单元转换成10节点四面体单元。以降低求解规模。显示单元的过渡区域：通过建立选择集来局部观察，UtilitySelectEntities弹出对话框，选择Elements和By Element name，输入Soli95，再显示单元UtilityplotElement，得到如图左上所示图形；再选择全部实体UtilitySelectEverything,UtilitySelectEntities弹出对话框，选择Elements和By Element name，输入Soli92，再显示单元UtilityplotElement，得到如图左下所示图形；可以清楚看到，六面体单元到四面体单元通过连接边界面附近的五面体单元过渡而成。保存单元网格数据：save as”Wheel_Anal_Mesh”（7）施加载荷并求解1）在面上施加对称约束：由于图中截面（A1+A39+A7+A20+A40+A27）及其反面（A15A41A3A34A41）均不能沿面的法向移动，面也不能沿两个方向摆动，因此应对这两个面施加对称约束。.DisplacementSymmetry B.C.On Areas在弹出的对话框中输入1,39,7,20,40,27，单击Apply，再输入15，21，3，34，41,单击OK注：这两个截面上的约束实际上是斜约束，即与总体坐标系成角度，因此无法用诸如UX，UY.等来表示，只能通过对称约束实现2）在关键点1上施加UY0的约束，原因是防止计算误差导致轴向合力近似0，而这一较小的轴向力将造成Y向刚体位移。因此只要约束一个角点就足够了。3）施加旋转角速度：ApplyStructuralIntertiaAngular VelocGlobal，在弹出的对话框OMEGY框中输入5254）保存：Save As“Wheel_Anal_Load”5)开始求解运算SolutionSolveCurrent LS(8)浏览分析结果浏览Von Mises应力：General postProcPlot ResultsContour PlotNodal Solu，弹出对话框，在Item to be contoured后面左栏选择Stress，右栏选择Von Mises SEQV，单击OK。结果如图所示。周向扩展浏览1）将工作平面平移到全局直角坐标系的原点，并在此工作平面的原点上定义局部柱坐标系；2）绕局部坐标柱系的Z轴（轴向）扩展结果：UtilityPlotCtrlsStyleSymmetry ExpansionUser Specified Expansion在对话框中分别输入和选择：NREPEAT16，TYPELocal PolarPATTERNAlternate Symm，DY22.5