有限元上机实验报告

1、有限元基础理论报告学院：机电与车辆工程学院班级：制造1班姓名：马龙学号：631224110114任课老师:张继祥二0一五年十一月4题目一：三维托架实体受力分析题目：1、三维托架实体受力分析：托架顶面承受50psi的均匀分布载荷。托架通过有孔的表面固定在墙上，托架是钢制的，弹性模量E=29Xl06psi,泊松比v=0.3.试通过ANSYS俞出其变形图及其托架的vonMises应力分布。题目1的分析：先进行建模，此建模的难点在对V3的构建（既图中的红色部分）。要想构建V3,首先应将A15做出来，然后执行MainMenuPreprocessorModelingOperateBooleansAddVo

2、lumes命令,将所有的实体合并为一个整体。建模后，就对模型进行网格的划分，实行MainMenuPreprocessorMeshingMeshTool,先对网格尺寸进彳T编辑，选0.1,然后点Meshing,Pickall进行网格划分。划分网格后，就可以对模型施加约束并进行加载求解了。施加约束时要注意，由于三维托架只是通过两个孔进行固定，故施加约束应该只是针对两孔的内表面，执行MainMenuSolutionDefineLoadsApplyStructrualDisplacementSymmetryB.COnAreas命令，然后拾取两孔的内表面，单击OKM行了。施加约束后，就可以对实体进行加载

3、求解了，载荷是施加在三维托架的最顶上的表面的，加载后求解运算，托架的变形图。1.指定分析标题：1 .选取菜单路径UtilityMenu|File|ChangeJobname,将弹出ChangeJobname（修改文件名）对话框。2 .在Enternewjobname（输入新文件名）文本框中输入文字“bracket”，为本分析实例的数据库文件名。单击对话框中的“OK按钮，完成文件名的修改。3 .选取菜单路径UtilityMenu|File|ChangeTitle,将弹出ChangeTitle（修改标题）对话框。4 .在Enternewtitle（输入新标题）文本框中输入文字“pressanaly

4、sisofbracketstructure”,为本分析实例的标题名。单击对话框中的“OK按钮，完成对标题名的指定。2.定义单元类型:1 .选取菜单路径MainMenu|Preprocessor|ElementType|Add/Edit/Delete,将弹出ElementTypes（单元类型定义）对话框。单击对话框中的按钮，将弹出LibraryofElementTypes（单元类型库）对话框。2 .在左边的滚动框中单击“StructuralSolid”,选择结构壳单元类型。在右边的滚动框中单击“brick8node45，在对话框中单击“OK按钮，完成对这种单元的定义。3 .设置单元尺寸：执行Ma

5、inMenuPreprocessorMeshingSizeCntrlsManualSizeGlobalSize3.指定材料特性:1 .选取菜单路径MainMenu|Preprocessor|MaterialProps|MaterialModels,将弹出DefineMaterialModelBehavior（材料模型定义）对话框。2 .依次双击Structural,Linear,Elastic和Isotropic,将弹出1号材料的弹性模量EX和泊松比PRXY的定义对话框。3 .在EX文本框中输入2.9E11,PRXY文本框中输入0.3。定义材料的弹性模量为2.9E7N/m2,泊松比为0.3。单

6、击“OK按钮，关闭对话框。ALinearIsotropicPropertiesforMaterialNumber1LinearIsotropicMaterialPropertiesforMaterialMmuber1TemperaturesEXPRXYOKAddTenperatureDeleteTeirperature4 .在De巾neMaterialModelBehavior（材料模型定义）对话框中，选取路径MaterialExit,完成对材料模型的定义。5 .单击ANSYS12.0的ANSYSToolbar（工具条）上的“SAVE按钮，保存数据库文件4.建立托架的有限元模型:1 .选取路径

7、MainMenu|Preprocessor|Modeling|Create|Volume|Block|ByDimensions,将弹出CreateBlockbyDimensions（根据坐标创建体）对话框。在对话框输入:X1,X2X-coordinates:-1,1Y1,Y2Y-coordinates:-1.5,1.5Z1,Z2Z-coordinates：0,1/8然后单击“APPLY”按钮，再次在对话框输入:X1,X2X-coordinates:-1,1Y1,Y2Y-coordinates:1.5,1.625Z1,Z2Z-coordinates:0,32 .选取路径MainMenu|Prep

8、rocessor|Modeling|Operate|Booleans|Add|Volumes将弹出AddVolumes（体相力口）对话框，在对t框单击“PICKALL”按钮完成体相加操作。3 .在菜单栏PlotCtrls选择Number,在PlotNumberingcontrol对话框中显示keypointNumber和lineNumber的操作，选取路径MainMenu|Preprocessor|Modeling|Create|Lines|Lines|StrainghtLine将弹出CreateStraightLine对话框，然后选择关键点5、13生成L13直线，单击“OK”按钮完成操作。4

9、 .选取路径MainMenu|Preprocessor|Modeling|Create|Area|Arbitrary|ByLines将弹出对话框CreateAreaByANSYSAREAS冲W 18 201513:37:47Lines然后选择直线L1,L9,L13,L20,L24,L25,单击“OK”即可生成面A4。（注：以上两步需连续进行，实验时点线不同时显示会造成L13丢失，如随后无法显示面，可在菜单栏中Plot选择Area操作显示。）5 .选取路径MainMenu|Preprocessor|Modeling|Operate|Extrude|Areas|AlongNormal弹出preaa

10、 ana 1 yis of bracket 3triictUEeExtrudeAreaby选择面A4单击OK按钮，将弹出另一个对话框，在Lengthextrusion项输入-1/8”,单击OK按钮。然后选取路径MainMenu|Preprocessor|Modeling|Operate|Booleans|Add|Volumes将弹出AddVolumes（体相加）对话框，在对t框单击“PICKALL”按钮完成体相加操作。6 .选取路径MainMenu|Preprocessor|Modeling|Create|Areas|Circle|SolidCircle将弹出SolidCircularArea

11、对话框，在WPx,WPy,Radius项分别输入0,-0.5,0.25单击Apply,再次输入0,0.5,0.25然后单击“OK”按钮。然后选取路径MainMenu|Preprocessor|Modeling|Operate|Extrude|Areas|AlongNormal弹出对话框ExtrudeAreaby选择其中一个圆面单击“Apply”按钮，将弹出另一个对话框，在Lengthextrusion项输入05,单击OK”按钮。对另一个圆面做相同操作。7 .选取路径MainMenu|Preprocessor|Modeling|Operate|Booleans|Subtract|Volumes将

12、弹出SubtractVolumes（体相减）对话框，选择两个圆柱体所在的体，单击“Apply”按钮，然后选择两个圆柱体，单击“OK”按钮，完成体相减操作。5.网格划分:1.选取路径MainMenu|Preprocessor|ANSYSAREASNUMpress analyslB of bracket striiat-ureSTOV IE 201S13:49 3:02Modeling|Meshing|MeshTool将弹出MeshTool对话框，单击“Mesh”按钮，弹出另一对话框，再次单击“PICKALL”按钮完成网格划分。6.施加约束，载荷并求解：1 .选取菜单路径MainMenu|Prep

13、rocessor|Loads|DefineLoads|Apply|Structural|Displacement|OnAreas,将会弹出拾取对话框，选择两圆孔，单击对话框中的OK按钮。完成施加约束操作。2 .选取菜单路径MainMenu|Solution|DefineLoads|Apply|Structural|Pressure|OnAreas将会弹出拾取对话框，选择面A26,A12单击“OK按钮，弹出另一对话框，在LoadPRESvalue项输入“50”单击“OK完成施加载荷操作。3 .选取菜单路径MainMenu|Solution|Solve|CurrentLS弹出SokeCurrent

14、LoadStep对话框，单击OK开始求解，求解结束后，关闭相应对话框。4 .选取菜单路径MainMenu|GeneralPostproc|PlotResults|ContourPlot|NodalSolu将弹出ContourNodalSolutionData对话框，选择NodalSolution|DOFSolution|Displacementvectorsum将得到位移图：33托架位移变形图分析：由图可看出，由于载荷的作用，托架上面板明显变形了，变形最严重的就是红色部分，这是因为该部分离托板较远而且悬空，没有任何物体支撑并分担载荷，所以比较容易变形甚至断裂。当托架在承受载荷时，应尽量避免该部

15、承受较大载荷。将弹出 Contour Nodal5选取菜单路径MainMenu|GeneralPostproc|PlotResults|ContourPlot|NodalSoluSolutionData对话框，选择NodalSolution|Stress|vonMisesStress将得到应力图:托架应力分布图分析：如图所示，应力主要集中在两孔，尤其是下面的一个孔，出现了最大应力。换句话说，两孔是该托板的薄弱环节，容易出现裂痕。我们在应用托架的时候，应当注意采取一些设施，以便减缓其应力集中。特别是在施加载荷时，绝对不能够超过托架所能承受的极限，否则必然导致托板断裂，甚至引发安全事故。题目二：圆

16、柱齿轮模态分析题目：对图1-1所示的一简化齿轮模型进行模态分析，要求确定齿轮的低阶固有频率。然后扩展模态，求出各阶模态的相对应力值、相对应变值和相对位移值。另外齿轮的齿根部分是该圆周的2/3,齿端部分是该圆周的2/7,齿数是24,齿顶与齿根看做圆弧，齿面看做直线，齿轮厚0.2m,弹性模量为2x1011Pa,泊松2比为0.3,密度为7.8x103Kg/m3。其中内圆孔R=1.0,齿根圆R=1.8,齿顶圆R=2.10齿轮简化模型1 .定义工作文件名和工作标题(1) 定义工作文件名：执行UtilityMenuChangeJobname命令，弹出【ChangeJobnme】对话框。输入Cylinder

17、gear”并选择【Newloganderrorfiles复选框，单击OK按钮。(2)定义工作标题：执行UtilityMenuFileChangeTitle命令，弓t出ChangeTitle对话框。输入ThemodelAnalysisofCylindergear,单击OK按钮。2 .定义单元类型和材料属性(1) 设置单元类型：执行MainMenuPreprocessorElementTypeAdd/Edit/Delete命令，弹出【ElementTypes对话框。单击Apply”按钮，弹出如图1-2所示的【LibraryofElementTypes对话框。选择Solid和10node92”选项。

18、(2)设置材料属T执行MainMenuPreprocessorMaterialPropsMaterialModels,弹出如图1-3的DefineMaterialModelsBehavior】对话框。双击【MaterialModelsAvailable】列表框中的StructuralLinearElasticIsotropic”选项，弓t出【LinearIsotropicMaterialPropertiesForMaterialNumber1对话框。输入EX=2e11,PRXY=0.3”,单击OK按钮。双击Density选项，在弹出的【DensityofMaterialNumber1对话框输入

19、密度7800”,单击“OK”按钮。执行MaterialExit命令，完成材料属性的设置。DefineMaterialModelsBehavior3 .创建圆柱齿轮模型(1)激活工作坐标系：执行UtilityMenuWorkplaneChangeActivecstoGlobalCylinder。(2)定义关键点：执行MainMenuPreprocessormodelingCreatekeypointsInActiveCs。弹出如图1-4所示的【CreateKeypointsinActiveCoordinationSystem对话框，在NPT和X,Y,ZLocationinactiveCS”中输入

20、“1”和“0,0,0”、“2”单击“Apply，而后，依次输入“2.1,90,0”、“3”和“2.1,87.8571,0”、“4”和“1.8,85,0”、“5”和“1.8,82.5,0，单击“OK”按钮。【CreateKeypointsinActiveCoordinationSystem(3)生成线：执行MainMenuPreprocessormodelingCreatelinesInActiveCoord,选取关键点2和3、4和5。再执行MainMenuPreprocessormodelingCreatelinesStraightlines,选取关键点3和4、5和1。(4)激活工作坐标系：执

21、行UtilityMenuWorkplaneChangeActivecstoGlobalCartesian。(5)镜像线：执行MainMenuPreprocessormodelingReflectlines,选取上面建立的直线，弹出【ReflectLines对话框，单击“OK”按钮。【ReflectLines(6)粘接线：执行MainMenuPreprocessormodelingOperateBooleansGlueLines,选取所有直线。(7)生成面：执行MainMenuPreprocessormodelingCreateAreasArbitraryBylines,选取所有直线即可。(8)

22、激活工作坐标系：执行UtilityMenuWorkplaneChangeActivecstoGlobalCylinder。(9)圆周整列齿轮廓线：执行MainMenuPreprocessormodelingcopyAreas,选取之前创建的面，弹出【CopyAreas对话框，在Numberofcopies和Y-offsetinactiveCS中输入24和15,单击OK按钮。ACopyArea5AGENCopyAreasITIMENumberofcopies-includingoriginaDXX-ottsetinactiveC$DYY-offsetinactiveCSDZZ-off?etina

23、ctiveCSIQNCKeypointincrementNOELEMItemstobecopiedQK|ApplyCaned IHelpCopyLinesMIEASTYPE NUMANSYSNOVIS201514:53:58themodelanalysisofcylinder齿轮面(10)面相加：执行MainMenuPreprocessormodelingOperateBooleansAddAreas,选取所有面。(11)面拉伸：执行MainMenuPreprocessormodelingOperateExtrudeAreasAlongNormal，在弹出的【ExtrudeAreasAlong

24、Normal对话框中输入拉伸长度“0.2并确定得齿轮图。齿轮图(12)创建圆柱并想减：执行MainMenuPreprocessormodelingCreateVolumesCylinderOK”。SolidCylinder,弹出【SolidCylinder对话框，在“Radius”和“Depth”中输入“1”和“-0.2”，单击执行MainMenuPreprocessormodelingOperateBooleansSubstractVolumes,先后选取齿轮和圆柱，想减后得取终模型。ANSYSTOLUME5TYPEr IS 2015 15:02:154 .网格划分执行MainMenuPre

25、processorMeshingSizeCntrlsManualSizeGlobalsize,在【GlobalElementSizes对话框中的SIZEElementedgelength”输入“0.1”。而后执行网格划分，得有限元模型。ri八GlobalElementSizesESIZEGlobalelementsizesanddivisions(appliesonlytaunsized11lines)SIZEElementedgelength01NDIVNo.ofelementdivisions0-(usedonlyifelementedgelengthBSIZErisblankorzero

26、)OKICancelHelp亘【GlobalElementSizesANSYS1ELEMENT 旨the model analysis of cylinderNOV 16 201515:C:15有限元模型5 .模态分析(1)创建分析：执行MainMenuSolutionAnalysisTypeNewAnalysis,在弹出【NewAnalysis】对话框中选取Modal”。【NewAnalysis(2)选择菜单:执行MainMenuSolutionAnalysisTypeAnalysisOptions,弹出【ModalAnalysis】对话框，在“No.ofmodestoextract中输入5

27、。弹出【BlockLanczosMethod对话框，在StareFreq和EndFreq输入“1”和“100000”。【Modal Analysis A Block Lanczos MethodMO DO PT Options for Block Uhczoe Modal AnalysisFREQB Start Freq (initial shift)FREQE End FrequencyN rm key Normalize mode shapes0SolutionDefineLoadsApplyStructuralDisplacementOnAreas,弹出【ApplyU,ROTonArea

28、s拾取菜单，选取内圆柱面。(4)将模态扩展并求解：执行MainMenuSolutionLoadsStepOptsExpansionPassSingleExpandModes,弹出【ExpandModes对话框，在No.ofmodestoexpand中输入5。S30.001CancelHelpAExpandModesMXPANDExpandModesNMODENqhofmodestoexpandFREQB.FREQEFrequencyrang*ElcalcCalculateelemresults?StGNlFSignificantThreshold-onlyvalidforSPRSandDDAM

29、OKI【ExpandModes(5)进行结果运算：执行MainMenuSolutionSolveCurrentLS,开始求解。6,结果处理(1)执行MainMenuGeneralPosprocResultsSummary，观察求得的5阶频率。其中低阶固有频率为：163.41HZ。ASET,LISTCommandfile*INDEXOFDATASETSONRESULTSFILE*SETTIME/FREQLOADSTEPSUBSTEPCUMULATIVE1163.411112165.351223165-371334176.261445176.29155(2)执行MainMenuGeneralPos

30、procReadResultsFirstSet。(3)执行MainMenuGeneralPosprocPlotResultsContourPlotNodalSolu,在弹出的图1-17ContourNodalSolutionData对话框中，分别选择DOFSolution”下面的Displacementvectorsum、Stress”下面的vonMisesstress和ElasticStrain”下面的vonMisesElasticStrain,分别得到在第1阶频率下的相对位移图、相对应力图、和相对应变图ANSYS009903themodelanalysisofcylinder1NODALS

31、OLUTIONHWIB201515:31:54,003961.007923.001981,0-05942.011884.015845.013g5.0178265TEP=1SUB=1TREQ-163.408USUM(fiWS)RSY5=0DMX=-0176265MX=-017826DisplacementvectorsumNOEUkLSOLCniOiiSIEP=1SUB=1FREQiea-iQSSEQV(fiVG)DHX=.017826SHN-.199E+07SMX=.134E+1GANSYSNOVIE201515i32i4e,199E+07,29能+09-595E+09-891E+09-119

32、E+10fc150E+09,446Eg,743E+09-104E+10.134E+10themodelanalysisofcylindervonMisesstressANSYSHODALSOLUTION.006621themodelanalysisofcylinder】；CVLS201515:34i38,102E-04-75IE-035TEP=1SUB=1FREgl石3J08EFELEQVAVIjJDMX=.017626SMIT=.102E-04SMX=,006681+0QT493B002375.004453.QOS944.002234003716,005159ivonMisesElastic

33、Strain(4)依次执行(1)(3)步骤，得各阶的相对位移、相对应力和相对应变，结果如下图所示。并将所有阶次的结果统计整理如表1-1所示。NODAL SOLUTIONEIEP=1 SUB =2 FREQ=165.35 OSUH (AVG) RSYSfl DMX =.024973 5MX =.024973ANSYSNOV IE 2015IS:57:050.00555.011099.016649.022190-002775,008324.013874,01M24,024973themodelanalysisofcylinderDisplacementvectorsum二阶vonMisesstre

34、ssMODAL SOLUTIONSTEP-1 SUB 2FREQ-165-35 EPELEQV (AVG) Di必=.024973 SMN =.352E-05 SMX =,009447ANSYSNOV 18 2015 15:58:29.3S2E-05.002102.004201.006299.006396.OO 10 S3,003151.00 52 5.00734S,009447themodelanalysisofcylinder二阶vonMisesElasticStrainWOnAL SOLUTIOHSTEP=1 SUB FREQsl5a 374USUM (AVG) RSYS-d 函=602

35、4915 SMX =.024915ANSYSNOT 18 2015 15:59:270.005537.011073.01661.022147,002768.008305.013842019379.024915themodelanalysisofcylinder三阶DisplacementvectorsumANSYSSTEP=1 SUB =3 FREQ=165.374SEQM (AVG)WK -.024915 SMN 654546 smx -.iaeE+ioNOV 16 201516:00:12654546.41fiE+09.它 35E+09.125E-F10.167E+10.209E+09-6

36、27E+09.104E+10.146E+10-1BSE+101LOCAL5OL1II0EIt-hemodelanalysisofcylinder三阶vonMisesstressNODALSOLUTION5TEP=1SUB=3FREQ=165*374EPELEQVAVE)DMC=.024915SMN=.345E-G5SMX=,009396ANSYSMWIB2015IfisOliDZ.34SE-05002091s.0D4178.D0C25B008352001047,003134.005221007305g939themodelanalysisofcylinder三阶vonMisesElasticS

37、trainANSYS.011229.01683，。口M22.01403E,019651themodelanalysisofcylinder1NCDftLSOLUTIONKOV1&2Q1516;Q115B,005614.0028075TEP=1SUB=4FKEQ=17hSUSCTMAVG)R5S=DIM=.025265SHX=,02525四阶DisplacementvectorsumNODALSOLUTION5TEP=1SOB=4FREQ-176+26SEQV(AVG)DMX=.02S25SMN=-127E+07SMX=a192E+10.127E+07215EM9ANSYSMOVIE201516

38、SO2：:5B.428E+09.55E+09+12SE+10.171E+10.641E+09.107E+10.150E+10*192E+10themodelanalysisocylinder四阶vonMisesstress四阶vonMisesElasticStrainNODALSOLUTION5TEP=1SUB=5FRQ=176+.2e6USUM(AVG)RSYS-QDMX=.025252SWX=-025252ANSYSKOVie201516:05:22.005612,011223.01B35,02246.002806.003417,014029.01954filechangejobname命

39、令，在弓t出的【changejobname对话框中输入pinandpinhole,选中【newloganderrorfiles复选框，单击OK按钮。（2）定义工作标题：执行utilitymenufilechangetitle命令，在弹出的【changetitle】对话框中输入“thecontactanalysisofpinandhole”,单击“OK”按钮。（3）重新显示：执行utilitymenuplotreplot命令。2、定义单元类型及材料属性（1）定义单元类型：执行mainmenupreprocessorelementtypeadd/eidt/delete命令，弹出【elementty

40、pe】对话框，单击add按钮，弹出【libraryofelementtypes对话框。在左右列表框中分别选择structuralsolid和brick8model85选项，单击OK按钮。（2）定义材料属性：执行mainmenupreprocessormaterialpropsmaterialmodels命令，在弹出的【definematerialmodelbehavior】窗口。双击structurallinearelasticisotropic,弹出【linearisotropicmaterialpropertiesformaterialnumber1】对话框。在【ex】和【prxy】文本框

41、中分别输入“36e6”和“0.3”，单击OK按钮。执行materialexit命令。3、建立几何模型（1）生成插座：执行mainmenupreprocessormodelingcreateblockbydimension命令，弹出【createblockbydimension对话框，输入参数“X2=Y2=2,Z1=2.5,Z2=3.5。单击OK按钮。（2）改变视图角度：打开pan-zoom-rotate工具栏：执行utilitymenuplotctrlspanzoomrotate命令，单击iso按钮。（3）生成圆柱体：执行mainmenupreprocessormodelingcreatevo

42、lumescylinderbydimensions命令，弹出【createcylinderbydimension】对话框，输入参数“RAD=0.49,Z1=2.5,Z2=3.5”.单击OK按钮。（4）生成插孔：执行mainmenupreprocessormodelingoperatebooleanssubtractvolumes命令，弹出拾取对话框，拾取编号为v1的长方体，单击OK按钮，然后拾取编号为v2的圆柱体，单击OK按钮，得到带插孔的插座。（5）生成插销：执行mainmenupreprocessormodelingcreatevolumescylinderbydimensions命令，弹

43、出【createcylinderbydimension】对话框，输入参数“rad1=0.5,z1=2,z2=4.5”。单击OK生成一个圆柱体。（6）打开体积编号：执行utilitymenupoltctrlsnumbering命令，弹出【plotnumberingcontrol】对话框，选中【volumenumbers】复选框，单击OK按钮。（7）切分模型的1/4：由于问题的对称条件，只需要完整插销和底座模型1/4来进行分析。设置工作平面，以备切分模型。执行utilitymenuworkplanewpsetting命令，弹出【wpsetting】对话框，选中单选按钮”catesian”和”gri

44、dandtriad”,单击OK按钮。显示工作平面，执行utilitymenuworkplanedisplayworkingplane命令。旋转工作平面，绕Y轴正方向旋转90.切分模型的1/2，删除1/2模型。还原工作平面位置。旋转工作平面，绕x轴正向旋转90.切分模型的1/4，删除模型的1/4.AREASANSYiTYPE NUMWOV ie 2015 15:50:16nn andan hole(8)保存几何模型：单击工具条中的save-db,保存几何模型。4、划分有限元网格(1)执行utilitymenupreprocessormeshingmeshtool命令，寸T开【meshtool对话

45、框。(2)设置插销单元网格密度：在【meshtool对话框中的sizecontrol域单击linesset按钮，弹出拾取对话框拾取插销前端的水平和垂直直线，单击OK按钮，弓t出【elementsizesonpicklines对话框；输入“ndiv=3，取消复选框SIZE,NDIVcanbechanged,单击OK按钮确定。(3)设置插座单元网格密度：在【meshtool】对话框中sizecontrol域单击linesset按钮，弹出拾取对话框,拾取插座前端的曲线，单击OK按钮，在弹出的【elementsizesonpickedlines对话框，输入NDIV=4,取消复选框【SIZE.NDIVc

46、anbechanged，单击OK按钮确定。(4)设置单元形状和网格划分方法；在【meshtool对话框中mesh列表中选择volumes选项，单击“hex/wedge和sweep，选择autosrc/trg选项(5)单击“sweep”按钮，弹出拾取对话框，单击pickall,进行网格划分。(6)单击close,关闭【meshtool对话框。(7)执行utilitymenuplotctrlsstylesizeandshape命令，弹出sizeandshapd对话框，在【facets/elementedge列表中选择2facets/edge选项，单击OK按钮。5、建立接触单元(1) 执行mainm

47、enumodelingcreatecontactpair命令，弓t出【contactmanager,对话框。(2)单击最左边的按钮，启动【contactwizard(3)在“targetsurface”域选择areas选项，在“targettype”域选择flexible选项，单击“picktarget按钮，弹出拾取对话框，拾取插座上与插销接触的曲面，单击OK按钮。(4)单击next按钮，接触向导，定义接触曲面。(5)在“contactsurface”域选择areas选项，在“contactelementtype域选择surface-to-surface选项。单击“pickcontact”按钮

48、，弹出拾取对话框，拾取插销上与插座接触的曲面，单击OK。(6)单击create”按钮，单击finish按钮。(7)关闭【contactmanager】对话框。6、定义位移约束(1)执行utilitymenuplotareas命令，重新绘制面积。(2)在面上施加对称约束：执行mainmenusolutiondefinedefineloadsapplystructuraldisplacementsymmetricB.COnArea命令，弹出拾取对话框，拾取插座和插销被切出来的四个面。单击OK。(3)执行mainmenusolutiondefinedefineloadsapplystructural

49、displacementonareas命令，弹出拾取对话框，拾取插座白左侧面，单击OK按钮，弹出【applyU,rotonareas对话框，固定该面所有自由度。7、求解装配预应力(1)设置求解选项：执行mainmenusolutionanalysistypesolcontrol命令，弓t出【solutioncontrols对话框，设置求解选项。(2) 在“analysistype”歹U表中选择largedisplacementstatic选项，在“timeatendofloadstep域输入100,在automatictimestepping列表中选择off选项，在numberofsubsteps”域中车入“1,单击OK确定。(3) 求解：执行mainmenusolutionsolvecurrentls命令。(4) 绘制装配应力图：执行mainmenugeneralpostprocplotresultscontourplotnodalsolution命令，在弹出的对话框中选择stressvonmisesSEQV,单击OK按钮，得到应力分