ANSYS例子 (2)

1 来自 simwe 多物理场 热结构 参考一下这些例子 1 2 Sample Induction Heating Analysis of a Circular Billet help 3 5 1 3 induction heating of a solid cylinder billet 间接顺序法 5 8 5 4 University of Alberta ANSYS Tutorials 包括一些 ansys 分析例题 9 5 两个稳态热分析实例 9 6 接触分析 钢板卷制成圆筒实例 13 1 1 来自来自 simwesimwe 多物理场 热结构 参考一下这些例子 多物理场 热结构 参考一下这些例子 原创例子 热形变分析 关于热结构耦合的计算 在 WB12 中流 热 固分析 原创 2 2 SampleSample Induction HeatingInduction Heating AnalysisAnalysis ofof a a CircularCircular Billet helpBillet help 3 5 3 5 物理环境法 batch list title Induction heating of a solid cylinder billet prep7 shpp off com com row 015 outer radius of workpiece ric 0175 inner radius of coil roc 0200 outer radius of coil ro 05 outer radius of model t 001 model thickness freq 150000 frequency Hz pi 4 atan 1 pi cond 392e7 maximum conductivity muzero 4e 7 pi free space permeability mur 200 maximum relative permeability skind sqrt 1 pi freq cond muzero mur skin depth ftime 3 final time tinc 05 time increment for harmonic analysis time 0 initialize time delt 01 maximum delta time step Electromagnetic model et 1 53 1 PLANE53 axisymmetric AZ dof et 2 53 1 emunit mks set magnetic units mp murx 1 1 air relative permeability mp murx 3 1 coil relative permeability mptemp 1 25 5 160 291 5 477 6 635 698 temps for relative permeability mptemp 7 709 720 3 742 761 1000 mpdata murx 2 1 200 190 182 161 135 104 steel relative permeability 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 rectng 0 row 0 t billet rectng row ric 0 t air gap rectng ric roc 0 t coil rectng roc ro 0 t outer air aglue all numcmp area ksel s loc x row select keypoints at outer radius of workpiece kesize all skind 2 set meshing size to 1 2 skin depth ksel s loc x 0 select keypoints at center kesize all 40 skind set meshing size lsel s loc y t 2 select vertical lines lesize all 1 set 1 division through thickness lsel all asel s area 1 aatt 2 1 1 set attributes for billet region asel s area 3 aatt 3 1 2 set attributes for coil region asel s area 2 4 2 aatt 1 1 2 set attributes for air region asel all mshape 0 2d mshk 1 amesh 1 mesh billet area lsel s loc y 0 lsel a loc y t lsel u loc x row 2 lesize all 001 lsel all amesh all mesh remaining areas nsel s loc x d all az 0 apply flux normal b c nsel all esel s mat 3 bfe all js 15e6 apply current density to coil alls Thermal model et 4 55 1 PLANE55 thermal element axisymmetric agen 2 1 1 aatt 2 1 4 ksel s loc x row select keypoints at outer radius of workpiece kesize all skind 2 set meshing size to 1 2 skin depth ksel s loc x 0 select keypoints at center kesize all 40 skind set meshing size lsel s loc y t 2 select vertical lines lesize all 1 set 1 division through thickness lsel all mshape 0 2d mshk 1 amesh 5 asel 5 allsel below area nsel r loc x row sf all rdsf 0 68 1 Radiation alls spctemp 1 25 Space temperature v2dopt 1 radopt 0 01 stef 5 67e 8 Stefan Boltzman constant esel s mat 2 select billet material bfe all fvin 1 define volumetric interface finish solu mfan on Activate ANSYS Multi field solver analysis mfel 1 1 2 Field 1 ET s Emag mfel 2 4 Field 2 ET s Thermal mfor 1 2 Field solution order mfti ftime Final time mfdt tinc Stagger time increment mfco all 1e 3 Convergence criteria antyp harm Emag analysis options harfrq 150000 outres all all tunif 100 mfcm 1 Write Emag analysis options mfclear solu Clear analysis options antype trans Thermal analysis options toffst 273 tunif 100 initial uniform temperature kbc 1 step loads trnopt full autos on auto time stepping deltim 01 005 01 on time step control mfcm 2 Write Thermal analysis options mfvo 1 1 hgen 2 Transfer hgen from Emag to Thermal mfvo 1 2 temp 1 Transfer Temp from Thermal to Emag solve save finish post26 file field2 rth esel s type 4 nsle s nsel r loc x row get node at outer radius nsel r loc y 0 get nor node num min nsle s nsel r loc x 0 nsel r loc y 0 get node at centerline get nir node num min nsol 2 nor temp outerR Outer radius nsol 3 nir temp inner Inner radius centerline axlab y Temperature plvar 2 3 plot temperature prvar 2 3 finish post1 file field2 rth set last Solution at 3 seconds esel s type 4 select thermal elements plns temp plot temperature finish 3 3 inductioninduction heatingheating ofof a a solidsolid cylindercylinder billet billet 间接顺序法间接顺序法 5 8 5 8 batch list filenam induc config nres 100000 将载荷步最大设置为 100000 prep7 shpp off title induction heating of a solid cylinder billet com et 1 13 1 PLANE13 axisymmetric AZ dof et 2 13 1 et 3 151 1 1 1 SURF151 thermal radiation r 3 0 Real constant set for SURF151 row 015 outer radius of workpiece ric 0175 inner radius of coil roc 0200 outer radius of coil ro 05 outer radius of model t 001 model thickness freq 150000 frequency Hz pi 4 atan 1 pi cond 392e7 maximum conductivity muzero 4e 7 pi free space permeability mur 200 maximum relative permeability skind sqrt 1 pi freq cond muzero mur skin depth ftime 3 final time tinc 05 time increment for harmonic analysis time 0 initialize time delt 01 maximum delta time step emunit mks set magnetic units mp murx 1 1 air relative permeability mp murx 3 1 coil relative permeability mptemp 1 25 5 160 291 5 477 6 635 698 temps for relative permeability mptemp 7 709 720 3 742 761 1000 mpdata murx 2 1 200 190 182 161 135 104 steel relative permeability 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 rectng 0 row 0 t billet rectng row ric 0 t air gap rectng ric roc 0 t coil rectng roc ro 0 t outer air aglue all numcmp area 对面重新编号 ksel s loc x row select keypoints at outer radius of workpiece kesize all skind 2 set meshing size to 1 2 skin depth 为什么是 半个集肤深度 ksel s loc x 0 select keypoints at center kesize all 40 skind set meshing size lsel s loc y t 2 select vertical lines lesize all 1 set 1 division through thickness lsel all asel s area 1 aatt 2 1 1 set attributes for billet region asel s area 3 aatt 3 1 2 set attributes for coil region asel s area 2 4 2 aatt 1 1 2 set attributes for air region asel all mshape 0 2d mshk 1 amesh 1 mesh billet area lsel s loc y 0 lsel a loc y t also select lsel u loc x row 2 lesize all 001 lsel all amesh all mesh remaining areas n create space node for SURF151 get nmax node num max lsel s loc x row type 3 real 3 REAL NSET Sets the element real constant set attribute pointer mat 2 lmesh all mesh billet outer radius with SURF151 get emax elem num max emodif emax 3 nmax modify element to add space node for radiation et 3 0 reset type 3 to null element nsel s loc x d all az 0 apply flux normal b c az 向量磁位能 nsel all esel s mat 3 bfe all js 15e6 apply current density to coil 施加体力载荷 esel all finish solu antyp harm harfrq 150000 physics write emag write emag physics file finish prep7 lsclear all clear all b c s and options et 1 55 1 PLANE55 thermal element axisymmetric et 2 0 null element type for coil and air region et 3 151 1 1 1 SURF151 element for radiation keyopt 3 9 1 r 3 1 5 67e 8 form factor Stefan Boltzmann constant 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 finish solu antype trans toffst 273 tunif 100 initial uniform temperature d nmax temp 25 ambient temperature cnvtol heat 1 convergence tolerance kbc 1 step loads trnopt full autos on auto time stepping deltim 1e 5 1e 6 delt on time step control outres basic all save all load step information physics write thermal write thermal physics file finish do i 1 ftime tinc solution do loop 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 restart assign emat therm emat solu if i gt 1 then antype trans rest thermal restart endif time time time at end of thermal run esel s mat 2 select billet region ldread hgen 2 rmg apply coupled joule heating load from emag esel all solve finish assign esav reassign files to default assign emat enddo end of solution looping finish save save database post26 time history postprocessor show nsol 2 1 temp tempcl store temperature at billet centerline nsol 3 2 temp tempsurf store temperature at billet outer diameter plvar 2 3 plot temperature rise over time prvar 2 3 print temperature rise over time finish n 定义节点 lsel s loc x row 选择线 type 3 指定元素类型 real 3 指定单元实常数 mat 2 材料属性 lmesh all 对线进行划分网格 emodif emax 3 nmax 改变选中的单元类型为映射 et 3 0 重新设置类型 3 为空单元 4 4 UniversityUniversity ofof AlbertaAlberta ANSYSANSYS TutorialsTutorials 包括一些包括一些 ansysansys 分析例题分析例题 ANSYS is a general purpose finite element modeling package for numerically solving a wide variety of mechanical problems These problems include static dynamic structural analysis both linear and non linear heat transfer and fluid problems as well as acoustic and electromagnetic problems Most of these tutorials have been created using ANSYS 7 0 therefore make note of small changes in the menu structure if you are using an older or newer version google http www mece ualberta ca tutorials ansys index html 5 5 两个稳态热分析实例两个稳态热分析实例 一个轴对称模型的热应力分析一个轴对称模型的热应力分析 finish clear filname ex3 1 prep7 title thermal stress in long cylinder PNUM NODE 1 PNUM ELEM 1 et 1 plane55 1 定义 1 号单元 et 2 plane55 1 定义 2 号单元 mp kxx 1 3 对流系数 n 1 1875 n 8 625 fill 2 ngen 2 10 1 8 1 1 e 11 1 2 12 由节点形成 1 类型单元 type 2 e 12 2 3 13 由节点形成 2 类型单元 egen 5 1 2 由 2 号单元复制 5 个 type 1 e 8 18 17 7 节点形成 1 类型单元 finish solu antype static 指定为稳态热分析 形成 sda 文件 D 1 TEMP 10 11 10 定义内表面 1 和 11 节点温度为 10 度 D 8 TEMP 100 18 10 outpR BASIC ALL 输出所有子步的基本解 SOLVE 可求出结构温度分布 形成 rth 文件 finish post1 LFT NODE NODE 0 1875 0 0 IN NODE NODE 0 2788 0 0 RT NODE NODE 0 625 0 0 GET LFT TEMP NODE LFT NODE TEMP GET IN TEMP NODE IN NODE TEMP GET RT TEMP NODE RT NODE TEMP FINISH PREP7 ETCHG TTS KEYOPT 1 3 1 KEYOPT 2 3 1 KEYOPT 1 6 1 element formulation mixed U P MP EX 1 30E6 MP ALPX 1 435E 5 热膨胀系数 MP NUXY 1 3 CPNGEN 7 UY 11 18 Defines modifies or adds to a set of coupled degrees of freedom CPNGEN NSET Lab NODE1 NODE2 NINC 其中 lab 为耦合自由度的方向 NSET 为耦合集的编号 NINC 默认为 1 CP 8 UX 1 11 Defines or modifies a set of coupled degrees of freedom CP NSET Lab NODE1 NODE17 可选 17 个节点 NSET 和 Lab 含义同上 CPSGEN 8 1 8 Generates sets of coupled nodes from existing sets S 代表 sets CPSGEN ITIME INC NSET1 NSET2 NINC 耦合约束问题 FINISH SOLU ANTYPE STATIC D 1 UY 8 施加 1 到 8 节点在 UY 的位移为 0 LDREAD TEMP rth Reads results from the results file and applies them as loads LDREAD Lab LSTEP SBSTEP TIME KIMG Fname Ext 其中 Ext 为文件扩 展名 solve FINISH POST1 LFT NODE NODE 0 1875 0 0 RT NODE NODE 0 625 0 0 GET LFT AXST NODE LFT NODE S Y GET LFT TST NODE LFT NODE S Z GET RT AXST NODE RT NODE S Y GET RT TST NODE RT NODE S Z 一圆筒形的罐有一接管一圆筒形的罐有一接管 罐外径为罐外径为 3 3 英尺 壁厚为英尺 壁厚为 0 20 2 英尺 英尺 接管外径为接管外径为 0 50 5 英尺 壁厚为英尺 壁厚为 0 10 1 英尺 英尺 罐与接管的轴线垂直且接管远离罐的端部罐与接管的轴线垂直且接管远离罐的端部 70 200 300 400 500 70 200 300 400 500 分别对应不同的导热系数分别对应不同的导热系数 比热比热 对流系数对流系数 罐内流体温度为华氏罐内流体温度为华氏 450450 度度 与罐壁的对流换热系数与罐壁的对流换热系数 prep7 title Steady state thermal analysis of pipe junction units bin 使用英制单位 et 1 90 定义热单元 mp dens 1 285 密度 mptemp 70 200 300 400 500 建立温度表 mpdata kxx 1 8 35 12 8 90 12 9 35 12 9 80 12 10 23 12 导热系数 mpdata c 1 0 133 0 177 0 119 0 122 0 125 比热 specyfic heat mpdata hf 2 426 144 405 144 352 144 275 144 221 144 接管对流系数 heat flow 上面只用一条命令 MPTEMP 定义了温度 随后三条命令就可以定义与温度对应的 参数了 MPDATA Lab MAT STLOC C1 C2 C3 C4 C5 C6 其中 MAT 为材料编号 STLOC 默认 定义几何模型参数 ri1 1 3 罐内半径 ro1 1 5 罐外半径 z1 2 罐长 ri2 0 4 接管内半径 ro2 0 5 接管外半径 z2 2 接管长 建立几何模型 cylind ri1 ro1 z1 90 1 4 罐体 wprota 0 90 将工作平面旋转到垂直于接管轴线 cylind ri2 ro2 z2 90 1 4 接管 wpstyl defa 将工作平面恢复到默认状态 vovlap 1 2 进行 OVERLAP 布尔操作 会生成更多的体 pnum volu 1 打开实体编号 view 3 1 1 定义显示角度 type 4 title Volumes used in building pipe tank junction vplot 显示实体 vdele 3 4 1 删除多余实体 1 表示删除体及其以下的所有图元 划分网格 asel loc z z1 选择罐上 Z Z1 的面 asel a loc y 0 添加选择罐上 Y 0 的面 cm aremote area 创建名为 AREMOTE 的面组 pnum area 1 pnum line 1 title lines showing the portion being modeled aplot noerase 多窗口显示问题 lplot erase accat all 组合罐远端的面及线 为映射划分网格做准备 lccat 12 7 lccat 10 5 lesize 20 4 在接管壁厚方向分 4 等分 lesize 40 6 在接管长度方向分 6 等分 lesize 6 4 在罐壁厚方向分 4 等分 allsel 选择 EVERYTHING esize 0 4 设定默认的单元大小 mshape 0 3d 选择 3D 映射网格 mshkey 1 save 保存数据文件 vmesh all 划分网格 产生节点与单元 pnum defa title elements in portion being modeled eplot 显示单元 finish 加载求解 solu antype static 定义为稳态分析 nropt auto 设置求解选项为 Program chosen Newton Raphson tunif 450 设定初始所有节点温度 csys 1 变为柱坐标 nsel s loc x ri1 选择罐内表面的节点 sf all conv 250 144 450 定义对流边界条件 cmsel aremote 选择 AREMOTE 面组 Utility Menu Select Comp Assembly Select Comp Assembly nsla 1 选择属于 AREMOTE 面组的节点 Selects those nodes associated with the selected areas d all temp 450 定义节点温度 wprota 0 90 将工作平面旋转到垂直于接管轴线 cswpla 11 1 创建局部柱坐标 nsel s loc x ri2 选择接管内壁的节点 sf all conv 2 100 定义对流边界条件 allsel 选择 EVERYTHING pbc temp 1 显示所有温度约束 psf conv 2 显示所有对流边界 title Boundary conditions nplot 显示节点 wpstyle defa 工作平面恢复默认状态 csys 0 变为直角坐标 autots on 打开自动步长 nsubst 50 设定子步数量 kbc 0 设定为阶越 outpr nsol last 设置输出 solve 进行求解 finish 进入后处理 post1 title Temperature contrours at pipe tank junction plnsol temp 显示温度彩色云图 finish exit all 参考其它文章 稳态热分析综述 6 6 接触分析接触分析 钢板卷制成圆筒实例钢板卷制成圆筒实例 Finish clear config nres 2000 设置最大子步数 prep7 Et 1 shell43 用于划分上下辊 Et 2 solid186 Mp ex 1 2e11 Mp dens 1 7800 Mp prxy 1 0 3 Mp ex 2 2e11 Mp dens 2 7800 Mp prxy 2 0 3 Tb bkin 2 1 Tbtemp 0 Tbdata 240e6 2e8 定义材料的屈服极限 切变模量 R 1 0 02 壳单元的厚度 Cylind 0 19 0 3 2 3 360 创建上辊 Cylind 0 095 0 3 360 Cylind 0 095 2 6 2 3 360 Wpoff 0 27 0 4 Cylind 0 15 0 3 2 3 360 Cylind 0 075 0 3 360 Cylind 0 075 2 6 2 3 360 Wpoff 0 54 Cylind 0 15 0 3 2 3 360 Cylind 0 075 0 3 360 Cylind 0 075 2 6 2 3 360 Vglue all Allsel Vdele all Csys 0 Hptcreate area 5 coord 在辊端面中心处创建硬点 Hptcreate area 10 coord 2 6 Hptcreate area 34 coord 0 27 0 4 2 6 Hptcreate area 29 coord 0 27 0 4 Hptcreate area 22 coord 0 27 0 4 2 6 Hptcreate area 17 coord 0 27 0 4 Mat 1 Type 1 Real 1 Esize 0 05 Smrtsize 8 Mshape 1 Amesh all Csys 4 Block 0 7 1 0 17 0 2 0 5 2 1 Smrtsize off Lesize 96 1 Mat 2 Type 2 Mshkey 1 Mshape 0 Esize 0 05 Vmesh all view 1 1 1 1 eplot pnum area 1 aplot Et 4 targe170 Et 5 conta174 Keyopt 5 12 1 Keyopt 5 5 3 R 10 0 01 0 1 R 11 0 01 0 1 R 12 0 01 0 1 Asel s 14 Nsla s 1 Nsel u loc z 0 8 1 8 Mat 3 Real 10 Type 4 Esurf Allsel Asel s 3 4 1 Nsla s 1 Nsel u loc z 0 9 1 7 Nsel u loc z 0 3 0 4 Nsel u loc z 2 2 2 3 Mat 3 Real 10 Type 5 Esurf Allsel Asel s 13 Nsla s 1 Nsel u loc z 0 8 1 8 Mat 4 Real 11 Type 4 Esurf Allsel Asel s 27 28 1 Nsla s 1 Nsel u loc z 0 9 1 7 Nsel u loc z 0 3 0 4 Nsel u loc z 2 2 2 3 Mat 4 Real 11 Type 5 Esurf Allsel Asel s 13 Nsla s 1 Nsel u loc z 0 8 1 8 Mat 4 Real 12 Type 4 Esurf Allsel Asel s 15 16 1 Nsla s 1 Nsel u loc z 0 9 1 7 Nsel u loc z 0 3 0 4 Nsel u loc z 2 2 2 3 Mat 4 Real 12 Type 5 Esurf Allsel Finish solu Csys 0 D all uz D all rotx D all roty Nsel s loc x Nsel r loc y D all uy 0 08 D all ux D all uz Allsel Csys 4 Asel s loc y 2 0 16 Nsla s 1 D all all Allsel Antype trans Neqit 100 Lnsrch on Nlgeom on Time 1 Autot on Nsubst 30 80 20 Kbc 0 outres all all Solve Finish clear config nres 2000 设置最大子步数 prep7 Et 1 shell43 用于划分上下辊 Et 2 solid186 Mp ex 1 2e11 Mp dens 1 7800 Mp prxy 1 0 3 Mp ex 2 2e11 Mp dens 2 7800 Mp prxy 2 0 3 Tb bkin 2 1 Tbtemp 0 Tbdata 240e6 2e8 定义材料的屈服极限 切变模量 R 1 0 02 壳单元的厚度 Cylind 0 19 0 3 2 3 360 创建上辊 Cylind 0 095 0 3 360 Cylind 0 095 2 6 2