! 2.5 端部有间隙的杆的热膨胀
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH ! 退出以前模块
/CLEAR, START ! 清除系统中的所有数据, 读入启动文件设置
! （1）定义分析类型、单元类型、平面单元厚度和材料的力学、热学参数
/FILNAME, EX2.5 ! 指定当前工程的文件名
/PREP7 ! 进入前处理模块
/TITLE, EX2.5, THERMAL EXPANSION TO CLOSE A GAP AT A RIGID SURFACE.
ANTYPE, STATIC ! 定义分析类型为静力分析
ET, 1, CONTAC26 ! 定义第1类单元为二维接触单元CONTACT26
R, 1, 2.625E15 ! 定义接触单元的接触刚度 K = 2.625E15
ET, 2, PLANE42, , , 3 ! 定义第2类单元为二维带厚度的平面应力单元PLANE42
R, 2, 1 ! 定义第2类实常数, 平面单元的厚度1 mm
MP, EX, 1, 72.4E3 ! 定义第1类材料的弹性模量 EX = 72.4 GPa
MP, ALPX, 1, 2.25E-5 ! 定义第1类材料的热膨胀系数 ALPX = 2.25E-5
MP, NUXY, 1, 0 ! 定义第1类材料的泊松比 NUXY = 0
! （2）定义节点信息
N, 1, 2000, 1000 ! 在（2000 mm, 1000 mm）定义1号节点
N, 2, 3000, 1000 ! 在（3000 mm, 1000 mm）定义2号节点
N, 3, 3000, 4000 ! 在（3000 mm, 4000 mm）定义3号节点
N, 4, 2000, 4000 ! 在（2000 mm, 4000 mm）定义4号节点
N, 11, 1000, 998 ! 在（1000 mm, 998 mm）定义11号节点
N, 12, 4000, 998 ! 在（4000 mm, 998 mm）定义12号节点
! （3）定义单元信息
E, 1, 11, 12 ! 过11、12号节点定义接触单元
E, 2, 11, 12
TYPE, 2 ! 选择单元类型号为2
REAL, 2 ! 选择实常数类型号为2
E, 1, 2, 3, 4 ! 过节点1、2、3、4定义平面单元
TREF, 21 ! 设定起始参考温度21度
BFUNIF, TEMP, 77 ! 将所有节点的温度均匀提升到77度
! （4）定义位移约束
D, 3, ALL, , , 4 ! 约束3号节点的所有位移, 并循环到4号节点
D, 1, UX, , , 2, 1 ! 约束1号节点的UX位移, 并按增量1循环到2号节点
OUTPR, BASIC, LAST ! 选择基本输出选项, 直到最后一个荷载步
FINISH ! 退出前处理模块PREP7
! （5）在求解模块开始求解
/SOLU ! 进入求解模块
SOLVE ! 开始求解
FINISH ! 退出求解模块SOLU
! （6）在一般后处理模块, 定义水平应力和铅直应力单元表, 并提取3号单元的应力结果
/POST1 ! 进入一般后处理模块
ETABLE, STRSX, S, X ! 定义X方向的应力单元表STRSX
ETABLE, STRSY, S, Y ! 定义Y方向的应力单元表STRSY
*GET, STRSSX, ELEM, 3, ETAB, STRSX ! 从单元表STRSX提取3号单元X方向的应力, 存入STRSSX
*GET, STRSSY, ELEM, 3, ETAB, STRSY ! 从单元表STRSY提取3号单元Y方向的应力, 存入STRSSY
*STATUS, PARM ! 显示所有内存变量的值
FINISH ! 退出后处理模块
