1、5.1 结构静力学分析中用到的单元,ANSYS软件中结构静力分析用来分析由于稳态外载荷引起的系统或部件的位移、应变、应力和力。稳态外载荷包括稳定的惯性力（如重力、旋转件所受的离心力）和能够等效为静载荷的随时间变化的载荷。这种分析类型有很广泛的应用，如确定确定结构应力集中程度，预测结构最大应力等。,5.2 结构静力学分析的类型,静力分析,线性静力分析,非线性静力分析,非线性静力分析允许有大变形、蠕变、应力刚化、接触单元、超弹性单元等。,线性静力分析根据结构特征和所受外载荷的形式通常归结为如下几个问题：平面问题、轴对称问题、周期对称问题、三维问题等。,5.2.1 平面问题,平面问题,平面应变,平面
2、应力,特点：物体沿某坐标轴（如Z轴）的尺寸远大于其它两个坐标轴的尺寸；垂直于Z轴各截面的形状和尺寸均相同；所有外力与Z轴垂直且不随Z坐标变化；物体的约束条件不随Z坐标而变化。在这种情况下，可以认为物体沿Z轴方向各截面没有z向位移，而沿x和y向的位移对各截面均相同（与z坐标无关），各截面内将产生平面应变。,任何机械零部件，一般说来都是空间结构。但是在某些条件下，它们可以简化为平面问题来处理。平面问题包括平面应变问题和平面应力问题两类。,特点：物体沿某坐标轴（如Z轴）的尺寸远小于其它两个坐标轴的尺寸；外力沿周边作用且与XY平面平行，且体积力也垂直于z轴；由于物体在z方向厚度很小，故外载的表面力和体
3、积力都可看成是沿z向不变化的。约束条件在XY平面内。在这种情况下，可以认为物体沿Z轴方向无应力，所有应力都发生在XY平面内。,5.2.2 轴对称结构问题与周期对称结构问题,轴对称结构问题,周期对称结构问题,特点： （1）结构为回转体（截面绕它的回转中心轴旋转而形成）；（2）载荷关于轴心线对称。ANSYS软件中提供了专门的分析方法对这类问题进行求解，与普通方法相比可以节约大量的人力和计算机资源，大大提高求解问题的效率。,特点： 如果结构绕其轴旋转一个角度，结构（包括材料常数）与旋转前完全相同，则将这种结构称为周期对称结构（循环对称结构）。符合这一条件的最小旋转角,5.3 杆系结构分析实例,问题描
4、述：图5-1所示为一平面桁架，长度L=0.1m，各杆横截面面积均为A=1 10-4m2 ，力P=2000N，计算各杆的轴向力Fa、轴向应力a。,图5-1 平面桁架,实例1 平面桁架分析,思考: 根据静力平衡条件，很容易计算出轴向力Fa、轴向应力a，如表5-1所示。,表5-1 各杆的轴向力和轴向应力,求解步骤：, 创建单元类型 选择Structusral Beam类的2D elastic 3。 定义实常数 定义AREA=1E-4。 定义材料特性 输入Ex=2e11，PRXY=0.3。 创建节点 GUI ：Main MenuPreprocessorModelingCreate NodesIn Ac
5、tive 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 MenuPreprocessorModelingCreateEle MentsAuto NumberedThru Nodes。 弹出拾取窗口，拾取节点1和2，单击拾取窗口的“Apply”按钮，于是在节点
6、1和2之间创建了一个单元。重复以上过程，在节点2和3、1和4、3和4间分别创建单元，建模图形如左图5-2所示。,图5-2 桁架建模图, 施加约束 GUI ：Main Menu SolutionDefine LoadesApply StructuralDisplacementOn Node。 弹出拾取窗口，拾取节点1，单击ok按钮，在弹出的对话框列表中选择“UX、UY”，单击Apply按钮。再次弹出拾取窗口，拾取节点3，单击 ok按钮，在弹出的对话框列表中选择“UY”，单击ok按钮，完成对模型的约束施加。 施加载荷 GUI ：Main Menu SolutionDefine LoadesAppl
7、y StructuralForce/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 Table Define Table。 弹出“Element Tab
8、le 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 Ta
9、bleList Elem Table。 在弹出的对话框的列表中选择FA、 SA，单击ok按钮，即显示出求解结果。与表5-1对照，二者完全一致。,5.4 梁结构分析实例,实例1简支梁载荷工况组合实例,图5-3 简支梁,问题描述及解析解： 图5-3（a）所示为一圆截面简支梁，跨度L=1m，圆截面直径 D=30mm，作用在梁上的集中力为P=1000N,作用点距支座A距a=0.2m，已知弹性模量E=2e11.,由材料力学知识可得： 梁截面的惯性距为,最大挠度,思路分析:当进行线性分析时，简支梁的应力、应变和变形等于如图5-3（b）、（c）所示两个简支梁的结果叠加。如图5-3（b)所示的简支梁结构和载荷
10、均对称于梁的中点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.0
11、3。 定义材料特性 输入Ex=2e11（弹性模量），PRXY=0.3（泊松比）。 创建关键点 GUI：MainMenuPreprocessorModelingCreatKeypoints In 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：MainMenuPrepr
12、ocessorModelingCreateLines LinesStraight Line 弹出拾取窗口，拾取关键点1和2，单击ok按钮。 创建硬点 GUI:MainMenuPreprocessorModelingCreateKeypointsHard PT On line Hard PT by ratio。 弹出拾取窗口，拾取直线，单击ok按钮，在文本框中输入0.4，单击ok按钮。（为了施加集中载荷，作用点处创建了一个硬点，该硬点属于直线）, 划分单元 GUI ： Main MenuPreprocessorMeshingMeshTool。弹出MeshTool对话框，单击“Size Contr
13、ols”区域中的 “Line”后的set按钮，弹出拾取窗口，拾取直线，单击ok按 钮，在NDIV文本框中输入50，单击Mesh按钮，弹出拾取窗口， 拾取直线，单击ok按钮。 显示点、线、单元 GUI：UtilityPlotMulti-Plots。 施加载荷 施加第一个载荷步 施加第一个载荷步的位移载荷 GUI:MainMenuSolutionDefineLoadsApplyStructuralDisplacementOn Keypoint。,对关键点1加载UX、UY方向的载荷，关键点2加载UX、ROTZ方向的载荷。（注意：约束关键点2的转动自由度，相当于在该处时加一个弯矩） 施加第一个载荷步的
14、力载荷 GUI:MainMenuSolutionDefineLoadsApply StructuralForce/MomentOn Keyopints。 弹出拾取窗口，拾取关键点3，单击ok按钮，选择Lab为FY,在VALUE文本框中输入-500，单击ok按钮。 写第一个载荷步文件 GUI ： MainMenuGeneral PostprocLoad CaseRead Load Case。在LCNO文本框中输入1，单击ok按钮。 删除位移载荷 GUI ：MainMenuSolutionDefine LoadsDelete StructuralDisplacementOn Keypoint。 弹
15、出拾取窗口，拾取关键点2，单击ok按钮，弹出一个对话框，然后单击ok按钮。,施加第二个载荷步 施加第二个载荷步的位移载荷 GUI:MainMenuSolutionDefineLoadsApply StructuralDisplacementOn Keypoint。 对关键点2加载UX、UY方向的载荷。 写第二个载荷步文件 GUI ：MainMenuSolutionLoad Step OptsWrite LS File. 在LSNUM文本框中输入2，单击ok按钮。 求解 GUI ： Main Menu SolutionSolveFrom LS Files。 在LSMIN文本框中输入1,在LSMA
16、X文本框中输入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 ： MainMenuGe
17、neral 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 PostprocWr
18、ite Results。 在弹出的对话框的LSTEP文本框中输入9999，单击ok按钮。 查看简支梁有半部分的结果 重复以上- 的过程，在步骤中，拾取菜单MainMenu General PostprocLoad CaseSubtract，对载荷工况进行减运算，可得到简支梁右半部分的结果。（注意：与简支梁的实际情况呈对称关系）。,5.5 平面问题分析实例,问题描述： 图5-2所示为一厚壁圆筒，其内径r1=50mm,外径r2=100mm,作用在内孔上的压力p=10MPa,无轴向压 力，轴向长度视为无穷。要求计算后壁圆筒的径向应力r和切向应力t沿半径r方向的分布。,图5-2厚壁圆筒问题,实例1厚壁
19、圆筒问题,思路:根据材料力学的知识，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
20、.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 Menu
21、SolutionDefine LoadsApply Structural DisplacementOn Lines。 弹出拾取窗口，拾取水平直线边，单击ok按钮，在列表中选择UY，单击Apply按钮，再次弹出拾取窗口，拾取面的垂直直线边，单击ok按钮，在列表中选择UX，单击ok按钮。, 施加载荷 GUI ：MainMenuSolutionDefineLoadsApplyStructuralPressure On Lines。 弹出拾取窗口，拾取面的内侧圆弧边（较短的一条圆弧），单击ok按钮，在VALUE文本框中输入10e6，单击ok按钮。 求解 GUI ：Main Menu SolutionS
22、olveCurrent LS。 单击“Solution Current Load Step”对话框ok按钮。出现“solution is done！”提示时，求解结束，即可查看求解结果。,结果显示： 显示节点 GUI ： Utility MenuPlotNodes。 定义路径 GUI ： MainMenuGeneralPostprocPathOperations Define PathBy Nodes。 弹出拾取窗口，由内向外依次拾取与X轴平行的水平直线边上的所有节点，单击ok按钮，在Name文本框中输入p1，单击ok按钮。 将数据映射到路径上 GUI ： Main MenuGeneral P
23、ostprocPath Operations Map 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 Operations Plot Path ItemOn Graph。 在弹出的对话框中的列表中选择
24、SR、ST，单击ok按钮。即可显示径向应力r和切向应力t 随半径的分布情况。,5.6 空间问题实例分析,实例1 各种坐标系下的扭转问题,设等直圆轴的圆截面直径D=50mm，长度L=120mm，作用在圆轴两端上的转矩Mn=1.5103N M ,求圆柱的最大剪切应力和转角 。,问题描述及解析解： 由材料力学知识可得： 圆截面对圆心的极惯性矩为,圆截面的抗扭截面模量为,圆截面上任意一点的剪应力与该点半径成正比，在圆截面的边缘上有最大值,等直圆轴距离为0.09m的两截面间的相对转角,求解步骤： 创建单元类型 选择Structusral Solid类的Quad 4node 42和Brick 8node
25、45 单元。 定义材料特性 输入Ex=2e11（弹性模量），PRXY=0.3(泊松比）。 创建矩形面 GUI：MainMenuPreprocessorModelingCreatAreas By Dimensions。 在X1、X2文本框中输入0，0.025；在Y1、Y2文本框中输入0，0.12., 划分单元 GUI ：Main MenuPreprocessorMeshingMeshTool 弹出MeshTool对话框，单击“Size Controls”区域中的“Line” 后的set按钮，弹出拾取窗口，拾取矩形面的任一短边，单击ok按钮，在NDIV文本框中输入4，单击Apply按钮，再次弹出拾
26、取窗口，拾取矩形面的任一长边，再次弹出单元尺寸对话框，在NDIV文本框中输入8，单击ok按钮。在Mesh区域，选择单元形状Quad（四边形），选择划分单元的方法为Mapped，单击Mesh按钮，弹出拾取窗口，拾取面，单击ok按钮。 设定挤出选项 GUI：MainMenuPreprocessorModelingOperateExtrude Elem Ext opts。,在VAL1文本框中输入3（挤出段数），选定ACLEAR为Yes（清除矩形面上的单元），单击ok按钮。 由面旋转挤出体 GUI：MainMenuPreprocessorModelingOperateExtrude AreasAbou
27、t 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：U
28、tility 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：MainMenuP
29、reprocessorModelingMove/ModifyRoate Node CSTo Active CS。弹出拾取窗口，单击Pick All按钮。 施加约束 GUI ：MainMenuSolutionDefineLoadsApplyStructuralDisplacementOn Nodes。 弹出拾取窗口，选择Pick All按钮。在弹出的对话框中的Lab2列表框中选择UX，单ok按钮。 选中圆柱面上最上端的所有节点 GUI：Utility MenuSlectEntity。 在弹出的对话框的下拉列表框、文本框、单选按钮中依次选择或输入Nodes、By Location、Z Coordi
30、nates、-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 ：
31、MainMenuSolutionDefineLoadsApplyStructuralDisplacementOn Areas。,弹出拾取窗口，拾取圆柱体下侧底面（由四部分组成），单击ok按钮。在弹出的对话框中的Lab2列表框中选择All DOF，单ok按钮。 求解 GUI ： Main Menu SolutionSolveCurrent LS。 单击“Solution Current Load Step”对话框ok按钮。出现 “solution is done！”提示时，求解结束，即可查看求解结 果。,结果显示： 显示变形 GUI ： Main MenuGeneral PostprocPlot
32、 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 C
33、oordinates、-0.09、From Full，单击Apply按钮。再在下拉列表框、文本框、单选按钮中依次选择或输入Nodes、By Location、YCoordinates、0、Reselect，单击Apply按钮。 列表显示节点位移 GUI ： Main MenuGeneral PostprocList Results Nodal Solution。,在左侧列表中选择DOF Solution，在右侧列表中选择UY，单击ok按钮。 选择单元 GUI：Utility MenuSlectEntity。 在弹出的对话框的下拉列表框、文本框、单选按钮中依次选择或输入Nodes、By Loca
34、tion、Z Coordinates、-0.1，0、From Full，单击Apply按钮。再在下拉列表框、文本框、单选按钮中依次选择或输入Elements、Attachedto、Nodes all、 Reselect，单击Apply按钮。这样做的目的是，在下一步显示应力时，不包含集中力作用点附近的单元，以得到更好的计算结果。 查看结果，用等高线显示剪应力 GUI ：Main MenuGeneral PostprocPlot Results Contour PlotElements Solu。 在Item,Comp两个列表中分别选择Stress、“YZ-shear SYZ”,即可查看结果。,如
35、图所示轴承座结构，尺寸单位为mm。已知沉孔径向内柱面承受外推压力0.0069MPa，轴承孔柱面下部分承受向下压力0.0345MPa，计算轴承的最大应力及其位置。,实例2:轴承座分析,分析： 约束方式： 四个通孔为轴承座安装固定孔，必须给予对称约束； 除此之外，还必须限制轴承座垂直方向的移动，避免该方向的刚体位移。 几何建模 （1）建立底座：VolumesBlockBy Dimension,轮子二维结构如图所示（尺寸单位为英寸），现要分析该轮仅承受绕其中心轴旋转角速度的作用下，轮的受力及其变形情况。,已知角速度：w=525rad/s, 材料属性：弹性模量E30106psi泊松比为0.3密度为0.
36、000731bf-2/in4,实例3轮子分析,思路分析：这是一个空间问题的静力学分析，由于结构的复杂性，不能采用由2D网格拖拉生成3D网格的方式，宜采用自由网格和映射网格相结合的网格生成方式，根据该轮的对称性，在分析式只要分析其中的一部分即可。1）先建立三维几何模型；2）根据三维模型的特点，对三维模型进行分割，使其中一部分模型能够采用映射网格方式划分，其余部分则可以用自由网格划分。因此必须选用20节点六面体单元（Solid95)，它可以进行映射网格划分，同时也可以变成20节点的四面体网格进行自由网格划分。因此能够实现不同网格形状的连接。,3）对划分网格的连接处进行网格转换。即将20节点的四面体
37、自由网格转换为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
38、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
39、.475 X=3.5,Y=3.505 生成圆弧线(如图所示） ArcBy End Kps&Rad R=0.35 生成由圆弧线围成的面（共两处） 面相加（Pick All） 线相加以减少线的数量（可不作）：L17L12L13；L22L10L20,压缩编号操作，并重新显示线 保存结果数据：存为Wheel_Anal_2D （4）通过拖拉生成3D模型 生成轴线的关键点： (0,0)-(0,5) 2D拖拉成3D PreprocessorModelingOperateExtrudeAbout Axis 在Arc length in degree框中输入22.5度，输入NSEG1，即生成的实体由一块体积组成
40、 关闭线号显示 改变视图方向：ISO 结果如下页所示。,（5）生成一个圆柱孔 激活工作平面 旋转工作平面：WorkplaneOffset by increments，在”XY，YZ，ZX Angles”下面的输入栏中输入”0，90“,工作平面的”WZ”将指向上方； 生成一个圆柱体：CreateCylinderSolid Cylinder，出现一个对话框，输入WP X=2.375 WP Y=0，Radius0.45，Depth2.6(注：2.5) 体相减 得到的3D几何模型如图所示 另存盘为：Save as “Wheel_Anal_3D”,（6）生成网格 平移工作平面到关键点9上。用工作平面切分
41、体：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工具条上，
42、选择Shape下的“Tet”四面体网格类型和“Free”自由网格生成方式，单击Mesh，拾取V6，则完成自由网格划分。,注：20节点六面体网格到10节点四面体网格之间的连接，必须在会合面上通过五面体（20节点的退化形式）过渡，即20节点六面体20节点五面体20节点四面体，退化是通过节点重合实现的。,转变的单元类型：MeshingModify MeshChange Tets， 出现一个对话框，在“Change from”后面的框中选择”95 to 92”。即20节点95退化四面体单元转换成10节点四面体单元。以降低求解规模。,显示单元的过渡区域：通过建立选择集来局部观察，UtilitySelec
43、tEntities 弹出对话框，选择Elements和By Element name，输入Soli95，再显示单元UtilityplotElement，得到如图左上所示图形； 再选择全部实体UtilitySelectEverything, UtilitySelectEntities 弹出对话框，选择Elements和By Element name，输入Soli92，再显示单元UtilityplotElement，得到如图左下所示图形； 可以清楚看到，六面体单元到四面体单元通过连接边界面附近的五面体单元过渡而成。 保存单元网格数据： save as”Wheel_Anal_Mesh”,（7）施加载
44、荷并求解 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，而这一较小的轴向力将造
45、成Y向刚体位移。因此只要约束一个角点就足够了。,3）施加旋转角速度： ApplyStructuralIntertiaAngular Veloc Global，在弹出的对话框OMEGY框中输入525 4）保存：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 Polar PATTERNAlternate Symm，DY22.5,
