声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1560|回复: 11

[PAM-VAOne] va one 同时存在壳单元和体单元的隔声量计算问题

[复制链接]
发表于 2016-10-8 17:43 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x

图片1 (2).png

图片2.png

图片3.png

图片4.png

就是上图的一个由体单元和壳单元组成的结构,在算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

回复
分享到:

使用道具 举报

发表于 2016-10-9 08:55 | 显示全部楼层
接触的问题吧
 楼主| 发表于 2016-10-9 09:48 | 显示全部楼层

没明白,壳单元和体单元接触了啊

点评

接触设置的不对就算不了  详情 回复 发表于 2016-10-9 13:39
发表于 2016-10-9 13:39 | 显示全部楼层
csuluhao 发表于 2016-10-9 09:48
没明白,壳单元和体单元接触了啊

接触设置的不对就算不了
 楼主| 发表于 2016-10-9 17:19 | 显示全部楼层
jiangwj 发表于 2016-10-9 13:39
接触设置的不对就算不了

我是新手,怎么去设置这个接触呢?

点评

这个我也不清楚 因为我用ansys的时候 如果接触设置不好 就算不了 所以猜测你这个也是  详情 回复 发表于 2016-10-10 08:49
发表于 2016-10-10 08:49 | 显示全部楼层
csuluhao 发表于 2016-10-9 17:19
我是新手,怎么去设置这个接触呢?

这个我也不清楚  因为我用ansys的时候 如果接触设置不好  就算不了  所以猜测你这个也是
 楼主| 发表于 2016-10-10 09:00 | 显示全部楼层
jiangwj 发表于 2016-10-10 08:49
这个我也不清楚  因为我用ansys的时候 如果接触设置不好  就算不了  所以猜测你这个也是

ansys的结构模态可以算,本身这个模型就是在ansys中导出来的

点评

那我也不知道了 你现在解决了吗  详情 回复 发表于 2016-10-14 10:44
发表于 2016-10-14 10:44 | 显示全部楼层
csuluhao 发表于 2016-10-10 09:00
ansys的结构模态可以算,本身这个模型就是在ansys中导出来的

那我也不知道了  你现在解决了吗
 楼主| 发表于 2016-10-15 16:29 | 显示全部楼层
jiangwj 发表于 2016-10-14 10:44
那我也不知道了  你现在解决了吗

问题解决了,不过用的是其他方法,应该还是网格质量的问题,

点评

哦 啥方法啊 方便说吗  详情 回复 发表于 2016-10-17 08:58
发表于 2016-10-17 08:58 | 显示全部楼层
csuluhao 发表于 2016-10-15 16:29
问题解决了,不过用的是其他方法,应该还是网格质量的问题,

哦  啥方法啊  方便说吗
 楼主| 发表于 2016-10-21 16:56 | 显示全部楼层
jiangwj 发表于 2016-10-17 08:58
哦  啥方法啊  方便说吗

我把体单元的下表面的fe face作为双层加筋板的上下板,然后就可以算了,也不知道为啥。
发表于 2016-10-25 09:01 | 显示全部楼层
csuluhao 发表于 2016-10-21 16:56
我把体单元的下表面的fe face作为双层加筋板的上下板,然后就可以算了,也不知道为啥。

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-9 01:28 , Processed in 0.072862 second(s), 28 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表