一种非连续介质中热传导过程数值模拟方法.PDF

第33 卷第1 期 Vol.33 No.1 工 程 力 学 2016 年 1 月 Jan. 2016 ENGINEERING MECHANICS 25 文章编号：1000-4750(2016)01-0025-07 一种非连续介质中热传导过程的数值模拟方法 1,2 3 李腊梅 ，冯 春 (1. 北京科技大学土木与环境工程学院，北京 100083 ；2. 中国冶金地质总局矿产资源研究院，北京 100025 ； 3. 中国科学院力学研究所，流固耦合系统力学重点实验室，北京 100190) 摘 要：传统热传导的分析基于连续模型，无法刻画热量在两个接触体之间的传递。该文提出了一种非连续介质 中热传导过程的数值计算方法，并编制了相应的C++计算程序。该方法首先将计算域离散为一系列的块体，块体 内部划分若干连续介质单元，块体边界设定为潜在接触界面，并利用半弹簧-半棱联合接触模型进行接触对的快速 检索及标记。每个块体内部的热传导采用传统连续模型进行计算(该文采用有限体积法) ，每个接触界面采用点面 接触型及棱棱接触型热传导模型进行描述。通过调整接触界面热传导系数中的刚度因子，可以实现接触界面对热 传导过程不同的抵抗效应。数值算例表明，该文所述方法可以较为准确地模拟热量在非连续介质中的传递过程； 接触界面上的刚度因子越大，界面对热传导过程的抵抗效应越小；当刚度因子大于100，界面抵抗效应基本消失， 非连续介质的计算结果与连续介质的计算结果完全一致；此外，接触界面上的刚度因子仅影响热传导的瞬态过程， 而不影响其稳态解。 关键词：非连续介质；散体；接触；热传导；刚度因子 中图分类号：TK124 文献标志码：A doi: 10.6052/j.issn.1000-4750.2014.05.0444 A NUMERICAL SIMULATION METHOD FOR HEAT CONDUCTION IN DISCONTINUOUS MEDIA LI La-mei1,2 , FENG Chun3 (1. School of Civil and Environmental Engineering, Beijing University of Science and Technology, Beijing 100083, China; 2. Institute of Mineral Resources Research, China Metallurgical Geology Bureau, Beijing 100025, China; 3. Key Laboratory for Mechanics in Fluid Solid Coupling Systems, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China) Abstract: Conventional heat conduction analysis is based on continuous model, so it cannot describe heat transport between two bodies with contact. A numerical simulation method for heat conduction in discontinuous media is proposed, and a C++ computing code is developed.