是的
就是上图的一个由体单元和壳单元组成的结构,在算structure FE的时候报错 如图示。但是单一把体单元和壳单元分析的时候又不会出错。
Apdl程序:/CLEAR,START
LENGTH=2.25, !长度
WIDTH=1.800, !宽度
CHICUNE=0.03, !单元尺寸
FLOOR=60, !
!*************单元材料属性定义**********************
/PREP7
ET,1,181 !选取单元类型1
MP,EX,1,7.1e10 !定义材料的弹性模量
MP,DENS,1,2700 !定义材料的密度
MP,PRXY,1,0.33 !定义材料的泊松比
ET,2,SOLID45 !选取单元类型1
MP,EX,2,7.1e10 !定义材料的弹性模量
MP,DENS,2,2700 !定义材料的密度
MP,PRXY,2,0.33 !定义材料的泊松比
MP,EX,3,340000 !定义材料的弹性模量
MP,DENS,3,1300 !定义材料的密度
MP,PRXY,3,0.47 !定义材料的泊松比
MP,EX,4,2.3e9 !定义材料的弹性模量
MP,DENS,4,1100 !定义材料的密度
MP,PRXY,4,0.49 !定义材料的泊松比
MP,EX,5,24832 !定义材料的弹性模量
MP,DENS,5,48 !定义材料的密度
MP,PRXY,5,0.2 !定义材料的泊松比
R,1,0.0018, ! 定义厚度(1.8)肋板
R,2,0.0025, ! 定义厚度(2.5)上层板
R,3,0.003, ! 定义厚度(3) 下层板
!**************清空上一次数据*********************
/PREP7
DATE=LSEL(1), ! 查询否存在线 不存在返回0
*IF,DATE,EQ,0,THEN
*ELSE
ALLSEL,ALL, ! 选择所有实体
DDELE,ALL ! 删除所有约束
ACLEAR,ALL ! 清除所有面网格
ADELEL,ALL ! 删除所有面
LDELEL,ALL ! 删除所有线
KDELEL,ALL ! 删除所有点
*ENDIF
*set,FLNumP(1) ! 清空数组FLNumP 控制点编号
*set,FLNumberA(1) ! 清空数组FLNumberA 面编号
!**************绘制点、线、面************************
P1=0, ! 定义下表面Z值
P2=0.03, ! 定义上表面Z值
*DIM,FLNumP,,(FLOOR+3),1,1 ! 定义数组FLNumP
FLNumP(1)=0,
FLNumP(2)=0,
*DO,i1,1,FLOOR,1,
a1=WIDTH/FLOOR*i1,
FLNumP(i1+2)=a1,
*ENDDO
FLNumP(FLOOR+3)=FLNumP(i1+2),
*DO,i2,1,(FLOOR+3),1
C=MOD(i2,2)
*IF,C,EQ,1,THEN
K,i2,FLNumP(i2,1,1),LENGTH,P2
*ELSE
K,i2,FLNumP(i2,1,1),LENGTH,P1
*ENDIF
*ENDDO
*GET, KPNUM1, KP, 0, NUM, MAXD, !获取最大的控制点编号 存入变量KPNUM1
*DO,i1,1,(KPNUM1-1),1
a=i1,
b=a+1,
lstr,a,b, !依次连接关键点形成肋板截面线 例如 1-2 2-3 3-4........
*ENDDO
*DO,i1,1,(KPNUM1-2),2
a=i1,
b=a+2,
lstr,a,b, !依次连接关键点形成上表面截面线 例如 1-3 3-5 5-7........
*ENDDO
*DO,i1,2,(KPNUM1-2),2
a=i1,
b=a+2,
lstr,a,b, !依次连接关键点形成下表面截面线 例如 2-4 4-6 6-8........
*ENDDO
*GET, Lnum1, LINE, 0, NUM, MAXD,!获取最大的线编号Lnum1
K,(KPNUM1+1),0,0,0.03 !绘制引导线所需另一个点。
lstr,1,(KPNUM1+1), !绘制拉伸引导线 1号点与KPNUM1+1号点连接
*GET, Lnum2, LINE, 0, NUM, MAXD,!获取最大的线编号Lnum2,及引导线编号
*DO,i2,1,Lnum1,1 !除了引导线共有Lnum1根线,依次拉伸线形成面
b=i2,
ADRAG,b,,,,,,Lnum2,,,,,, !将既有线沿一定路径拖拉成面。
*ENDDO
NUMCMP,ALL
NUMMEG,ALL
wpof,,,-LENGTH/2,
wprot,0,0,90,
SHULIANG=30,
cc2=WIDTH/FLOOR !30
SHULIANG2=WIDTH/cc2 !60
*do,i1,1,SHULIANG2-1,1,
asel,all
asel,s,loc,z,-0.001,0.001,
asel,a,loc,z,29.9/1000,30.1/1000,
wpof,,,cc2
ASBW,ALL
*enddo
CSYS,0
WPAVE,0,0,0
CSYS,0
ALLSEL,ALL
asel,s,loc,z,30/1000,30.1/1000, !1
asel,r,loc,x,0 ,1.8,
EXTOPT,ESIZE,1,0,
VEXT,ALL, , ,0,0,0.003,,,,
asel,all
*do,i1,1,FLOOR/6,1,
ALLSEL,ALL
asel,s,loc,z,33/1000,33.1/1000, !2
asel,r,loc,x,(0.18*i1-0.15) ,(0.18*i1),
EXTOPT,ESIZE,1,0,
VEXT,ALL, , ,0,0,0.003,,,,
*enddo
asel,all
*do,i1,1,FLOOR/6,1, !3
ALLSEL,ALL
asel,s,loc,z,33/1000,33.1/1000,
asel,r,loc,x,(0.18*(i1-1)) ,(0.18*i1-0.15),
EXTOPT,ESIZE,1,0,
VEXT,ALL, , ,0,0,0.003,,,,
*enddo
asel,all
ALLSEL,ALL !4
asel,s,loc,z,36/1000,36.1/1000,
asel,r,loc,x,0 ,1.8,
EXTOPT,ESIZE,1,0,
VEXT,ALL, , ,0,0,0.003,,,,
asel,all
ALLSEL,ALL !5
asel,s,loc,z,39/1000,39.1/1000,
asel,r,loc,x,0 ,1.8,
EXTOPT,ESIZE,1,0,
VEXT,ALL, , ,0,0,0.015,,,,
asel,all
ALLSEL,ALL
vsel,s,loc,z,30/1000,33/1000,
vsel,r,loc,x,0 ,1.8,
VATT,2,,2,
MSHAPE,0,
ESIZE,CHICUNE
VMESH,all
NUMCMP,ALL
NUMMEG,ALL
asel,all
*do,i1,1,FLOOR/6,1,
ALLSEL,ALL
vsel,s,loc,z,33/1000,36/1000,
vsel,r,loc,x,(0.18*i1-0.15) ,(0.18*i1),
VATT,2,,2,
MSHAPE,0,
ESIZE,CHICUNE
VMESH,all
*enddo
NUMCMP,ALL
NUMMEG,ALL
asel,all
*do,i1,1,FLOOR/6,1,
ALLSEL,ALL
vsel,s,loc,z,33/1000,36/1000,
vsel,r,loc,x,(0.18*(i1-1)) ,(0.18*i1-0.15),
VATT,3,,2,
MSHAPE,0,
ESIZE,CHICUNE
VMESH,all
*enddo
NUMCMP,ALL
NUMMEG,ALL
asel,all
ALLSEL,ALL
vsel,s,loc,z,36/1000,39/1000,
vsel,r,loc,x,0 ,1.8,
VATT,4,,2,
MSHAPE,0,
ESIZE,CHICUNE
VMESH,all
NUMCMP,ALL
NUMMEG,ALL
asel,all
ALLSEL,ALL
vsel,s,loc,z,39/1000,54/1000,
vsel,r,loc,x,0 ,1.8,
VATT,5,,2,
MSHAPE,0,
ESIZE,CHICUNE
VMESH,all
NUMCMP,ALL
NUMMEG,ALL
asel,all
ALLSEL,ALL
asel,s,loc,z,0.001,0.029,
AATT,1,1,1,
ESIZE,CHICUNE
MSHKEY,0
AMESH,all
NUMCMP,ALL
NUMMEG,ALL
ALLSEL,ALL
asel,s,loc,z,29.9/1000,30.1/1000,
AATT,1,2,1,
ESIZE,CHICUNE
MSHKEY,0
AMESH,all
NUMCMP,ALL
NUMMEG,ALL
ALLSEL,ALL
asel,s,loc,z,-0.001,0.001,
AATT,1,3,1,
ESIZE,CHICUNE
MSHKEY,0
AMESH,all
NUMCMP,ALL
NUMMEG,ALL
asel,all
|