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.

Structural ensemble dynamics based closure model for wall-bounded turbulent flow

Wall-bounded turbulent flow involves the development of multi-scale turbulent eddies, as well as a sharply varying boundary layer. Its theoretical descriptions are yet phenomenological. We present here a new framework called structural ensemble dynamics （SED）, which aims at using systematically all relevant statistical properties of turbulent structures for a quantitative description of ensemble means. A new set of closure equations based on the SED approach for a turbulent channel flow is presented. SED order functions are defined, and numerically determined from data of direct numerical simulations （DNS）. Computational results show that the new closure model reproduces accurately the solution of the original Navier-Stokes simulation, including the mean velocity profile, the kinetic energy of the streamwise velocity component, and every term in the energy budget equation. It is suggested that the SED-based studies of turbulent structure builds a bridge between the studies of physical mechanisms of turbulence and the development of accurate model equations for engineering predictions.

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.

PGA and structural dynamics input motion at a given site

The computation of the representative ground motions,to be used as input for the dynamic analyses of a struc- ture at a particular site,can be approached by several methods.The choice of the approach depends on two factors:the da- ta available and the type of problem to be solved.This paper reports the experience of the authors in approaching a specific case study:the Southern Memnon Colossus,located in Luxor,Egypt.The results are of interest when the hazard analysis estimation in developing countries and the safeguard of cultural heritage are concerned.Monuments have to be treated as important structures,due to their historical and economical value.Hence,standard procedures of probabilistic seismic haz- ard analysis for the seismic classification of common buildings have to be disregarded.On the other hand,the consequences of the collapse of a monument are not comparable to those related to structures such as nuclear power plants and large dams, for which the deterministic seismic hazard analysis provides a straightforward framework for evaluation of the worst case ground motions.An'intermediate'approach,which requites a lower amount of input data with respect to the deterministic one,is adopted.Its stochastic component can eapture significant eharacteristics of earthquakes,primarily the frequency contents which depend on the magnitude(often referred to as the earthquake scaling law).

A time integral formulation and algorithm for structural dynamics with nonlinear stiffness

A newly-developed numerical algorithm, which is called the new Generalized-α （G-α） method, is presented for solving structural dynamics problems with nonlinear stiffness. The traditional G-α method has undesired overshoot properties as for a class of α-method. In the present work, seven independent parameters are introduced into the single-step three-stage algorithmic formulations and the nonlinear internal force at every time interval is approximated by means of the generalized trapezoidal rule, and then the algorithm is implemented based on the finite difference theory. An analysis on the stability, accuracy, energy and overshoot properties of the proposed scheme is performed in the nonlinear regime. The values or the ranges of values of the seven independent parameters are determined in the analysis process. The computational results obtained by the new algorithm show that the displacement accuracy is of order two, and the acceleration can also be improved to a second order accuracy by a suitable choice of parameters. Obviously, the present algorithm is zero- stable, and the energy conservation or energy decay can be realized in the high-frequency range, which can be regarded as stable in an energy sense. The algorithmic overshoot can be completely avoided by using the new algorithm without any constraints with respect to the damping force and initial conditions.

HOMOTOPY SOLUTION OF THE INVERSE GENERALIZED EIGENVALUE PROBLEMS IN STRUCTURAL DYNAMICS

The structural dynamics problems,such as structural design,parameter identification and model correction,are considered as a kind of the inverse generalized eigenvalue problems mathematically.The inverse eigenvalue problems are nonlinear.In general,they could be transformed into nonlinear equations to solve.The structural dynamics inverse problems were treated as quasi multiplicative inverse eigenalue problems which were solved by homotopy method for nonlinear equations.This method had no requirements for initial value essentially because of the homotopy path to solution.Numerical examples were presented to illustrate the homotopy method.

Structural,Electrical,and Lithium Ion Dynamics of Li2MnO3 from Density Functional Theory

The layered Li2MnO3 is investigated by using the first-principles calculations within the GGA and GGA-t-U scheme, respectively. Within the GGA4-U approach, the calculated intercalation voltage （ranges from 4,5 V to 4.9 V） is found to be in good agreement with experiments. From the analysis of electronic structure, the pure phase Li2MnO3 is insulating, which is indicative of poor electronic-conduction properties. However, further studies of lithium ion diffusion in bulk Li2MnO3 show that unlike the two-dimensional diffusion pathways in rock salt structure layered cathode materials, lithium can diffuse in a three-dimensional pathway in Li2MnO3, with moderate lithium migration energy barrier ranges from 0.57 to 0.63 e V.

A strategy for lightweight designing of a railway vehicle car body including composite material and dynamic structural optimization

Rolling stock manufacturers are finding structural solutions to reduce power required by the vehicles,and the lightweight design of the car body represents a possible solution.Optimization processes and innovative materials can be combined in order to achieve this goal.In this framework,we propose the redesign and optimization process of the car body roof for a light rail vehicle,introducing a sandwich structure.Bonded joint was used as a fastening system.The project was carried out on a single car of a modern tram platform.This preliminary numerical work was developed in two main steps:redesign of the car body structure and optimization of the innovated system.Objective of the process was the mass reduction of the whole metallic structure,while the constraint condition was imposed on the first frequency of vibration of the system.The effect of introducing a sandwich panel within the roof assembly was evaluated,focusing on the mechanical and dynamic performances of the whole car body.A mass saving of 63%on the optimized components was achieved,corresponding to a 7.6%if compared to the complete car body shell.In addition,a positive increasing of 17.7%on the first frequency of vibration was observed.Encouraging results have been achieved in terms of weight reduction and mechanical behaviour of the innovated car body.

Solid-state NMR spectroscopy at ultrahigh resolution for structural and dynamical studies of MOFs

To characterize the structure and dynamics of metal--organic frameworks(MOFs)indepth at the molecular level,it is necessary to pursue high-resolution solid-state magic angle spinning(MAS)nuclear magnetic resonance(NMR)spectroscopy.Spectral resolution is usually affected by the quality of materials and various experimental conditions,of which magic angle(MA)accuracy is a crucial determinant.The current industrial criteria for MA calibration based on the common standard of KBr were found insufficient in guaranteeing optimal resolution MAS NMR for highly ordered MOFs.To drive towards higher-resolution MAS NMR spectroscopy,we propose_a calibration protocol for more accurate MA with a higher-precision criterion based on 79Br MAS NMR of KBr,where the linewidth ratio of the fifth-order spinning sideband to the central band of KBr should be less than 1.00.As a result,ultrahigh-resolution 13C cross-polarization(CP)MAS NMR of MOF-5 is achieved with minimal linewidths as low as 4 Hz,and therefore MOF-5 can be used as a new standard convenient for verifying MA accuracy and also optimizing 13c CP conditions.Maintaining high-precision MA under variable temperature(VT)was found challenging on certain commercial MAS NMR probes,as was systematically investigated by VT NMR using KBr and MOF-5.Nevertheless,ultrahigh-resolution MAS NMR spectroscopy with stable MA under VT is employed to reveal fine structures and linker dynamics of a series of Zn-based MOFs with highly regulated structures.The ultrahigh-resolution NMR methodcan be generally applied to study a broad range of MOFs and other materials.

An Educational GUI-Based Software for Dynamic Analysis of Framed Structural Models

This paper introduces a new version of the open-source educational software, LESM (Linear Elements Structure Model), developed in MATLAB for structural analysis of one-dimensional models such as frames, trusses, and grillages. The updated program includes dynamic analysis, which incorporates inertial and damping effects, time-dependent load conditions, and a transient solver with multiple time integration schemes. The software assumes small displacements and linear-elastic material behavior. The paper briefly explains the theoretical basis for these developments and the reorganization of the source code using Object-Oriented Programming (OOP). The updated Graphical User Interface (GUI) allows interactive use of dynamic analysis features and displays new results such as animations, envelope diagrams of internal forces, phase portraits, and the response of degrees-of-freedom in time and frequency domain. The new version was used in a structural dynamics course, and new assignments were elaborated to improve students’ understanding of the subject.

Capturing the non-equilibrium state in light–matter–free-electron interactions through ultrafast transmission electron microscopy

Ultrafast transmission electron microscope(UTEM) with the multimodality of time-resolved diffraction, imaging,and spectroscopy provides a unique platform to reveal the fundamental features associated with the interaction between free electrons and matter. In this review, we summarize the principles, instrumentation, and recent developments of the UTEM and its applications in capturing dynamic processes and non-equilibrium transient states. The combination of the transmission electron microscope with a femtosecond laser via the pump–probe method guarantees the high spatiotemporal resolution, allowing the investigation of the transient process in real, reciprocal and energy spaces. Ultrafast structural dynamics can be studied by diffraction and imaging methods, revealing the coherent acoustic phonon generation and photoinduced phase transition process. In the energy dimension, time-resolved electron energy-loss spectroscopy enables the examination of the intrinsic electronic dynamics of materials, while the photon-induced near-field electron microscopy extends the application of the UTEM to the imaging of optical near fields with high real-space resolution. It is noted that light–free-electron interactions have the ability to shape electron wave packets in both longitudinal and transverse directions, showing the potential application in the generation of attosecond electron pulses and vortex electron beams.

Seismic response of a mid-story isolated structure considering SSI in mountainous areas under long-period earthquakes

At present,there is not much research on mid-story isolated structures in mountainous areas.In this study,a model of a mid-story isolated structure considering soil-structure interaction(SSI)in mountainous areas is established along with a model that does not consider SSI.Eight long-period earthquake waves and two ordinary earthquake waves are selected as inputs for the dynamic time history analysis of the structure.The results show that the seismic response of a mid-story isolated structure considering SSI in mountainous areas can be amplified when compared with a structure that does not consider SSI.The structure response under long-period earthquakes is larger than that of ordinary earthquakes.The structure response under far-field harmonic-like earthquakes is larger than that of near-fault pulse-type earthquakes.The structure response under near-fault pulse-type earthquakes is larger than that of far-field non-harmonic earthquakes.When subjected to long-period earthquakes,the displacement of the isolated bearings exceeded the limit value,which led to instability and overturning of the structure.The structure with dampers in the isolated story could adequately control the nonlinear response of the structure,effectively reduce the displacement of the isolated bearings,and provide a convenient,efficient and economic method not only for new construction but also to retrofit existing structures.

Evaluation on Spot Weld Models in Structural Dynamic Analysis of Automotive Body in White

Spot weld models are widely used in finite element analysis（FEA） of automotive body in white（BIW） to predict static,dynamic,durability and other characteristics of automotive BIW.However,few researches are done on evaluation of the validity of these spot weld models in structural dynamic analysis of BIW.To evaluate the validity and accuracy of spot weld models in structural dynamic analysis of BIW,two object functions,error function and deviation function,are introduced innovatively.Modal analysis of Two-panel and Double-hat structures,which are the dominated structures in BIW,is conducted,and the values of these two object functions are obtained.Based on the values of object functions,the validity of these spot weld models are evaluated.It is found that the area contact method（ACM2） and weld element connection（CWELD） can give more precise prediction in modal analysis of these two classical structures,thus are more applicable to structural dynamic analysis of automotive BIW.Modal analysis of a classical BIW is performed,which further confirms this evaluation.The error function and deviation function proposed in this research can give guidance on the adaptability of spot weld models in structural dynamic analysis of BIW.And this evaluation method can also be adopted in evaluation of other finite element models in static,dynamic and other kinds of analysis for automotive structures.

A SINGULAR VALUE DECOMPOSITION BASED TRUNCATION ALGORITHM IN SOLVING THE STRUCTURAL DAMAGE EQUATIONS

The structural damage identification through modal data often leads to solving a set of linear equations. Special numerical treatment is sometimes required for an accurate and stable solution owing to the ill conditioning of the equations. Based on the singular value decomposition (SVD) of the coefficient matrix, an error based truncation algorithm is proposed in this paper. By rejection of selected small singular values, the influence of noise can be reduced. A simply-supported beam is used as a simulation example to compare the results to other methods. Illustrative numerical examples demonstrate the good efficiency and stability of the algorithm in the nondestructive identification of structural damage through modal data.

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.

Global interpolating meshless shape function based on generalized moving least-square for structural dynamic analysis

A global interpolating meshless shape function based on the generalized moving least-square （GMLS） is formulated by the transformation technique. Both the shape function and its derivatives meet the Kronecker delta function property. With the interpolating GMLS （IGMLS） shape function, an improved element-free Galerkin （EFG） method is proposed for the structural dynamic analysis. Compared with the conven- tional EFG method, the obvious advantage of the proposed method is that the essential boundary conditions including both displacements and derivatives can be imposed by the straightforward way. Meanwhile, it can greatly improve the ill-condition feature of the standard GMLS approximation, and provide good accuracy at low cost. The dynamic analyses of the Euler beam and Kirchhoff plate are performed to demonstrate the feasi- bility and effectiveness of the improved method. The comparison between the numerical results of the conventional method and the improved method shows that the proposed method has better stability, higher accuracy, and less time consumption.

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.

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.

Mass-Stiffness Templates for Cubic Structural Elements Dedicated to Professor Karl Stark Pister for his 95th birthday

1 This paper considers Lagrangian finite elements for structural dynamics constructed with cubic displacement shape functions.The method of templates is used to investigate the construction of accurate mass-stiffness pairs.This method introduces free parameters that can be adjusted to customize elements according to accuracy and rank-sufficiency criteria.One-and two-dimensional Lagrangian cubic elements with only translational degrees of freedom(DOF)carry two additional nodes on each side,herein called side nodes or SN.Although usually placed at the third-points,the SN location may be adjusted within geometric limits.The adjustment effect is studied in detail using symbolic computations for a bar element.The best SN location is taken to be that producing accurate approximation to the lowest natural frequencies of the continuum model.Optimality is investigated through Fourier analysis of the propagation of plane waves over a regular infinite lattice of bar elements.Focus is placed on the acoustic branch of the frequency-vs.-wavenumber dispersion diagram.It is found that dispersion results using the fully integrated consistent mass matrix(CMM)are independent of the SN location whereas its lowfrequency accuracy order is O(κ8),whereκis the dimensionless wave number.For the diagonally lumped mass matrix(DLMM)constructed through the HRZ scheme,two optimal SN locations are identified,both away from third-points and of accuracy order O(κ8).That with the smallest error coefficient corresponds to the Lobatto 4-point integration rule.A special linear combination of CMM and DLMM with nodes at the Lobatto points yields an accuracy of O(κ10)without any increase in the computational effort over CMM.The effect of reduced integration(RI)on both mass and stiffness matrices is also studied.It is shown that singular mass matrices can be constructed with 2-and 3-point RI rules that display the same optimal accuracy of the exactly integrated case,at the cost of introducing spurious modes.The optimal SN location in two-dimensional,bicubic,isoparametric plane stress quadrilateral elements is briefly investigated by numerical experiments.The frequency accuracy of flexural modes is found to be fairly insensitive to that position,whereas for bar-like modes it agrees with the one-dimensional results.

Structural Parameter Analyses on Rotor Airloads with New Type Blade-Tip Based on CFD/CSD Coupling Method

For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics(CSD)module is introduced to handle blade flexibility.In CFD module,dual time-stepping algorithm is employed in temporal discretization,Jameson two-order central difference(JST)scheme is adopted in spatial discretization and B-L turbulent model is used to illustrate the viscous effect.The CSD module is developed based on Hamilton′s variational principles and moderate deflection beam theory.Grid deformation is implemented using algebraic method through coordinate transformations to achieve deflections with high quality and efficiency.A CFD/CSD loose coupling strategy is developed to transfer information between rotor flowfield and blade structure.The CFD and the CSD modules are verified seperately.Then the CFD/CSD loose coupling is adopted in airloads prediction of UH-60A rotor under high speed forward flight condition.The calculated results agree well with test data.Finally,effects of torsional stiffness properties on airloads of rotors with different tip swept angles(from 10° forward to 30° backward)are investigated.The results are evaluated through pressure distribution and airloads variation,and some meaningful conclusions are drawn the moderated shock wave strength and pressure gradient caused by varied tip swept angle and structural properties.