supericon 发表于 2011-12-19 08:57:31

为什么solid117的速度效应的实常数加不上去

   一个线圈通过感应涡流加热钢板的模型,这个问题已经困惑我很久了。。也做了很多次的模拟,但是却一直弄不出速度效应来,命令流文件如下。神曲了划分网格的部分。该程序完全可以运算,但算完之后却没有速度效应出来,和静止的居然一摸一样。。
   
shpp,off
vovlap,all
VGLUE,ALL   
numcmp,all
EMUNIT,MKS                  !
*AFUN,DEG                     !
ET,1,SOLID117            !          空气
ET,2,SOLID117,1,1               !   钢板属性
r,2
rmore, , ,10
ET,3,SOLID117,1               !       线圈属性
mp,murx,1,1 ! air relative permeability
mp,murx,3,1 ! coil relative permeability
mp,rsvx,1,1e-9
mptemp,1,25.5,160,291.5,477.6,635,698 ! temps for relative permea   
mptemp,7,709,720.3,742,761,1000
mpdata,murx,2,1,200,190,182,161,135,104 ! steel relative permeabili
mpdata,murx,2,7,84,35,17,1,1
mptemp
mptemp,1,0,125,250,375,500,625 ! temps for resistivity
mptemp,7,750,875,1000   
mpdata,rsvx,2,1,.184e-6,.272e-6,.384e-6,.512e-6,.656e-6,.824e-6
mpdata,rsvx,2,7,1.032e-6,1.152e-6,1.2e-6 ! steel resistivity
mptemp
mptemp,1,0,730,930,1000 ! temps for conductivity
mpdata,kxx,2,1,60.64,29.5,28,28
mptemp ! temps for enthalpy
mptemp,1,0,27,127,327,527,727   
mptemp,7,765,765.001,927
mpdata,enth,2,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e9
mpdata,enth,2,7,3.548547e9,3.548556e9,4.3520e9
mp,emis,2,.68 ! emissivity               !
*SET,FTIME,0.2                  !!!
*SET,TINC,0.05                                 !
*SET,TIME,0                                  !
*SET,DELT,0.01
VSEL,S, , ,      4,5,1
vatt,2,2,2          !钢的属性
VSEL,S, , ,       1
vatt,1,1,3       !线圈属性
VSEL,S,, ,2,3,1
VSEL,A, , ,6,7,1
vatt,3,1,1       !空气属性

网格划分部分


allsel
eplot   
NSEL,EXT
D,ALL,AZ,0
FINISH
/SOL
ANTYP,HARM                              !
HARFR,50000
ASEL,S, , ,      5
NSLA,S,1
d,all,volt,0
ASEL,S, , ,   9
nsla,s,1
cp,1,volt,all   
*get,n1,node, ,num,max
F,n1,AMPS,14000
allsel
PHYSICS,WRITE,EMAG                        !   写入emag
FINISH
/PREP7
LSCLEAR,ALL                               !
ET,1,0                                    !   将1号材料转换成0
ET,2,SOLID90                              !   将2号材料转换为SOLID90   
ET,3,0                                    !   将3号单元转换成0
MP,KXX,2,30             !!!!!传导率各向同性Thermal conductivities (also KYY, KZZ
!MP,KYY,2,30
!MP,KZZ,2,30
MP,C,2,500                              !
MP,DENS,2,7900                            !
FINISH
/SOLU   
ANTYPE,TRANS                              !
TOFFST,273
TUNIF,100                                  !
CNVTOL,HEAT,1                           !
KBC,1                                     !
TRNOPT,FULL
AUTOS,ON                                  !   自动时间-   
DELTIM,5E-3,5E-4,DELT,ON                  !
OUTRES,BASIC,ALL                           !
PHYSICS,WRITE,THERMAL                      !   写入THERMAL
FINISH
*do,i,1,ftime/tinc ! solution *do loop
*SET,time,time+tinc ! increment time
physics,read,emag ! read emag physics file
/solu   
*if,i,eq,1,then
tunif,100 ! initial temperature
*else   
ldread,temp,last,,,,,rth ! read thermal analysis temperatures   
*endif
solve ! solve harmonic analysis
finish
physics,read,thermal ! read thermal physics file
/assign,esav,therm,esav ! redirect files for use in thermal res
/assign,emat,therm,emat
/solu   
parsav,scalar,parameter,sav !save parameters before multiframe rest
*if,i,gt,1,then
antype,trans,rest ! thermal restart
*endif
parres,new,parameter,sav !restore parameters after multiframe re
time,time ! time at end of thermal run
esel,s,mat,,2 ! select billet region
ldread,hgen,,,,2,,rst ! apply coupled joule heating load from   
esel,all
solve   
finish
/assign,esav ! reassign files to default
/assign,emat
*enddo ! end of solution looping

supericon 发表于 2011-12-19 08:59:28

红色线圈 紫色钢板。

supericon 发表于 2011-12-19 11:12:18

难道我是在玩单机版么。。

lwx27 发表于 2011-12-20 18:42:42

赞LZ,如何用SOLID117的速度效应也是我所关注的。
LZ试着把问题简化去验算SOLID117的速度效应了吗?LZ确定会使用SOLID117了吗?
这方面我是初学者,希望LZ能来继续探讨。

yezj361 发表于 2012-5-25 10:49:18

SOLID117单元 是基于棱边单元法的,棱边单元法是不能有运动效应和电路耦合的,所以结果里面何来运动效应?
页: [1]
查看完整版本: 为什么solid117的速度效应的实常数加不上去