共查询到19条相似文献，搜索用时 421 毫秒
1.
在实际工程结构动力反应分析中,往往由于结构型式十分复杂,常用的两种直接积分方法,即显式积分方法和隐式积分方法,在使用中都存在着一定的局限性,如何将这两种积分方法合理有效地结合起来,是一个十分有意义的研究课题。针对实际工程问题中整体结构计算时间步长的选择往往受局部区域的材料特性、尺寸大小等因素影响的这一现象,提出了一种对结构局部区域进行隐式积分、对其余区域进行显式积分的显隐式积分方法,这种积分格式相对于显式积分格式而言,能显著提高整体结构的计算速度。最后采用两个数值计算实例对这一方法进行验证。 相似文献
2.
显-隐式组合数值积分算法结合了显式算法无需迭代和隐式算法无条件稳定的各自优点,是结构抗震拟动力试验顺利运行的关键.在对传统显式中央差分法和隐式Newmark β组合算法进行参数修正的基础上,建立了修正CD-Newmark算法,考虑阻尼的影响分析了组合算法的稳定性条件、周期失真率和数值阻尼比,分别得到了试验子结构的稳定性条件和计算子结构无条件稳定的参数合理取值范围,并对计算精度进行了分析.通过算例分析验证了算法的数值特性,从而初步解决了CD-Newmark算法存在稳定性界限过严的问题,为结构抗震拟动力混合试验提供了研究参考. 相似文献
3.
加权近似解析离散化(WNAD) 方法是近年发展的一种在粗网格步长条件下能有效压制数值频散的数值模拟技术. 在地震勘探的实际应用中, 不是所有情况都适合使用空间大网格步长. 为适应波场模拟的实际需要, 本文给出了求解波动方程的非一致网格上的WNAD算法. 这种方法在低速区、介质复杂区域使用细网格, 在其他区域采用粗网格计算. 在网格过渡区域, 根据近似解析离散化方法的特点, 采用了新的插值公式, 使用较少的网格点得到较高的插值精度. 数值算例表明, 非一致网格上的WNAD方法能够有效压制数值频散, 显著减少计算内存需求量和计算时间, 进一步提高了地震波场的数值模拟效率. 相似文献
4.
如何更精确高效地模拟大型、复杂系统内域的波动是发展和完善近场波动数值模拟技术的一个重要研究课题。内域波动的数值模拟通常采用计算量较小的显式方法,但现有的时空解耦显式有限元方法的精度只有二阶;低精度不仅影响数值模拟的精度,而且制约着计算效率的提高。鉴于此,本文旨在探索内域波动数值模拟具有更高精度且稳定的显式方法。作者发展了内域波动数值模拟的现有显式有限元解耦技术,提出了一种高精度且稳定的显式数值模拟方法。 相似文献
5.
《地球物理学进展》2017,(3)
复杂地表条件下的有限差分地震波场的数值模拟,由于受到低速层和地表起伏的限制,模型速度分布范围变大,一般使用精细的差分网格来抑制频散,提高模拟分辨率,但精细网格会显著增加计算成本.为了能有效地解决这一问题,本文提出一种步长自适应有限差分波动方程数值模拟方法.(1)新方法根据模型中的介质速度分布,对不同的速度区域采用与该速度匹配的空间步长,实现对模型空间网格的步长自适应精细划分.对于速度分布范围大的复杂地表模型,新方法不仅能够极大地减少模型的网格节点数,同时又能提高波场的时间采样步长,减少时间采样数,提高计算效率.(2)推导了不同步长边界网格节点Laplace算子的二阶有限差分表达式,避免了在这些结点进行插值计算产生的假扰动和数值不稳定问题.(3)为了降低有限差分产生的数值频散,本文在常规的差分方程中增加了一频散校正项,能有效地衰减了高波数成分,抑制了数值频散.对复杂近地表的波场数值模拟结果表明,本文提出的步长自适应新方法能够有效减少网格节点数和时间采样数,极大地提高计算效率,计算量比常规粗网格增加一些,但效果能够达到了常规精细网格的模拟结果. 相似文献
6.
本文较为系统地给出了线性单自由度振动系统不等时间步长振动反应数值计算的连锁公式,并且讨论了该算法的某些数值分析特性,证明了该算法是无条件稳定的,且满足收敛性要求.文中给出了同时计算相对位移、相对速度和绝对加速度反应谱的计算程序. 相似文献
7.
海洋电磁法是一种有效的钻前储层评价手段,可识别出海底构造是否储油气,从而减少干井率降低勘探成本.近年来,频率域三维海洋电磁反演得到了快速发展,但受到空气波的影响,其在浅海环境中无法取得很好的效果.为解决这一问题,本文研究时间域海洋电磁数据三维正反演.正演模拟和计算中,我们选择基于非结构网格的矢量有限元方法.该算法中空间离散采用非结构四面体网格,可很好地拟合复杂海底地形条件和地下复杂结构;而对于时间离散,我们采用无条件稳定的后推欧拉方法,以确保任意时间步长数值计算的稳定性.反演计算中,灵敏度信息采用伴随正演隐式进行计算;同时,依据时间域反演方法的特点采用L-BFGS方法计算模型修正量.通过利用合成数据反演结果证明了本文提出的三维时间域反演方法可用于复杂海底环境,特别是在浅海环境下的有效性. 相似文献
8.
瑞利阻尼介质有限元离散模型动力分析的数值稳定性 总被引：12，自引：0，他引：12
本文针对几种有一般阻尼的动力系数数值积分的显式方法，讨论了阻尼对稳定性的影响，并建议了瑞利阻尼介质有限元离散模型中动力分析数值稳定性的实用稳定判别方法。 相似文献
9.
交错网格方法（SG）和最优近似解析离散化方法（ONADM）是两类典型的地震波场数值模拟方法.两类方法各有其优势,相对于ONADM方法,SG方法在单个时间层内的计算更为简单;相对于SG方法,ONADM方法可以在较大空间步长条件下有效压制数值频散.结合两种方法的优势,本文提出了一种新的地震波场模拟方法（SG-ONADM）.该方法对控制方程中的一阶偏导数采用SG方法给出的一阶偏导数近似公式,高阶偏导数采用ONADM方法给出的高阶偏导数逼近公式.理论分析及数值算例表明,SG-ONADM方法保留了两种方法的优势,不仅能在较大空间步长条件下有效压制数值频散,同时具有较低的内存需求量;在对同一计算区域进行波场模拟时,SG-ONADM方法的计算效率要高于SG方法和ONADM方法.最后,我们使用SG-ONADM方法进行黏滞声波波场模拟,研究了黏滞声波在复杂介质中的传播. 相似文献
10.
声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.有限差分格式分为显式有限差分和隐式有限差分.隐式有限差分能够进一步压制数值频散效应.因此本文提出了给定频率范围满足时间-空间域隐式有限差分频散关系的方法,并根据震源频率、波速和网格间距确定波数范围,在此基础上建立方程确定了相应的隐式有限差分系数,使得差分系数能在更大频率范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性. 相似文献
11.
Introduction In linear elastic medium,motion equation for lumped-mass finite element simulation of wave motion is expressed as(LIAO,2002)∑=+lililii ttGtM)()()(FUU&&(1)where Mi is lumped-mass at node i,Gil is stiffness coefficient of node i with respect to node l,üi(t)is acceleration vectors at node i,Ul(t)is displacement vectors at node l,Fi(t)is the external nodal force vectors acting at node i.If acceleration vectorsüi(t),displacement vectors Ul(t)and the external nodal force vectors… 相似文献
12.
Viscoelastic artificial boundaries are widely adopted in numerical simulations of wave propagation problems. When explicit time-domain integration algorithms are used, the stability condition of the boundary domain is stricter than that of the internal region due to the influence of the damping and stiffness of an viscoelastic artificial boundary. The lack of a clear and practical stability criterion for this problem, however, affects the reasonable selection of an integral time step when using viscoelastic artificial boundaries. In this study, we investigate the stability conditions of explicit integration algorithms when using three-dimensional (3D) viscoelastic artificial boundaries through an analysis method based on a local subsystem. Several boundary subsystems that can represent localized characteristics of a complete numerical model are established, and their analytical stability conditions are derived from and further compared to one another. The stability of the complete model is controlled by the corner regions, and thus, the global stability criterion for the numerical model with viscoelastic artificial boundaries is obtained. Next, by analyzing the impact of different factors on stability conditions, we recommend a stability coefficient for practically estimating the maximum stable integral time step in the dynamic analysis when using 3D viscoelastic artificial boundaries.
相似文献13.
Implementation and application of the unconditionally stable explicit parametrically dissipative KR‐α method for real‐time hybrid simulation 
下载免费PDF全文
下载免费PDF全文 Chinmoy Kolay James M. Ricles Thomas M. Marullo Akbar Mahvashmohammadi Richard Sause 《地震工程与结构动力学》2015,44(5):735-755
In real‐time hybrid simulations (RTHS) that utilize explicit integration algorithms, the inherent damping in the analytical substructure is generally defined using mass and initial stiffness proportional damping. This type of damping model is known to produce inaccurate results when the structure undergoes significant inelastic deformations. To alleviate the problem, a form of a nonproportional damping model often used in numerical simulations involving implicit integration algorithms can be considered. This type of damping model, however, when used with explicit integration algorithms can require a small time step to achieve the desired accuracy in an RTHS involving a structure with a large number of degrees of freedom. Restrictions on the minimum time step exist in an RTHS that are associated with the computational demand. Integrating the equations of motion for an RTHS with too large of a time step can result in spurious high‐frequency oscillations in the member forces for elements of the structural model that undergo inelastic deformations. The problem is circumvented by introducing the parametrically controllable numerical energy dissipation available in the recently developed unconditionally stable explicit KR‐α method. This paper reviews the formulation of the KR‐α method and presents an efficient implementation for RTHS. Using the method, RTHS of a three‐story 0.6‐scale prototype steel building with nonlinear elastomeric dampers are conducted with a ground motion scaled to the design basis and maximum considered earthquake hazard levels. The results show that controllable numerical energy dissipation can significantly eliminate spurious participation of higher modes and produce exceptional RTHS results. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
14.
Real‐time hybrid testing combines experimental testing and numerical simulation, and provides a viable alternative for the dynamic testing of structural systems. An integration algorithm is used in real‐time hybrid testing to compute the structural response based on feedback restoring forces from experimental and analytical substructures. Explicit integration algorithms are usually preferred over implicit algorithms as they do not require iteration and are therefore computationally efficient. The time step size for explicit integration algorithms, which are typically conditionally stable, can be extremely small in order to avoid numerical stability when the number of degree‐of‐freedom of the structure becomes large. This paper presents the implementation and application of a newly developed unconditionally stable explicit integration algorithm for real‐time hybrid testing. The development of the integration algorithm is briefly reviewed. An extrapolation procedure is introduced in the implementation of the algorithm for real‐time testing to ensure the continuous movement of the servo‐hydraulic actuator. The stability of the implemented integration algorithm is investigated using control theory. Real‐time hybrid test results of single‐degree‐of‐freedom and multi‐degree‐of‐freedom structures with a passive elastomeric damper subjected to earthquake ground motion are presented. The explicit integration algorithm is shown to enable the exceptional real‐time hybrid test results to be achieved. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
15.
Two primary concerns in performing watershed overland flow routing are the numerical instability and computational efficiency. The stability of executing an explicit scheme has to be maintained by observing the Courant–Friedrich–Lewy criterion, which is adopted to confirm that the numerical marching speed is larger than the wave celerity. Moreover, there is another criterion of time step devised in previous studies to avoid back‐and‐forth refluxing between adjacent grids. The situation of refluxing usually occurs on flat regions. In light of this, the selection of a small time increment to honor both restrictions simultaneously is believed to decrease the computational efficiency in performing overland flow routing. This study aims at creating a robust algorithm to relax both restrictions. The proposed algorithm was first implemented on a one‐dimensional overland plane to evaluate the accuracy of the numerical result by comparing it with an analytical solution. Then, the algorithm was further applied to a watershed for 2D runoff simulations. The results show that the proposed integrated algorithm can provide an accurate runoff simulation and achieve satisfactory performance in terms of computational speed. 相似文献
16.
两相介质波动问题显式有限元方法稳定性研究 总被引：2，自引：0，他引：2
基于动力反应递推计算格式传递矩阵的性质,进行了饱和两相介质波动问题时域显式有限元方法稳定性问题的研究。定义了综合考虑各种影响因素的稳定性判别指标——传递因子;研究了稳定性影响因素,包括时间步长、空间步距和渗透系数取值的作用规律;给出了方法稳定性的实用判别准则。研究结果表明:时间步长、空间步距和渗透系数取值都对饱和两相介质波动问题的时域显式有限元方法的稳定性有较为显著的影响。随着时间步长的增大,方法的稳定性降低;随着空间步距的增大,方法的稳定性增加;当渗透系数取值增大时,方法的稳定性增加。 相似文献
17.
湖泊三维风生流隐式差分模型的研究 总被引：1，自引：0，他引：1
建立了湖泊三维风生流的隐格式差分模型,通过对该模型计算稳定性的分析发现,它的时间步长的取值具有较大的灵活性．该模型波用于武汉市墨水湖的民生流模拟,实例计算表明该模型计算结果合理,能较好地反映湖泊风生流的流态特征． 相似文献
18.
19.
多次透射公式的一种稳定实现措施 总被引：3，自引：2，他引：3
基于对双曲偏微分议程数值解稳定性的GKS（Gustafsson,Kreiss and Sundstrom）准则的物理解释，本文提出了在波动模拟中多次透射公式飘移失稳的一个简单措施，并用波源问题和散射问题的数值试验检验了建议的稳定措施的有效性。 相似文献