该单元是solid45单元怎么才能划成六面体网格呢?
我做的一个桩土分析的模型命令流如下:
finish$/clear$/prep7
k,1,0,0$k,2,0,-1$k,3,-0.2,-1.3$k,4,0,-1.5$k,5,0,-2$k,6,-0.2,-2.3$k,7,0,-2.5$k,8,,-3.5$k,9,0.25,-3.5$k,10,0.25,0
*do,i,1,9$l,i,i+1$*enddo
l,10,1
al,all
vrotat,1,,,,,,10,9,90
et,1,solid45$et,2,170$et,3,173
keyopt,3,9,0$keyopt,3,12,2
mp,ex,1,2.6e10$mp,prxy,1,0.167$mp,dens,1,2500
mp,ex,2,5e6$mp,prxy,2,0.42$mp,dens,2,1900
tb,dp,2$tbdata,,19000,31,29
mp,mu,3,0.2$r,1$r,2$r,3
blc4,0.25,,-2,-7,-2$vsbv,2,1
wpoff,-0.3$wprota,,,90$vsbw,all$wpcsys,-1
wpoff,,,-0.55$vsbw,all$wpcsys,-1
wpoff,,-3.6$wprota,,90$vsbw,all$wpcsys,-1$numcmp,all
vatt,2,2,1$esize,0.2$mshape,0,3d$mshkey,1$vmesh,1,7
vatt,2,2,1$smrtsize,5$mshape,1,3d$mshkey,0$vmesh,8
k,50,,0$k,51,,-1$k,52,-0.2,-1.3$k,53,,-1.5$k,54,,-2$k,55,-0.2,-2.3$k,56,,-2.5$k,57,,-3.5$k,58,0.25,-3.5$k,59,0.25,0
*do,i,50,58$l,i,i+1$*enddo
l,59,50
lsel,s,,,79,88
!$lsel,a,,,43,46$lsel,a,,,53,54$lsel,a,,,19,38,19
al,all
vrotat,45,,,,,,58,59,-90
vatt,1,1,1$smrtsize,5$vmesh,all
asel,s,,,46,52$nsla,s,1
real,3$type,2$mat,3$esurf,all$allsel,all
asel,s,,,1,7$nsla,s,1$real,3,type,3$esurf,all
asel,s,loc,z,0$asel,a,loc,x,0.25$da,all,symm
asel,s,loc,x,-1.75$asel,a,loc,z,-2$da,all,all
asel,s,loc,y,-7$da,all,uy
asel,s,,,54$sfa,all,1,pres,50e4$allse,all
/solu$antype,static$nlgeom,on$autots,on
acel,,,9.8$time,1$nsubst,100$outres,all,all
solve$finish
模型做好之后算时总是超限,所以想全部划成六面体单元,但是由于形状有点怪异现在不知道怎么才能运行!
|