How can I simulate joule heating in a stranded conductor?
The coil has temperature-dependent resistivity.
The input below is the simplest way, but it does use an undocumented plane13 keyopt(1)=5. The results can be checked by hand (compare the etable item jheat with element volu x RSVZ(temp) x JSZ^2). /prep7 et,1,13,5,,1 ! axisymmetric az,temp mp,kxx,1,.005 mptemp,1,0,500 mpdata,RSVZ,1,1,1e-7,1e-6 mp,murx,1,1 ! to prevent a warning message rect,0.1,0.12,,0.02 esize,.002 amesh,1 d,all,az,0 ! constrain unused dof bf,all,js,,,1e5 ! known current density of a stranded coil nsel,s,ext sf,all,conv,.5,25 nsel,all fini /solu solve fini /post1 set,last plns,temp etab,jh,jheat pletab,jh ! higher at center where resistivity is greater |
||
|