ansys面与面接触分析实例

面与面接触实例：插销拨拉问题分析定义单元类型Element/add/edit/delete定义材料属性Material Props/Material ModelsStructural/Linear/Elastic/Isotropic定义材料的摩擦系数建立几何模型Modeling/Create/Volumes/Block/By DimensionsX1=Y1=0，X2=Y2=2，Z1=2.5，Z2=3.5Modeling/Create/Volumes/Cylinder/By DimensionsModeling/Operate/Booleans/Subtract/Volumes先拾取长方体，再拾取圆柱体。Modeling/Create/Volumes/Cylinder/By Dimensions划分掠扫网格Meshing/Size Cntrls/ManualSize/Lines/Picked Lines拾取插销前端的水平和垂直直线，输入NDIV=3再拾取插座前端的曲线，输入NDIV=4PlotCtrls/Style/Size and Shape，在Facets/element edge列表中选择2 facets/edge建立接触单元Modeling/Create/Contact pair，弹出Contact Manager对话框，如图所示。单击最左边的按钮，启动Contact Wizard（接触向导），如图所示。单击Pick Target，选择目标面。选择接触面定义位移约束施加对称约束，Define Loads/Apply/Structural/Displacement/Symmetric B.C/On Areas，选择对称面。再固定插座的左侧面。设置求解选项Analysis Type/Sols Control求解：Solve/Current LS绘制装配应力图General Postproc/Plot Results/Contour Plot/Nodal Solution，选择Stress/von Mises stress求解拨拉过程选择Z=4.5处的所有节点。Define Loads/Apply/Structural/Displacement/On Nodes，弹出Apply U，ROT on Nodes拾取框，单击Pick All按钮，选择UZ，在Displacement value输入1.7Select/EverythingAnalysis Type/Sols ControlSolve/Current LS结果后处理扩展模型：Style/Symmetry Expansion/Priodic/Cyclic Smmetry，在弹出的对话框中选择1/4 Dihedral Sym选择General Postproc/Read Results/By time/frequency，在TIME域输入120。选择插销中与插座接触的单元，在Select Entities中选择Element，在列表中选择By Element name，再Element Name域输入174Plot/ElementsGeneral Postproc/Plot Results/Contour Plot/Nodal Solution，在对话框中选择Contact / Contact Pressure 读入载荷步2结果。Read Results/By Load Step绘制拨拉过程的应变变化动画PlotCtrls/Animate/Over Results，弹出如图所示的对话框。命令流操作：(1)建立几何模型/filename,bolt/title,bolt_pulling analysis/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,contact pair creation-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.e-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,contact pair creation-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/pos