混凝土坝瞬态热传导分析的快速边界面法.pdf
《混凝土坝瞬态热传导分析的快速边界面法.pdf》由会员分享,可在线阅读,更多相关《混凝土坝瞬态热传导分析的快速边界面法.pdf(6页珍藏版)》请在咨信网上搜索。
1、水利水电科技进展,2016,36(1)摇 Tel:02583786335摇 E鄄mail:jz hhu. edu. cn摇 http:/ / www. hehaiqikan. cn 第 36 卷第 1 期 Vol. 36 No. 1 水 利 水 电 科 技 进 展 Advances in Science and Technology of Water Resources 2016 年 1 月 Jan. 2016 基金项目:国家自然科学基金(11172098);国家重点基础研究项目(973 计划)(2010CB328005) 作者简介:何建平(1966),男,高级工程师,博士,主要从事水电、机械
2、计算机仿真研究。 E鄄mail:jianping_he163. com DOI:10. 3880/ j. issn. 10067647. 2016. 01. 013 混凝土坝瞬态热传导分析的快速边界面法 何建平1,2,张见明1,李光耀1,郭帅平1 (1. 湖南大学汽车车身先进设计制造国家重点实验室,湖南 长沙摇 410082; 2. 中国电建集团中南勘测设计研究院有限公司,湖南 长沙摇 410014) 摘要:采用与快速算法相结合的边界面法,对大坝浇筑过程中的瞬态热传导问题进行仿真分析,预 测混凝土大坝施工过程中温度的分布及随时间的变化。 采用时间卷积法计算时间积分,通过Taylor 级数展开将
3、基本解中的空间变量和时间变量分离,一次性计算并存储空间变量的积分,提高了卷积 积分的计算效率。 算例仿真分析表明:边界面法的计算温度分布与有限元法计算结果高度吻合,将 边界面法应用于真实大坝浇筑过程瞬态热传导分析是可行的。 关键词:边界面法;有限元法;时间卷积法;瞬态热分析;混凝土坝 中图分类号:O241; O414摇 摇 摇 文献标志码:A摇 摇 摇 文章编号:10067647(2016)01007106 Fast boundary face method for transient heat conduction analysis of concrete dams/ / HE Jianpi
4、ng1, 2, ZHANG Jianming1, LI Guangyao1, GUO Shuaiping1(1. State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan University, Changsha 410082, China; 2. Power China Zhongnan Engineering Corporation Limited, Changsha 410014, China) Abstract: The boundary face method (BFM), co
5、mbined with a fast algorithm, was applied to transient heat conduction analysis of concrete dams during the pouring process, so as to predict the temperature distribution and its time history in concrete dams during the construction process. Using the time convolution method to compute the time inte
6、gration and the Taylor series expansion method to separate the spatial variables and time variables in fundamental solutions, the spatial variable integrals can be calculated and stored at the same time, improving the calculation efficiency of the convolution integral. A case study shows that the te
7、mperature distribution obtained by the fast BFM agrees with the results of the finite element method, demonstrating the feasibility of application of the fast BFM to transient heat conduction analysis during the actual dam construction process. Key words: boundary face method; finite element method;
8、 time convolution method; transient heat conduction analysis; concrete dam 摇 摇计算机辅助工程(CAE)对于推动产品研发和 指导工程建设具有重要意义。 计算机技术的快速发 展推动了 CAE 的发展。 但 CAE 仍面临许多难 题1,比如如何对复杂几何形体进行离散才能进行 有效计算,如何处理大规模工程问题的数值计算。 陆续涌现出的有限差分法(FDM)、有限体积法 (FVM)、有限元法(FEM)和边界元法2(BEM)都未 能很好地解决这些问题。 目前大多数 CAE 软件采 用的是有限元法。 有限元法需对整个求解域进行离 散
9、处理,其应力计算精度较位移精度低。 相比有限元法,另一种有吸引力的数值方法是 以边界积分方程为基础的边界元法3,该法不仅可 以减轻有限元法的网格离散困难,还具有高精度和 便于处理无限域问题、边界及界面问题、奇异性问题 等突出优点,因此在工程数值分析领域被公认是有 限元法的一个重要补充。 近年来一些关键技术的突 破使边界元法又发生了革命性的变化:淤稠密矩阵 稀疏化技术,包括快速多极算法4、分级矩阵和小 波算法(其中快速多极算法是 20 世纪十大算法之 一),可以将边界元满阵方程的求解复杂度降为线 性,从而使得边界元法在解题规模上优于有限元法; 于对偶互易法、径向积分法5和奇异边界法6的提 出大大
10、扩展了边界元法的应用范围,使得边界元法 不仅能够求解基本解存在的线弹性问题,还可以求 解非均质问题、非线性问题和时域问题;盂奇异积 分7鄄8、近奇异积分处理技术的发展9,使得边界元 17 水利水电科技进展,2016,36(1)摇 Tel:02583786335摇 E鄄mail:jz hhu. edu. cn摇 http:/ / www. hehaiqikan. cn 法可以基于三维理论求解超薄型结构。 然而,在传 统的边界元法中,三维 CAD 几何模型被离散成边界 元分析模型后,CAD 模型的原始几何信息基本被丢 掉,导致设计和分析变成了两个相互独立的过程。 由于 CAD 模型和分析模型的分离
11、,在边界元自适应 网格细分过程中,需要反复地与 CAD 系统进行交 互,而每个交互过程都是很烦琐的。 为了克服上述缺点,张见明等10在边界元法和 边界点法11的基础上提出了边界面法(BFM)。 在 边界面法中,对边界的数值积分和场变量的插值都 在边界曲面的二维参数空间中进行,CAE 分析直接 在 CAD 模型上进行,实现了复杂结构的 CAE 分析 自动化。 由于 CAE 模型与 CAD 几何模型融为一 体,不管网格离散有多么粗糙,分析模型在几何上都 是精确的,并且在自适应网格细分过程中,不再需要 与 CAD 系统反复地进行交互,使自适应分析变得简 单。 目前边界面法已经取得了许多研究成果,如参
12、 数曲面内的网格自动生成算法12、多域稳态热传导 分析13、弹性静力分析、温度应力分析、声场分 析14以及近奇异积分计算9。 张见明等15通过快 速多级子算法、自适应交叉拟合(ACA)和分级矩阵 (H鄄matrix)技术,成功地将计算量级从 O(N2)降至 O( NlgN), 并且基于 GPU 高性能并行计算, 在 CUDA 编程环境中实现了边界面法正则积分的并行 加速,使边界面法的大规模工程应用成为可能。 混凝土大坝施工过程的温度控制对于防止大坝 混凝土开裂进而保证大坝安全运行具有重要意义。 温度控制的目的是保证基础温差、上下层温差、内外 温差在许可范围内,以保证温度应力在可控范围内, 对大
13、坝施工期及运行期瞬时温度场的预测是温度控 制的关键。 笔者运用边界面法对大坝逐层浇筑施工 期及运行期的瞬态温度进行数值仿真,结合自适应 交叉拟合算法和基本解时间积分与空间积分分离方 法解决瞬态热传导问题中存在的计算时间长、需要 存储空间大等问题,对大坝自然冷却过程中经典瞬 态热传导过程16及真实大坝冷却过程进行仿真,以 验证边界面法计算大坝瞬态热传导问题的有效性。 1摇 瞬态热传导问题的时域卷积法 1. 1摇边界积分方程 大坝浇筑及其温度场演变过程是一个包含水化 热生成的瞬态热传导过程,坝体几何形状、内部温 度、热流密度随时间而变化。 瞬态热传导问题的数 学模型为 k 籽c 2u + Q 籽c
14、 = 觶 u摇 (x 沂 赘) u = 軈u 摇 (x 沂 Su) - k 鄣u 鄣n 以 q = 軃q摇 (x 沂 Sq) q = 茁(u - ua)摇 (x 沂 SR ) (1) 式中:赘 为目标结构所占区域,其边界 祝=Su胰Sq胰 SR;Su、Sq、SR分别为已知温度、热流密度、对流关系 的边界;n 为外法向量;u、q 分别为温度和热流密 度;軈u、軃q 为温度和热流密度的边界已知量; 觶 u 为温度 对时间的导数;Q 为热源密度函数;ua为对流表面 流体的平均温度; k、籽、c 分别为材料的传热系数、密 度、比热容;茁 为对流系数。 瞬态热传导问题的基本解17为 u*(y,x,t,子
15、) = 1 (4仔K(t - 子)1郾 5 e - r2 4K(t-子) q*(y,x,子,t) =- k 鄣 鄣n(x)u *(y,x,子,t ) (2) 其满足 K 2u*+ 觶 u*=- 啄(y - x)啄(t - 子)(3) 式中:u*为位移基本解;K 为导温系数;q*为位移基 本解的法向导数;y、x 分别为源点和场点;t、子 分别 为时间变量和待分析时刻;r 为源点和场点间的距 离函数;啄 为狄利克雷函数。 运用式(2)中的基本解 可得到式(1)的等效积分形式17: 1 籽c乙 子 0乙祝u *(y,x,子,t) k鄣u(x,t) 鄣n(x () ) d祝(x)dt - 1 籽c乙
16、子 0乙祝u(x,t) k 鄣u*(y,x,子,t) 鄣n(x () ) d祝(x)dt + 乙 赘u *(y,x,子,0)u(x,0)d赘(x) + 1 籽c乙 子 0乙赘u *(y,x,子,t)Q(x,t)d赘(x)dt = C(y)u(y,子) (4) 其中C(y) = 0摇 y 埸 赘,y 埸 祝 1摇 y 沂 赘 兹 4仔 摇 y 沂 祝 式中 兹 为 y 点处的立体角。 式(4)中的边界积分方 程为实现瞬态热传导数值仿真的基础。 为简单起 见,暂不考虑热源及初始条件,亦即假设结构内部无 热源并且结构的初始温度为零。 实际上,当热源密 度分布、初始温度分布满足 Laplace 方程时
17、,式(4) 中相应的体积分项可以通过互易原理转化为边界积 分。 边界积分方程可表示为 C(y)u(y,子) = 1 籽c乙 子 0乙祝u(x,t)q *(y,x,子,t)d祝(x)dt - 27 水利水电科技进展,2016,36(1)摇 Tel:02583786335摇 E鄄mail:jz hhu. edu. cn摇 http:/ / www. hehaiqikan. cn 1 籽c乙 子 0乙祝u *(y,x,子,t)q(x,t)d祝(x)dt (5) 1. 2摇 时域卷积法 运用边界面法分析瞬态热传导问题时,时域方 法计算时间积分的方案有拟初始条件方法和时间卷 积法。 拟初始条件法的基本原
18、理是每一个时间步的 计算结果都作为下一步计算的初始条件,因此需要 计算对初始条件的体积分;时间卷积法的基本原理 是每一步的计算都需要用到初始时刻到当前时间步 之间所有时间步的边界温度及热流密度。 采用时间 卷积法,结合边界面法离散技术,将方程(5)最终转 化为线性方程组进行求解。 将时间区间0,子均匀分成 M 个长度为 驻t 的 时间步,在每个时间步内假设温度和热流密度不随 时间改变,于是式(5)通过时间离散为 C(y)u(y,tM) = 1 籽c移 M m = 1乙祝u(x,t m)乙 tm tm-1q *(y,x,t M,t) dtd祝(x) - 1 籽c移 M m = 1乙祝 q(x,t
19、m)乙 tm tm-1u *(y,x,t M,t)dtd祝(x) (6) 其中 tm=m驻t,记 TM+1-m(y,x) =乙 tm tm-1u *(y,x,t M,t)dt (7) FM+1-m(y,x) =乙 tm tm-1q *(y,x,t M,t)dt (8) 则式(6)可化简为 C(y)u(y,tM) = 1 籽c移 M m = 1乙祝 u(x,tm)FM+1-m(y,x)d祝(x) - 1 籽c移 M m = 1乙祝q(x,t m)TM+1-m(y,x)d祝(x) (9) 摇 摇 使用边界面法将上述边界积分方程进行空间离 散,具体过程可以参见稳态热问题的离散过程13, 以下只列出离
20、散之后的边界积分方程: C(y)u(y,tM) = 1 籽c移 M m = 1移 N j = 1 u(xj,tm)移 NP k = 1乙祝kN j(x) FM+1-m(y,x)d祝(x) - 1 籽c移 M m = 1移 N j = 1 q(xj,tm) 移 NP k = 1乙祝kN j(x)TM+1-m(y,x)d祝(x) (10) 方程(10)组织成矩阵形式为 移 M m = 1 HmMum=移 M m = 1 GmMqm(11) 其中摇 摇 HmM,ij= 1 籽c移 NP k = 1乙祝kN j(x)FM+1-m(yi,x)d祝(x) - C(yi)啄mM(12) GmM,ij= 1
21、籽c移 NP k = 1乙祝kN j(x)TM+1-m(yi,x)d祝(x) (13) um,j= u(xj,tm)(14) qm,j= q(xj,tm)(15) 式中:N 为单元形函数;下标 i、j 分别为域边界号和 单元号。 式(11)在考虑边界条件后重新组织可以得到 AMMXM= bM- 移 M-1 m = 1 HmMum-移 M-1 m = 1 GmMq() m (16) 式中 XM为 tM时边界上的未知量。 1. 3摇 基本解的展开 在式(16)右端积分的计算中,对于每一个 M, 前面时间步对应的矩阵 HmM、GmM都需要进行更新, 并且分别与节点温度向量、节点热流密度向量相乘。 当
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 混凝土 瞬态 热传导 分析 快速 边界
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【ho****x】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【ho****x】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。