波动问题的时空高精度模拟研究报告.ppt

波动问题的时空高精度模拟 内容提纲 引言—时空高精度模拟 引言—时空高精度模拟 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 结构弹塑性地震反应分析的微分求积混合方法 一维谱元—微分求积单元格式 一维谱元—微分求积单元格式 一维谱元—微分求积单元格式 一维谱元—微分求积单元格式 一维谱元—微分求积单元格式 一维谱元—微分求积单元格式 一维谱元—微分求积单元格式 一维谱元—微分求积单元格式 一维谱元—微分求积单元格式 一维谱元—微分求积单元格式 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 精度计算量分析 谱元格式的质量矩阵 谱元格式的质量矩阵 谱元格式的质量矩阵 谱元格式的质量矩阵 谱元格式的质量矩阵 谱元格式的质量矩阵 谱元格式的质量矩阵 谱元格式的质量矩阵 谱元格式的质量矩阵 谱元格式的质量矩阵 谱元格式的质量矩阵 数值算例1： 一维空间取100m，在坐标50m处，施加一个初始位移函数，位移函数为三次B样条函数，初始位移函数的施加区域为(50m, 60m),此问题是一个柯西初值波动问题。两端边界取为固定边界。初始速度皆取为0.波速取为1m/s,分别取空间坐标为60m，70m，80m,90m的四个空间点A,B,C,D作为记录点。求取四个点的时程的最大值与解析解的误差，在计算中采取多种参数组合。 位移最大值误差对比（Legendre谱元） 位移最大值误差对比（Chebyshev谱元） 由此可见高精度格式的精度处于一个较高的水平 考虑进一步放宽时间与空间步距：采用参数（5，7；0.2，5），分别代表（ △x， N； △t， n） 不同空间记录点计算时程对比 Legendre谱元 Chebyshev谱元 不同时间记录点空间快照对比 Legendre谱元 Chebyshev谱元 采用有限元-中心差分法进行对比：空间网格取0.5，时间步距取0.04 不同空间记录点计算时程对比 不同时间记录点空间快照对比 数值算例2： 一维空间取100m，在坐标50m处，施加一个初始位移函数，位移函数为三次B样条函数，初始位移函数的施加区域分为三种情况(1).(50m,60m),(2).(50m,70m), (3).(50m,80m)。对应的空间波数为(1),0.1。 (2), 0.05。(3),0.033。初始速度皆取为0.波速取为1m/s。两端边界取为固定边界。取空间坐标80m的空间点A点作为记录点，记取最大位移值并且和解析解进行对比 精度计算量对比—空间波数（1） Legendre谱元 Chebyshev谱元 精度计算量对比—空间波数（2） Legendre谱元 Chebyshev谱元 精度计算量对比—空间波数（3） Legendre谱元 Chebyshev谱元 一些初步结论与展望： 1.时空高精度波动模拟具有高精度的特性。 2.在合理控制精度的情况下，高精度波动模拟方法可以取得比有限元—中心差分方法大的空间与时间步距。 3.在精度控制在同一水平线上时，时空高精度模拟的计算效率优于传统的有限元—中心差分方法。 4.对于地震动荷载，高精度模拟方法受到地震动特性的影响，其不同于光滑的初值扰动问题，需要进行深入的工作与分析。 谱元格式的质量矩阵 Legendre谱元格式 Chebyshev谱元格式 对角阵 空间解耦 符合波速有限概念 满阵 空间耦联 不符合波速有限概念 Chebyshev质量阵 一致质量阵 集中质量阵 空间解耦 符合波速有限基本概念 空间耦连 不符合波速有限基本概念 集中质量阵的计算效果： 一维空间取100m，在坐标50m处，施加一个初始位移函数，位移函数为三次B样条函数，初始位移函数的施加区域为(50m, 60m),此问题是一个柯西初值波动问题。两端边界取为固定边界。初始速度皆取为0.波速取为1m/s,分别取空间坐标为60m，70m，80m的三个空间点A,B,C作为记录点。基本参数取为（2，5；0.1，5），分别代表（ △x， N； △t， n） 不同空间记