PLANE75 单元 说明PLANE75用作轴对称环单元,具有 3 维热传导能力。本单元有 4 个节点,每个节点只有一个自由度 – 温度 TEMP。本单元是轴对称单元 PLANE55 的推广 – 允许非轴对称载荷。不
PLANE75 单元 说明
PLANE75用作轴对称环单元,具有 3 维热传导能力。本单元有 4 个节点,每个节点只有一个自由度 – 温度 TEMP。本单元是轴对称单元 PLANE55 的推广 – 允许非轴对称载荷。不同载荷情况在 剪切变形 (有非对称载荷的轴对称单元 ?) 中说明。
本单元可以用于 2 维、轴对称、稳态或瞬态热分析。关于本单元的更多细节见 ANSYS 公司理论手册 中的 PLANE75。如果含有本单元的模型也进行结构分析,需要用等价的结构单元 (如 PLANE25) 替代本单元。类似的,有中间节点的单元是 PLANE78。
图 75.1 PLANE75 单元几何
PLANE75 输入数据
在图 75.1: "PLANE75 单元几何" 中给出了这一轴对称热实体单元的几何形状,节点位置和坐标系。输入数据与 PLANE55 单元基本相同,见 "PLANE55 Input Data"。本单元的输入数据还包括谐波数 (MODE) 和对称条件 (ISYM),用 MODE 命令输入。如果 MODE = 0 且 ISYM = 1,本单元与轴对称情况的 PLANE55类似。MODE 和 ISYM 参数说明了温度分布的类型,并在剪切变形 (有非对称载荷的轴对称单元 ?) 中讨论。
单元载荷在 "节点和单元载荷" 中说明。谐波变化的对流换热或热流密度 (二者不能同时) 可以作为单元边界上的面载荷输入,如图 75.1: "PLANE75 单元几何" 中带圆圈数字所示。谐波变化的生热率可以作为单元节点上的体载荷输入。如果输入了节点 I 处的生热率HG(I),但未给出其它节点处的生热率,则默认等于 HG(I)。
在 "PLANE75 输入汇总"中给出了本单元输入数据的汇总。关于本单元输入数据的一般说明,见 "单元输入"。
PLANE75 输入汇总
节点
I, J, K, L
自由度
TEMP
实常数
无
材料性能
KXX, KYY, KZZ, DENS, C, ENTH
面载荷
对流换热 -
边 1 (J-I), 边2 (K-J), 边3 (L-K), 边4 (I-L)
热流密度 -
边 1 (J-I), 边2 (K-J), 边3 (L-K), 边4 (I-L)
体载荷
生热率 - HG(I), HG(J), HG(K), HG(L)
谐波数 – 由 MODE 命令输入
载荷条件
由 MODE 命令 ISYM 选项输入
1 – 对称载荷
-1 – 反对称载荷
求解能力
生死单元
关键选项
无
PLANE75 输出数据
与单元有关的结果输出有两种形式:
• 包括在整个节点解中的节点温度。
• 附加的单元输出,见 表 75.1 "PLANE75 单元输出定义"。
对流热流以流出单元为正。施加的热流以流入单元为正。单元输出的方向平行于单元坐标系。侧面面积和热流密度都是对整个 360°圆周的。关于谐波单元的更多信息见 有非轴对称载荷的轴对称单元。在 "结果输出" 中给出了对于结果输出的一般说明。查看的方法见 "ANSYS 基本分析指南"。
单元输出定义表使用如下标记:
在名称列中的冒号 (:) 表示该项可以用分量名方法 [ETABLE, ESOL] 处理;O 列表示该项可用于 Jobname.OUT 文件;R 列表示该项可用于结果文件。
无论 O列或 R 列,Y 表示该项总是可用的,数字表示表的一个注解,其中说明了使用该项的条件;而减号 "-" 表示该项不可用。
表 75.1 PLANE75 单元输出定义
名称
定 义
O
R
EL
单元号
Y
Y
NODES
节点- I, J, K, L
Y
Y
MAT
材料号
Y
Y
VOLU
体积
Y
Y
XC,YC
结果输出点位置
Y
3
HGEN
生热率HG(I), HG(J), HG(K), HG(L)
Y
-
MODE
载荷谐波数
Y
-
TG:X, Y, SUM, Z
单元中心处的温度梯度分量和矢量和 (X 和 Y)
1
1
TF:X, Y, SUM
单元中心处的热流密度 (热流率/截面积) 分量和矢量和 (X 和 Y)
1
1
FACE
侧面标记
2
-
NODES
侧面节点
2
-
AREA
侧面面积
2
2
TAVG, TBULK
两端节点温度平均值,液体介质温度,都在峰值点计算
2
2
HEAT RATE
由对流产生的穿过边界的热流率
2
2
HEAT RATE/AREA
由对流产生的穿过边界的单位面积热流率
2
-
HFAVG
边界的平均膜系数
-
2
TBAVG
边界的平均介质温度
-
2
HFLXAVG
由输入热流产生的穿过边界的单位面积热流率
-
2
HFLUX
边界各节点的热流密度
2
-
注解:
1 梯度和热流峰值所在角度为 THETA = 0 和 THETA = 90 ÷ MODE;
2 仅在有面载荷输入时输出;
3 仅用于 *GET 命令,给出单元中心处结果。
表 75.2 "PLANE75 输出项和序列号" 列出了可以通过 ETABLE 命令,采用序列号方法输出的内容列表。更多信息见"ANSYS 基本分析指南" 中 一般后处理 (POST1) 部分和本手册中关于 "输出项和序列号表"部分。在 表 75.2 "PLANE75 输出项和序列号" 中使用如下标识符:
Name
与表 75.1: "PLANE75 单元输出定义" 中相同定义的输出量;
Item
用于 ETABLE 命令的预先定义的输出项;
FCn
单元第 n 条边输出项的序列号。
表 75.2 PLANE75 输出项和序列号
输出量名称
ETABLE 和 ESOL 命令输入
Item
FC1
FC2
FC3
FC4
AREA
NMISC
1
7
13
19
HFAVG
NMISC
2
8
14
20
TAVG
NMISC
3
9
15
21
TBAVG
NMISC
4
10
16
22
HEAT RATE
NMISC
5
11
17
23
HFLXAVG
NMISC
6
12
18
24
PLANE75 假设和限制
· 单元的面积必须非负和非零;
· 对于轴对称问题,本单元必须位于总体坐标的X-Y 平面中,总体 Y 轴必须是对称轴,如 图 75.1 "PLANE75 单元几何" 所示。
· 轴对称结构必须在 +X 象限中建模;
· 如果定义节点号 K 和 L 相同,可以形成三角形单元 (见 "三角形、金字塔形和四面体单元");
·如果将这一热单元替换为有表面应力的结构单元 (PLANE25),本热单元的方向应该使 IJ 边 ( 以及KL 边 – 如果节点 K 和 L 不相同) 为自由表面。
·单元的自由表面(既不与其它单元相邻,又没有边界约束) 假设为绝热的;
·具有较小的积分时间步长和严重的表面热梯度的瞬态热问题要求该表面划分精细的网格。
·假定温度相关材料特性(包括膜系数)是轴对称的,即使温度是谐波分布。
·如果 MODE = 0,材料特性按上一子步得到的温度进行计算 (或对第一步,按 TUNIF 计算);
·如果 MODE > 0,材料特性按上一MODE = 0 子步得到的温度进行计算;如果不存在 MODE = 0 的子步,则按零角度的温度计算。 |