ansys面与面接触分析实例

面与面接触实例：插销拨拉问题分析定义单元类型Element/add/edit/delete定义材料属性MaterialProps/MaterialModelsStructural/Linear/Elastic/Isotropic定义材料的摩擦系数建立几何模型Modeling/Create/Volumes/Block/ByDimensionsX1=Y1=0，X2=Y2=2，Z1=2.5，Z2=3.5Modeling/Create/Volumes/Cylinder/ByDimensionsModeling/Operate/Booleans/Subtract/Volumes先拾取长方体，再拾取圆柱体。Modeling/Create/Volumes/Cylinder/ByDimensions划分掠扫网格Meshing/SizeCntrls/ManualSize/Lines/PickedLines拾取插销前端的水平和垂直直线，输入NDIV=3再拾取插座前端的曲线，输入NDIV=4PlotCtrls/Style/SizeandShape，在Facets/elementedge列表中选择2facets/edge建立接触单元Modeling/Create/Contactpair，弹出ContactManager对话框，如图所示。单击最左边的按钮，启动ContactWizard（接触向导），如图所示。单击PickTarget，选择目标面。选择接触面定义位移约束施加对称约束，DefineLoads/Apply/Structural/Displacement/SymmetricB.C/OnAreas，选择对称面。再固定插座的左侧面。设置求解选项AnalysisType/Sol’sControl求解：Solve/CurrentLS绘制装配应力图GeneralPostproc/PlotResults/ContourPlot/NodalSolution，选择Stress/vonMisesstress求解拨拉过程选择Z=4.5处的所有节点。DefineLoads/Apply/Structural/Displacement/OnNodes，弹出ApplyU，ROTonNodes拾取框，单击PickAll按钮，选择UZ，在Displacementvalue输入1.7Select/EverythingAnalysisType/Sol’sControlSolve/CurrentLS结果后处理扩展模型：Style/SymmetryExpansion/Priodic/CyclicSmmetry，在弹出的对话框中选择1/4DihedralSym选择GeneralPostproc/ReadResults/Bytime/frequency，在TIME域输入120。选择插销中与插座接触的单元，在SelectEntities中选择Element，在列表中选择ByElementname，再ElementName域输入174Plot/ElementsGeneralPostproc/PlotResults/ContourPlot/NodalSolution，在对话框中选择Contact/ContactPressure读入载荷步2结果。ReadResults/ByLoadStep绘制拨拉过程的应变变化动画PlotCtrls/Animate/OverResults，弹出如图所示的对话框。命令流操作：(1)建立几何模型/filename,bolt/title,bolt_pullinganalysis/PREP7Block,-2,2,-2,2,2.5,3.5/view,1,1,1,1/ang,1/rep,fastCylind,0.49,,2.5,3.5,0,360Vsbv,1,2Cylind,0.5,,2,4.5,0,360/pnum,volu,1Wpstyle,0.05,0.1,-1,1,0.003,0,0,,5Wpstyle,,,,,,,,1Wpro,,,90Wsbw,allVdele,4,,,1Vdele,6,,,1Wpcsys,-1,0Wpro,,90Vsbw,allVdele,p51x,,,1Wpcsys,-1,0Wpstyle,,,,,,,,0（2）定义单元类型、材料模型和网格划分Et,1,solid185Mptemp,,,,,,,,Mptemp,1,0Mpdata,ex,1,,36e6Mpdata,prxy,1,,0.3Lesize,4,,,3,,,,,0Lesize,10,,,3,,,,,0Lesize,18,,,4,,,,,0Vsweep,all/shrink,0/eshape,0.0/Efacet,2/ratio,1,1,1/cformat,32,0（3）定义接触单元/com,contactpaircreation-startCm,_nodecm,nodeCm,_elemcm,elemCm,_kpcm,kpCm,_linecm,lineCm,_areacm,areaCm,_volucm,volu/gsav,cwz,gsav,,tempMp,mu,1,0.2Mat,1Mp,emis,1,7.88860905221e-031R,3Real,3Et,2,170Et,3,174R,3,,,0.1,0.1,0Rmore,,,1.0e20,0.0,1.0Rmore,0.0,0,1.0,0.5Rmore,0,1.0,1.0,0.0,,1.0Keyopt,3,4,0Keyopt,3,5,0Nropt,unsymKeyopt,3,7,0Keyopt,3,8,0Keyopt,3,10,1Keyopt,3,11,0Keyopt,3,12,0Keyopt,3,2,0Keyopt,3,5,0Asel,s,,,23Cm,_target,areaType,2Nsla,s,1Esln,s,0EsurfCmsel,s,_elemcmAsel,s,,,27Cm,_contact,areaType,3Nsla,s,1Esln,s,0EsurfAllselEsel,allEsel,s,type,,2Esel,a,type,,3Esel,r,real,,3/psymb,esys,1/pnum,type,1/num,1EplotEsel,allEsel,s,type,,2Esel,a,type,,3Esel,r,real,,3Cmsel,a,_nodecmCmdel,_nodecmCmsel,a,_elemcmCmdel,_elemcmCmsel,s,_kpcmCmdel,_kpcmCmsel,s,_linecmCmdel,_linecmCmsel,s_areacmCmdel,_areacmCmsel,s,_volucmCmdel,_volucm/gres,cwz,gsavCmdel,_targetCmdel,_contact/com,contactpaircreation-end(4)定义位移约束Finishaplot/soluFlst,2,4,5,orde,4Fitem,2,3Fitem,2,7Fitem,2,11Fitem,2,14Da,p51x,symmFlst,2,1,5,orde,1Fitem,2,19Da,p51x,all,0(5)求解装配预应力Antype,0Nlgeom,1Nsubst,1,0,0Autots,0Time,100/status,soluSolveFinish/post1/efacet,1Avprin,0Plnsol,s,eqv,0,1.0Save(6)求解拨拉过程AplotNsel,s,loc,z,4.5Finish/solAntype,restD,all,,1.7,,,,uz,,,,,Allsel,allNsrbst,100,10000,10Outres,eraseOutres,all,allOutots,1Time,200/status,soluSolveFinish(7)结果后处理/expand,4,polar,half,,90Eplot/post1Set,,,1,,120Esel,s,ename,,174Eplot/efacet,1Avprin,0Plnsol,cont,pres,0,1.0Allsel,allSet,2,last,1/efacet,1Avprin,0Plnsol,s,eqv,0,1.0Plns,s,eqvAndata,0.5,,0,0,0,1,1,1finish/exit,all