共查询到19条相似文献，搜索用时 156 毫秒
1.
有阻尼体系动力问题的一种显式差分解法 总被引：15，自引：7，他引：15
本文以中心差分方法为基础,结合Newmark常平均加速度法的基本假定,推导出了一种求解有限自由度有阻尼体系动力方程的自起步显式差分格式。此格式的稳定性条件与一般中心差分格式的相当,其计算精度不低于二阶精度。 相似文献
2.
水平成层场地地震反应非线性分析 总被引：12，自引：0，他引：12
本文首先推导了李小军积分格式(中心差分与Newmark平均加速度法相结合)的增量形式,并据此离散动力平衡方程,同时,采用Pyke提出的土动力本构模型以及多次透射人工边界条件,提出了一种水平成层场地地震反应非线性分析的显式有限元方法,并据此编制了计算机程序。数值实验表明,这种方法能较好地模拟土层在强地震作用下的非线性特性。 相似文献
3.
在频率-空间域显式叠前深度偏移中，波场深度延拓是通过显 式差分短算子与波场的空间褶积完成的. 基于对显式差分短算子的设计方法的研究，提出了 一种基于相移算子约束的离散光滑插值的构造一维显式短算子的方法. 通过离散光滑插值法 ，在频率-波数域中，以传播区内的相移算子为约束，在传播区外的算子两端处以零点为约 束，进行离散光滑插值，使得所得算子具有二阶光滑可导性，则其对应的频率-空间域中的 算子就可以取得很短. 该方法设计简单，精度高，能够满足波场深度延拓的需要. 相似文献
4.
求解加速度反应的显式积分格式研究 总被引：1，自引：0，他引：1
目前波动显式有限元分析多以位移或速度作为输入,而加速度记录更直接地保留了地震动的原始信息。为此,本文发展了一种直接以地震加速度作为输入的显式算法,该算法以中心差分和New-m ark-β法结合,并通过平衡方程约简得到。以水塔及框架结构为例,与现有振型分解联合Duham el积分方法、Newm ark-β隐式解法及五种显式算法进行了对比分析。结果表明:该算法与文献[3]的方法,可以在保证计算效率的前提下,得到与振型分解联合Duham el积分方法、Newm ark-β隐式解法相吻合的加速度反应。 相似文献
5.
6.
在深度偏移方法中，把二维隐式方法推广到三维，就会面对一个分块对角矩阵求逆问题. 通常，这种矩阵的求逆将耗费大量计算时间，严重制约了三维隐式方法偏移在实际资料处理中的广泛应用. 在螺旋边界条件下，该矩阵H具有Toeplitz结构的正定厄密矩阵，其快速求逆可由谱法LU分解或直解法快速实现. 本文结合谱法LU分解和直接解法方法的优点，提出了一种混合算法. 文中采用谱分解方法建立起矩阵列元素的谱分解表，并采用直解法的递推公式，可以快速给出矩阵的分解. 通过与谱法分解和直解法在分解精度和分解速度两方面的比较表明，本文方法与谱法相比，在非均匀介质中亥姆霍兹算子矩阵分解时的精度提高10倍；在计算速度方面，混合方法比简化后的直解法快. 因此，该方法的提出，在计算精度许可的条件下，最大限度地减少三维隐式差分偏移中矩阵求逆占用的时间，从而使得该方法能真正用于实际地震资料的处理. 相似文献
7.
显-隐式组合数值积分算法结合了显式算法无需迭代和隐式算法无条件稳定的各自优点,是结构抗震拟动力试验顺利运行的关键.在对传统显式中央差分法和隐式Newmark β组合算法进行参数修正的基础上,建立了修正CD-Newmark算法,考虑阻尼的影响分析了组合算法的稳定性条件、周期失真率和数值阻尼比,分别得到了试验子结构的稳定性条件和计算子结构无条件稳定的参数合理取值范围,并对计算精度进行了分析.通过算例分析验证了算法的数值特性,从而初步解决了CD-Newmark算法存在稳定性界限过严的问题,为结构抗震拟动力混合试验提供了研究参考. 相似文献
8.
结构动力反应分析的三阶显式方法 总被引：14，自引：6，他引：14
本通过对传统动力反应分析方法的总结，阐明了建立隐式和显式方法的一般思路及数学本质，提出了使用系统位移反向向量三阶导数的隐工和实用显式积分方法－3阶显式方法，分析了该显式方法的精度和稳定性，并对建立更高阶隐式和显式方法以及方法的精度和稳定性作了初步讨论。最后，通过算例对本方法、献[1]方法和经典的常平均加速度法（隐式方法、视为精确解）的精度和稳定性进行了比较分析。结果表明，本方法具有明显的优点。 相似文献
9.
提出一种新的资料同化方法——显式变分四维同化方法, 该方法将奇异值分解(SVD)技术用于四维空间的预报集合提取正交基向量, 这些基向量不但能够表现分析变量的空间结构, 也能反映它的时间演变特征. 将分析变量依截断的基向量展开后, 控制变量会显式地出现在代价函数中, 避免了传统的变分四维同化方法所必需的伴随模式的运用, 使同化过程变得简单. 用浅水方程模式和人造资料进行的一系列数值试验对所提方法的有效性作了检验并和传统的变分四维同化方法进行比较. 结果表明, 在观测点很密集, 观测和模式都没有误差的情况下, 它不如传统的变分四维同化方法好. 但是当观测点稀疏时显式方法会好于传统的方法, 它对模式误差及观测误差的敏感性也远远小于传统的方法. 相似文献
10.
关于Newmark-β法机理的一种解释 总被引：1，自引：0，他引：1
一般认为,Newmark-β法属于积分类型的动力数值分析方法,和基于荷载分段插值类型的数值方法不是相同类型的方法.在本文中,研究了这两类方法之间的关系,以最常使用的两种Newmark方法-平均常加速度法和线性加速度法为例,从Newmark基本假定出发推导出这两种方法所具有的荷载分布模式.结果发现:平均常加速度法和线性加... 相似文献
11.
We present a comparison of methods for the analysis of the numerical substructure in a real‐time hybrid test. A multi‐tasking strategy is described, which satisfies the various control and numerical requirements. Within this strategy a variety of explicit and implicit time‐integration algorithms have been evaluated. Fully implicit schemes can be used in fast hybrid testing via a digital sub‐step feedback technique, but it is shown that this approach requires a large amount of computation at each sub‐step, making real‐time execution difficult for all but the simplest models. In cases where the numerical substructure poses no harsh stability condition, it is shown that the Newmark explicit method offers advantages of speed and accuracy. Where the stability limit of an explicit method cannot be met, one of the several alternatives may be used, such as Chang's modified Newmark scheme or the α‐operator splitting method. Appropriate methods of actuator delay compensation are also discussed. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
12.
This study examines the performance of integration methods for hybrid simulation of large and complex structural systems in the context of structural collapse due to seismic excitations. The target application is not necessarily for real-time testing, but rather for models that involve large-scale physical sub-structures and highly nonlinear numerical models. Four case studies are presented and discussed. In the first case study, the accuracy of integration schemes including two widely used methods, namely, modified version of the implicit Newmark with fixed-number of iteration (iterative) and the operator-splitting (non-iterative) is examined through pure numerical simulations. The second case study presents the results of 10 hybrid simulations repeated with the two aforementioned integration methods considering various time steps and fixed-number of iterations for the iterative integration method. The physical sub-structure in these tests consists of a single-degree-of-freedom (SDOF) cantilever column with replaceable steel coupons that provides repeatable highlynonlinear behavior including fracture-type strength and stiffness degradations. In case study three, the implicit Newmark with fixed-number of iterations is applied for hybrid simulations of a 1:2 scale steel moment frame that includes a relatively complex nonlinear numerical substructure. Lastly, a more complex numerical substructure is considered by constructing a nonlinear computational model of a moment frame coupled to a hybrid model of a 1:2 scale steel gravity frame. The last two case studies are conducted on the same porotype structure and the selection of time steps and fixed number of iterations are closely examined in pre-test simulations. The generated unbalance forces is used as an index to track the equilibrium error and predict the accuracy and stability of the simulations. 相似文献
13.
The existing on‐line numerical integration algorithms are derived from the Newmark method, which is based on an approximation of derivatives in the differential equation. The state–space procedure (SSP), based on an interpolation of the discrete excitation signals for piecewise convolution integral, has been confirmed as more reliable than the Newmark method in terms of numerical accuracy and stability. In an attempt to enhance the pseudodynamic test, this study presents an on‐line integration algorithm (referred to as the OS–SSP method) via an integration of the state–space procedure with Nakashima's operator‐splitting concept. Numerical stability and accuracy assessment of the proposed algorithm in addition to the explicit Newmark method and the OS method were investigated via an eigenvalue, frequency‐domain and time‐domain analysis. Of the on‐line integration algorithms investigated, the OS–SSP method is demonstrated as the most accurate method with an acceptable stability (although not unconditionally stable) characteristic. Therefore, the OS–SSP method is the most desirable method for pseudodynamic testing if the numerical stability criterion (Δt/T⩽0.5) is ensured for every vibration mode involved. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
14.
Finite element(FE) is a powerful tool and has been applied by investigators to real-time hybrid simulations(RTHSs). This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time(TET) decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method(CDM), the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ(λ = 4) method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay. 相似文献
15.
Converting the second-order differential equation to a first-order equation by integrating it with respect to time once as the governing equation of motion for a structural system can be very promising in the pseudodynamic testing. This was originally found and developed by Chang. The application of this time-integration technique to the Newmark explicit method is implimented and investigated in this paper. The main advantages of using the integral form of Newmark explicit method instead of the commonly used Newmark explicit method in a pseudodynamic test are: a less-error propagation effect, a better capability in capturing the rapid changes of dynamic loading and in eliminating the adverse linearization errors. All these improvements have been verified by theoretical studies and experimental tests. Consequently, for a same time step this time-integration technique may result in less-error propagation and achieve more accurate test results than applying the original form of Newmark explicit method in a pseudodynamic test due to these significant improvements. Thus, the incorporation of this proposed time-integration technique into the direct integration method for pseudodynamic testings is strongly recommended. © 1998 John Wiley & Sons, Ltd. 相似文献
16.
It has been shown that the operator‐splitting method (OSM) provides explicit and unconditionally stable solutions for quasi‐static pseudo‐dynamic substructure testing. However, the OSM provides only an explicit target displacement but not an explicit target velocity, so that it is essentially an implicit method for real‐time substructure testing (RST) when the velocity‐dependent restoring force is considered. This paper proposes a target velocity formulation based on the forward difference of the predicted displacements so as to render the OSM explicit for RST. The stability and accuracy of the resulting OSM‐RST algorithm are investigated. It is shown that the OSM‐RST is unconditionally stable so long as the non‐linear stiffness and damping are of the softening type (i.e. the tangent stiffness and damping never exceed the initial values). The stability of the OSM‐RST for structures with infinite tangent damping coefficient or stiffness is also proved, and the stability of the method for MDOF structures with a non‐classical damping matrix is demonstrated by an energy criterion. The effects of actuator delay and compensation are analysed based on the bilinear approximation of the actuator step response. Experiments on damped SDOF and MDOF structures verify that the stability of the OSM‐RST is preserved when the experimental substructure generates velocity‐dependent reaction forces, whereas the stability of real‐time substructure tests based on the central difference method is worsened by the damping of the specimen. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
17.
The central difference method (CDM) that is explicit for pseudo‐dynamic testing is also believed to be explicit for real‐time substructure testing (RST). However, to obtain the correct velocity dependent restoring force of the physical substructure being tested, the target velocity is required to be calculated as well as the displacement. The standard CDM provides only explicit target displacement but not explicit target velocity. This paper investigates the required modification of the standard central difference method when applied to RST and analyzes the stability and accuracy of the modified CDM for RST. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
18.
Retracted: Discussion of paper ‘Development of a family of unconditionally stable explicit direct integration algorithms with controllable numerical energy dissipation’ by Chinmoy Kolay and James M. Ricles,Earthquake Engineering and Structural Dynamics 2014; 43:1361–1380 
下载免费PDF全文
下载免费PDF全文 Shuenn‐Yih Chang 《地震工程与结构动力学》2015,44(8):1325-1328
It seems that the explicit KR‐α method (KRM) is promising for the step‐by‐step integration because it simultaneously integrates unconditional stability, explicit formulation, and numerical dissipation together. It was shown that KRM can inherit the numerical dispersion and energy dissipation properties of the generalized‐α method (GM) for a linear elastic system, and it reduces to CR method (CRM) if ρ∞ = 1is adopted, where ρ∞ is the spectral radius of the amplification matrix of KRM as the product of the natural frequency and the step size tends to infinity. However, two unusual properties were found for KRM and CRM, and they might limit their application to solve either linear elastic or nonlinear systems. One is the lack of capability to capture the structural nonlinearity, and the other is that it is unable to realistically reflect the dynamic loading. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
19.
应用地震危险性分析理论和地震动人工合成技术，给出Newmark法中所需的地震动时程，解决了斜坡稳定性分析Newmark法中难以选取合适地震动时程的难题。通过对黄土斜坡实例计算，给出了坡体中地震动峰值加速度与深度的关系；在计算坡体位移时，提出了等效峰值加速度的概念；对比了使用地面地震动时程和使用坡体内等效地震动时程的计算结果。 相似文献