24小时热门版块排行榜    

查看: 3859  |  回复: 0

绝无葬爱

新虫 (初入文坛)

[求助] ansys 电场仿真

*DIM,load,CHAR,24,1,1, , ,
*SET,load(1),'a11','a12','a13','a14','a21','a22','a23','a24','a31','a32','a33','a34','a41','a42','a43','a44','a51','a52','a53','a54','a61','a62','a63','a64'
FINISH  
/PREP7  
*set,tij,0.0225
*do,i,1,24!循环求解12次
FINISH  
/SOL
DA,pingbi,VOLT,0 !给屏蔽施加0V电压   
allsel,all  
*do,j,1,24  
DA,load(j),VOLT,0 !给12个电极施加0V电压
*enddo  
DA,load(i),VOLT,5 !给第i个电极施加5V电压   
solve!求解
finish  
/post1!进入后处理  
ALLSEL,ALL  
FLST,5,1,6,ORDE,1   
FITEM,5,31
VSEL,S, , ,P51X
ALLSEL,BELOW,VOLU   
EPLOT   
etable,%load(i)%x,ef,x !提取被测材料的电场强度X值单元表  
etable,%load(i)%y,ef,y !提取被测材料的电场强度Y值单元表  
etable,%load(i)%z,ef,z !提取被测材料的电场强度Z值单元表  
*VGET,%load(i)%x,ELEM, ,ETAB,%load(i)%x, , ,2  !将电场强度表储存为向量   
*VGET,%load(i)%y,ELEM, ,ETAB,%load(i)%y, , ,2   
*VGET,%load(i)%z,ELEM, ,ETAB,%load(i)%z, , ,2   
FINISH  
/SOL
allsel,all  
DADELE,all,VOLT!删除所施加的电压
*enddo  
*do,i,1,23!计算S矩阵的66次循环  
*do,j,i+1,24
*voper,%load(i)%%load(j)%x,%load(i)%x,mult,%load(j)%x!计算S值
*voper,%load(i)%%load(j)%y,%load(i)%y,mult,%load(j)%y!计算S值
*voper,%load(i)%%load(j)%z,%load(i)%z,mult,%load(j)%z!计算S值
*voper,%load(i)%%load(j)%xy,%load(i)%%load(j)%x,add,%load(i)%%load(j)%y!计算S值  
*voper,%load(i)%%load(j)%,%load(i)%%load(j)%xy,add,%load(i)%%load(j)%z!计算S值   
*vfact,-tij/25!计算S值   
*vfun,%load(i)%%load(j)%,copy,%load(i)%%load(j)%!计算S值
*dim,%load(i)%%load(j)%new,,8192,1  
*vfact,1
*vfun,%load(i)%%load(j)%new,copy,%load(i)%%load(j)%!将前5120行复制到新向量名字为load(i)load(j)new   
*enddo  
*enddo  
SAVE
/REPLOT,RESIZE  
write   
FINISH  
! /EXIT
这是24块电容极板循环接5v电压激励求电场线,然后计算灵敏度S的程序,前面程序我省略了,但是运行完求取的灵敏度不完整,下面这段程序是12极板的求取过程,这个没有错,请大家帮我找找错哪里了
*do,i,1,12!循环求解12次
FINISH  
/SOL
DA,pingbi,VOLT,0 !给屏蔽施加0V电压   
allsel,all  
*do,j,1,12  
DA,load(j),VOLT,0 !给12个电极施加0V电压
*enddo  
DA,load(i),VOLT,5 !给第i个电极施加5V电压   
solve!求解
finish  
/post1!进入后处理  
ALLSEL,ALL  
FLST,5,1,6,ORDE,1   
FITEM,5,14  
VSEL,S, , ,P51X
ALLSEL,BELOW,VOLU   
EPLOT   
etable,%load(i)%x,ef,x !提取被测材料的电场强度X值单元表  
etable,%load(i)%y,ef,y !提取被测材料的电场强度Y值单元表  
etable,%load(i)%z,ef,z !提取被测材料的电场强度Z值单元表  
*VGET,%load(i)%x,ELEM, ,ETAB,%load(i)%x, , ,2  !将电场强度表储存为向量   
*VGET,%load(i)%y,ELEM, ,ETAB,%load(i)%y, , ,2   
*VGET,%load(i)%z,ELEM, ,ETAB,%load(i)%z, , ,2   
FINISH  
/SOL
allsel,all  
DADELE,all,VOLT!删除所施加的电压
*enddo  
*do,i,1,12!计算S矩阵的66次循环  
*do,j,i+1,12
*voper,%load(i)%%load(j)%x,%load(i)%x,mult,%load(j)%x!计算S值
*voper,%load(i)%%load(j)%y,%load(i)%y,mult,%load(j)%y!计算S值
*voper,%load(i)%%load(j)%z,%load(i)%z,mult,%load(j)%z!计算S值
*voper,%load(i)%%load(j)%xy,%load(i)%%load(j)%x,add,%load(i)%%load(j)%y!计算S值  
*voper,%load(i)%%load(j)%,%load(i)%%load(j)%xy,add,%load(i)%%load(j)%z!计算S值   
*vfact,-tij/25!计算S值   
*vfun,%load(i)%%load(j)%,copy,%load(i)%%load(j)%!计算S值
*dim,%load(i)%%load(j)%new,,5120,1  
*vfact,1
*vfun,%load(i)%%load(j)%new,copy,%load(i)%%load(j)%!将前5120行复制到新向量名字为load(i)load(j)new   
*enddo  
*enddo  
SAVE
/REPLOT,RESIZE  
write   
FINISH  
! /EXIT,ALL
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 绝无葬爱 的主题更新
信息提示
请填处理意见