有关电弧模拟的
大家好，我初来乍到，我想有没有做过电弧模拟的，涉及到电磁场和流畅耦合，看过一些论文，，说是先要计算电场，再将电场的电流密度导入磁场计算焦耳热，将生热率导入流场进行顺序耦合，具体怎么实现？ 我代码如下：/units,si
pi=4*atan(1)
i=300
r=1.6e-3 !钨极半径
/PREP7
k,1,0,8.3e-3,,,,,
k,2,8.3e-3,8.3e-3,,,,,
k,3,8.3e-3,0,,,,,
k,4,0,0,,,,,
k,5,0,5.5e-3,,,,,
k,6,0.84e-3,5.5e-3,,,,,
k,7,1.6e-3,8.3e-3,,,,,
k,8,0.84e-3,0,,,,,
ET,1,PLANE53
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,5,0
KEYOPT,1,7,0
ET,2,0
MPTEMP,,,,,,,,!密度kg/m3
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,DENS,1,,1.67
MPDATA,DENS,1,,0.2
MPDATA,DENS,1,,0.1
MPDATA,DENS,1,,0.07
MPDATA,DENS,1,,0.04
MPDATA,DENS,1,,0.03
MPDATA,DENS,1,,0.02
MPDATA,DENS,1,,0.015
MPDATA,DENS,1,,0.014
MPTEMP,,,,,,,,!热导率J/m*s*K
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,KXX,1,,0
MPDATA,KXX,1,,0.1
MPDATA,KXX,1,,0.1
MPDATA,KXX,1,,0.3
MPDATA,KXX,1,,0.6
MPDATA,KXX,1,,1.6
MPDATA,KXX,1,,2.4
MPDATA,KXX,1,,2
MPDATA,KXX,1,,2.2
MPTEMP,,,,,,,,!电阻率 欧姆.m
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,RSVX,1,,1000
MPDATA,RSVX,1,,100
MPDATA,RSVX,1,,1
MPDATA,RSVX,1,,1.3e-3
MPDATA,RSVX,1,,3.3e-4
MPDATA,RSVX,1,,1.6e-4
MPDATA,RSVX,1,,1.2e-4
MPDATA,RSVX,1,,9.8e-5
MPDATA,RSVX,1,,8.5e-5
MPTEMP,,,,,,,,!相对磁导率
MPTEMP,1,0
MPDATA,MURX,1,,0.9
MPTEMP,,,,,,,,!比热J/kg.K
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,c,1,,521
MPDATA,c,1,,521
MPDATA,c,1,,521
MPDATA,c,1,,600
MPDATA,c,1,,1250
MPDATA,c,1,,6000
MPDATA,c,1,,10000
MPDATA,c,1,,4000
MPDATA,c,1,,2500
MPTEMP,,,,,,,,!粘度kg/m.s
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,visc,1,,2e-5
MPDATA,visc,1,,7e-5
MPDATA,visc,1,,1e-4
MPDATA,visc,1,,2.2e-4
MPDATA,visc,1,,2.7e-4
MPDATA,visc,1,,2e-4
MPDATA,visc,1,,6e-5
MPDATA,visc,1,,3e-5
MPDATA,visc,1,,2.6e-5
!钨极材料性能
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,19300
MPDATA,RSVX,2,,0.0532
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,180
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,130
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,2,,1.26e-6
!建模及网格划分
LSTR, 1, 5
LSTR, 5, 4
LSTR, 3, 2
LSTR, 2, 7
LSTR, 7, 1
LSTR, 5, 6
LSTR, 6, 7
LSTR,6,8
LSTR, 4, 8
LSTR, 8, 3
al,1,6,7,5
al,6,2,9,8
al,7,8,10,3,4
aglue,all
numcmp,area
lesize,6,,,20
lesize,1,,,10
lesize,2,,,60
lccat,4,3
mshape,0,2d
mshk,1
aatt,2,,1,,,
amesh,1
aatt,1,,1,,,
asel,s,area,,2,3,1
amesh,all
allsel,all
!加载及边界条件
ri=i/(pi*r**2)
BFA,1,JS, , ,ri,0
lsel,s,,,10
lsel,a,,,3
lsel,a,,,4
lsel,a,,,5
lsel,a,,,9
sfl,all,inf
lsel,s,,,5,7,1
lsel,a,,,1
lsel,a,,,9,10
BFL,all,TEMP,10000
lsel,s,,,3,4
bfl,all,temp,6000
/solu
ANTYPE,0
/SOL
MAGSOLV,0,3,0.001, ,25, 谢谢分享，好资料啊
页:
[1]