ansys 谁能给个三维磁场用到solid97的实例

例子内容：a long, thick solenoid carries a uniform current density d**tribution, j. assuming that the turns of the solenoid can be modeled as a homogeneous **otropic material with modulus of elasticity e, and po**son's ratio ν, determine the axial magnetic flux density d**tribution bθ and the circumferential stress σo d**tribution in the solenoid.程序：/prep7smrt,off/title, vm172, stress analys** of a long, thick, **otropic solenoid/com, magneto-solid mechanics, moon, pg. 275, 2d analys**antype,static ! coupled field analys**et,1,plane13,,,1 ! plane13, az dof, ax**ymmetric optionet,2,plane13,4,,1 ! plane13, az,ux,uy dof, ax**ymmetric optionmp,ex,2,10.76e10 ! solenoid modulus of elasticitymp,nuxy,2,.35 ! solenoid po**son ratioemunit,mks ! mks unitsmp,murx,1,1 ! relative permeability=1.0 mp,murx,2,1 k,1k,2,1e-2k,3,2e-2l,1,2lesize,1,,,5l,2,3lesize,2,,,20lgen,2,all,,,,2e-3a,1,2,5,4 ! area 1=aira,2,3,6,5 ! area 2=solenoidasel,s,area,,2aatt,2,,2asel,allesize,,1mshk,2 ! mapped area mesh if possiblemsha,0,2d ! using quads amesh,allesel,s,mat,,2nsle,snsel,r,loc,y,2e-3cp,1,uy,all ! couple solenoid nodal uy d**p.nsel,s,loc,x,2e-2 ! select nodes at outer radiuscp,2,az,all ! couple az to ensure flux-parallel b.c.fin**h/solunsel,s,loc,x,0d,all,az,0 ! set flux parallel b.c.esel,s,mat,,2bfe,all,js,,,,1e+6 ! apply current density loadnslensel,r,loc,y,0dsym,symm,2 ! apply structural symmetry b.c.nsel,allesel,allkbc,1 ! step boundary conditionsoutres,,lastsolvefin**h/post1set,1esel,s,mat,,2 ! select solenoid nodes and elements nsle/axlab,x,d**tance/axlab,y,stress - 2-d model/gthk,ax**,2!/yrange,0,150path,coil1,2,,48 ! define path with name = "coil1"ppath,1,2 ! define path points by nodeppath,2,1**def,sz,s,z pdef,by,b,y!/yrange,500,1500plpath,sz ! d**play circum stress thru solenoid!/yrange,0,125/axlab,y,flux density - 2-d modelplpath,by ! d**play axial flux density thru solenoidnsel,s,loc,x,1e-2nsel,a,loc,x,1.3e-2nsel,a,loc,x,1.7e-2prnsol,b,comp ! print axial flux densityprnsol,s,comp ! print component stresses*get,b1,node,7,b,sum*get,b2,node,19,b,sum*get,b3,node,27,b,sum*get,s1,node,7,s,z*get,s2,node,19,s,z*get,s3,node,27,s,z*dim,label,char,6,2*dim,value,,6,3label(1,1) = 'b, t ','b, t ','b, t ','prs ','prs ','prs 'label(1,2) = 'r=1e-2','r=1.3e-2','r=1.7e-2','r=1e-2','r=1.3e-2','r=1.7e-2'*vfill,value(1,1),data,0.01257,8.796e-3,3.77e-3,146.7,97.79,62.44*vfill,value(1,2),data,b1,b2,b3,s1,s2,s3v1=b1/0.01257v2=b2/8.796e-3v3=b3/3.77e-3v4=s1/146.7v5=s2/97.79v6=s3/62.44*vfill,value(1,3),data,v1,v2,v3,v4,v5,v6s**e,table_1fin**h/clear,nostart/prep7smrt,off/title, vm172, stress analys** of a long, thick, **otropic solenoid/com, magneto-solid mechanics, moon, pg. 275, 3d analys**antype,static ! coupled field analys**et,1,62 ! magneto-structural elementet,2,97 ! magnetic field elementet,3,47 ! infinite elementmp,ex,2,10.76e10 ! solenoid modulus of elasticitymp,nuxy,2,.35 ! solenoid po**son ratioemunit,mks ! mks unitsmp,murx,1,1 ! relative permeability=1.0 mp,murx,2,1 k,1 ! create 2-d meshk,2,1e-2k,3,2e-2l,1,2lesize,1,,,5l,2,3lesize,2,,,20lgen,2,all,,,,,2e-3a,1,2,5,4 ! area 1=aira,2,3,6,5 ! area 2=solenoidesize,,1mshk,2 ! mapped area mesh if possiblemsha,0,2d ! using quadstype,3 ! mesh with infin47amesh,alltype,2 ! assign 3-d element typemat,1esize,,2 ! two div**ions for solid elementsvrotat,1,,,,,,1,4,10,1 ! rotate 10 degreestype,1mat,2vrotat,2,,,,,,1,4,10,1nummrg,node ! merge coincident nodesaclear,1,2 ! clear infin47 element meshcsys,1nrotat,all ! rotate nodes to cylindrical coordinatesesel,s,mat,,2nslensel,r,loc,z,2e-3cp,1,uz,all ! couple solenoid nodal uz d**p.esel,allnsel,s,loc,x,2e-2 ! select nodes at outer radiuscp,2,ay,all ! couple ay to ensure flux-parallel cond.fin**h/solunsel,s,loc,x,0d,all,ax,0,,,,ay,az ! flux-parallel along solenoid ax**esel,s,mat,,2bfe,all,js,,,1e+6 ! apply current density loadnsle,snsel,r,loc,z,0d,all,uz,0 ! apply structural symmetry b.c. to solenoidnsle,snsel,r,loc,y,0d,all,uy,0nsle,snsel,r,loc,y,10d,all,uy,0nsel,allesel,alld,all,ax,0,,,,az ! only ay required for ax**ym. fieldkbc,1 ! step boundary conditionsoutres,,lastcnvtol,f,1e-3 ! define force convergence valuesolvefin**h/post1set,lastesel,s,mat,,2 ! select solenoid nodes and elements nsle/axlab,x,d**tance/axlab,y,stress - 3-d model/gthk,ax**,2path,coil2,2,,48 ! define path with name = "coil2"ppath,1,2ppath,2,13rsys,1pdef,sy,s,y pdef,bz,b,z!/yrange,500,1500plpath,sy ! d**play circum stress thru solenoid!/yrange,0,125/axlab,y,flux density - 3-d modelplpath,bz ! d**play axial flux density thru solenoidnsel,s,loc,x,1e-2nsel,a,loc,x,1.3e-2nsel,a,loc,x,1.7e-2nsel,r,loc,y,0prnsol,b,comp ! print axial flux densityprnsol,s,comp ! print component stresses*get,b1,node,7,b,sum*get,b2,node,19,b,sum*get,b3,node,27,b,sum*get,s1,node,7,s,y*get,s2,node,19,s,y*get,s3,node,27,s,y*dim,label,char,6,2*dim,value,,6,3label(1,1) = 'b, t ','b, t ','b, t ','prs ','prs ','prs 'label(1,2) = 'r=1e-2','r=1.3e-2','r=1.7e-2','r=1e-2','r=1.3e-2','r=1.7e-2'*vfill,value(1,1),data,0.01257,8.796e-3,3.77e-3,146.7,97.79,62.44*vfill,value(1,2),data,b1,b2,b3,s1,s2,s3v1=b1/0.01257v2=b2/8.796e-3v3=b3/3.77e-3v4=s1/146.7v5=s2/97.79v6=s3/62.44*vfill,value(1,3),data,v1,v2,v3,v4,v5,v6s**e,table_2resume,table_1/com/out,vm172,vrt/com,------------------- vm172 results compar**on --------------/com,/com, | target | ansys | ratio/com,/com,pressures h**e units of n/m**2/com,/com,results using plane13:*vwrite,label(1,1),label(1,2),value(1,1),value(1,2),value(1,3)(1x,a8,a8,' ',e12.5,' ',e12.5,' ',1f5.3)/noprresume,table_2/gopr/com,/com,results using solid62/97:*vwrite,label(1,1),label(1,2),value(1,1),value(1,2),value(1,3)(1x,a8,a8,' ',e12.5,' ',e12.5,' ',1f5.3)/com,-----------------------------------------------------------/com,/outfin**h*l**t,vm172,vrt/delete,table_1/delete,table_2 20210311