土动力学第5章地震动数值模拟.ppt

问题提出 大跨结构抗震不适合一致地震动输入 多点地震动输入数值模拟研究 1.地震波在传播过程中进行复杂的反射和散射；同时，地震动场的不同位置地震波的叠加方式也不同，存在相干性，即部分相干效应(incoherence effect)； 2.在地震动场的不同位置，地震波到达的时间上存在差异，即行波效应(traveling wave effect)； 3.在地震动场不同位置，地质条件不同，影响地震波的振幅和频率，即局部场地效应(site effect)； 4.地震波在传播过程中，随着能量的耗散，振幅会减小，即衰减效应(attenuation effect)。 基于开方分解的谱表示法： 相干函数矩阵的开方分解 模型合理性的证明: 自功率谱模型(Clough and Penzien)： 相干函数模型(Harichandran and Vanmarcke )： 各点平稳的地震动时程： 各点非平稳的地震动时程： 第一点与第二点的互功率谱 第二点与第三点的互功率谱 第一点与第三点的互功率谱 原型谱表示法是地震动模拟中比较普遍的一种方法，该方法分为基于功率谱矩阵Cholesky分解的谱表示法和基于功率谱矩阵特征正交分解的谱表示法。 自功率谱模型(Clough and Penzien)： 由于互功率谱矩阵三角分解后下三角矩阵的元素随频率的变化较陡，使得基于互功率谱矩阵Cholesky分解的谱表示法精度较低。而互功率谱矩阵三角分解后的元素随频率的变化较陡，主要是由自功率谱函数随频率变化较陡引起。改进的方法将自功率谱函数分离了出去且可精确计算，而迟滞相干矩阵Cholesky分解后的元素随频率变化较为平缓，因此，只需要对平缓变化的下三角矩阵元素进行插值，即可实现谱表示法的插值模拟。 用谱表示法模拟随机动力场时，模拟的若干组随机序列的功率谱的平均值与目标值是一致的，但是其中任何一组随机序列的功率谱一般与目标值是不同的。因此，有必要分析模拟方法的每一组序列功率谱的离散程度。 偏度误差： 3. 多点地震动的模拟方法---方法2 改进的谱表示方法： 基于功率谱矩阵Cholesky分解 原型谱表示法 基于功率谱矩阵特征正交分解(模态分解) 特征正交分解谱表示法 谱表示法 Hao(1989);屈铁军(1997) Di Paola(2000) 3. 多点地震动的模拟方法---方法2 3. 多点地震动的模拟----方法2 原型谱表示法（Hao,1989）： 改进原型谱表示法： 3. 多点地震动的模拟----方法2 Cholesky分解 3. 多点地震动的模拟----算例2 渐进功率谱 采用Harichandran and Vanmarcke相干函数模型 3. 多点地震动的模拟----算例2 计算量的主要组成 三角函数和的求解(FFT技术) 功率谱矩阵的生产与分解 系数的存储占用内存较大 功率谱矩阵维数增加 插值模拟 当模拟点数较多时 3. 多点地震动的模拟方法---方法2 3. 多点地震动的模拟----方法2 插值模拟前后计算效率比较（酷睿2处理器，主频2.8） ： 4.5 131.1 插值 83.3 3976.3 未插值 500 1.8 45.6 插值 33.6 1446.3 未插值 300 0.1 7.1 插值 1.5 164.7 未插值 100 计算时间(min) 占用内存(Mb) 方法 模拟点数 3. 多点地震动的模拟---方法2 误差函数的定义： 最大自谱误差： 最大互谱误差： 最大自谱误差： 平均互谱误差： 算例： 1000m范围内等间距分布的101个地震动模拟点 v1 v2 v3 v4 v5 v100 v101 1000m 10m 3. 多点地震动的模拟---算例3 引自Der Kiureghian and Neuenhofer (EESD,1992) 截断频率取为120（rad/s） 仍采用Harichandran and Vanmarcke相干函数模型 3. 多点地震动的模拟---算例3 2 4 8 16 32 64 128 1E-3 0.01 0.1 1 10 100 (c) 最大互谱误差 (%) 频率插值基点数 谱表示法 改进方法 误差比较 功率谱比较 0.011 0.020 0.021 0.037 0.057 0.060 Ratio 0.01 0.06 0.25 0.92 4.83 5.77 Scheme2(%) 0.93 2.96 11.9 24.8 84.7 97.1 Scheme1(%) Error4 0.021 0.027 0.027 0.060 0.073 0.088 Ratio 0.02 0.08 0.33 1.54 6.30 8.60 Scheme2