瞬态热弹性问题的混合边界元法自适应网格式.pdfVIP

第九届全国建筑工程计算机应用学术套议论文集 黑龙江·哈尔滨1998 瞬态热弹性问题的混合边界元法自适应网格式 王建省 (北方工业大学建筑学院) 【摘要】结合建筑工程中瞬态热弹性／热载荷问题，提出r利用混合边界元新的处理方法 自适应网提 炼格式。茸先建立瞬态热弹性问题的数学模掣，并刘不同类型边界的适应条件进行讨论，分离了混合边 界单元各子项，给出了网提炼的形式及分析结果并和已知解进行，比较．指…了本文分析的意义及有效 性。 一、引言 建筑工程中由于冲击及温度的变化，瞬i蠡热载荷造成的危害现象是普遍存在的。早在 本世纪七十年代就引起人们的重视∽法国CETIM在这方面进行了一系列的工作[2,31，并 研制出了当时先进的cA．sT．OR软件用以解决工程中的有关问题，其计算的精度和有限元 秆j比有了较大提高Pj。以后在三维弹性及热传导方面，瞬态热载荷的有关研究均取得了重 要成就【=”。在这些研究中，均使用常规单7c格式。zienkjewicz【。J提出了自适应网格提炼 方法，用其中P提炼法进行了有关汁算，RankVl总结了常规边界单元法的网点函数关系， 作者分别对热弹性及弹塑性问题进行门_开究，给出了类似的结果189J。本文将基于目前计 算机的发展及Fortr．PowerStation4．xFor nin95的曰益广泛应用，提出分析瞬态热 载荷In1N-的混合边界元方法，并结合问题的特点给出自适应网格的形式。 二、瞬态热弹性问题数学模型 对各向同性、力学性能均匀的材料，若材料性能为定常，基本控制方程为[21 量，y为Poisson比，上标e，0分别表示与弹性、温度有关的位移项，口为与位移和时间 有关的温度，记包含弹性应变及热应变的伞应变为Ph，材料服从Hooke’定律，有关系 38 苎垄壁全璺壅苎兰堡苎墨垫皇璺兰查全垡堡圭叁 墨垄兰：堡鱼苎!塑! (2) e8=05@“：+％u。e)．C0=0．5(4‰0+瓦坼0) ‘‘e e。。右k。+ JPMkm 嗡(3) okm七i。——■勺k。=艽j‰+黑酝：∥％屯+22(2 温度0与热通量≯共同影响的混台边界(记为F“)，此种边界存在以下形式的关系 叫譬一T4)+妒(g，功=0 (4) 式中，T为绝对温度，e下标表示边界外部温度，q表示边界位置。事实上，还存在一种情 况，即仅有位移和面力的部分，现记为厂“，r，，其关系为 “(g，们=疗，吼，心％=7 (5) 任意两种边界的井集即ar构成不超过全空I刚维数的相应维空间．整个边界的构成也可 以并集的形式表示，作以卜J标记 r=r。UF9Ur“UF“UFj (6) 7=甲 且有关系 F‘Ur (7) 式中，，，J=只以厅，Ⅳ，．厂掣为v维窄间(v≤rt)，其中全空问为n维。 三、混合边界单元及自适应网格式 记所考虑的热弹性混合控制体为Q。 边界L，应力函数在F7上满足平衡条件及边 界条件。位