基于有限差分的微分变换法在混凝土溶质运移问题中的.doc

混合微分变换法 在混凝土溶质运移问题中的应用收稿日期：2015-06-23 收稿日期：2015-06-23 基金项目：国家“九七三”重点基础研究发展计划(2011CB013800) 第一作者：钟新然(1992-)，男，博士生，主要研究方向为混凝土多场耦合，E-mail: 101839@tongji.edu.cn 通讯作者：戴瑛(1966-)，女，教授，博士生导师，工学博士，主要研究方向为计算力学，E-mail: ydai@tongji.edu.cn 钟新然，戴 瑛，陈熹 （同济大学 航空航天与力学学院，上海 200092） 摘要：微分变换法是一种基于泰勒公式的求解方法，在求解包含间断条件的溶质运移偏微分方程（组）时会失效。可以采用将有限差分法与微分变换法相结合的求解思路，使问题得到解决。利用该方法得到混凝土SO42-的反应性溶质运移方程以及CO2-SO42-和CO2-Cl-的全耦合溶质运移方程组的近似解表达式并将结果与有限元结果进行比较，证明了该方法的有效性。 关键词：混凝土；溶质运移；微分变换法；有限差分法 中图分类号：TU37 Application of Hybrid DTM and FDM in the Solute Transport Problem of Concrete Abstract: Differential transform method (DTM) ceases to be effective when solving the solute transport equation with discontinuous condition. In order to solve the problem, hybrid DTM and FDM is adopted to solve the solute transport equation of SO42-, coupled CO2- SO42- and coupled CO2-Cl- of concrete, and the approximate solutions are obtained. Finally, the results are compared with those obtained by FEM, and it is shown that the proposed method is effective to obtain the approximate solution of solute transport problem. Key words: concrete; solute transport; finite difference method; differential transform method 混凝土是一种多孔介质，当它作为地下结构的材料时，土壤或地下水中的CO2、SO42-、Cl-等化学物质的侵蚀会导致混凝土力学性能劣化[1-4]，从而影响混凝土的承载力。因此，为了保证地下结构的使用安全，有必要对混凝土中的溶质运移过程进行分析。 1997年，Gospodinov P N和Kazandjiev R[4]通过相关实验标定SO42-的扩散系数，采用不稳定差分格式计算时间步长，获得二维硫酸盐溶质运移方程的数值解。2001年，窦红[5]将有限体积和有限元进行结合，求解非线性对流扩散问题并证明其数值解的稳定性。2009年，左晓宝和孙伟[6]对一维硫酸盐侵蚀混凝土的化学-力学(M-C)耦合方程进行Matlab自主编程，得到数值解。2012年，Bonakdar A和Mobasher B等[7]结合闭合形式的误差函数及有限元对氯离子扩散方程进行计算，分别获得考虑反应项为1阶以及2阶的近似解析解。 混凝土的溶质运移过程通常是由偏微分方程（组）描述的[4-7]，而微分变换法(Differential Transform Method，DTM)是近二十年发展起来的一种求解偏微分方程的方法，它由赵家奎于1986年首次提出[8]，此后， Arikoglu A[9]、 ?zdemir ? [10]、 Kangalgil F[11]和Shahed M E[12]等将其用于求解多种偏微分方程。在此，本文亦将其运用于求解混凝土中的溶质运移问题。 微分变换法是一种基于泰勒级数的求解方法，但在处理某些含间断条件的偏微分方程时，会出现级数的待定系数值不确定的问题，如混凝土溶质运移问题中初始条件间断时，级数的待定系数取值不唯一，从而使后续的求解无法进行。为解决该问题，本文利用Chu H P[13]的混合微分变换法(hybrid differential transform and finite difference method, hybrid DTM