APDL裂纹的建模

ban_length=0.15ban_width=0.1crack_width=0.0008crack_length=0.01pi=2*asin(1)theta=pi/4EX_ban=2.1e11NUXE_ban=0.3DENS_ban=7800pre_ban=-1.4e8kdes=12/prep7/units,siet,1,plane82KEYOPT,1,3,3r,1,0.08mp,ex,1,EX_banmp,nuxy,1,NUXY_banmp,dens,1,DENS_bank,1,0,-(ban_length)/2,0k,2,0,-(crack_width)/2,0 k,3,crack_length*(cos(theta),crack_length*(sin(theta),0 k,4,0,(crack_width)/2,0k,5,0,(ban_length)/2,0 k,6,(ban_width)/2,(ban_length)/2,0 k,7,(ban_width)/2,-(ban_length)/2,0 k,8,-(crack_length*(cos(theta),-crack_length*(sin(theta),0 k,9,-(ban_width)/2,(ban_length)/2,0 k,10,-(ban_width)/2,-(ban_length)/2,0l,1,2l,2,3l,3,4l,4,5l,5,6l,6,7l,7,1al,alllsel,invelsel,s,line,1l,2,8l,4,8lsel,a,line,4l,5,9l,9,10l,10,1al,allasel,s,area,1aatt,1,1,1 kscon,3,0.002,1,kdes lsel,s,line,5 lesize,all,10 lsel,s,line,7lplolesize,all,10allslsel,s,line,1 lesize,all,10,0.2 lsel,s,line,4 lesize,all,10,5 allslsel,s,line,6 lesize,all,20 allslesize,3,7lesize,2,9allsamesh,1allsasel,s,area,2aatt,1,1,1 kscon,8,0.002,1,kdes lsel,s,line,10 lesize,all,10 lsel,s,line,12lplolesize,all,10allslsel,s,line,11lesize,all,20allslesize,8,7lesize,9,9allsamesh,2alls/soluantype,0nsel,s,loc,y,-ban_length/2 nsel,a,loc,y,banength/2 nsel,a,loc,x,banength/2 nsel,a,loc,x,-ban_length/2 sf,all,pres,pre_banallssolvesave