面与面接触实例:插销拨拉问题分析定义单元类型Element/add/edit/delete定义材料属性Material Props/Material ModelsStructural/Linear/Elastic/Isotropic定义材料的摩擦系数建立几何模型Modeling/Create/Volumes/Block/By Dimensions X1=Y1=0,X2=Y2=2,Z1=2.5,Z2=3.5Modeling/Create/Volumes/Cylinder/By DimensionsModeling/Operate/Booleans/Subtract/V olumes先拾取长方体,再拾取圆柱体。
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/Sol’s 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/Sol’s 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,,5 Wpstyle,,,,,,,,1Wpro,,,90Wsbw,allVdele,4,,,1Vdele,6,,,1Wpcsys,-1,0Wpro,,90Vsbw,allVdele,p51x,,,1Wpcsys,-1,0Wpstyle,,,,,,,,0(2)定义单元类型、材料模型和网格划分Et,1,solid185 Mptemp,,,,,,,,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.88860905221e-031 R,3Real,3Et,2,170Et,3,174R,3,,,0.1,0.1,0Rmore,,,1.0e20,0.0,1.0 Rmore,0.0,0,1.0,0.5Rmore,0,1.0,1.0,0.0,,1.0 Keyopt,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,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,1Plnsol,s,eqv,0,1.0 Save(6)求解拨拉过程AplotNsel,s,loc,z,4.5Finish/solAntype,restD,all,,1.7,,,,uz,,,,, Allsel,allNsrbst,100,10000,10 Outres,erase Outres,all,all Outots,1Time,200/status,soluSolveFinish(7)结果后处理/expand,4,polar,half,,90 Eplot/post1Set,,,1,,120Esel,s,ename,,174 Eplot/efacet,1Avprin,0Plnsol,cont,pres,0,1.0 Allsel,allSet,2,last,1/efacet,1Avprin,0Plnsol,s,eqv,0,1.0 Plns,s,eqvAndata,0.5,,0,0,0,1,1,1 finish/exit,all。