例10连铸过程的热力学分析过程.doc

PAGE  PAGE 1 例10 连铸过程的热力学分析过程 10.1 问题描述 利用例9中温度场分析网格，在右端网格的单元边界施加压力为89300N/m2的边界条件。网格左端轴向位移和径向位移均设为零。并给定无初应力的参考温度为1162℃。 计算热应力时所需要的各时刻的温度分布从前面温度分析的后处理文件中读入。 分析500s的时间历程期间的热应力变化。设置加载历程的时间步长控制参数时，允许每个增量步的最大温度变化为100℃。 10.2 网格模型  MAIN FILES NEW OK MODEL: OPEN e10x5a5.mud OK SAVE AS e10x5b5.mud OK RETURN PLOT （控制几何显示的格式） ELEMENTS: SOLID(off) （关闭实体显示） REGEN  10.3定义边界条件 MAIN BOUNDARY CONDITIONS NEW (设定轴向的位移为0) NAME fix_axis MECHANICAL FIXED DISPLACEMENT （定义位移边界条件，并将其施加在节点上） ON Y DISPLACEMENT OK NODES:ADD 1 12 23 34 45 56 67 78 89 100 111 122 133 144 155 166 177 188 199 210 221 232 243 254 265 276 287 298 309 320 331 342 353 364 END LIST(#) (施加Y向的位移约束如下图所示) NEW （设定网格左端的位移为0） NAME fix_front FIXED DISPLACEMENT ON X DISPLACEMENT OK NODES:ADD 1 2 3 4 5 6 7 8 9 10 11 END LIST(#) NEW （在右边界上施加压力如下图） NAME pressure EDGE LOAD PRESSURE(TOP) 89300 OK EDGES:ADD 321:1 322:1 323:1 324:1 325:1 326:1 327:1 328:1 329:1 330:1 NEW （计算温度载荷时的边界条件） NAME thermal MORE STATE VARIABLE （定义状态变量将其施加在单元上） ON METHOD: POST FILE(on) POST FILE e10x5a5_heatransf.t16 OK FIRST INCREMENT 1 # STEP 413 OK ELEMENTS:ADD （施加在单元体上） ALL:EXIST. 10.4定义初始条件 MAIN （将无应力的初始温度定为1162℃） INITIAL CONDITIONS NEW NAME inistate MECHANICAL STATE VARIABLE ON (on) VALUE 1162 OK ELEMENTS:ADD ALL:EXIST. 10.5定义材料参数 MAIN MATERIAL PROPERTIES NEW NAME mechanical ISOTROPIC OK TABLES NEW （给定杨氏模量与温度的变化曲线） NAME youngsm TABLE TYPE temperature OK ADD POINT 0 130.2 100 126.7 200 121.8 300 116.2 400 107.1 500 98.0 600 91 700 84 800 77 900 68.5 1000 60 1083 53 1500 1 FIT SCALE （将表格中的数据点沿X/Y方向放大/缩小，实际上是Y方向扩大了1e9倍） 1 1.e9 FIT （将所有数据点显示出来） NEW （指定柏松比-温度关系） NAME poissr TABLE TYPE temperature OK ADD POINT 0 0.35 20 0.35 1000 0.45 1500 0.45 F