声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2129|回复: 14

[噪声治理] vaone隔声量的计算

[复制链接]
发表于 2016-9-29 20:53 | 显示全部楼层 |阅读模式

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

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

x
在vaone里只有壳单元的隔声量能算,只有体单元的能算,但是碰到同时有体单元和壳单元的模型就报错是怎么回事啊?
回复
分享到:

使用道具 举报

发表于 2016-9-30 08:28 | 显示全部楼层
你说的体单元是指结构模型是体单元吗?VAOne支持体单元的计算吗?
发表于 2016-9-30 08:46 | 显示全部楼层
你把错误贴出来
回复 支持 1 反对 0

使用道具 举报

 楼主| 发表于 2016-10-8 16:05 | 显示全部楼层
倾冷月 发表于 2016-9-30 08:28
你说的体单元是指结构模型是体单元吗?VAOne支持体单元的计算吗?

是的


就是上图的一个由体单元和壳单元组成的结构,在算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-8 16:06 | 显示全部楼层
Agoni 发表于 2016-9-30 08:46
你把错误贴出来



就是上图的一个由体单元和壳单元组成的结构,在算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-8 17:00 | 显示全部楼层
支持体网格计算,关键是不是体单元与壳单元连接的问题?
 楼主| 发表于 2016-10-8 18:41 | 显示全部楼层
sujunshou 发表于 2016-10-8 17:00
支持体网格计算,关键是不是体单元与壳单元连接的问题?

va one 里体单元与壳单元是怎么连接的呢?
发表于 2016-10-9 08:52 | 显示全部楼层
接触没设置好?
 楼主| 发表于 2016-10-9 10:26 | 显示全部楼层

啥意思?不是本身就是在一起的吗,用了autoconnect也不行
发表于 2016-10-10 09:00 | 显示全部楼层
错误提示是什么
 楼主| 发表于 2016-10-10 09:03 | 显示全部楼层

the cousmic running failed
 楼主| 发表于 2016-10-10 09:10 | 显示全部楼层

the cosmic running failed
发表于 2016-10-10 11:11 | 显示全部楼层
共用节点?

点评

啥意思 模型还不是一个整体?  详情 回复 发表于 2016-10-11 08:42
 楼主| 发表于 2016-10-10 17:28 | 显示全部楼层
发表于 2016-10-11 08:42 | 显示全部楼层

啥意思  模型还不是一个整体?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-23 10:42 , Processed in 0.061767 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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