共查询到19条相似文献,搜索用时 140 毫秒
1.
本通过对传统动力反应分析方法的总结,阐明了建立隐式和显式方法的一般思路及数学本质,提出了使用系统位移反向向量三阶导数的隐工和实用显式积分方法-3阶显式方法,分析了该显式方法的精度和稳定性,并对建立更高阶隐式和显式方法以及方法的精度和稳定性作了初步讨论。最后,通过算例对本方法、献[1]方法和经典的常平均加速度法(隐式方法、视为精确解)的精度和稳定性进行了比较分析。结果表明,本方法具有明显的优点。 相似文献
2.
有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程,高阶显式有限差分方法通常用来提高求解精度,已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数,发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明:忽略这些小的差分系数能够近似维持有限差分的精度,但是显著减小了计算量。然后,引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性,采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后,给出了针对均匀和非均匀介质模型的弹性波模拟例子,表明了本文方法的优点。 相似文献
3.
Newm ark-更新精细积分法是动力方程求解的隐式的时域逐步积分法,其稳定性条件非常容易满足。与隐式方法相比较,显式积分方法不需要求解耦联的方程组,可以有效地减少内存占用和机时耗费。因此,根据显式积分方法的特点和优点,基于Newm ark-更新精细积分法的基本思想,提出其显式积分格式。对显式积分方法的精度与稳定性进行了初步的分析,指出该显式积分方法具有极好的稳定性,其精度比隐式积分方法的精度稍低。随着时间步长的增加,其精度优于传统的方法。 相似文献
4.
有限差分方法因其操作简单、计算消耗低而成为地震勘探领域中最为常用的数值模拟方法之一, 然而用离散的显式差分算子数值逼近地震波动方程中的连续导数容易导致数值频散, 并且基于正方形网格离散形式的有限差分方法对不同地质模型的适应性较低.针对一阶变密度声波方程的数值模拟, 本文发展了一种适用于矩形网格离散形式的时间高阶空间隐式有限差分格式, 可以有效压制时间和空间频散, 同时灵活的网格剖分增强了其应用的广泛性.基于本文矩形交错网格时间高阶空间隐式有限差分格式的时空域频散关系和变量替换的思想, 首先采用泰勒级数展开方法求解不同方向的非轴上时间差分系数及轴上空间差分系数, 使本文差分格式可以获得任意偶数阶时间和空间精度.为了进一步提高本文差分格式在更大波数区域的空间模拟精度, 我们采用线性优化方法来求取新的轴上空间差分系数用于一阶变密度声波方程的波场迭代求解中.频散、稳定性分析及数值模拟算例表明: 相比于传统十字形空间域隐式有限差分格式, 本文矩形交错网格时间高阶空间隐式有限差分格式在精度、稳定性和效率方面均具有优势. 相似文献
5.
目前隐式方法是动力弹塑性时程分析最常用的分析方法,然而隐式方法在强非线性分析中常常存在迭代不收敛的问题,并且刚度矩阵求解消耗的存储空间随结构自由度增加呈几何级数增长。因此,在求解高层建筑这种大规模问题中,极易遭遇计算瓶颈。显式分析方法直接求解解耦的方程组,不需要迭代。本文对隐式方法和显式方法进行了对比分析,研究了显式分析方法在高层建筑弹塑性地震反应分析中适用性。实例分析表明,从计算精度来讲,隐式方法和显式方法在稳定条件下都能得到较好的精度。从计算效率来讲,对于自由度较少的结构,隐式方法计算效率较高;对于自由度庞大的结构,显式方法计算效率较高。建议在进行自由度多的高层建筑弹塑性地震反应分析时,采用显式分析方法。 相似文献
6.
介绍了一种求解有阻尼体系振动方程的显式积分格式及其逐步求解的过程,并对其计算精度和稳定性进行了分析。该方法不但能同时求得体系的位移、速度和加速度反应,而且所得到的加速度反应的精度能满足工程需要。 相似文献
7.
本文主要介绍了用于求解有阻尼振动方程的四种常用显式积分方法,并针对其稳定性和精度进行了分析对比,讨论了其适用范围及阻尼对稳定性的影响。 相似文献
8.
针对弹性波数值模拟, 现有的时间高阶有限差分方法主要基于显式差分和正方形网格单元建立的, 因此容易引起较强数值频散且缺乏灵活性.基于上述问题, 本文通过结合十字形和菱形模板设计了一种改进的矩形组合差分模板, 在此基础上发展了隐式交错网格差分方法求解弹性波方程.该方案采用改进差分模板和二阶离散格式共同近似时间导数, 采用隐式差分格式求解空间导数.进一步给出了泰勒级数展开和最小二乘优化两种算法求取的高阶差分系数.联合应用高阶差分方案和波场分离技术产生高精度的弹性波场.将本文方法与几种现有差分方法进行了对比, 包括: 显式和隐式时间二阶差分法, 显式时间高阶差分法.数值分析和模型算例表明, 本文方法比其他方法具有更高的精度和灵活性. 相似文献
9.
显-隐式组合数值积分算法结合了显式算法无需迭代和隐式算法无条件稳定的各自优点,是结构抗震拟动力试验顺利运行的关键.在对传统显式中央差分法和隐式Newmark β组合算法进行参数修正的基础上,建立了修正CD-Newmark算法,考虑阻尼的影响分析了组合算法的稳定性条件、周期失真率和数值阻尼比,分别得到了试验子结构的稳定性条件和计算子结构无条件稳定的参数合理取值范围,并对计算精度进行了分析.通过算例分析验证了算法的数值特性,从而初步解决了CD-Newmark算法存在稳定性界限过严的问题,为结构抗震拟动力混合试验提供了研究参考. 相似文献
10.
本文依据波速有限物理概念,从一维波动方程柯西初值问题的解析解出发,推广应用一种在空间域上采用Lagrange多项式内插、时间域内精确积分的显式方法;并基于包含7个节点的3层网格离散模型,通过权衡精度和稳定性的要求,构建出时空离散精度皆为6阶的稳定递推格式,且具有显式和时空解耦特性。最后,通过算例进行与本文同源的低阶格式(2阶显式格式、4阶显式格式)的对比分析,验证该6阶显式格式的精度和稳定性等理论结果,指出高阶公式对提高计算效率的价值。 相似文献
11.
主要研究简体结构动力分析的新方法,先利用QR法建立了简体结构的新动力模型,然后利用5次样条加权残数法建立了求解动力模型的新算法,最后提出了高层与超高层建筑结构动力分析的新方法,文中附有算例。 相似文献
12.
传统的拟静力法和安全系数时程分析法在评价边坡地震稳定性时存在一定的局限性。在提出准确的评价边坡地震稳定性必需因素的基础上,建议对边坡地震稳定性分析方法重新进行分类。根据动力分析得到的边坡在地震作用下的破坏机制和破裂面的性质和位置,提出基于拉-剪破坏的动力时程分析法和强度折减动力分析法。第一种方法将FLAC计算得到破坏时刻的动应力施加到静力情况下边坡上,采用动力分析得到的拉-剪破裂面,结合极限平衡法求解边坡地震安全系数,是一种改进的动力有限元时程分析法;第二种方法考虑了拉-剪破坏的FLAC强度折减动力分析法,是完全动力的方法。最后通过算例分析验证了新方法的可行性,为边坡地震安全系数计算提供了一种新的思路。 相似文献
13.
通过对多自由度体系结构在地震作用下的能量分析,考虑地震动因素和结构自身特性,提出基于能量概念的二阶段设计方法,取得了较有意义的研究成果,对于基于能量的抗震设计方法走向实用化,提供一条思路。 相似文献
14.
本文简要总结,评述了结构动力反应分析中两类输入反演方法(脉动函数法和直接时域法)的优缺点,在此基础上通过应用文献[1]中提出的3阶显式方法,获得了一个关于结构输入反演问题的新解法-显式直接时域法,新方法集合了脉冲函数法和直接时域法的优点,既可按显式求解动力微方程,又适用于非线性系统的动力反演问题,最后,通过算例对该方法进行了数值验证。 相似文献
15.
为实现地震作用下土-结构动力相互作用问题的有限元模拟,需要在人工边界上完成地震动的有效输入,目前工程和科研中常用的地震动输入方法有两种:波动输入方法和振动输入方法。波动输入方法的模拟精度高,但实施上相对复杂且耗时,而振动输入方法处理简单,但模拟精度较低。针对应力型人工边界提出一种在人工边界上实现地震动输入的新方法,该方法通过对土-结构有限元模型中由人工边界节点及相邻节点组成的局部子结构施加自由波场位移时程并进行动力分析,从而直接获得可实现地震波动有效输入的等效地震荷载,然后在土-结构有限元模型的人工边界节点上施加等效输入地震荷载并完成动力计算,由此完成土-结构动力相互作用问题的地震动输入和地震反应计算。与原有波动输入方法相比,新方法避免了原方法需分别计算人工边界上自由场应力和由引入人工边界条件引起的附加力,以及需要根据不同人工边界面的外法线方向确定荷载作用方向等较为复杂的处理过程,具有等效地震荷载计算简便、地震动输入过程更易于实施的特点。采用均匀弹性半空间和成层弹性半空间一维地震反应算例初步验证新方法的正确性和可靠性。 相似文献
16.
土-结构相互作用分析的关键是建立以土-结构界面定义的无限半空间的动力刚度矩阵。本文介绍了一种求解半无限地基动力刚度的新方法,通过两个算例验证了该方法的精度,并给出了一种利用频域刚性基础动力刚度计算基础时域荷载响应的实用方法,该研究为刚性基础设计提供了一种新的,可靠的理论方法。 相似文献
17.
随着地下空间大规模开发利用,地震灾害对其造成的潜在威胁不容忽视。基于《城市轨道交通结构抗震设计规范(GB50909-2014)》和《地下铁道建筑结构抗震设计规范(DG/TJ08-2064-2009)》建议的分析方法,选取惯性力法、反应位移法(国家规范法、上海规范法)、动力时程方法(线弹性方法、等效线性化方法)三类共5种计算方法,以典型两层双柱三跨地铁车站结构为分析对象进行地震反应的对比验算,对上述计算方法的适应性进行评价。分析结果表明,与动力时程方法相比较,惯性力法计算得到的侧墙剪力值偏大,中柱结果较为接近;对于反应位移法,国家规范方法和上海规范方法的计算模型略有不同,但两者计算结果基本相近,其中土体强制位移、集中地基弹簧、土体动剪切模量等参数取值对计算结果影响显著;对于动力时程方法,线弹性方法和等效线性化方法的结果较为接近,且变化趋势相同。 相似文献
18.
本文重点关注对盾构隧道的纵向地震响应分析方法的研究,提出了一种简化的动力有限元分析方法。这种方法采纳了响应位移法的一些基本思想,采用和响应位移法一样的环梁和土弹簧模型,但采用动力有限元法来评估地基地震动,而摒弃了响应位移法中采用的地基振动为谐波形式的假定。与响应位移法相比,简化动力有限元法能够提高地基震动的计算精度,从而提高结构响应的计算精度;而与三维连续模型相比,新方法的计算规模小得多,从而能够大幅提高计算效率。然后,本文分别采用相应位移法和新提出的简化动力有限元法,对武汉长江隧道盾构段隧道进行了纵向地震响应分析,并对这2种方法计算结果的差异进行了分析。 相似文献
19.
This paper introduces two new types of precise integration methods based on Chebyshev polynomial of the first kind for dynamic response analysis of structures, namely the integral formula method (IFM) and the homogenized initial system method (HISM). In both methods, nonlinear variable loadings within time intervals are simulated using Chebyshev polynomials of the first kind before a direct integration is performed. Developed on the basis of the integral formula, the recurrence relationship of the integral computation suggested in this paper is combined with the Crout decomposed method to solve linear algebraic equations. In this way, the IFM based on Chebyshev polynomial of the first kind is constructed. Transforming the non-homogenous initial system to the homogeneous dynamic system, and developing a special scheme without dimensional expansion, the HISM based on Chebyshev polynomial of the first kind is able to avoid the matrix inversion operation. The accuracy of the time integration schemes is examined and compared with other commonly used schemes, and it is shown that a greater accuracy as well as less time consuming can be achieved. Two numerical examples are presented to demonstrate the applicability of these new methods. 相似文献
|