联系我们 Contact Us
技术热线:(+86)13592550753
电话:(0371)86528660
电子邮件:fealab@163.com
如何提升创新和竞争能力?
技术应用文章>>>
- CAE咨询--有限元分析技术服...
- 有限元培训_ANSYS培训_W...
- ANSYS新旗舰版本ANSYS...
- 基于CFX的涡流燃烧器气体流场...
- LINK180—三维有限应变杆...
- 3维的显式的结构实体单元
- PLANE121--2 维 8...
- 8 节点轴对称-谐波热单元-P...
- 2 维 8 节点热实体单元-P...
- 4 节点轴对称谐波热单元
- 2 维热-电耦合实体单元
- SOLID65单元性质
- SOLID65单元性质
- 3-D 各向异性结构实体单元
- 2 维 8 节点磁实体单元-P...
- 3-D结构实体单元
- 电流源SOURC36
- LINK32--二维热传导杆单...
- 弹簧-阻尼器Spring-Da...
- LINK8单元描述
- SOLID5--三维耦合场体单...
- Beam4 单元描述
- Beam3二维弹性单元特性
- PLANE2
- 液压支架参数的确定
PLANE78用作轴对称环单元,具有 3 维热传导能力。本单元有 8 个节点,每个节点只有一个自由度 – 温度 TEMP。本单元是轴对称单元 PLANE77 的推广 – 允许非轴对称载荷。不同载荷情况在 有非对称载荷的轴对称单元 中说明。
这一 8 节点单元具有一致的温度形函数,可以较好地适应具有曲线边界的模型。
本单元适用于 2 维、轴对称,稳态或瞬态热分析。关于本单元的更多细节见 ANSYS 公司理论手册 中的PLANE78。如果包含本单元的模型还要用于结构分析,应该用等价的结构单元 (如 PLANE83) 替换本单元。
![]() |
图 78.1 PLANE78 单元几何
在图 78.1: "PLANE78 单元几何" 中给出了本单元的几何形状,节点位置和坐标系。输入数据基本与 PLANE77 相同,可见 "PLANE77 输入数据"。本单元的输入数据还包括谐波数 (MODE 命令中的 MODE 项) 和对称条件 (MODE 命令中的 ISYM项)。如果 MODE = 0 且 ISYM = 1,本单元与轴对称情况的PLANE77 类似。如果 MODE = 1,假定沿着某直径的温度为零。MODE 和 ISYM 参数说明了温度分布的类型,其细节在 有非对称载荷的轴对称单元 中讨论。
单元载荷在 "节点和单元载荷" 中说明。谐波变化的对流换热或热流密度 (二者不能同时) 可以作为单元边界上的面载荷输入,如图 78.1: "PLANE78 单元几何" 中带圆圈数字所示。谐波变化的生热率可以作为单元节点上的体载荷输入。如果输入了节点 I 处的生热率HG(I),但未给出其它节点处的生热率,则默认等于 HG(I)。如果输入了所有角节点的生热率,各中间节点的生热率默认为其相邻角节点生热率的平均值。
在 "PLANE78 输入汇总"中给出了本单元输入数据的汇总。关于本单元输入数据的一般说明,见 "单元输入"。
PLANE78 输入汇总
节点
I, J, K, L, M, N, O, P
自由度
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), HG(M), HG(N), HG(O), HG(P)
谐波数 – 由 MODE 命令输入
载荷条件
由 MODE 命令 ISYM 选项输入
1 – 对称载荷
-1 – 反对称载荷
求解能力
生死单元
关键选项
KEYOPT(1) – 比热矩阵
0 – 一致比热矩阵
1 – 对角线比热矩阵
与单元有关的结果输出有两种形式:
• 包括在整个节点解中的节点温度。
• 附加的单元输出,见表 78.1 "PLANE78 单元输出定义"。
对流热流以流出单元为正。施加的热流以流入单元为正。单元输出的方向平行于单元坐标系。侧面面积和热流密度都是对整个 360°圆周的。关于谐波单元的更多信息见 有非轴对称载荷的轴对称单元。在 "结果输出" 中给出了对于结果输出的一般说明。查看的方法见 "ANSYS 基本分析指南"。
单元输出定义表使用如下标记:
在名称列中的冒号 (:) 表示该项可以用分量名方法 [ETABLE, ESOL] 处理;O 列表示该项可用于 Jobname.OUT 文件;R 列表示该项可用于结果文件。
无论 O列或 R 列,Y 表示该项总是可用的,数字表示表的一个注解,其中说明了使用该项的条件;而减号 "-" 表示该项不可用。
表 78.1 PLANE78 单元输出定义
名称 | 定 义 | O | R |
EL | 单元号 | Y | Y |
NODES | 节点- I, J, K, L, M, N, O, P | Y | Y |
MAT | 材料号 | Y | Y |
MODE | 载荷谐波数 | Y | Y |
VOLU | 体积 | Y | Y |
XC,YC | 结果输出点位置 | Y | 3 |
HGEN | 生热率HG(I), HG(J), HG(K), HG(L), HG(M), HG(N), HG(O), HG(P) | Y | - |
TG:X, Y, SUM, Z | 单元中心处的温度梯度分量和矢量和 (X 和 Y) | 1 | 1 |
TF:X, Y, SUM | 单元中心处的热流密度 (热流率/截面积) 分量和矢量和 (X 和 Y) | 1 | 1 |
FACE | 侧面标记 | 2 | 2 |
NODES | 侧面节点 | 2 | 2 |
AREA | 侧面面积 | 2 | 2 |
HFILM | 膜系数 | 2 | 2 |
TAVG, TBULK | 两端节点温度平均值,液体介质温度,都在峰值点计算 | 2 | 2 |
HEAT RATE | 由对流产生的穿过边界的热流率 | 2 | 2 |
HEAT RATE/AREA | 由对流产生的穿过边界单位面积的热流率 | 2 | 2 |
HFAVG | 边界的平均膜系数 | - | 2 |
TBAVG | 边界的平均介质温度 | - | 2 |
HFLXAVG | 由输入热流产生的穿过边界单位面积的热流率 | - | 2 |
HFLUX | 边界各节点的热流密度 | 2 | 2 |
1 梯度和热流峰值所在角度为 THETA = 0 和 THETA = 90 ÷ MODE;
2 仅在有面载荷输入时输出;
3 仅用于 *GET 命令,给出单元中心处结果。
表 78.2 "PLANE78 输出项和序列号" 列出了可以通过 ETABLE 命令,采用序列号方法输出的内容列表。更多信息见"ANSYS 基本分析指南" 中 一般后处理 (POST1) 部分和本手册中关于"输出项和序列号表"部分。在表 78.2 "PLANE78 输出项和序列号" 中使用如下标识符:
Name
与表 78.1: "PLANE78 单元输出定义" 中相同定义的输出量;
Item
用于 ETABLE 命令的预先定义的输出项;
FCn
单元第 n 条边输出项的序列号。
表 78.2 PLANE78 输出项和序列号
输出量名称 | 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 |
PLANE78 假设和限制
· 单元的面积必须非负和非零;
· 对于轴对称问题,本单元必须位于总体坐标的X-Y 平面中,总体 Y 轴必须是对称轴,如图 78.1 "PLANE78 单元几何" 所示。
· 轴对称结构必须在 +X 象限中建模;
·删除某边的中间节点,意味着沿该边的温度呈线性变化,而不再是抛物线变化的。关于中间节点的使用,见 ANSYS 建模和分网指南 中的 二次单元 (中间节点);
·如果将这一热单元替换为有表面应力的结构单元 (PLANE83),本热单元的方向应该使 IJ 边 ( 以及KL 边 – 如果节点 K 和 L 不相同) 为自由表面。单元的自由表面(既不与其它单元相邻,又没有边界约束) 假设为绝热的;
·具有较小的积分时间步长和严重的表面热梯度的瞬态热问题要求该表面划分精细的网格。
·假定温度相关材料特性(包括膜系数)是轴对称的,即使温度是谐波分布。
·如果 MODE = 0,材料特性按上一子步得到的温度进行计算 (或对第一步,按 TUNIF 计算);
·如果 MODE > 0,材料特性按上一MODE = 0 子步得到的温度进行计算;如果不存在 MODE = 0 的子步,则按零角度的温度计算。