• Structural Engineering and Mechanics
• /
• v.60 no.5
• /
• pp.891-902
• /
• 2016

The beam structure with breathing cracks subjected to harmonic excitations was modeled by FEM based on Euler-Bernoulli theory, and a piecewise dynamical system was deduced. The precise integration method (PIM) was employed to propose an algorithm for analyzing the dynamic responses of the deduced system. This system was first divided into linear sub-systems, between which there are switching points resulted from the breathing cracks. The inhomogeneous terms due to the external excitations were tackled by introducing auxiliary variables to express the harmonic functions, hence the sub-systems are homogeneous. The PIM was then applied to solve the homogeneous sub-systems one by one. During the procedures, a predictor-corrector algorithm was presented to determine the switching points accurately. The presented method can provide solutions with an accuracy to a magnitude of $10^{-12}$ compared with exact solutions obtained by the theories of ordinary differential equations. The PIM results are much more accurate than Newmark ones with the same time step. Moreover, it is found that the PIM can maintain a high level of accuracy even when the time step increases within a relatively wide range.

• Structural Engineering and Mechanics
• /
• v.54 no.4
• /
• pp.623-648
• /
• 2015

A numerical approach is presented for the analysis of the forced vibration of a rigid surface foundation with arbitrary shape. In the analysis, the foundation is discretized into a number of sub squaree-lements. The dynamic response within each sub-element is described by the Green's function, which is obtained by the Fourier-Bessel transform and Precise Integration Method (PIM). Incorporating the displacement boundary condition and force equilibrium of the foundation, it obtains a system of linear algebraic equation in terms of the contact forces within each sub-element. Solving the equation leads to the desired dynamic impedance functions of the foundation. Numerical results are obtained for foundation not only with simple geometrical configurations, such as rectangular and circular foundation, but also the case of irregularly shaped foundation. Several comparisons between the proposed approach and other methods are made. Very good agreement is reached. Also, parametric studies are carried out on the dynamic response of foundation. Addressed in this study are the effects of Poisson's ratio, material damping and contact condition of soil-foundation interface. Several conclusions are drawn the significance of the factors.

• Structural Engineering and Mechanics
• /
• v.43 no.1
• /
• pp.1-13
• /
• 2012

A fast precise integration method (FPIM) is proposed for solving structural dynamics problems. It is based on the original precise integration method (PIM) that utilizes the sparse nature of the system matrices and especially the physical features found in structural dynamics problems. A physical interpretation of the matrix exponential is given, which leads to an efficient algorithm for both its evaluation and subsequently the solution of large-scale structural dynamics problems. The proposed algorithm is accurate, efficient and requires less computer storage than previous techniques.

• Structural Engineering and Mechanics
• /
• v.57 no.2
• /
• pp.327-355
• /
• 2016

A new analytical derivation of the elastodynamic point load solutions for an isotropic multi-layered half-space is presented by means of the precise integration method (PIM) and the approach of dual vector. The time-harmonic external load is prescribed either on the external boundary or in the interior of the solid medium. Starting with the axisymmetric governing motion equations in a cylindrical coordinate system, a second order ordinary differential matrix equation can be gained by making use of the Hankel integral transform. Employing the technique of dual vector, the second order ordinary differential matrix equation can be simplified into a first-order one. The approach of PIM is implemented to obtain the solutions of the ordinary differential matrix equation in the Hankel integral transform domain. The PIM is a highly accurate algorithm to solve sets of first-order ordinary differential equations and any desired accuracy of the dynamic point load solutions can be achieved. The numerical simulation is based on algebraic matrix operation. As a result, the computational effort is reduced to a great extent and the computation is unconditionally stable. Selected numerical trials are given to validate the accuracy and applicability of the proposed approach. More examples are discussed to portray the dependence of the load-displacement response on the isotropic parameters of the multi-layered media, the depth of external load and the frequency of excitation.