24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1078  |  回复: 3

魔石乱枫

新虫 (小有名气)

[求助] 求大佬指点迷津,下面已建模,求大佬给写个能显示电压的后处理已有1人参与

求大佬指点迷津,下面已建模,但是压电片的电压读取始终不行,求大佬给写个能显示电压的后处理
/prep7                                                    !进入前处理器
/title,Transient and harmonic analyses of a piezoelectric circuit
/nopr
L1=0.068                                                  !金属基片的尺寸
W1=0.010
H1=0.00064

L2=0.028                                                  !上压电片尺寸
W2=0.007
H2=0.0003

L3=0.028                                                  !片尺寸
W3=0.0015
H3=0.00064

L4=0.02969
W4=0.007
H4=0.00064

L5=0.01031                                                 !shang质量块尺寸
W5=0.007
H5=0.0049

BLOCK,-W3,0,-H3/2,0,0,L3                     !建立金属基片的模型
BLOCK,0,W2/2,-H3/2,0,0,L3
BLOCK,W2/2,W2,-H3/2,0,0,L3
BLOCK,W2,W2+W3,-H3/2,0,0,L3
BLOCK,-W3,0,-H3/2,0,L3,L3+L4
BLOCK,0,W2/2,-H3/2,0,L3,L3+L4
BLOCK,W2/2,W2,-H3/2,0,L3,L3+L4
BLOCK,W2,W2+W3,-H3/2,0,L3,L3+L4
BLOCK,-W3,0,-H3/2,-H3,0,L3                     !建立金属基片的模型
BLOCK,0,W2/2,-H3/2,-H3,0,L3
BLOCK,W2/2,W2,-H3/2,-H3,0,L3
BLOCK,W2,W2+W3,-H3/2,-H3,0,L3
BLOCK,-W3,0,-H3/2,-H3,L3,L3+L4
BLOCK,0,W2/2,-H3/2,-H3,L3,L3+L4
BLOCK,W2/2,W2,-H3/2,-H3,L3,L3+L4
BLOCK,W2,W2+W3,-H3/2,-H3,L3,L3+L4
BLOCK,0,W2/2,0,H2+H2,0,L3         !建立上压电片模型
BLOCK,W2/2,W2,0,H2+H2,0,L3
BLOCK,0,W2/2,-H3,-H3-H2,0,L3         !建立下压电片模型
BLOCK,W2/2,W2,-H3,-H3-H2,0,L3
BLOCK,-W3,0,0,H2,L3+L4,L3+L4+L5                    !建立shangya质量块模型
BLOCK,0,W2/2,0,H2,L3+L4,L3+L4+L5
BLOCK,W2/2,W2,0,H2,L3+L4,L3+L4+L5
BLOCK,W2,W2+W3,0,H2,L3+L4,L3+L4+L5
BLOCK,-W3,0,H2,H2+0.0046,L3+L4,L3+L4+L5                    !建立shangyashang质量块模型
BLOCK,0,W2/2,H2,H2+0.0046,L3+L4,L3+L4+L5
BLOCK,W2/2,W2,H2,H2+0.0046,L3+L4,L3+L4+L5
BLOCK,W2,W2+W3,H2,H2+0.0046,L3+L4,L3+L4+L5
BLOCK,-W3,0,-H3/2,0,L3+L4,L3+L4+L5                             !建立中质量块模型
BLOCK,0,W2/2,-H3/2,0,L3+L4,L3+L4+L5
BLOCK,W2/2,W2,-H3/2,0,L3+L4,L3+L4+L5
BLOCK,W2,W2+W3,-H3/2,0,L3+L4,L3+L4+L5
BLOCK,-W3,0,-H3/2,-H3,L3+L4,L3+L4+L5                            !建立中质量块模型
BLOCK,0,W2/2,-H3/2,-H3,L3+L4,L3+L4+L5
BLOCK,W2/2,W2,-H3/2,-H3,L3+L4,L3+L4+L5
BLOCK,W2,W2+W3,-H3/2,-H3,L3+L4,L3+L4+L5
BLOCK,-W3,0,-H3,-H3-H2,L3+L4,L3+L4+L5                   !建立中下质量块模型
BLOCK,0,W2/2,-H3,-H3-H2,L3+L4,L3+L4+L5
BLOCK,W2/2,W2,-H3,-H3-H2,L3+L4,L3+L4+L5
BLOCK,W2,W2+W3,-H3,-H3-H2,L3+L4,L3+L4+L5
BLOCK,-W3,0,-H3-H2,-H3-H2-H5,L3+L4,L3+L4+L5                   !建立xia质量块模型
BLOCK,0,W2/2,-H3-H2,-H3-H2-H5,L3+L4,L3+L4+L5
BLOCK,W2/2,W2,-H3-H2,-H3-H2-H5,L3+L4,L3+L4+L5
BLOCK,W2,W2+W3,-H3-H2,-H3-H2-H5,L3+L4,L3+L4+L5
BLOCK,0,W2/2,-H3/2,0,L3+L4+L5+L5,L3+L4+L5+L5+L5/2      !外部固定
BLOCK,W2/2,W2,-H3/2,0,L3+L4+L5+L5,L3+L4+L5+L5+L5/2
BLOCK,0,W2/2,-H3/2,-H3,L3+L4+L5+L5,L3+L4+L5+L5+L5/2
BLOCK,W2/2,W2,-H3/2,-H3,L3+L4+L5+L5,L3+L4+L5+L5+L5/2
vglue,all

!载荷参数
V=100  
Uy=10.e-3

!   材料属性
E1=2.0e9  
NU12=0.29  
G12=0.775e9  
d31=2.2e-11  
d32=0.3e-11
d33=-3.0e-11
ept33=12

!PZT上材料特性参数
mp,dens,1,7750                                     !定义密度
emunit,epzro,8.85e-12                         !自由空间介电常数设置为8.85 e-12 F/m。
MP,PERX,1,729                                          !介电常数
MP,PERY,1,635
MP,PERZ,1,729   
tb,piez,1
tbdata,2,-5.2                  ! e31
tbdata,5,15.1                  ! e33
tbdata,8,-5.2                  ! e31
tbdata,10,12.7                 ! e15
tbdata,15,12.7                 ! e15
TB,ANEL,1                                                            !弹性劲度系数
tbdata,1,13.9E10,7.43E10,7.78E10  ! c11,c13,c12
tbdata,7,11.5E10,7.43E10       ! c33,c13
tbdata,12,13.9E10              ! c11
tbdata,16,2.56E10              ! c44
tbdata,19,2.56E10              ! c44
tbdata,21,3.06E10              ! c66
MP,DENS,1,7750                          !PZT参数
MP,NUXY,1,0.29
mp,dens,2,8300            !AL 层材料参数
mp,ex,2,110e9
mp,nuxy,2,0.324
MP,DENS,3,8900               !质量块参数
MP,EX,3,180E9
MP,NUXY,3,0.3

ET,1,SOLID5,3
ET,2,SOLID45
ET,3,combin39

MAT,1$TYPE,1 !PZT66上进行网格划分
VSEL,S,,,66
VSEL,A,,,55
ESIZE,5E-4
MSHKEY,1
MSHAPE,0,3D
numstr,node,14
VMESH,all
!*get,Epz,elem,,count           ! Get the number of solid elements得到固体元素的数量
MAT,1$TYPE,1 !PZT74上进行网格划分
VSEL,S,,,74
VSEL,A,,,56
ESIZE,5E-4
MSHKEY,1
MSHAPE,0,3D
VMESH,all

MAT,2$TYPE,2 !AL53层网格划分
VSEL,S,,,53
VSEL,A,,,52
VSEL,A,,,51
VSEL,A,,,1
VSEL,A,,,68
VSEL,A,,,65
VSEL,A,,,63
VSEL,A,,,50
VSEL,A,,,69
VSEL,A,,,67
VSEL,A,,,64
VSEL,A,,,49
VSEL,A,,,91
VSEL,A,,,90
VSEL,A,,,89
VSEL,A,,,73
VSEL,A,,,94
VSEL,A,,,72
VSEL,A,,,70
VSEL,A,,,71
VSEL,A,,,93
VSEL,A,,,92
VSEL,A,,,54
VSEL,A,,,75
ESIZE,5E-4
MSHKEY,1
MSHAPE,0,3D
VMESH,all
MAT,3$TYPE,2 !上、下质量块80网格划分
VSEL,S,,,80
VSEL,A,,,81
VSEL,A,,,78
VSEL,A,,,79
VSEL,A,,,76
VSEL,A,,,77
VSEL,A,,,58
VSEL,A,,,57
VSEL,A,,,86
VSEL,A,,,87
VSEL,A,,,84
VSEL,A,,,85
VSEL,A,,,82
VSEL,A,,,83
VSEL,A,,,60
VSEL,A,,,59
ESIZE,5E-4
MSHKEY,1
MSHAPE,0,3D
VMESH,all

MAT,2$TYPE,2 !AL62层网格划分
VSEL,S,,,62
VSEL,A,,,88
VSEL,A,,,61
VSEL,A,,,45
ESIZE,5E-4
MSHKEY,1
MSHAPE,0,3D
VMESH,all

asel,s,loc,z,,,
da,all,all,,
allsel,all
ACEL,0,9.8,0
NSEL,ALL
FINISH
/SOLU   
MODOPT,LANB,8,400,1200   !采用LANCZOS算法,设置计算8阶模态以及频率范围
nsel,s,loc,z,0,L3                                              !定义上压电片下层电极
nsel,r,loc,x,0,W2
nsel,r,loc,y,0
cp,1,volt,all
*get,n_bot,node,0,num,min
nsel,s,loc,z,0,L3                                       !定义上压电片上层电极
nsel,r,loc,x,0,W2
nsel,r,loc,y,H2+H2
cp,2,volt,all
*get,n_top,node,0,num,min@wuming524
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
你使用了 solid5  耦合单元,做的是模态分析,但是看帮助文件,solid5 单元好像没有模态分析的功能。
  而且,你定义了两个节点耦合:
    cp,1,volt,all
    cp,2,volt,all
    节点 1 和 2 互为从节点;其它节点在两个 cp 命令中都是从节点,这两条都是错误的,将不能求解。
  而且,你的命令流本身也没有进行求解。没有结果,如何显示电压。

» 本帖已获得的红花(最新10朵)

2楼2019-03-09 10:06:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

魔石乱枫

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by shgao20 at 2019-03-09 10:06:53
你使用了 solid5  耦合单元,做的是模态分析,但是看帮助文件,solid5 单元好像没有模态分析的功能。
  而且,你定义了两个节点耦合:
    cp,1,volt,all
    cp,2,volt,all
    节点 1 和 2 互为从节点;其它节 ...

哦哦,谢谢大神指点,我再修改下试试

发自小木虫Android客户端
3楼2019-03-10 16:28:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyaobo0509

专家顾问 (著名写手)

mark一下,改天运行看看;1 solid5可以模态分析;2 看你应该是想获得输出电压信号,谐响应分析可以用完全法试试(更简单),模态叠加法用于大模型计算可节省时间,对新手来说相对麻烦

发自小木虫Android客户端
4楼2019-04-18 09:38:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 魔石乱枫 的主题更新
信息提示
请填处理意见