Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.     

On precise time integration method for non-classically damped MDOF systems

In the complex mode superposition method, the equations of motion for non-classically damped multiple-degree-of-freedom （MDOF） discrete systems can be transferred into a combination of some generalized SDOF complex oscillators. Based on the state space theory, a precise recurrence relationship for these complex oscillators is set up; then a delicate general solution of non-classically damped MDOF systems, completely in real value form, is presented in this paper. In the proposed method, no calculation of the matrix exponential function is needed and the algorithm is unconditionally stable. A numerical example is given to demonstrate the validity and efficiency of the proposed method.     

基于离散时间传递矩阵法的变截面梁地震时程分析方法研究

为提高变截面梁地震动力求解的计算效率,提出了基于离散时间传递矩阵法的时程分析方法。首先,从欧拉梁的偏微分振动方程出发,基于逐步时间积分法的线性化方法并结合张量变换原理,建立了变截面梁的动力时程计算方法;其次,考虑地震动激励的非一致输入效应,采用数值迭代求解的方式建立了变截面梁地震动力时程分析的离散时间传递矩阵算法;最后,编制了数值仿真计算程序,并结合具体算例进行了算法的有效性和高效性验证。算例结果表明:在采用相同计算模型的前提下,离散时间传递矩阵法不仅能够在计算精度方面与有限元法保持一致,同时还拥有更高的计算效率。     

频率波数域内层状地基动力柔度的精细求解

精细积分法既可得到在计算机精度意义下的精确解,又能够保持哈密顿体系的辛结构。其是求解一阶线性常微分方程组的精确数值方法,既可以用于时间域的初值问题,又可以应用于空间域的两点边值问题。运用精细积分法求解微层区段矩阵,并对微层区段矩阵合并得到整个层状地基的区段矩阵,最终得到层状地基的动力柔度值。运用数值算例验证了本文方法的计算精度。     

基于修正拟牛顿公式的全波形反演

波形反演是一种利用全波场信息,通过最小化预测波场和实际波场的残差来揭示地下岩性和构造信息的方法.本文首先简述了常规拟牛顿算法的原理,之后利用一种新的拟牛顿公式对Davidon-Fletcher-Powell(DFP)和Broyden-Fletcher-Goldfarb-Shanno(BFGS)算法进行了修正,改进后的BFGS算法在近似Hessian矩阵逆矩阵时,不仅考虑了梯度和模型信息,还加入了目标函数本身的信息,而且对于每次迭代,基本没有增加计算量.数值试验表明,相对常规拟牛顿方法,修正BFGS算法在保证反演精度的同时,明显提高了反演效率.     

重力异常场空间-波数混合域三维数值模拟

重力勘探中复杂条件下的三维正演计算量大存储要求高,使得这种条件下重力勘探高效、精细正反演变得困难.针对这一问题,提出一种空间-波数混合域数值模拟方法,该方法将空间域引力位积分进行水平方向二维傅里叶变换,将三维空间域卷积问题转换为多个不同波数之间相互独立的空间垂向一维积分问题,一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征密度变化,可得出单元积分的解析表达式.该方法计算量和存储需求少,算法高度并行;保留垂向为空间域,优势之一在于可根据实际情况合理调整单元疏密程度,准确模拟任意复杂地形和密度异常体的重力异常,兼顾计算精度与计算效率;优势之二在于用形函数拟合求得积分的解析解,计算精度和效率高;充分利用一维形函数积分的高效和高精度,不同波数之间一维积分高度并行性及快速傅里叶变换的高效性,实现重力异常场三维数值模拟.设计棱柱体模型,通过数值解和解析解对比验证了该方法的正确性、适用性和高效性.针对任意复杂地形条件下的重力场及其张量的模拟问题,提出一种快速算法,对其有效性进行了验证.探究标准FFT法的截断效应对计算精度的影响,对比分析Gauss-FFT法和标准FFT扩边法两种方法的计算精度和效率,总结了二者的选取策略,结果表明选用标准FFT扩边法计算效率更高.实际地形的数值模拟表明本文算法适用于任意复杂地形的高效计算.     

General formulation and solution procedure for harmonic response of rigid foundation on isotropic as well as anisotropic multilayered half-space

A general formulation and solution procedure are proposed for harmonic response of rigid foundation on multilayered half-space. It is suitable for isotropic as well as anisotropic soil medium. The wave motion equation is formulated in frequency wave-number domain in the state space. A hybrid approach is proposed for its solution, where the precise integration algorithm (PIA) is employed to carry out the integration. Very high accuracy can be achieved. The mixed variable form of wave motion equation enables the assembly of layers simple and convenient. The surface Green׳s function is regarded as rigorous, because it is free from approximations and discretization errors. The algorithm is unconditionally stable. The numerical implementation is based on algebraic matrix operation. Numerical examples of vibration of rigid foundation validate the efficiency and accuracy of the proposed approach.     

结构地震反应的频域精细传递矩阵法

传统的传递矩阵法需要对控制微分方程进行求解,获得相应的传递矩阵。公式繁琐、复杂。文中提出将传递矩阵法与精细积分法中的指数矩阵运算技巧结合起来,在频域内对结构进行动力分析。与传统的传递矩阵法相比,无需对微分方程进行求解,只需按照迭代公式进行计算,就可以得到所需要的传递矩阵。这种方法公式简单,理论上可实现任意精度,而且计算效率较高,能够快速、高精度的进行结构的地震反应分析。算例显示了精细传递矩阵法的有效性。     

Axisymmetric dynamic response of the multi-layered transversely isotropic medium

By virtue of the precise integration method (PIM) and the technique of mixed variable formulations, solutions for the dynamic response of the multi-layered transversely isotropic medium subjected to the axisymmetric time-harmonic forces are presented. The planes of cross anisotropy are assumed to be parallel to the horizontal surface of the stratified media. Four kinds of vertically acting axisymmetric loads are prescribed either at the external surface or in the interior of the soil system. Thicknesses and number of the medium strata are not limited. Employing the Hankel integral transform in cylindrical coordinate, the axisymmetric governing equations in terms of displacements of the multi-layered media are uncoupled. Applying mixed variable formulations, more concise first-order ordinary differential matrix equations from the uncoupled motion equations can be obtained. Solutions of the ordinary differential matrix equations in the transformed domain are acquired by utilizing the approach of PIM. Since PIM is highly accurate to solve the sets of first-order ordinary differential equations, any desired accuracy of the solutions can be achieved. All calculations are based on the corresponding algebraic operations and computational efforts can be reduced to a great extent. Comparisons with the existing numerical solutions are made to confirm the accuracy of the present solutions proposed by this procedure. Several examples are illustrated to explore the influences of the type and degree of material anisotropy, the frequency of excitation and loading positions on the dynamic response of the stratified medium.     

A least-squares penalty method algorithm for inverse problems of steady-state aquifer models

Based on the generalized Gauss–Newton method, a new algorithm to minimize the objective function of the penalty method in (Bentley LR. Adv Wat Res 1993;14:137–48) for inverse problems of steady-state aquifer models is proposed. Through detailed analysis of the “built-in” but irregular weighting effects of the coefficient matrix on the residuals on the discrete governing equations, a so-called scaling matrix is introduced to improve the great irregular weighting effects of these residuals adaptively in every Gauss–Newton iteration. Numerical results demonstrate that if the scaling matrix equals the identity matrix (i.e., the irregular weighting effects of the coefficient matrix are not balanced), our algorithm does not perform well, e.g., the computation cost is higher than that of the traditional method, and what is worse is the calculations fail to converge for some initial values of the unknown parameters. This poor situation takes a favourable turn dramatically if the scaling matrix is slightly improved and a simple preconditioning technique is adopted: For naturally chosen simple diagonal forms of the scaling matrix and the preconditioner, the method performs well and gives accurate results with low computational cost just like the traditional methods, and improvements are obtained on: (1) widening the range of the initial values of the unknown parameters within which the minimizing iterations can converge, (2) reducing the computational cost in every Gauss–Newton iteration, (3) improving the irregular weighting effects of the coefficient matrix of the discrete governing equations. Consequently, the example inverse problem in Bentley (loc. cit.) is solved with the same accuracy, less computational effort and without the regularization term containing prior information on the unknown parameters. Moreover, numerical example shows that this method can solve the inverse problem of the quasilinear Boussinesq equation almost as fast as the linear one.In every Gauss–Newton iteration of our algorithm, one needs to solve a linear least-squares system about the corrections of both the parameters and the groundwater heads on all the discrete nodes only once. In comparison, every Gauss–Newton iteration of the traditional method has to solve the discrete governing equations as many times as one plus the number of unknown parameters or head observation wells (Yeh WW-G. Wat Resour Res 1986;22:95–108).All these facts demonstrate the potential of the algorithm to solve inverse problems of more complicated non-linear aquifer models naturally and quickly on the basis of finding suitable forms of the scaling matrix and the preconditioner.     

Numerical technique for dynamic analysis of structures with friction devices

In recent years a number of studies on employing friction elements for the seismic protection of buildings has demonstrated conclusively that such devices can markedly reduce earthquake-induced vibrations. Any numerical estimate of the effectiveness of such isolation systems implies a correct solution of the pertinent nonlinear equations of motion. In direct integration algorithms, the phase transitions between adherence and sliding, or the sliding phase may be accompanied by marked high-frequency oscillation of the relative velocity difference. The paper presents a numerical technique for overcoming these problems, thus leading to increased accuracy of the solutions of equations of motions with Coulomb damping. Since only the damping matrix and the loading vector are involved, the procedure is also computationally efficient. In order to validate the proposed numerical technique, an experimental study of a friction system has been carried out. The dynamic response of a four-storey braced frame with friction devices is presented as an example for the practical application of the proposed numerical technique.     

Preconditioned non‐linear conjugate gradient method for frequency domain full‐waveform seismic inversion

We present preconditioned non‐linear conjugate gradient algorithms as alternatives to the Gauss‐Newton method for frequency domain full‐waveform seismic inversion. We designed two preconditioning operators. For the first preconditioner, we introduce the inverse of an approximate sparse Hessian matrix. The approximate Hessian matrix, which is highly sparse, is constructed by judiciously truncating the Gauss‐Newton Hessian matrix based on examining the auto‐correlation and cross‐correlation of the Jacobian matrix. As the second preconditioner, we employ the approximation of the inverse of the Gauss‐Newton Hessian matrix. This preconditioner is constructed by terminating the iteration process of the conjugate gradient least‐squares method, which is used for inverting the Hessian matrix before it converges. In our preconditioned non‐linear conjugate gradient algorithms, the step‐length along the search direction, which is a crucial factor for the convergence, is carefully chosen to maximize the reduction of the cost function after each iteration. The numerical simulation results show that by including a very limited number of non‐zero elements in the approximate Hessian, the first preconditioned non‐linear conjugate gradient algorithm is able to yield comparable inversion results to the Gauss‐Newton method while maintaining the efficiency of the un‐preconditioned non‐linear conjugate gradient method. The only extra cost is the computation of the inverse of the approximate sparse Hessian matrix, which is less expensive than the computation of a forward simulation of one source at one frequency of operation. The second preconditioned non‐linear conjugate gradient algorithm also significantly saves the computational expense in comparison with the Gauss‐Newton method while maintaining the Gauss‐Newton reconstruction quality. However, this second preconditioned non‐linear conjugate gradient algorithm is more expensive than the first one.     

时间域航空电磁2.5维非线性共轭梯度反演

对于时间域航空电磁法二维和三维反演来说,最大的困难在于有效的算法和大的计算量需求.本文利用非线性共轭梯度法实现了时间域航空电磁法2.5维反演方法,着重解决了迭代反演过程中灵敏度矩阵计算、最佳迭代步长计算、初始模型选取等问题.在正演计算中,我们采用有限元法求解拉式傅氏域中的电磁场偏微分方程,再通过逆拉氏和逆傅氏变换高精度数值算法得到时间域电磁响应.在灵敏度矩阵计算中,采用了基于拉式傅氏双变换的伴随方程法,时间消耗只需计算两次正演,从而节约了大量计算时间.对于最佳步长计算,二次插值向后追踪法能够保证反演迭代的稳定性.设计两个理论模型,检验反演算法的有效性,并讨论了选择不同初始模型对反演结果的影响.模型算例表明:非线性共轭梯度方法应用于时间域航空电磁2.5维反演中稳定可靠,反演结果能够有效地反映地下真实电性结构.当选择的初始模型电阻率值与真实背景电阻率值接近时,能得到较好的反演结果,当初始模型电阻率远大于或远小于真实背景电阻率值时反演效果就会变差.     

A mixed method for transient analysis of soil–structure interaction under SH-motion

A Fourier transform approach is applied to the transient analysis of dynamic soil–structure interaction under SH-motion. The governing equations are formulated in the frequency domain using a Finite Element–Boundary Element (FE–BE) coupling method. After solving the transformed problem, the transient solution is obtained using the discrete inverse Fourier transform with a fast Fourier transform algorithm. Two examples are presented in order to show the numerical performance of the proposed technique.     

An efficient numerical model for multicomponent compressible flow in fractured porous media

An efficient and accurate numerical model for multicomponent compressible single-phase flow in fractured media is presented. The discrete-fracture approach is used to model the fractures where the fracture entities are described explicitly in the computational domain. We use the concept of cross flow equilibrium in the fractures. This will allow large matrix elements in the neighborhood of the fractures and considerable speed up of the algorithm. We use an implicit finite volume (FV) scheme to solve the species mass balance equation in the fractures. This step avoids the use of Courant–Freidricks–Levy (CFL) condition and contributes to significant speed up of the code. The hybrid mixed finite element method (MFE) is used to solve for the velocity in both the matrix and the fractures coupled with the discontinuous Galerkin (DG) method to solve the species transport equations in the matrix. Four numerical examples are presented to demonstrate the robustness and efficiency of the proposed model. We show that the combination of the fracture cross-flow equilibrium and the implicit composition calculation in the fractures increase the computational speed 20–130 times in 2D. In 3D, one may expect even a higher computational efficiency.     

半拉格朗日模式风的半解析解及与差分解的比较试验

质点的轨迹计算是半拉格朗日模式的重要基础,传统的数值计算方法由于采用时间差分代替微分,只能得到质点运动轨迹终点的速度,因此质点的移动轨迹（位移）只能靠风速外推的方法计算,导致了模式计算不稳定等问题.借鉴精细积分法中使用半解析解的思路,利用正压原始方程研究了用运动方程的半解析解构建数值模式的可能性.求解了运动方程的一阶和二阶微分方程组的半解析解,通过时间积分半解析解计算质点运动轨迹.数值试验表明,一阶微分方程组的半解析解比差分解略有优势.二阶微分方程组的半解析解在时间步长增大时优势非常明显,而且在保证计算精度的前提下,节省计算时间,这对提高模式性能有重要作用.     

曲线箱梁桥地震反应的频域精细传递矩阵法

本文将高精度的精细积分法和力学概念清晰的传递矩阵法结合起来,以微分方程和矩阵分析理论为基础,提出了一种新的精细传递矩阵形式,在频域内对曲线箱梁桥地震反应进行分析. 与传统的传递矩阵法相比,无需对微分方程进行求解,只需按照迭代公式进行计算,就可以得到所需要的传递矩阵.这种方法公式简单,理论上可实现任意精度要求,而且计算效率较高.能够快速、高精度地进行曲线梁桥的地震反应分析.算例显示了精细传递矩阵法的有效性..     

An efficient approach for dynamic impedance of surface footing on layered half-space

A semi-analytical model for the evaluation of dynamic impedance of rigid surface footing bonded to multi-layered subsoil is proposed. The technique is based on the dual vector form of wave motion equation and Green's influence function of subdisk for horizontally layered half-space. The multi-layered half-space is divided into a quite large number of mini-layers and the precise integration method (PIM) is introduced for the numerical implementation. The PIM is highly accurate for solving sets of first-order ordinary differential equations with specified two-end boundary conditions. It can produce numerical results of Green's influence functions up to the precision of computer used. The dual vector form of wave motion equation makes the combination of two adjacent mini-layers/layers very easy. As a result, the computational effort for the evaluation of Green's influence function of the multi-layered half-space is reduced to a great extent. In order to satisfy the mixed boundary condition at the surface, the footing–soil interface is discretized into a number of uniformly spaced subdisk-elements. Comparisons illustrating the efficiency and accuracy of the proposed approach are made with a number of solutions available in the literature.     

塔里木地区勘探地震正演模拟研究

以塔里木前陆盆地地质背景为模型,采用三维任意差分精细积分方法并行算法实现了库车地区三维正演模拟.三维任意差分精细积分方法通过时间域采用局部积分半解析方法求得波动方程的递推算子,与常规的差分法相比计算精度有较大提高;文中计算稳定性根据实际算例进行分析,采用稳定因子约束,得到较好的稳定性;边界条件采用改进的自适应吸收边界,并通过串行程序并行化,大大减少了三维正演模拟的耗时,完成了大数据量三维正演模拟.     

Flexibility-based linear dynamic analysis of complex structures with curved-3D members

A flexibility-based formulation of a new mass matrix for the dynamic analysis of spatial frames consisting of curved elements with variable cross-sections is presented. The main characteristic of such formulations is the exact equilibrium of forces at any interior point, with no additional hypotheses about the distribution of displacements, strains or stresses. Accordingly, the derived element mass matrix takes into account the exact stiffness and mass distribution throughout each element. In validation tests, results obtained with this method are compared with those obtained by other numerical or analytical formulations, showing the accuracy of the proposed method. The comparison of experimental results for a multispan arch bridge subjected to a dynamic load with those achieved by means of the proposed method are finally included to illustrate its efficiency in the treatment of complex structures. © 1998 John Wiley & Sons, Ltd.