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

[电磁与结构流体等的耦合] 有关电弧模拟的

[复制链接]
发表于 2011-6-16 17:35:02 | 显示全部楼层 |阅读模式 来自 清华大学
大家好,我初来乍到,我想有没有做过电弧模拟的,涉及到电磁场和流畅耦合,看过一些论文,,说是先要计算电场,再将电场的电流密度导入磁场计算焦耳热,将生热率导入流场进行顺序耦合,具体怎么实现?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
 楼主| 发表于 2011-6-16 17:35:44 | 显示全部楼层 来自 清华大学
Simdroid开发平台
我代码如下:
/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,
回复 不支持

使用道具 举报

发表于 2019-2-28 20:20:27 | 显示全部楼层 来自 上海嘉定区
谢谢分享,好资料啊
回复 不支持

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-3-28 18:37 , Processed in 0.033307 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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