LS Dyna - Response Spectra
LS Dyna - Response Spectra
LS Dyna - Response Spectra
Abstract
A set of new features for frequency domain dynamic and acoustic computation have been implemented in LS
DYNA . They include random vibration and high-cycle fatigue analysis, frequency response functions, steady state
dynamics, response spectrum analysis, and acoustic analysis based on boundary element methods or finite element
methods. The objective of introducing these new features is to add capabilities to LS-DYNA to solve frequency
domain vibration and noise radiation problems. This class of problems is very common in auto and aerospace
industries and many other industries.
The paper provides a brief introduction of the new features. Keywords for the features are introduced. Areas of
applications are discussed. Some examples are given to illustrate how to use these features.
1. Introduction
A set of new features for frequency domain dynamic and acoustic computation have been
implemented in LS-DYNA. The features include random vibration and fatigue analysis,
frequency response function (FRF), steady state dynamics (SSD), response spectrum analysis
and BEM/FEM acoustics (see Table 1). They are important numerical tools for dynamic
simulation and safety and reliability analysis of structures and parts. They have been extensively
used in various fields, such as NVH (Noise, Vibration and Harshness) analysis in auto industry,
aerodynamics in aircraft and aerospace industries, fatigue and lifespan prediction in machinery
manufacturing, seismic analysis, and acoustic design for theatres and auditoriums.
Frequency
domain
Features
Random vibration
Applications
Structural vibration
Random fatigue
FRF
SSD
Response spectrum
Acoustics (BEM/FEM)
Page 1
problems, such as impact and crashworthiness, in some cases, frequency domain analysis is more
advantageous than the traditional time domain analysis. For example, for structures subjected to
periodic or harmonic loading, it may be more efficient to approach the problem in frequency
domain as there might be only a few frequencies involved. For some problems such as wheels
running on rough road, ocean wave loads on offshore platforms, the nature of the loads is
random. The loads are usually described in non-deterministic sense such as frequency spectrum.
One has to resort to random analysis in frequency domain, which solves the problem in a
statistical manner.
Frequency domain analysis is also suitable for problems with extremely long load history, such
as wind turbine exposed to long-term wind load and parts operated in a high-cycle fatigue
environment, where it is too expensive to run a time domain analysis. Meanwhile the frequency
domain analysis can reveal the dependency of the structural response to the excitation
frequencies so that the resonant frequencies can be avoided in the early design phase for new
structures or systems.
With the addition of the new frequency domain features, LS-DYNA, which has been renowned
for its advanced capabilities in simulating complex transient real world problems, becomes more
versatile and can serve as a competitive candidate for tackling even more challenging
engineering problems in wider areas.
The frequency domain dynamic features (random vibration, random fatigue, FRF, SSD and
response spectrum) are based on the results of modal analysis of the structures, e.g. the natural
frequencies and modal shapes. They use either mode superposition method, or mode acceleration
method or some other modal combination methods. Thus one needs to run implicit eigenvalue
analysis preceding the frequency domain dynamic analysis. As usual, double precision
executables are suggested to run the implicit analysis. The unknown variables in the equation
system are modal coefficients instead of the physical variables. As the number of modes
involved in the solution is usually much less than the number of physical degree of freedom, the
frequency domain approach is more efficient than the direct time domain approach.
If stress results are requested in the analysis, modal stress output is required when running the
implicit eigenvalue analysis. Linear element formulation is suggested to retrieve the modal stress
output.
The frequency domain analysis results are given as binary plot files and ASCII database files,
which are accessible by the post-processing software LS-PREPOST.
The new keywords for the frequency domain dynamic and acoustic computation include
*FREQUENCY_DOMAIN_RANDOM_VIBRATION_{FATIGUE}
*FREQUENCY_DOMAIN_FRF
*FREQUENCY_DOMAIN_SSD
*FREQUENCY_DOMAIN_RESPONSE_SPECTRUM
*FREQUENCY_DOMAIN_ACOUSTIC_BEM_{PANEL_CONTRIBUTION}
*FREQUENCY_DOMAIN_ACOUSTIC_FEM
They can be found in the Section *FREQUENCY_DOMAIN in the new release of LS-DYNA
Keyword Users Manual [1].
Page 2
The random vibration capability of LS-DYNA originated from Boeings in-house vibroacoustic code N-FEARA. This feature computes the dynamic response of structures exposed to
vibration or structural-acoustic coupling based on a known source. Various excitations and
acoustic environments can be considered, including base excitations, correlated or non-correlated
acoustic waves such as plane wave, progressive wave, reverberant wave, turbulent boundary
layer, etc. The coupling between the structures (represented by the modal shapes) and the
acoustic excitations is expressed through the concept of joint acceptance. The input spectrum can
be given in terms of PSD for mechanical or acoustic pressure waves, or SPL for acoustic
pressure waves. Here, PSD stands for Power Spectral Density and it displays the mean square
amplitude of each sinusoidal wave with respect to its frequency. It is a function to describe how
the energy or power of a signal is distributed across frequency. SPL stands for Sound Pressure
Level and is the sound pressure measured in a decibel scale. The results include PSD of the nodal
displacements, velocities, accelerations and element stresses, which are printed in the binary file
d3psd, and RMS (Root Mean Square) of those variables, which are printed in the binary files
d3rms. Both d3psd and d3rms are accessible by LS-PREPOST.
The solution is based on modal superposition method or modal acceleration method. If a preload
due to static pressure or temperature change is presumed, intermittent modal analysis is
performed to incorporate the consideration of these preload effects on structural response. A
restart option is provided in case the user has run modal analysis elsewhere and has the d3eigv
family files in the working directory. In this case, the modal analysis is skipped and LS-DYNA
reads the d3eigv files directly to get the modal information. For the case of base acceleration, e.g.
in the shaker table experiment, user can choose if absolute or relative results are requested. A
unit conversion option is also provided as the acceleration load for many industry problems are
defined in terms of G (gravity acceleration), instead of the compatible units ([length]/[time]2)
used by LS-DYNA by default.
This feature has been extended to random vibration or high-cycle fatigue analysis by LSTC,
based on knowledge of materials S-N fatigue curve and Miners cumulative fatigue damage
rule. For random fatigue analysis, a binary plot file d3ftg, which is accessible by LS-PREPOST,
is created to show fringe level of cumulative damage ratio, or expected remaining life-cycle for
the structure. User can choose parts or set of elements for running fatigue analysis, e.g. user may
run fatigue analysis for only the dangerous elements which have high risk of failure.
The keyword to activate the random vibration (and fatigue) analysis is
* FREQUENCY_DOMAIN_RANDOM_VIBRATION_{FATIGUE}.
A sample problem can be found in Figure 1. The structure is composed of 1972 nodes and 1865
shell elements. The material property of 1045 steel is adopted. The structure is subjected to base
acceleration 2G2/Hz for the range of frequency 100-2000 Hz. Constant modal damping ratio 0.03
is adopted. The first 15 natural modes are employed in the modal superposition. Figure 1 (a)
provides contour plot of the RMS of the Von-Mises stress; Figure 1 (b) provides the contour plot
of the cumulative damage ratio for the structure with the exposure time of 4 hours.
Page 3
Page 4
Page 5
FRF is a characteristic of a structure that has a measured or computed response resulting from a
known applied harmonic input. In other words, it is a transfer function between the excitation
and response of a system. The basic formula for a FRF is
X( f )
(1)
F( f )
Where F ( f ) is the input to the system in frequency domain and X ( f ) is the output of the system
in frequency domain.
H( f ) =
With the implementation of the FRF feature, LS-DYNA provides user the opportunity to
acquire a spectrum of structural response (displacement, velocity and acceleration) for applied
unit harmonic excitations. The harmonic excitations are given in the form of nodal force,
pressure or base acceleration, for a range of frequencies. The response can be given in terms of
displacement, velocity, or acceleration. FRFs are complex functions, with real and imaginary
components. They can also be expressed in the form of magnitude and phase angle pairs.
The FRF feature is activated by the keyword *FREQUENCY_DOMAIN_FRF. Through the
keyword, user provides information about the location, direction, range of frequencies for the
harmonic excitation, and the location where the response is requested. Damping information is
also provided in the keyword input. Damping can be given as constant or mode dependent
damping ratio, or Rayleigh damping. The location of the excitation and response area can be
given as node, or set of nodes, or set of segments. The direction of load can be in any of the x, y,
z directions or given as a vector by using *DEFINE_VECTOR. The results are given as ASCII
files frf_amplitude and frf_angle, which can be plotted as xy-curves with LS-PREPOST.
For more information about the FRF and its implementation in LS-DYNA, please refer to the
papers [4], [5] and [6].
Page 6
(a) Amplitude
Page 7
of Sum of Squares (SRSS) method, the United States Nuclear Regulatory Commission (NRC)
grouping method, and the Complete Quadratic Combination (CQC) method. The results are
written to a binary plot file d3spcm, which is accessible by LS-PREPOST. Both nodal and
elemental results are available. The input spectrum can be defined as a group of curves for
displacement, velocity, acceleration, nodal force, or pressure spectrums, for different damping
ratios.
This feature can be activated by the keyword
*FREQUENCY_DOMAIN_RESPONSE_SPECTRUM.
This feature can find applications in earthquake engineering and safety evaluation of buildings
and constructions.
An example of response spectrum analysis of dam-foundation system is given below. The damfoundation system is given in Figure 7. The foundation is assumed to be rigid for simplicity. The
arch dam is 464.88 ft high. Elastic material is assumed for the concrete dam.
Page 8
Page 9
Figure 12. Sound Pressure Level at a field point above the plate
One can notice that the frequencies corresponding to the peak SPL (17 Hz, 105 Hz and 295 Hz)
are also the frequencies corresponding to the peak acceleration response given by SSD (see
Figure 6 (a)).
Page 10
The BEM acoustic feature has been extended to provide the capability to study the acoustic panel
contribution. This helps the designers to understand the noise contribution from each of the
panels of a structure. Then remedies can be adopted to optimize the acoustic performance, to
meet the requirements from clients, or from environmental regulations.
To study the acoustic panel contribution, the following keyword should be used
*FREQUENCY_DOMAIN_ACOUSTIC_BEM_PANEL_CONTRIBUTION.
For more information about the BEM acoustic feature, please refer to the conference papers [9]
and [10].
8. Conclusion
This paper introduces briefly a set of new features of LS-DYNA for solving frequency domain
vibration and acoustic problems. They include random vibration and fatigue analysis, frequency
response function, steady state dynamics, response spectrum analysis and BEM/FEM acoustics.
The keyword commands and application areas of the features are introduced. The input
information and output files are discussed. Some examples are given to illustrate how to use
these features. References are listed to provide more information about the features.
The objective of implementing these features is to provide users the capabilities to deal with
frequency domain vibration and acoustic problems, which are very common in auto and
aerospace industries and many other industries. The implemented features as well as upcoming
developments establish LS-DYNA as an attractive tool for this class of problems.
Acknowledgements
The authors are grateful to Dr. George Laird from Predictive Engineering, Inc. for providing the
satellite model for the random vibration analysis example.
Page 11
References
1.
LS-DYNA Keyword User's Manual, Version 971, Livermore Software Technology Corporation,
Livermore, California, 2010.
2.
Rassaian M., Huang Y., Lee J., Structural analysis with vibro-acoustic loads in LS-DYNA. Proceedings of
the 10th International LS-DYNA Users Conference, June 8-10, 2008, Dearborn, Michigan.
3.
Shor O., Lev Y., Huang Y., Simulation of a thin walled aluminum tube subjected to base acceleration using
LS-DYNA's vibro-acoustic solver. Proceedings of the 11th International LS-DYNA Users Conference,
June 6-8, 2010, Dearborn, Michigan.
4.
Wang B., Tsao W., Application of FEA and EMA to Structural Model Verification. Proceedings of the 10th
CSSV Conference Taiwan, 2002; 131-138.
5.
Huang Y., Wang B., Implementation and validation of frequency response function in LS-DYNA. FEA
information, July 2009.
6.
Huang Y., Wang B., Mode-based frequency response function and steady state dynamics in LS-DYNA.
Proceedings of the 11th International LS-DYNA Users Conference, June 6-8, 2010, Dearborn, Michigan.
7.
Chopra A. K., Dynamics of Structures: Theory and Applications to Earthquake Engineering (2nd Edition),
Prentice Hall, 2000.
8.
Huang Y., Ashcraft C., L'Eplattenier P., A fast collocation boundary element method for vibro-acoustic
problems in LS-DYNA. Proceedings of the Inaugural International Conference of the Engineering
Mechanics Institute (EM08), May 18-21, 2008, Minneapolis, Minnesota.
9.
Huang Y., Souli M., Development of acoustic and vibro acoustic solvers in LS-DYNA. Proceedings of
the 10th International LS-DYNA Users Conference, June 8-10, 2008, Dearborn, Michigan.
10. Huang Y., Souli M., Liu R., New developments of frequency domain acoustic methods in LS-DYNA.
Proceedings of the 11th International LS-DYNA Users Conference, June 6-8, 2010, Dearborn, Michigan.
Page 12