| 查看: 3951 | 回复: 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 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有141人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
