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.

Development of Non-Dissipative Direct Time Integration Method for Structural Dynamics Application

A direct time integration scheme based on Gauss-Legendre quadrature is proposed to solve problems in linear structural dynamics.The proposed method is a oneparameter non-dissipative scheme.Improved stability,accuracy,and dispersion characteristics are achieved using appropriate values of the parameter.The proposed scheme has second-order accuracy with and without physical damping.Moreover,its stability,accuracy,and dispersion are analyzed.In addition,its performance is demonstrated by the two-dimensional scalar wave problem,the single-degree-of-freedom problem,two degrees-of-freedom spring system,and beam with boundary constraints.The wave propagation problem is solved in the high frequency wave regime to demonstrate the advantage of the proposed scheme.When the proposed scheme is applied to solve the wave problem,more accurate solutions than those of other methods are obtained by using the appropriate value of the parameter.For the single-degree-offreedom system,two degrees-of-freedom system,and the time responses of beam,the proposed scheme can be used effectively owing to its high accuracy and lower computational cost.

Adaptive Reduced Basis Methods Applied to Structural Dynamic Analysis

The reduced basis methods (RBM) have been demonstrated as a promising numerical technique for statics problems and are extended to structural dynamic problems in this paper. Direct step-by-step integration and mode superposition are the most widely used methods in the field of the finite element analysis of structural dynamic response and solid mechanics. Herein these two methods are both transformed into reduced forms according to the proposed reduced basis methods. To generate a reduced surrogate model with small size, a greedy algorithm is suggested to construct sample set and reduced basis space adaptively in a prescribed training parameter space. For mode superposition method, the reduced basis space comprises the truncated eigenvectors from generalized eigenvalue problem associated with selected sample parameters. The reduced generalized eigenvalue problem is obtained by the projection of original generalized eigenvalue problem onto the reduced basis space. In the situation of direct integration, the solutions of the original increment formulation corresponding to the sample set are extracted to construct the reduced basis space. The reduced increment formulation is formed by the same method as mode superposition method. Numerical example is given in Section 5 to validate the efficiency of the presented reduced basis methods for structural dynamic problems.

Towards a Unified Single Analysis Framework Embedded with Multiple Spatial and Time Discretized Methods for Linear Structural Dynamics

We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.

Precise integration method without inverse matrix calculation for structural dynamic equations

The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.

Time Collocation Method for Structural Dynamic Problems

In order to achieve highly accurate and efficient numerical calculations of structural dynamics, time collocation method is presented. For a given time interval, the numerical solution of the method is approximated by a polynomial. The polynomial coefficients are evaluated by solving algebraic equation. Once the polynomial coefficients are evaluated, the numerical solutions at any time in the interval can be easily calculated. New formulae are derived for the polynomial coefficients,which are more practical and succinct than those previously given. Two structural dynamic equations are calculated by the proposed method. The numerical solutions are compared with the traditional fourth-order Runge-Kutta method. The results show that the method proposed is highly accurate and computationally efficient. In addition, an important advantage of the method is the simplicity in software programming.

A lumped mass Chebyshev spectral element method and its application to structural dynamic problems

A diagonal or lumped mass matrix is of great value for time-domain analysis of structural dynamic and wave propagation problems,as the computational efforts can be greatly reduced in the process of mass matrix inversion.In this study,the nodal quadrature method is employed to construct a lumped mass matrix for the Chebyshev spectral element method(CSEM).A Gauss-Lobatto type quadrature,based on Gauss-Lobatto-Chebyshev points with a weighting function of unity,is thus derived.With the aid of this quadrature,the CSEM can take advantage of explicit time-marching schemes and provide an efficient new tool for solving structural dynamic problems.Several types of lumped mass Chebyshev spectral elements are designed,including rod,beam and plate elements.The performance of the developed method is examined via some numerical examples of natural vibration and elastic wave propagation,accompanied by their comparison to that of traditional consistent-mass CSEM or the classical finite element method(FEM).Numerical results indicate that the proposed method displays comparable accuracy as its consistent-mass counterpart,and is more accurate than classical FEM.For the simulation of elastic wave propagation in structures induced by high-frequency loading,this method achieves satisfactory performance in accuracy and efficiency.

Application of strength reduction method to dynamic anti-sliding stability analysis of high gravity dam with complex dam foundation

Considering that there are some limitations in analyzing the anti-sliding seismic stability of dam-foundation systems with the traditional pseudo-static method and response spectrum method, the dynamic strength reduction method was used to study the deep anti-sliding stability of a high gravity dam with a complex dam foundation in response to strong earthquake-induced ground action. Based on static anti-sliding stability analysis of the dam foundation undertaken by decreasing the shear strength parameters of the rock mass in equal proportion, the seismic time history analysis was carried out. The proposed instability criterion for the dynamic strength reduction method was that the peak values of dynamic displacements and plastic strain energy change suddenly with the increase of the strength reduction factor. The elasto-plastic behavior of the dam foundation was idealized using the Drucker-Prager yield criterion based on the associated flow rule assumption. The result of elasto-plastic time history analysis of an overflow dam monolith based on the dynamic strength reduction method was compared with that of the dynamic linear elastic analysis, and the reliability of elasto-plastic time history analysis was confirmed. The results also show that the safety factors of the dam-foundation system in the static and dynamic cases are 3.25 and 3.0, respectively, and that the F2 fault has a significant influence on the anti-sliding stability of the high gravity dam. It is also concluded that the proposed instability criterion for the dynamic strength reduction method is feasible.

SENSITIVITY ANALYSIS BASED ON LANCZOS ALGORITHM IN STRUCTURAL DYNAMICS

The sensitivity calculating formulas in structural dynamics was developed by utilizing the mathematical theorem and new definitions of sensitivities. So the singularity problem of sensitivity with repeated eigenvalues is solved completely. To improve the computational efficiency, the reduction system is obtained based on Lanczos vectors. After incorporating the mathematical theory with the Lanczos algorithm, the approximate sensitivity solution can be obtained. A numerical example is presented to illustrate the performance of the method.

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.

Nonlinear Time History Analysis for the Different Column Orientations under Seismic Wave Synthetic Approach

The significant impact of earthquakes on human lives and the built environment underscores the extensive human and economic losses caused by structural collapses. Over the years, researchers have focused on improving seismic design to mitigate earthquake-induced damages and enhance structural performance. In this study, a specific reinforced concrete (RC) frame structure at Kyungpook National University, designed for educational purposes, is analyzed as a representative case. Utilizing SAP 2000, the research conducts a nonlinear time history analysis to assess the structural performance under seismic conditions. The primary objective is to evaluate the influence of different column section designs, while maintaining identical column section areas, on structural behavior. The study employs two distinct seismic waves from Abeno (ABN) and Takatori (TKT) for the analysis, comparing the structural performance under varying seismic conditions. Key aspects examined include displacement, base shear force, base moment, joint radians, and layer displacement angle. This research is anticipated to serve as a valuable reference for seismic restraint reinforcement work on RC buildings, enriching the methods used for evaluating structures through nonlinear time history analysis based on the synthetic seismic wave approach.

Dynamic Analysis of a Structure Using Taylor Expansion Stochastic Finite Element Method

Due to spatial variability of material property, Young＇ s modulus is assumed to be a stochastic process. This paper proposes a new method of calculating stochastic field, which includes the information of nodes and the midpoint. The method is of high accuracy, and is easy to program. It is an improvement of the midpoint method of stochastic field. In the case of material properties, geometry parameters and applied loads are assumed to be stochastic, the vibration equation of a structure is transformed into a static problem using Newmark method. The Taylor expansion stochastic finite element method（TSFEM） is extended for the structure vibration analysis. An example is given and the calculated results are compared to validate the proposed methods.

An Integrated Structural Strength Analysis Method for Spar Type Floating Wind Turbine

An integrated structural strength analysis method for a Spar type floating wind turbine is proposed in this paper,and technical issues related to turbine structure modeling and stress combination are also addressed.The NREL-5MW ＂Hywind＂ Spar type wind turbine is adopted as study object.Time-domain dynamic coupled simulations are performed by a fully-coupled aero-hydro-servo-elastic tool,FAST,on the purpose of obtaining the dynamic characteristics of the floating wind turbine,and determining parameters for design load cases of finite element calculation.Then design load cases are identified,and finite element analyses are performed for these design load cases.The structural stresses due to wave-induced loads and wind-induced loads are calculated,and then combined to assess the structural strength of the floating wind turbine.The feasibility of the proposed structural strength analysis method for floating wind turbines is then validated.

An Explicit-Implicit Mixed Staggered Asynchronous Step Integration Algorithm in Structural Dynamics

Many engineering applications need to analyse the system dynamics on the macro and micro level,which results in a larger computational effort.An explicit-implicit asynchronous step algorithm is introduced to solve the structural dynamics in multi-scale both the space domain and time domain.The discrete FEA model is partitioned into explicit and implicit parts using the nodal partition method.Multiple boundary node method is adopted to handle the interface coupled problem.In coupled region,the implicit Newmark coupled with an explicit predictor corrector Newmark whose predictive wave propagates into the implicit mesh.During the explicit subcycling process,the variables of boundary nodes are solved directly by dynamics equilibrium equation.The dissipation energy is dynamically determined in accordance with the energy balance checking.A cantilever beam and a building two numerical examples are proposed to verify that the method can greatly reduce the computing time while maintaining a high accuracy.

Unascertained Factor Method of Dynamic Characteristic Analysis for Antenna Structures

The dynamic characteristic analysis model of antenna structures is built,in which the structural physical parameters and geometrical dimensions are all considered as unascertained variables.And a structure dynamic characteristic analysis method based on the unascertained factor method is given.The computational expression of structural characteristic is developed by the mathematics expression of unascertained factor and the principles of unascertained rational numbers arithmetic.An example is given,in which the possible values and confidence degrees of the unascertained structure characteristics are obtained.The calculated results show that the method is feasible and effective.

Probabilistic model and analysis of coupled train-ballasted track-subgrade system with uncertain structural parameters

Random dynamic responses caused by the uncertainty of structural parameters of the coupled train-ballasted track-subgrade system under train loading can pose safety concerns to the train operation.This paper introduced a computational model for analyzing probabilistic dynamic responses of three-dimensional(3D)coupled train-ballasted track-subgrade system(TBTSS),where the coupling effects of uncertain rail irregularities,stiffness and damping properties of ballast and subgrade layers were simultaneously considered.The number theoretical method(NTM)was employed to design discrete points for the multi-dimensional stochastic parameters.The time-histories of stochastic dynamic vibrations of the TBSS with systematically uncertain structural parameters were calculated accurately and efficiently by employing the probability density evolution method(PDEM).The model-predicted results were consistent with those by the Monte Carlo simulation method.A sensitivity study was performed to assess the relative importance of those uncertain structural parameters,based on which a case study was presented to explore the stochastic probability evolution mechanism of such train-ballasted track-subgrade system.

Precise integration methods based on the Chebyshev polynomial of the first kind

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.

Perturbation Analysis of the Attached Water Mass Effect on Offshore Platform Dynamic Response

Theoretically speaking, it is impossible to make the differential equation of motion uncoupled for the natural modes of a system in consideration of the attached water. The hydro-elastic structure is equal to the non-proportional damping system. In this paper a perturbation analysis method is put forward. The structure motion equation is strictly solved mathematically, and the non-proportional damping problem is transformed into a series of proportional damping ones in the superposition form. The paper also presents the calculation formula of the dynamic response of the structure being subjected to harmonic and arbitrary load. The convergence of the proposed method is also studied in this paper, and the corresponding convergence conditions are given. Finally, the proposed method is used to analyze the displacement response of a real offshore platform. The calculation results show that this method has the characteristics of high accuracy and fast convergence.

Analysis of Office-Teaching Comprehensive Buildings Using a Modified Seismic Performance Evaluation Method

Current building design codes allow the appearance of structural and nonstructural damage under design basis earthquakes.The research regarding probabilistic seismic loss estimation in domestic building structure is urgent.The evaluation in this paper is based on a 11-story reinforced concrete office building,incremental dynamic analysis(IDA)is conducted in Perform 3D program using models capable to simulate all possible limit states up to collapse.Next,the probability distribution of post-earthquake casualties,rebuild costs repair costs and business downtime loss are calculated in PACT software for the studied building considering the modified component vulnerability groups and population flow models.The evaluation procedure can also shed light on other types of buildings in China.For non-typical functional building structures,this article proposes to build a finite element model of structural components and to classify the vulnerability groups based on the construction drawings,and to supply and improve the vulnerability library of appendages in FEMA P-58 according to the actual situation.In this way,the application scope of building seismic performance evaluation can be expanded.

Tactile Response Characterization of a Dynamic System Using Craig-Bampton Method

Vibrational characteristics in small horizontal axis wind turbine system are presented in this study with a system concept called tactile response and substructuring.The main focus is on managing the dynamic properties like vibration,noise,and harshness that occur during the operational mode.Tactile response is defined as the response of subsystem which is induced when a human body touches a vibrating system.Sub structuring is a computational method used to reduce the dynamic behavior of a large complex system with a smaller number of degrees of freedom without disturbing the mesh size of the model.Sub structuring has the ability to combine numerical results along with the experimental results.Combination of substructuring and tactile response is applied in this study.Mode shapes of the system are analyzed and results obtained are correlated within this study to provide better optimization of the results.Wind turbine excited with wind energy depends on wind speed.Torsional vibration has a significant role in determining dynamic properties.Torsional vibration is caused as a result of the rotation of the turbine blade and depends on wind speed.The study gives importance to investigating the ability to simulate the numerical method and tactile response to predict and improve dynamic properties.