《有限元基本理论及应用》第七章 温度场的有限元分析.ppt

* * * * * * * 7.2 瞬态温度场分析 （4）瞬态热分析实例 TIME,3600 ! 设置瞬态分析的总时间 NSUBST,36,100,10 ! 设置瞬态分析的子步数 AUTOTS,1 ! 打开自动时间步长跟踪 TIMINT,1 ! 打开瞬态响应分析 KBC,1.0 ! 设置载荷为阶跃式 OUTRES,ALL,ALL ! 输出所有子步的结果 DDELE,ALL,TEMP ! 删除施加在所有节点上的温度载荷 SOLVE ! 开始瞬态分析求解 FINISH ! 退出求解器 ? /POST26 ! 进入时间历程后处理器 SOLU,2,DTIME,,DTIME ! 将每个子步的求解汇总结果赋予变量DTIME PLVAR,2, , , , , , , , , , ! 显示DTIME与时间的关系曲线 N1=NODE(0.2125,0.2075,0) ! 读取铜块中心点坐标的节点编号 N2=NODE(0.371,0.247,0) ! 读取铁块中心点坐标的节点编号 NSOL,3,N1,TEMP, ,T_COPPER ! 将铜块的温度赋予变量T_COPPer NSOL,4,N2,TEMP, ,T_IRON ! 将铁块的温度赋予变量T_IRON PLVAR,3, , , , , , , , , , ! 显示铜块的温度与时间的曲线 PLVAR,4, , , , , , , , , , ! 显示铁块的温度与时间的曲线 FINISH ! 退出后处理器 7.3 热应力的有限元分析 物体的热问题研究包括两个部分内容：(1) 传热问题研究，以确定温度场；(2) 热应力问题研究，即在已知温度场的情况下确定热应力与热变形。 传热问题所确定的温度将直接影响到物体的热应力分布，而后者对前者的耦合影响不大。 （1）基本理论 物体的物理方程由于增加了热膨胀量(正方向上的温度应变)而变为： 7.3 热应力的有限元分析 （1）基本理论 矩阵形式： 弹性问题的一般虚功原理为： 代入 热应力问题的虚功原理： 设节点的位移列阵为： 将单元内的力学参量都表达为节点位移的关系式： （1）基本理论 7.3 热应力的有限元分析 代入 （2）热应力分析过程与步骤 7.3 热应力的有限元分析 在ANSYS的多物理场(Multiphysics)模块中，包含有热—结构耦合的热应力计算。耦合场的计算方法有：载荷转换法(Load transfer)和直接法(direct)。 1）直接法 一般采用耦合场单元，其中单元的节点包含有所需耦合场的所有自由度。 如在热—结构耦合分析中，单元的节点所包含的自由度则由位移与温度组成。 通过计算包含所有必须项的单元矩阵或单元载荷向量来进行耦合。 在强耦合(计算单元矩阵)的线性分析中，耦合场界面采用一次迭次计算，对于弱耦合(计算单元载荷向量)的线性分析，则至少需要迭代2次才能取得一个耦合响应。 结构—热耦合既具有强耦合，也可以是弱耦合，但如果涉及到接触单元，则一定是强耦合。 （2）热应力分析过程与步骤 7.3 热应力的有限元分析 1）直接法 为了加速一个耦合瞬态分析的收敛，可以关闭所不考虑自由度的时间积分效应，比如如果惯性和阻尼不考虑，则可以关闭结构自由度的时间积分效应。 结构—热分析的单元主要有：PLANE13、SOLID5、SOLID98、PLANE223、SOLID226和SOLID227。