动力有限元法清华大学于玉贞.pptVIP

* 缺点：仅适用于较简单边界条件； 没有考虑结构的影响； 仅适用于竖直传播的SH波。 二、特殊处理 3．边界条件 施加一维自由场 优点：直观、简单 基岩 地基 地震波 结构物 边界A 边界B 一维自由场A 一维自由场B 方法： （1）在距结构物较远处，假定结构物对地基运动影响较小，则可认为地基是自由场； （2）利用剪切层法或集中质量法等分析两侧边界的运动，将其时程运动作为边界条件施加于两侧。 2-4 土体动力分析有限单元法 * 对于P波、S波，无量纲参数a、b均取1－标准粘性边界 二、特殊处理 3．边界条件 粘性边界 i j ui, fi vi, gi uj, fj vj, gj θ l x y r,vS,vP ρ: 土的密度 vS: 剪切波速 vP: 压缩波速 2-4 土体动力分析有限单元法 * 优点：概念明确，处理方法简单； 与截断边界相比，可减小计算域的尺寸； 与吸收层相比，可定量设定阻尼的大小 二、特殊处理 3．边界条件 粘性边界 i j ui, fi vi, gi uj, fj vj, gj θ l x y r,vS,vP 缺点：边界距离仍较大，波长的3/4或1/2； 入射波与边界近乎平行时，精度较低，网格也不宜过小 2-4 土体动力分析有限单元法 * 二、特殊处理 3．边界条件 在边界上施加位移约束 位移约束等于入射波位移 相当于地震波可自由通过边界：只有入射波，没有反射波 位移约束用已知的位移表示 2-4 土体动力分析有限单元法 地基 地震波 建筑物 透射边界 * 透射边界 3．边界条件 计算区域 c cΔt ε θ caΔt cΔt/cosθ Q b x y a 人 工 边 界 入射波沿x轴的视传播 假定波动沿y轴方向不变化 f 为某函数 Dt为时间离散步长 如果点b在t+Dt时刻的位移用点Q在t时刻的位移代替，则入射波将穿过人工边界而不发生任何反射。 导出与入射角θ无关的透射公式 在包含不同波速的波动问题中导出统一的透射公式 采用人工波速ca 二、特殊处理 2-4 土体动力分析有限单元法 * 3．边界条件 计算区域 c cΔt ε θ caΔt cΔt/cosθ Q b x y a 人 工 边 界 误差为： Δ1u(·)是与原入射波传播方向和传播速度相同的波，称为一阶误差波 —透射边界 二、特殊处理 2-4 土体动力分析有限单元法 * 3．边界条件 由此可导出二阶误差波，重复该处理m次，可推广为 物理意义：任意阶误差波Δmu(·)在任一点x和任一时刻t的值可以通过将x和t分别倒退caΔt和Δt来近似。 近似误差为： 计算区域 c cΔt ε θ caΔt cΔt/cosθ Q b x y a 人 工 边 界 —透射边界 二、特殊处理 2-4 土体动力分析有限单元法 * 3．边界条件 计算区域 c cΔt ε θ caΔt cΔt/cosθ Q b x y a 人 工 边 界 在边界点x=xb处，可得N 阶透视公式 高阶误差波 暂态波统一透射公式 —透射边界 二、特殊处理 2-4 土体动力分析有限单元法 * 3．边界条件 计算区域 c cΔt ε θ caΔt cΔt/cosθ Q b x y a 人 工 边 界 暂态波统一透射公式的基本特征： （1）它对传播过程的模拟与入射角无关，任一人工边界点b(xb)在t+Δt时刻的位移可以用内部点(xb- ncaΔt)和t以前若干时刻的位移来计算 （2）对于给定的数值计算精度，人工波速ca允许在一定的范围内变动，一般情况下，ca可以调整为与多种物理波速中的一种相近或相等， ca≤2cmin，但不宜过小 （3）它对波动传播过程的模拟是直接以离散形式完成的， 误差量级：O[(Δx /λ)N] －其精度可与有限元的精度相匹配 —透射边界 二、特殊处理 2-4 土体动力分析有限单元法 * 3．边界条件 透射公式的的有限元格式 计算点和离散结点往往不重合 用插值法改写成以离散结点运动表示的形式 N：透射阶数，一般2－3阶可满足工程精度要求 Δx 人 工 边 界 J J-1 J-2 J-2N 人工边界点J在t=(p+1)Δt时刻的位移 Un：结点J及与结点J相邻的若干结点在(p-n+1)Δt时刻的位移 Ti(i=1,2,……,2n+1)：人工波速ca、Δt 和Δx的函数 —透射边界 （沈聚敏等，抗震工程学，p234-235） 二、特殊处理 2-4 土体动力分析有限单元法 * * * * * 土动力学与 土工抗震工程 于玉贞 清华大学水利水电工程系 岩土工程研究所 * 基岩 地基 地震波 土坝 结构 特性 地基 特性 地震 特性 地基及土工结构