csuluhao 发表于 2016-9-29 20:53

vaone隔声量的计算

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

倾冷月 发表于 2016-9-30 08:28

你说的体单元是指结构模型是体单元吗?VAOne支持体单元的计算吗?

Agoni 发表于 2016-9-30 08:46

你把错误贴出来

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-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,,,,,,    !将既有线沿一定路径拖拉成面。
*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-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,,,,,,    !将既有线沿一定路径拖拉成面。
*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

共用节点?

csuluhao 发表于 2016-10-10 17:28

billions1943 发表于 2016-10-10 11:11
共用节点?

???

Catsayer 发表于 2016-10-11 08:42

billions1943 发表于 2016-10-10 11:11
共用节点?

啥意思模型还不是一个整体?
页: [1]
查看完整版本: vaone隔声量的计算