*if,i,eq,1,then
save! save original node locations at the first run *endif fini /post1
cmsel,s,gasket nsort,u,sum,1,1
*get,dismax(i),sort,0,max! Get the maximum displacement value strcri(i)=toler*dismax(i) alls fini /prep7
mkey=2! Select level of mesh morphing for fluid
damorph,area2, ,mkey! Perform morphing of \!----------------
!!!!!Create element plot and write it in file gasket.grph fini /prep7 et,1,42 asel,s,,,1,3 esla,s
/Title, EPLOT ater DAMORPH,area2, ,%mkey% step number %i% eplot alls
!----------------- cmsel,s,gasket
nlist! List updated coordinates of bottom of gasket for comparison.. /com, +++++++++ UPDATED gasket coordinates -------- alls fini
/assign,esav /assign,emat
!!!!Checking convergence criteria imax= i index(i)=i *if,i,gt,1,then
strcri(i)=abs(dismax(i)-dismax(i-1))-toler*dismax(i-1) *if,strcri(i),le,0,then strcri(i)=0
*exit! Stop looping if convergence is reached *endif *endif *enddo
!!!!! End of the Computational loop
save! Nodal coordinates of deformed geometry are saved !!!!! Convergence printout *vwrite
(/\
/nopr *vlen,imax
*vwrite,index(1),dismax(1),strcri(1) (f7.0,2e17.4) finish
!!!!! Postprocessing of the results !!!1. Flotran results. physics,read,fluid /post1 set,last
/Title, Flotran: Streamlines Near Gasket plnsol,strm
/Title, Flotran: Pressure Contours plnsol,pres fini
!!!2. Structural results. physics,read,struc /post1 set,last
upcoord,-1! Return original node positions changed by morphing /Title, Structural results: Von Mises Stress plnsol,s,eqv,1,1 fini
!/exit,nosave
2.10使用物理环境方法进行热—传导分析的实例
本例示例了一个瞬态感应加热问题。该问题演示了顺序求解谐波电磁分析及有重启动的瞬态传热分析。
2.10.1问题描述
长钢坯通过感应线圈经历表面加热导致快速升温的过程。线圈紧贴钢坯表面放置并通过高频交变大电流。AC交流电感应加热钢坯，最显著的是表面，其温度快速升高。
简化的模型只考虑长钢坯的有限长度，如图2-16简化为一维问题进行研究。
图2-16 轴对称一维感应加热区域片
2.10.2步骤
钢坯将加热到700°C。对热问题及电磁场问题必须考虑温度相关的材料属性。必须顺序求解该问题。首先做一个交流谐波电磁场分析然后做一个瞬态热分析。而且，必须在不同时间间隔重复电磁分析来校正影响求解的温度相关属性及对钢坯的热载荷。如图2-17所示。
图2-17 求解流程图
感应加热问题求解步骤如下。
2.10.2.1 步骤1：建立属性关系
对模型区域定义如表2-6的属性关系
表2-6物理环境属性
区域 钢坯 线圈 空气 钢坯表面 2.10.2.2步骤2：建立模型
类型 1 2 2 3 材料 2 3 1 2 实常数 1 1 1 3 建立整个区域的模型。对不同区域分配属性。（钢坯表面用于定义热辐射用的表面效应单元。与实体区域的处理不同）。
2.10.2.3步骤3：创建电磁物理环境
通过定义如表2-7所示单元类型及材料属性定义电磁物理环境。