盾构隧道壁后注浆雷达检测数值模拟分析.pdfVIP

(日) (b， 图1隧道栩瑚质量问题的几何模型 Geomelry behmd F堙l nl¨dtImrgroutmgbuItof鼎k“Iunn“ne5egment 在对圈l表示的种不l司问题m行甫边波场∞数值模拟时，需要给几伺体赋r柏应的介 电参数．如介电常数5、磁导率∥、电导半q以及等效磁阻牢盯矿表l别j{l校型中所介 质的舟I乜参数值。值得注意的足，表1t}J管_It混凝土·I·的相对介电常数取值不是大部分资料 所说的625，m是95，返数值是}lI实验标定{i{到：¨lj：制筋埘于良导体，其电宁孳近似 敢无穷大，其相对介电常数取值的变化对模拟结果基本无影响。这种情况F介质的蝴埘碰导 率及等敬磁矾睾设置}日I表1所q、 表1琏遒壁后注浆模型中不同介质的电性参数列表 J 相对舟电常数s 电导率仃(s／m 十n村磁导平』z 相对磁阻事仃。 空’L 涅箍1． }f集体 {阵 俐箭 2时域有限差分模拟 电磁场问题样秩的最终要求楚求得满足实轴、条件的Maxwell方程的精确觯答，欹褂封|：{] 形式的觯柝解并给予正确的物蝌解释是研究^员所向钟的解次问逝的最佳结粜。。…。但作为 工程实际问趣，解析订池挣往光能为力，能够较厂泛发挥作用的唯有数值山往。FDTD被采用 Y∽氏差分网格的空间离散力式．把带时间变量的Maxwell旋度A程直接转化为差分格式．从 『Ⅲ实现在汁算机的数亨空间一}1对I也磁皱帕传播驶与介质作用进行模拟‘’。木文采用改进后迅 卜订耗彳r硪盟PMl．哦收边界条件的FDTI】沾，对隧道村6目模型，”电碰场的传播进行堆本质、最 肛＆“蚺数值模拟。 2 l差分计算式 在时域计算电磁场时．往{丰将二维问题苘化为二维问题，井假设入射波和放射体的形状 0结构均与z尢戈，从』町排导出TM(EJ或TE(M)波的差分解答。但列隧道衬砌结构这+具 体td趔来说，必颁将z打阳定望隧道轴向时，公式(1J才是竹台实“甫这波传播的二阶精 度的中心最商解粹；“r帕讨论以电型波TM(F)为倒； ￡?“(f，，)=一(，，／J￡；(L『)+口c(“，)【H?“2(，+】／2√)一 Ⅳ：“。“一1／王，)+Ⅳ■“(，l，一I／2)一H，”(，，，+】／2)】 (卜b H：㈠(i，+l／!)=，r。!(f_，+l，2)十B【￡：7(，√)一E：佤，+1门 上海市地球物理学会二十周年论文集 (1_c) 厶0+1坨(f+l／2，歹)=五誓_彪(f+1／2，_，)+B·【E；(f+l，歹)一E?(f，／)】 式中，各系数符号定义为： 彳(‘力=【l一盯(f，j『)·△f／2￡(‘歹)】·【1+仃(f，歹)·出／2s(元／)】一1 B=惝·再瓦、) C(f，歹)=占o·【占(L／)+仃(f，／)△f／2】一1 由于在公式(1)中出现了半个时间步及半个空间步，不方便计算机程序实现：但观察 Yee氏网格，对每个场分量来说，其取值的时间步及空间步均为一步，因此，公式(1)可 以迭代形式(2)表示，这样可以非常方便地由计算机实现计算： 厶乙(，，歹)=厶乙(f，／)+曰·【E(f，／一1)一E(f，／)】。 (2一a) ． (2_b) 以(f，歹)=q(f，歹)+曰‘【t(f，．，)一E：(f一1，．，)】 协c，、 ’ 【月I(f，力一月l(f，_，+1)+日，(f+l，歹)一日，(f，．，)】 Eo√卜篙