求助:我的模型怎么划分不了网格?

yxb345

!工作环境
fini
/cle
/filn,welding          !工作名
!进入前处理
/prep7
/units,si
et,1,shell63
et,2,shell63                   !焊条

!定义模型的各几何参数
boxb=0.4    !箱形梁宽度
boxh=0.7    !箱形梁高度
Ib=0.4        !工字形斜撑宽度
Ih=0.4      !工字形斜撑腹板高度
boxl=9.0    !箱形梁长度
IL=6.0      !工字形斜撑长度
Ith=3.14159/4      !工字形斜撑角度
Ibt=0.021        !工字形斜撑翼缘厚度
Iht=0.021        !工字形斜撑腹板厚度
boxbt=0.02    !箱形梁宽度方向厚度
boxht=0.02    !箱形梁高度方向厚度
Jc=0.7
h=0.2
!s=sin(Ith)
!c=cos(Ith)
!设置材料参数
r,1,0.021                                                        !实常数
r,2, 0.021                                                           !焊条
mptemp,1,30,500,1000,1500,2000
mpdata,ex,1,1,2.06e11,1.70e11,0.90e11,0.2e11,0.01e11             !弹性模量
mpdata,nuxy, 1,1, 0.3,0.3,0.3,0.3,0.3                            !泊松比
mpdata,dens,1,1,7800,7800,7800,7800,7800
mpdata,alpx, 1,1,1.06e-5,1.06e-5,1.06e-5,1.06e-5,1.06e-5         !热膨胀系数
mpdata,kxx,1,1,16.3,16.3,16.3,16.3,16.3                          !热传导系数
mpdata,c,1,1,502,502,502,502,502                                 !比热
tb,bskin,1,5                                       !双线性随动硬化弹塑性本构关系

!建立几何模型
k,1,0,0,0
k,2,0,-boxh,0
k,3,-boxl/2,0,0
k,4,-boxl/2,-boxh,0
k,5,-boxl/2+Ib*sin(Ith)/2,boxl*tan(Ith)/2+Ib*cos(Ith)/2,0
k,6,-boxl/2,boxl*tan(Ith)/2,0
k,7,-boxl/2-Ib*sin(Ith)/2,boxl*tan(Ith)/2-Ib*cos(Ith)/2,0
k,8,0.25*(-boxl/2)+Ib*sin(Ith)/2,0.25*boxl*tan(Ith)/2+Ib*cos(Ith)/2,0
k,9,0.25*(-boxl/2),0.25*(boxl*tan(Ith)/2),0
k,10,0.25*(-boxl/2)-Ib*sin(Ith)/2,0.25*boxl*tan(Ith)/2-Ib*cos(Ith)/2,0
k,11,-Jc*cos(Ith)+Ib*sin(Ith)/2,Jc*sin(Ith)+Ib*cos(Ith)/2,0
k,12,-Jc*cos(Ith),Jc*sin(Ith),0
k,13,-Jc*cos(Ith)-Ib*sin(Ith)/2,Jc*sin(Ith)-Ib*cos(Ith)/2,0
k,14,-boxl/2+Il*cos(Ith)+Ib*sin(Ith)/2,boxl*tan(Ith)/2+Ib*cos(Ith)/2-Il*sin(Ith),0
k,15,-boxl/2+Il*cos(Ith),boxl*tan(Ith)/2-Il*sin(Ith),0
k,16,-boxl/2+Il*cos(Ith)-Ib*sin(Ith)/2,boxl*tan(Ith)/2-Ib*cos(Ith)/2-Il*sin(Ith),0
!k,69,-boxl/2+Il*cos(Ith)-Ib*sin(Ith)/2-(h-boxl*tan(Ith)/2-Ib*cos(Ith)/2-Il*sin(Ith))/tan(Ith),h,0
!k,70,-boxl/2+Il*cos(Ith)-Ib*sin(Ith)/2-(h-boxl*tan(Ith)/2-Ib*cos(Ith)/2-Il*sin(Ith))/tan(Ith),h,0.4
k,17,0,(Ib+0.3)*cos(Ith)/2+Jc*sin(Ith),0
k,18,-Jc*cos(Ith)+(Ib+0.3)*sin(Ith)/2,Jc*sin(Ith)+(Ib+0.3)*cos(Ith)/2,0
k,19,-Jc*cos(Ith)-(Ib+0.3)*sin(Ith)/2,Jc*sin(Ith)-(Ib+0.3)*cos(Ith)/2,0
k,20,-Jc*cos(Ith)-(Ib+0.3)*sin(Ith)/2,0,0
k,21,-Jc*cos(Ith)+Ib*sin(Ith)/2,Jc*sin(Ith)+Ib*cos(Ith)/2,0
k,22,-Jc*cos(Ith),Jc*sin(Ith),0
k,23,-Jc*cos(Ith)-Ib*sin(Ith)/2,Jc*sin(Ith)-Ib*cos(Ith)/2,0
k,24,-boxl/2+Il*cos(Ith)+Ib*sin(Ith)/2,boxl*tan(Ith)/2+Ib*cos(Ith)/2-Il*sin(Ith),0
k,25,-boxl/2+Il*cos(Ith),boxl*tan(Ith)/2-Il*sin(Ith),0
k,26,-boxl/2+Il*cos(Ith)-Ib*sin(Ith)/2,boxl*tan(Ith)/2-Ib*cos(Ith)/2-Il*sin(Ith),0
k,27,-boxl/2+Il*cos(Ith)-Ib*sin(Ith)/2,0,0
k,28,0,boxl*tan(Ith)/2+Ib*cos(Ith)/2-Il*sin(Ith),0
k,29,-Jc*cos(Ith)-(Ib+0.3)*sin(Ith)/2,0,0
k,30,-boxl/2+Il*cos(Ith)-Ib*sin(Ith)/2,0,0
k,31,0,0,0
k,32,-Jc*cos(Ith)-(Ib+0.3)*sin(Ith)/2,-boxh,0
k,33,-boxl/8,-boxh,0
k,34,-boxl/8,0,0
kgen,2,1,34,1,0,0,boxb,34
a,3,34,20,27,1,2,32,33,4,3
a,5,6,7,10,13,16,15,14,11,8
a,23,19,29,30,26
a,17,18,21,24,28
a,24,25,26,30,31,28
a,23,22,21,24,25,26
a,6,40,43,46,49,15,12,9,6   !腹板
a,3+34,34+34,20+34,27+34,1+34,2+34,32+34,33+34,4+34,3+34
a,5+34,6+34,7+34,10+34,13+34,16+34,15+34,14+34,11+34,8+34
a,23+34,19+34,29+34,30+34,26+34
a,17+34,18+34,21+34,24+34,28+34
a,24+34,25+34,26+34,30+34,31+34,28+34
a,23+34,22+34,21+34,24+34,25+34,26+34

a,3,37,35,1
a,4,38,36,2
l,33,34

!划分网格
/pnum,line,1
allsel,all
lesize,1,0.05
lesize,8,0.05
lesize,9,0.05
lesize,87,0.05
type,1
mat,1
mshape,0,3d
mshkey,1

cxb305

我看看再说！