ANSYS Workbench 导热性随温度变化的杆的热应力分析
案例分析:
如下图所示,长杆的导热系数随着温度而变化。杆的两端受无摩擦表面的约束。在杆的一端(A端)施加 T ℃的温度。参考温度为5℃。在另一端,施加 h W/m2℃的恒定对流。环境温度为5℃。使用单元尺寸为2 m 的高级网格划分。
材料属性 | 几何特征 | 应用载荷 | |
E = 2 x 1011 Pa v = 0 α = 1.5 x 10-5 /℃ k = 0.038 x (1+0.00582*T) W/m ℃ | 长度 = 20 m 宽度 = 2 m 幅度 = 2 m | 背面温度 = 100 ℃ 换热系数 = 0.005 W/m2 ℃ 环境温度 = 5 ℃ 参考温度 = 5 ℃ | |
温度(℃) | 导热性(W/m ℃) | ||
5 | 3.91 x 10-2 | ||
800 | 0.215 |
1. 选择分析模块
此案例用于分析热应力问题问题,因此需要同时使用Steady-State Thermal和“Static Structural”模块进行耦合分析。
首先从Analysis System中拖拽Steady-State Thermal到Project Schematic中即可。
2. 定义材料属性
所用材料属性在材料库中无法找到,因此使用自定义属性。
— 右键Engineer Data并选择Edit/双击Engineer Data
— 新建材料并定义材料名称
— 从Toolbox中添加所需要的材料属性,材料的热膨胀系数以及弹性性能属性。
— 底部窗格中键入对应数值
其中热膨胀系数分为:割线热膨胀系数和瞬时(相切)热膨胀系数,这里使用前者。二者在此案例中对结果影响不大。
注意:热膨胀系数需要输入参考温度。
3. 建立几何数模
根据已知几何特征属性创建几何图形如下图所示。
右键Geometry并选择Edit Geometry in DesignModeler…
注意单位
4. 生成网格
— 双击Model/右键Model选择Edit
— 右键Mesh,添加Body Sizing,单元体尺寸设置为 2 m。网格类型设置为Nonlinear Mechanical。由于这里的单元尺寸和实体尺寸宽度一致。因此会生成一个方向上的单元体,故而需要将单元体的阶次降低,改为Linear。
注:在生成网格之前,可以先将材料属性赋予几何体。展开Geometry选中几何体 -->在相关详细信息视图中Material项关联菜单选中之前自定义的材料。
5. 添加边界条件和载荷
如下图所示,首先进行热力学分析,两端为Fixed Support;载荷已知
边界条件与载荷:Steady-State Thermal右键,插入“Temperature”和Convection约束条件并选择相应对象,设置对应换热系数。这里注意环境温度的设置。
可选:一般笔者习惯”Output Controls“均设置为”Yes“,即对象均输出至报告中,或用与后续求解。
求解(1):Solution右键,insert --> Thermal --> Temperature和Total Heat Flux添加温度分布。
温度分布求解结果
总热通量求解结果
6. 将热力学分析结果(温度分布)导入Static Structural中进行热应力分析
从Analysis System中拖拽Static Structural到Steady-State Thermal中的Solution,我们这里需要热力学模块的求解结果。
7. 进行相关分析设置
注意:这里需要激活Weak Spring,如下图。需要进行此步骤的原因是:弹簧的刚度比结构的刚度小很多个数量级,从而在不影响应力结果的情况下防止了刚体运动。如果不进行此设置,将会出现非期望的刚体运动报错。
Static Structural中的环境温度设置为5℃
杆的两端添加无摩擦支撑 Frictionless Support
8. 求解(2)
- Solution右键,insert --> Strain --> Thermal创建两个热应变求解并选中两个断端面.
- Solution右键,insert --> Deformation --> Directional并将坐标改为Z轴(与建模时模型的局部坐标系有关)
Z方向最大变形
最大热应变(z = 20) m/m
最大热应变(z = 0) m/m
ANSYS Mechanical结果比较
Results | Target | Mechanical | Error(%) |
---|---|---|---|
最低温度 ℃) | 38.02 | 38.015 | -0.013 |
最大热应变(z = 20) m/m | 0.000495 | 0.0004952 | 0.04 |
最大热应变(z = 0) m/m | 0.001425 | 0.001425 | 0 |
Z方向最大变形 | 0.00232 | 0.0023152 | -0.217 |
最大Z方向热通量(z = 20 m) (W/m2) | 0.165 | 0.16508 | 0.0485 |
注:上述Mechanical结果对于不同计算机或网格划分方式可能会有些许差别
Reference:
ANSYS Help 19.0