va one 同时存在壳单元和体单元的隔声量计算问题
就是上图的一个由体单元和壳单元组成的结构，在算structure FE的时候报错 如图示。但是单一把体单元和壳单元分析的时候又不会出错。Apdl程序：/CLEAR,START LENGTH=2.25, !长度WIDTH=1.800, !宽度CHICUNE=0.03, !单元尺寸FLOOR=60, !!*************单元材料属性定义**********************/PREP7ET,1,181 !选取单元类型1MP,EX,1,7.1e10 !定义材料的弹性模量MP,DENS,1,2700 !定义材料的密度MP,PRXY,1,0.33 !定义材料的泊松比 ET,2,SOLID45 !选取单元类型1MP,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) 下层板!**************清空上一次数据*********************/PREP7DATE=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 ! 定义数组FLNumPFLNumP(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-22-3 3-4........ *ENDDO *DO,i1,1,(KPNUM1-2),2 a=i1, b=a+2, lstr,a,b, !依次连接关键点形成上表面截面线 例如 1-33-5 5-7........*ENDDO *DO,i1,2,(KPNUM1-2),2 a=i1, b=a+2, lstr,a,b, !依次连接关键点形成下表面截面线 例如 2-44-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,,,,,, !将既有线沿一定路径拖拉成面。*ENDDONUMCMP,ALL NUMMEG,ALL wpof,,,-LENGTH/2,wprot,0,0,90, SHULIANG=30,cc2=WIDTH/FLOOR !30SHULIANG2=WIDTH/cc2!60 *do,i1,1,SHULIANG2-1,1,asel,allasel,s,loc,z,-0.001,0.001,asel,a,loc,z,29.9/1000,30.1/1000,wpof,,,cc2 ASBW,ALL*enddoCSYS,0WPAVE,0,0,0 CSYS,0 ALLSEL,ALLasel,s,loc,z,30/1000,30.1/1000, !1asel,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,ALLasel,s,loc,z,33/1000,33.1/1000, !2asel,r,loc,x,(0.18*i1-0.15) ,(0.18*i1),EXTOPT,ESIZE,1,0, VEXT,ALL, , ,0,0,0.003,,,, *enddoasel,all *do,i1,1,FLOOR/6,1, !3ALLSEL,ALLasel,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,,,,*enddoasel,all ALLSEL,ALL !4asel,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 !5asel,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,allNUMCMP,ALLNUMMEG,ALLasel,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*enddoNUMCMP,ALLNUMMEG,ALLasel,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*enddoNUMCMP,ALLNUMMEG,ALLasel,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,allNUMCMP,ALLNUMMEG,ALLasel,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,allNUMCMP,ALLNUMMEG,ALLasel,all ALLSEL,ALL asel,s,loc,z,0.001,0.029,AATT,1,1,1, ESIZE,CHICUNEMSHKEY,0 AMESH,allNUMCMP,ALLNUMMEG,ALL ALLSEL,ALL asel,s,loc,z,29.9/1000,30.1/1000,AATT,1,2,1, ESIZE,CHICUNEMSHKEY,0 AMESH,allNUMCMP,ALLNUMMEG,ALL ALLSEL,ALL asel,s,loc,z,-0.001,0.001,AATT,1,3,1, ESIZE,CHICUNEMSHKEY,0 AMESH,allNUMCMP,ALLNUMMEG,ALLasel,all
接触的问题吧 jiangwj 发表于 2016-10-9 08:55
接触的问题吧
没明白，壳单元和体单元接触了啊 csuluhao 发表于 2016-10-9 09:48
没明白，壳单元和体单元接触了啊
接触设置的不对就算不了 jiangwj 发表于 2016-10-9 13:39
接触设置的不对就算不了
我是新手，怎么去设置这个接触呢？ csuluhao 发表于 2016-10-9 17:19
我是新手，怎么去设置这个接触呢？
这个我也不清楚因为我用ansys的时候 如果接触设置不好就算不了所以猜测你这个也是 jiangwj 发表于 2016-10-10 08:49
这个我也不清楚因为我用ansys的时候 如果接触设置不好就算不了所以猜测你这个也是
ansys的结构模态可以算，本身这个模型就是在ansys中导出来的 csuluhao 发表于 2016-10-10 09:00
ansys的结构模态可以算，本身这个模型就是在ansys中导出来的
那我也不知道了你现在解决了吗 jiangwj 发表于 2016-10-14 10:44
那我也不知道了你现在解决了吗
问题解决了，不过用的是其他方法，应该还是网格质量的问题， csuluhao 发表于 2016-10-15 16:29
问题解决了，不过用的是其他方法，应该还是网格质量的问题，
哦啥方法啊方便说吗 jiangwj 发表于 2016-10-17 08:58
哦啥方法啊方便说吗
我把体单元的下表面的fe face作为双层加筋板的上下板，然后就可以算了，也不知道为啥。
csuluhao 发表于 2016-10-21 16:56
我把体单元的下表面的fe face作为双层加筋板的上下板，然后就可以算了，也不知道为啥。
{:3_61:}
页:
[1]