人工边界转换方法解析.docVIP

静-动力分析中人工边界转换方法的研究 摘要：通过将粘弹性动力人工边界应用于同时考虑静力效应和动力效应的工程算例，阐明了此类问题静-动力分析人工边界转换时保证模型为静力平衡体的必要性。通过将粘弹性静-动力统一人工边界应用于半无限空间体有限元模型的静力分析中，验证了静力计算中的误差将使模型动力分析的稳态反应出现相近的误差。在此基础上，系统阐述了适用于同时考虑静力效应和动力效应的工程问题的静-动力分析人工边界转换方法。 关键词：人工边界，静力分析，动力分析，边界转换 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 boundaries 1 前言 人工边界从广义上可分为静力人工边界和动力人工边界。静力人工边界由来已久，通常有固定边界、滚轴边界等。动力人工边界经过几十年的研究发展，动力人工边界发展到现在已有透射边界、粘性边界、粘弹性边界等几种类型。Deeks提出粘弹性人工边界[11]。1998年，刘晶波等人发展了二维的黏弹性人工边界[3]，又于2005年将其发展为三维时域黏弹性人工边界[4]。2006年，刘晶波等人再将二维黏弹性边界发展成一致粘弹性人工边界及粘弹性边界单元普遍做法是采用静力人工边界和动力人工边界分别对静力问题和动力问题进行计算，将计算结果进行叠加后得到完整的结果。但由于叠加原理仅在线弹性小变形范围内适用，原则上不能应用于涉及非线性或大变形问题的分析。涉及非线性或大变形问题 (1) (2) 式中，σV 为竖向应力，γ为围岩重度，σH为横向应力，λ为 (3) 对有限元计算中所取的有限区域，可以据此确定有限域边界条件。 2.2静力分析边界条件 根据上述静力分析方法，在静-动力共同作用问题的计算中，在进行动力分析之前须先确定地下结构模型的静应力场。地下结构开挖前可将大地假设为半无限空间体，其在重力作用下的静力计算，根据对称性，模型中任一处的水平位移。因此，计算该应力场时，有限区域模型两侧可用水平约束即法向约束，底部可用全约束或仅约束竖直方向，顶面即地面应为自由边界。 2.3 动力分析原理 在动力荷载作用下，有限元体系在t+(t时刻的运动平衡方程为： () 式中为体系的总质量矩阵；为体系的总阻尼矩阵；为体系的总刚度矩阵；为体系的节点加速度向量；为体系的节点速度向量；为体系的节点位移向量；为外荷载向量。体系的总阻尼矩阵采用瑞利阻尼[]： () 式中(、(为常数，可按两种不同的振动频率下测得的阻尼比ξ加以确定。计算中常数(、(可由α+β(2i=2(i(i和α+β(2j=2(j(j求得[]。则(和(可表示为： (6) (7) (i 、(i