找回密码
 注册
Simdroid-非首页
查看: 130|回复: 9

[压电] 请教:压电仿真,网格划分报错

[复制链接]
发表于 2010-3-20 19:16:12 | 显示全部楼层 |阅读模式 来自 北京海淀
各位老师,以下是我的命令流,运行报错,烦请各位老师帮帮我,不甚感激!
/prep7
! PZT-5H 材料特性参数
mp,DENS,1,7500                 ! 密度, kg/m**3
mp,perx,1,1100                 ! 介电常数
mp,pery,1,1100
mp,perz,1,827
tb,ANEL,1                      ! 弹性劲度系数, N/m^2
tbdata,1,13.2E10,7.1E10,7.31E10  ! c11,c12,c13
tbdata,7,13.2E10,7.3E10       ! c11,c33
tbdata,12,11.5E10              ! c33
tbdata,16,3.0E10              ! c44
tbdata,19,2.6E10              ! c44
tbdata,21,2.6E10              ! c66
tb,PIEZ,1                      ! 压电(应力)常数, C/m^2
tbdata,3,-5.2                  ! e13
tbdata,6,-5.2                 ! e23
tbdata,9,15.1                 ! e33
tbdata,14,12.7                ! e52
tbdata,16,12.7                ! e61
!定义主结构的材料参数
mp,dens,2,2300       !定义SI层材料参数  
MP,EX,2,130e9
MP,nuxy,2,0.27
mp,dens,3,8900       !定义Cu层材料参数  
MP,EX,3,110e9
MP,nuxy,3,0.324                  

mp,dens,4,8900       !定义NI层材料参数  
MP,EX,4,180e9
MP,nuxy,4,0.3

! 定义关键点
K,1,0,0,0             !SIO2层
k,2,4e-3,3e-4,0
K,3,4e-3,9e-4,0
k,4,0,12e-4,0
K,5,0,0,3e-6
k,6,4e-3,3e-4,3e-6
K,7,4e-3,9e-4,3e-6
k,8,0,12e-4,3e-6
K,9,0,0,13e-6             !Cu层
k,10,4e-3,3e-4,13e-6
K,11,4e-3,9e-4,13e-6
k,12,0,12e-4,13e-6
K,13,0,0,28e-6              !PZT层
k,14,4e-3,3e-4,28e-6
K,15,4e-3,9e-4,28e-6
k,16,0,12e-4,28e-6
K,17,3.5e-3,3e-4,28e-6              !NI层
k,18,4e-3,3e-4,28e-6
K,19,4e-3,9e-4,28e-6
k,20,3.5e-3,9e-4,28e-6
  
K,21,3.5e-3,3e-4,628e-6              
k,22,4e-3,3e-4,628e-6
K,23,4e-3,9e-4,628e-6
k,24,3.5e-3,9e-4,628e-6
! 定义体
V,1,2,3,4,5,6,7,8
V,5,6,7,8,9,10,11,12
V,9,10,11,12,13,14,15,16
V,17,18,19,20,21,22,23,24
vglue,all                             !将梁同压电片粘结
!定义有限元结构
et,1,solid5,3              !定义压电单元
et,2,solid45                  !定义主结构单元
!采用映射划分网格连接相邻面
asel,s,loc,z,28e-6
cm,CM_2,area
accat,CM_2

!进行网格划分
LESIZE,ALL,10e-5, , , ,1, , ,1,
mat,1 $ type,1   !压电片网格划分
Vmesh,6
mat,2 $ type,2   !SI层网格划分
Vmesh,1
mat,3 $ type,2   !Cu网格划分
vmesh,2
mat,4 $ type,2   !Ni网格划分
vmesh,5

nsel,s,loc,z,3e-6           !定义下层电极
cp,1,volt,all
*get,n_bot,node,0,num,min
d,n_bot,volt,0
nsel,s,loc,z,28e-6         !定义上层电极
cp,2,volt,all
*get,n_top,node,0,num,min
nsel,s,loc,x,0
d,all,ux,0,,,,uy,uz
   !下层电极短路
nsel,all
fini
/solu
antype,modal
modopt,lanb,3
mxpand,3
solve
fini
 楼主| 发表于 2010-3-20 20:07:50 | 显示全部楼层 来自 北京海淀
Simdroid开发平台
呼唤hy1966老师
回复 不支持

使用道具 举报

发表于 2010-3-21 13:32:01 | 显示全部楼层 来自 北京
把块分小,让每个块都是六面体,每个六面体的面都是单一完整的面(多个面时也可以合并)就可以了
回复 不支持

使用道具 举报

 楼主| 发表于 2010-3-21 23:06:38 | 显示全部楼层 来自 北京海淀
谢谢hy1966,不好意思,看了您的回答,我还是不知道怎么弄,请老师再具体指点下,再次感谢!
回复 不支持

使用道具 举报

 楼主| 发表于 2010-3-21 23:07:50 | 显示全部楼层 来自 北京海淀
把块变小是什么意思呢
回复 不支持

使用道具 举报

发表于 2010-3-22 12:54:18 | 显示全部楼层 来自 浙江杭州
你找找看,类似的问题好像讨论过,就在这个版块
回复 不支持

使用道具 举报

 楼主| 发表于 2010-3-22 16:57:31 | 显示全部楼层 来自 北京海淀
唉,版主,我找了个遍也没找到,还烦请hy1966老师指导,多谢多谢!
回复 不支持

使用道具 举报

发表于 2010-3-22 17:55:08 | 显示全部楼层 来自 北京
搜索了一下,上次也是你问的,呵呵
http://forum.simwe.com/viewthrea ... hlight=%D0%FC%B1%DB
回复 不支持

使用道具 举报

 楼主| 发表于 2010-3-22 21:53:03 | 显示全部楼层 来自 北京海淀
呵呵,上次估计也没搞懂,只是按hy1966老师的修改方式改正而已。我这里的命令流与上次的命令流只不过是模型从矩形梁改为梯形梁,其他都差不多,就是怎么样网格都变不过来,唉,还是烦请hy1966老师帮我彻底把这个问题解决,以后碰到这种问题就能自己解决了,谢谢老师,谢谢老师,在这里我也只能说谢谢两个字。
回复 不支持

使用道具 举报

 楼主| 发表于 2010-3-23 16:14:41 | 显示全部楼层 来自 北京海淀
厚着脸皮呼唤hy1966版主,谢谢。
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-9-30 11:25 , Processed in 0.034297 second(s), 14 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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