静-动力分析中人工边界转换方法的研究摘要：通过将粘弹性动力人工边界应用于同时考虑静力效应和动力效应的工程算例，阐明了此类问题静-动力分析人工边界转换时保证模型为静力平衡体的必要性。

通过将粘弹性静-动力统一人工边界应用于半无限空间体有限元模型的静力分析中，验证了静力计算中的误差将使模型动力分析的稳态反应出现相近的误差。

在此基础上，系统阐述了适用于同时考虑静力效应和动力效应的工程问题的静-动力分析人工边界转换方法。

关键词：人工边界，静力分析，动力分析，边界转换Abstract:Though the application of dynamic viscous-spring artificial boundary to an engineering case with a consideration of both static and dynamic effect, and the application of the unified viscous-spring boundary for static and dynamic analysis to static analysis of a finite modal of half space, the problems of the applications of viscous-spring artificial boundary to this kind of engineering calculation was pointed out, and its corresponding solving method was proposed. On the base, a systematic switching method of these artificial boundaries was specified.Keywords: artificial boundary, static analysis, dynamic analysis, switching of boundaries1 前言人工边界从广义上可分为静力人工边界和动力人工边界。

静力人工边界由来已久，通常有固定边界、滚轴边界等。

动力人工边界经过几十年的研究发展，已形成具有全局人工边界和局部人工边界的两大类别，并应用于各自适应的工程计算中[1]。

动力人工边界发展到现在已有透射边界、粘性边界、粘弹性边界等几种类型。

1994年，Deeks 提出粘弹性人工边界[11]。

1998年，刘晶波等人发展了二维的黏弹性人工边界[3]，又于2005年将其发展为三维时域黏弹性人工边界[4]。

2006年，刘晶波等人再将二维黏弹性边界发展成一致粘弹性人工边界及其对应的粘弹性边界单元[5]，并于2007年推导了三维一致粘弹性人工边界及等效粘弹性边界单元[6]。

目前对静-动力分析的普遍做法是采用静力人工边界和动力人工边界分别对静力问题和动力问题进行计算，将计算结果进行叠加后得到完整的结果[1]。

但由于叠加原理仅在线弹性小变形范围内适用，原则上不能应用于涉及非线性或大变形问题的分析。

目前对涉及非线性或大变形问题的静-动力分析，常用的人工边界转换方法主要有以下几种：（1）静力分析和动力分析都采用滚轴边界或固定边界；（2）静力分析中采用滚轴边界或固定边界，动力分析采用粘弹性边界、透射边界、粘性边界等人工边界；（3）静力分析和动力分析都采用静-动力统一边界，如粘弹性静-动力统一人工边界。

对第（1）种方法，由于固定边界使波动全部反射，已有许多文献证明其具有放大振动效应的作用，目前已经使用得不多。

刘晶波等人基于黏弹性动力人工边界和半无限空间中静力问题的基本解，建立了对动力问题和静力问题均适用的三维黏弹性静-动力统一人工边界，从而上述第（3）种方法得以解决[1]。

然而，在使用人工边界对地下结构进行动力分析时，还存在一些问题。

如第（2）种方法，由于在静-动力分析的人工边界转换时的方法存在问题，致使产生错误的结果。

在第（3）种方法中，将粘弹性静-动力统一人工边界应用于地下结构的静力分析时，其解与准确值存在误差。

本文将就此两问题进行论证和分析，并阐述合理的地下结构静-动力分析人工边界转换方法。

2 静力和动力有限元分析原理2.1 静力分析原理地铁等地下工程初始应力场的确定须先计算未开挖状态下围岩的自重应力场，进而根据施工步骤，采用释放荷载法，计算出衬砌结构和围岩的静应力场[13]。

许多地下结构的自重应力场模型可以假设为半无限空间体，根据经典围岩压力理论和弹性力学理论，半无限空间体中距地表面任一深度h 处的应力状态可定义为[7][8]：h V γσ= (1) h V H λγλσσ== (2)式中，σV 为竖向应力，γ为围岩重度，σH 为横向应力，λ为侧压力系数。

对浅层围岩，可假设其为各向同性介质，侧压力系数可用泊松比表示如下[7]：μμλ-=1 (3) 对有限元计算中所取的有限区域，可以据此确定有限域边界条件。

2.2静力分析边界条件根据上述静力分析方法，在静-动力共同作用问题的计算中，在进行动力分析之前须先确定地下结构模型的静应力场。

地下结构开挖前可将大地假设为半无限空间体，其在重力作用下的静力计算，根据对称性，模型中任一处的水平位移0=H u 。

因此，计算该应力场时，有限区域模型两侧可用水平约束即法向约束，底部可用全约束或仅约束竖直方向，顶面即地面应为自由边界。

2.3 动力分析原理在动力荷载作用下，有限元体系在t +∆t 时刻的运动平衡方程为：t t t t t t t t F u K u C uM ∆+∆+∆+∆+=++ (4) 式中M 为体系的总质量矩阵；C 为体系的总阻尼矩阵；K 为体系的总刚度矩阵； ut t +∆为体系的节点加速度向量； ut t +∆为体系的节点速度向量；u t t +∆为体系的节点位移向量；t t F ∆+为外荷载向量。

体系的总阻尼矩阵采用瑞利阻尼[9]：C M K =+αβ (5)式中α、β为常数，可按两种不同的振动频率下测得的阻尼比ξ加以确定。

计算中常数α、β可由α+βω2i =2ωi ξi 和α+βω2j =2ωj ξj 求得[9]。

则α和β可表示为：j i j i j i j i i j ωωωωωωωξωξα))(()(2-+-=(6)))(()(2j i j i j j i i ωωωωωξωξβ-+-=(7) ωi 、ξi 分别为振型向量φi 对应的自振圆频率和阻尼比。

根据振型分析结果可求得ωi 和ωj ，阻尼比在计算中取ξi =ξj =0.05。

2.4 粘弹性人工边界粘弹性人工边界从用途上可分为动力人工边界和静-动力统一人工边界；从具体实现方法上可分为弹簧-阻尼器边界单元和一致粘弹性边界单元。

下面先介绍动力人工边界。

具体计算中模型边界材料参数由其相邻的围岩介质材料决定。

则当人工边界采用等效的弹簧和阻尼器物理元件来模拟时，其弹簧系数和阻尼系数的计算分别如下[3][4]：法向边界：)1(2μαα+==R ER G K TTT ， s T c C ρ= (8) 切向边界：)1(2μαα+==R ER G K NNN ， p N c C ρ= (9) 式中K T 、K N 分别为法向与切向的弹簧刚度；R 为波源至人工边界点的距离；c s 和c p 分别为S波和P 波波速；E 和G 分别为介质弹性模量和剪切模量；μ为介质泊松比；ρ为介质质量密度；αT 与αN 分别为切向与法向粘弹性人工边界参数，具体取值情况见表1。

表1 粘弹性动力人工边界中参数α的取值模型类型 方向 α 二维人工边界平面内法向 2.0 平面内切向 1.5 出平面切向 0.5 三维人工边界法向 4.0 切向2.0若人工边界采用一致粘弹性边界单元来模拟，其边界单元等效剪切模量、等效弹性模量和阻尼系数可分别用如下几式计算[5][6]：)1(2~μαα+==R EhR G hG T T (10) ~~~~~~~)1()21)(1()1(2)1()21)(1(νννμανννα--++=--+=R E h R G hE N N (11)])1[()1(2])1[(~Np T s NpTsc c n nE R c c n nGRααμρααρη+-+=+-=(12)其中等效阻尼系数取的是法向和切向的平均值。

式中h 为等效边界单元厚度，~ν为等效泊松比，n 为计算模型维数。

其余符号意义同上。

表2 粘弹性静-动力统一人工边界中参数α的取值模型类型 人工边界位置 方向 α 二维人工边界模型底面法向 α*/3 切向 2/3 模型侧面 法向 α*/4 切向 1/8 三维人工边界模型底面法向 α* 切向 2 模型侧面法向 α* 切向1/2考虑到实施的方便，实际计算中边界材料常采用各向同性材料，这时上述等效剪切模量和等效弹性模量之间存在隐含关系式~E =2（1+~ν）~G 。

考虑到普通有限元材料泊松比应限制在0～0.5范围内，则该等效泊松比可按如下取值[5]。

⎪⎩⎪⎨⎧≥--=其它，，T N T N T N 02/)1/(22/~ααααααν (13) 即计算等效泊松比时应先确定T N αα/大小范围，再确定其计算式。

静-动力统一人工边界是在上述动力人工边界的基础上对人工边界参数α进行调整，其具体取值见表2[1]。

其中参数α*具体计算公式：])(1[])(1)[1(26222*RdR d +++-=μα (10)式中μ为围岩泊松比，d 为位置坐标，R 为荷载作用点到人工边界点的距离。

对于底面人工边界，式中D 取荷载作用点至边界单元的水平距离；对于侧面人工边界，式中D 取荷载作用点至边界单元的垂直距离。

3 静-动力分析中的人工边界转换的几种方法在引言中已经介绍了地下结构静-动力分析中常用的三种人工边界转换方法，下面采用算例对第二种和第三种方法存在的问题进行分析。

3.1算例一某双线铁路隧道，其净空跨度10.2m ，为三心圆断面；衬砌初期支护15cm ，二次衬砌30cm ；隧道埋深45m ，Ⅳ类围岩。

采用ANSYS 大型有限元分析软件对该隧道进行建模和计算，模型区域竖向取100m ，横向取100m 。

隧道及围岩物理力学参数见表4。

静力分析采用的边界条件如图1所示，其中g 为重力加速度。

动力分析采用如图2如示动力模型，其最外一层即为粘弹性动力人工边界，厚度2米。

图中P (t )为某实测的列车振动荷载，其时程曲线见图3。

第（2）方法在静力分析中采用滚轴边界或固定边界（可称之为传统静力边界条件），在动力分析中采用粘弹性动力人工边界，在静力分析和动力分析中的荷载和边界条件见表3。

图1隧道静力分析边界条件 图2 动力分析边界条件表3 计算步骤及条件表表4 围岩物理力学参数表图3 列车荷载时程曲线由于隧道结构是对称的，整个模型有限域以及荷载也是对称的，监测点可以取自其对称轴及其一侧即可；取隧道衬砌的拱顶、拱腰、拱脚以及仰拱底等几个点监测其动力响应情况。