[toc]

问题说明

样块上下两端固定，在室温20℃下进行夹紧，分析其升温到150℃时的热应力。 采用间接法进行分析，温度场单元选择278，应力场单元为185 首先进行稳态热分析，然后读取温度场的数据，进行结构分析，具体的编程代码如下

温度场分析APDL

1/CLEAR
2/FILENAME TEST
3/PREP7
4
5!定义温度场实体单元
6ET,1,SOLID278  
7!定义材料1的导热系数
8MPTEMP,,,,,,,,  
9MPTEMP,1,0  
10MPDATA,KXX,1,,20
11
12!建立样件模型
13BLC5,0,0,0.2,0.1,0.05  
14
15!设置体1的网格尺寸
16CM,_Y,VOLU  
17VSEL, , , ,       1
18CM,_Y1,VOLU
19CMSEL,S,_Y  
20CMSEL,S,_Y1
21VATT,       1, ,   1,       0   
22CMSEL,S,_Y  
23CMDELE,_Y   
24CMDELE,_Y1  
25
26FLST,5,1,4,ORDE,1   
27FITEM,5,5   
28CM,_Y,LINE  
29LSEL, , , ,P51X
30CM,_Y1,LINE
31CMSEL,,_Y   
32LESIZE,_Y1, , ,20, , , , ,1
33 
34FLST,5,1,4,ORDE,1   
35FITEM,5,6   
36CM,_Y,LINE  
37LSEL, , , ,P51X
38CM,_Y1,LINE
39CMSEL,,_Y   
40LESIZE,_Y1, , ,10, , , , ,1
41  
42FLST,5,1,4,ORDE,1   
43FITEM,5,10  
44CM,_Y,LINE  
45LSEL, , , ,P51X
46CM,_Y1,LINE
47CMSEL,,_Y   
48LESIZE,_Y1, , ,5, , , , ,1  
49
50CM,_Y,VOLU  
51VSEL, , , ,       1
52CM,_Y1,VOLU
53CHKMSH,'VOLU'   
54CMSEL,S,_Y  
55VCLEAR,_Y1  
56VMESH,_Y1   
57CMDELE,_Y   
58CMDELE,_Y1  
59CMDELE,_Y2  
60
61
62!在一侧施加温度150
63FLST,2,1,5,ORDE,1   
64FITEM,2,5   
65/GO
66DA,P51X,ALL,150
67
68!设置其他表面的对流换热系数
69FLST,2,5,5,ORDE,3   
70FITEM,2,1   
71FITEM,2,-4  
72FITEM,2,6   
73/GO
74SFA,P51X,1,CONV,50,
75
76!设置分析类型为稳态求解
77FINISH  
78/SOL
79ANTYPE,0
80
81!设置温度偏移
82TOFFST,273,
83
84!写温度场物理分析文件
85FINISH  
86/PREP7  
87PHYSICS,WRITE,Thermal1, , ,
88
89!清空物理环境数据
90PHYSICS,CLEAR  
91!删除边界和载荷 
92LSCLEAR,ALL
93
94!转换单元类型，thermal to struc
95ETCHG,TTS
96
97!设置材料参数
98MPTEMP,,,,,,,,  
99MPTEMP,1,0  
100MPDATA,EX,1,,2.1E+011   !弹性模量
101MPDATA,PRXY,1,,0.3       !泊松比
102UIMP,1,REFT,,,              !参考温度
103MPDATA,ALPX,1,,2.35E-5      !热膨胀系数
104
105!定义结构边界条件--左右固定
106
107FLST,2,2,5,ORDE,2   
108FITEM,2,5   
109FITEM,2,-6  
110/GO
111DA,P51X,ALL,0
112
113!设置参考温度
114TREF,20,
115
116!求解设置
117FINISH  
118/SOL
119ANTYPE,0
120NLGEOM,1
121NSUBST,8,10,6   
122TIME,1
123
124!写结构物理分析文件
125FINISH  
126/PREP7  
127PHYSICS,WRITE,struct1, , ,  
128
129!读温度场物理分析文件
130PHYSICS,READ,THERMAL1  
131
132!求解
133FINISH  
134/SOL
135/STATUS,SOLU
136SOLVE  
137
138!设置图像显示
139/PLOPTS,INFO,1  !inf0为legend on
140
141!显示温度云图
142FINISH  
143/POST1  
144/EFACET,1   
145PLNSOL, TEMP,, 0
146
147!读取结构场物理分析文件
148FINISH  
149/PREP7  
150PHYSICS,READ,STRUCT1
151
152!定义温度边界条件
153LDREAD,TEMP,,, , ,'test_temp1','rth',' '
154
155!求解
156FINISH  
157/SOL
158/STATUS,SOLU
159SOLVE 
160
161!查看结果
162FINISH  
163/POST1  
164SET,LIST,999
165SET,FIRST
166/EFACET,1   
167PLNSOL, S,EQV, 0,1.0  
168

#结果