唉,水平实在很菜,分配材料属性的时候又出问题了
我把程序简化了下,ICEMMAN兄弟,能再帮忙看看吗,真是不好意思,麻烦你这么多次,多谢了!!
!电场分析
/prep7
/title, model
emunit,mks !国际单位制
ET,1,SOLID5
KEYOPT,1,1,9
KEYOPT,1,3,0
KEYOPT,1,5,0
!solid5单元(六面体20节点线性元,电压自由度,传导电流分析)
i=1500
!建立模型(单位:m)
ma=0.580$mb=1.230$mc=0.580$md=0.930$me=0.580$mf=1.080 !镁液槽
block,ma,mb,mc,md,me,mf
yi1a=0.600$ yi1b=0.610$yic=0.730$yid=0.880$yie=0.790$yif=1.860 !阴极1-5
block, yi1a,yi1b,yic,yid,yie,yif
yi2a=0.750$yi2b=0.760$block,yi2a,yi2b,yic,yid,yie,yif
yi3a=0.900$yi3b=0.910$block,yi3a,yi3b,yic,yid,yie,yif
yi4a=1.050$yi4b=1.060$block,yi4a,yi4b,yic,yid,yie,yif
yi5a=1.200$yi5b=1.210$block,yi5a,yi5b,yic,yid,yie,yif
ya1a=0.660$ya1b=0.700$yac=0.730$yad=1.380$yae=0.820$yaf=0.990 !阳极1-4
block,ya1a,ya1b,yac,yad,yae,yaf
ya2a=0.810$ya2b=0.850$block,ya2a,ya2b,yac,yad,yae,yaf
ya3a=0.960$ya3b=1.000$block,ya3a,ya3b,yac,yad,yae,yaf
ya4a=1.110$ya4b=1.150$block,ya4a,ya4b,yac,yad,yae,yaf
ha=0.600$hb=1.210$hc=0.790$hd=0.820$he=0.780$hf=0.790 !焊条
block, ha,hb,hc,hd,he,hf
!overlap
vsel,all$vovlap,all
numcmp,volu
/pnum,volu,1
!材料参数(待修改)rsvx:电阻率; perx:介电常数;murx:磁导率
mp,rsvx,4,0.25$mp,perx,4,80$ mp,murx,4,1.2 !溶液(氯化镁)
mp,rsvx,5,1e-7$ mp,perx,5,15$mp,murx,5,2 !阴极(铁)
mp,rsvx,6,5e-7$mp,perx,6,15$ mp,murx,6,2 !阳极(铁)
mp,rsvx,7,5e-6$mp,perx,7,15$ mp,murx,7,2 !焊条(铁)
vsel,s,loc,x,ma,mb !镁槽
vsel,r,loc,y,mc,md
vsel,r,loc,z,me,mf
cm,mc,volu
vatt,4,1,1
VSEL,S,MAT,,4
/replot
vsel,s,loc,x,ha,hb !焊条
vsel,r,loc,y,hc,hd
vsel,r,loc,z,he,hf
cm,ht,volu
vatt,7,1,1$/replot
VSEL,S,MAT,,7
/replot
Asel,s,loc,y,yic,yid !阴极1~5
Asel,r,loc,z,yie,yif
Vsla,s,1
Cm,yi5,volu
vatt,5,1,1
/replot
VSEL,S,MAT,,5
/replot
Asel,s,loc,y,yac,yad !阳极1~4
Asel,r,loc,z,yae,yaf
Vsla,s,1
Cm,ya4,volu
vatt,6,1,1
/replot
VSEL,S,MAT,,6
/replot
Cmsel,s,mc !镁槽
Cmsel,u,yi5
Cmsel,u,ya4
Cmsel,u,ht
cm,mc,volu
vatt,4,1,1
VSEL,S,MAT,,4
/replot
VSEL,S,MAT,,4
/replot
VSEL,S,MAT,,7
/replot
VSEL,S,MAT,,5
/replot
VSEL,S,MAT,,6
/replot |