设为首页收藏本站联系我们
切换到宽版

声振论坛

 找回密码
 我要加入

QQ登录

只需一步，快速开始

声振论坛»论坛 › 基础理论 › 声学基础及其分析软件 › vaone隔声量的计算
返回列表 发新帖
查看: 1883|回复: 14

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

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

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

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

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

使用道具 举报

倾冷月
发表于 2016-9-30 08:28 | 显示全部楼层
你说的体单元是指结构模型是体单元吗？VAOne支持体单元的计算吗？
回复 支持 反对

使用道具 举报

Agoni
发表于 2016-9-30 08:46 | 显示全部楼层
你把错误贴出来
回复 支持 1 反对 0

使用道具 举报

csuluhao
 楼主| 发表于 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


回复 支持 反对

使用道具 举报

csuluhao
 楼主| 发表于 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


回复 支持 反对

使用道具 举报

sujunshou
发表于 2016-10-8 17:00 | 显示全部楼层
支持体网格计算，关键是不是体单元与壳单元连接的问题？
回复 支持 反对

使用道具 举报

csuluhao
 楼主| 发表于 2016-10-8 18:41 | 显示全部楼层
sujunshou 发表于 2016-10-8 17:00
支持体网格计算，关键是不是体单元与壳单元连接的问题？

va one 里体单元与壳单元是怎么连接的呢？
回复 支持 反对

使用道具 举报

Edinburgh
发表于 2016-10-9 08:52 | 显示全部楼层
接触没设置好？
回复 支持 反对

使用道具 举报

csuluhao
 楼主| 发表于 2016-10-9 10:26 | 显示全部楼层
Edinburgh 发表于 2016-10-9 08:52
接触没设置好？

啥意思？不是本身就是在一起的吗，用了autoconnect也不行
回复 支持 反对

使用道具 举报

Raspberry
发表于 2016-10-10 09:00 | 显示全部楼层
错误提示是什么
回复 支持 反对

使用道具 举报

csuluhao
 楼主| 发表于 2016-10-10 09:03 | 显示全部楼层
Raspberry 发表于 2016-10-10 09:00
错误提示是什么

the cousmic running failed
回复 支持 反对

使用道具 举报

csuluhao
 楼主| 发表于 2016-10-10 09:10 | 显示全部楼层
Raspberry 发表于 2016-10-10 09:00
错误提示是什么

the cosmic running failed
回复 支持 反对

使用道具 举报

billions1943
发表于 2016-10-10 11:11 | 显示全部楼层
共用节点？

点评

Catsayer
啥意思 模型还不是一个整体？  详情 回复 发表于 2016-10-11 08:42
回复 支持 反对

使用道具 举报

csuluhao
 楼主| 发表于 2016-10-10 17:28 | 显示全部楼层
billions1943 发表于 2016-10-10 11:11
共用节点？

???
回复 支持 反对

使用道具 举报

Catsayer
发表于 2016-10-11 08:42 | 显示全部楼层
billions1943 发表于 2016-10-10 11:11
共用节点？

啥意思  模型还不是一个整体？
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-4-24 15:36 , Processed in 0.097852 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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