# Dynamics

Evaluate your design in real-time using SimSolid dynamics analysis.

## Setup

- Specify the modal results to which the analysis is linked. The modal solution must exist in the current design study. In SimSolid, the time integration of the equations of motion is extremely fast and all modes are always included in the analysis.
- For Frequency and Random Response, specify the frequency span upper and lower limits. For Transient response, specify Time span.
- Specify damping using Rayleigh damping coefficients or Modal damping.
- Select the Evaluate peak responses during solving check box to evaluate peak responses during solving phase.

See Create Analysis for additional information.

## Damping

- Rayleigh Damping Coefficient
- Assumes the damping matrix is proportional to the mass and stiffness matrices. You need to specify values for Mass (F1) and Stiffness (F2) in the Dynamics creation dialog to use this method.
- Modal Damping
- Creates critical damping ratio for each mode. You can specify this value in the Dynamics analysis creation dialog.

## Notes for Dynamics Analysis

- When the base excitation type is displacement, the initial condition for displacement and velocity is always assumed to be zero.
- In SimSolid, the boundary compatibility is approximately met. The response at the constrained end is not going to be an absolute zero but is relatively small compared to the peak responses.
- Equivalent radiated power density is calculated as:
(1) $$\text{ERP}\text{\hspace{0.33em}}\text{Density}\text{\hspace{0.17em}}\text{=}\text{\hspace{0.17em}}\text{ERPRLF}\text{\hspace{0.17em}}\text{*}\text{\hspace{0.17em}}\text{(0}\text{.5}\text{\hspace{0.17em}}\text{*}\text{\hspace{0.17em}}\text{ERPC}\text{\hspace{0.17em}}\text{*}\text{\hspace{0.17em}}\text{ERPRHO)}\text{\hspace{0.17em}}*\text{\hspace{0.17em}}{v}^{2}$$Where:- $v$
- Normal velocity of the picked point
- ERPC (Speed of sound in air)
- 343 m/s
- ERPRHO (Density of air)
- 1.225 Kg/m3
- ERPRLF (Radiation loss factor)
- 1

Equivalent radiated power is calculated as an integral of ERP density over picked faces as:

(2) $$\text{ERP}\text{\hspace{0.17em}}\text{=}\text{\hspace{0.17em}}\text{ERPRLF}\text{\hspace{0.17em}}\text{*}\text{\hspace{0.17em}}\text{(0}\text{.5}\text{\hspace{0.17em}}\text{*}\text{\hspace{0.17em}}\text{ERPC}\text{\hspace{0.17em}}\text{*}\text{\hspace{0.17em}}\text{ERPRHO)}\text{\hspace{0.17em}}{\int}_{S}\text{\hspace{0.05em}}{v}^{2}\text{\hspace{0.05em}}ds$$ - Phase for Absolute displacement can be queried using Pick Info for frequency dynamics.
- Relationship between relative and absolute displacement in frequency
dynamics is as follows.Relative motion is calculated as:
(3) $$\ddot{x}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}2\zeta {\omega}_{n}\dot{x}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\omega}_{n}{}^{2}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}-\ddot{b}$$Where base excitation $b$ is:(4) $$-\ddot{b}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{Y}_{0}\left(\omega \right){e}^{i\left(\omega t\text{\hspace{0.17em}}+\text{\hspace{0.17em}}\varphi \right)}$$Solving this differential equation, the relative displacement can be calculated as:(5) $${x}_{r}\left(t\right)\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{Y}_{0}\left(\omega \right){e}^{i\varphi}\text{\hspace{0.17em}}\left(\frac{{\omega}_{n}{}^{2}\text{\hspace{0.17em}}-\text{\hspace{0.17em}}{\omega}^{2}}{{\left({\omega}_{n}{}^{2}\text{\hspace{0.17em}}-\text{\hspace{0.17em}}{\omega}^{2}\right)}^{2}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\left(2\zeta {\omega}_{n}\omega \right)}^{2}}\text{\hspace{0.17em}}-\text{\hspace{0.17em}}i\frac{2\zeta {\omega}_{n}\omega}{{\left({\omega}_{n}{}^{2}\text{\hspace{0.17em}}-\text{\hspace{0.17em}}{\omega}^{2}\right)}^{2}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\left(2\zeta {\omega}_{n}\omega \right)}^{2}}\right){e}^{i\omega t}$$Absolute motion is calculated as:(6) $$\ddot{x}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}2\zeta {\omega}_{n}\dot{x}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\omega}_{n}{}^{2}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}2\zeta {\omega}_{n}\dot{b}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\omega}_{n}{}^{2}b$$Where,(7) $$b\text{\hspace{0.17em}}=\text{\hspace{0.17em}}\frac{{Y}_{0}\left(\omega \right)}{{\omega}^{2}}{e}^{i\left(\omega t\text{\hspace{0.17em}}+\text{\hspace{0.17em}}\varphi \right)}$$Solving this equation, the absolute displacement can be calculated as:(8) $${x}_{a}\left(t\right)\text{\hspace{0.17em}}=\text{\hspace{0.17em}}\left({\left(\frac{{\omega}_{n}}{\omega}\right)}^{2}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}i2\zeta \frac{{\omega}_{n}}{\omega}\right){x}_{r}\left(t\right)$$ - In Frequency and Random dynamics, the Complex Function Method is used to
solve differential equations. Displacement, velocity, and acceleration
results have complex components.Given complex values of displacement as:
(9) ${D}_{x}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{a}_{x}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}i{b}_{x}$${D}_{y}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{a}_{y}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}i{b}_{y}$${D}_{z}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{a}_{z}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}i{b}_{z}$Displacement magnitude is calculated as:(10) $$\left|\sqrt{{D}_{x}{}^{2}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{D}_{x}{}^{2}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{D}_{x}{}^{2}}\right|\text{\hspace{0.17em}}=\text{\hspace{0.17em}}\sqrt{\left|{D}_{x}{}^{2}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{D}_{x}{}^{2}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{D}_{x}{}^{2}\right|}\text{\hspace{0.17em}}$$ - When a transient dynamics analysis is linked to a prestressed modal analysis, a new result type is offered called Total Displacement. Total Displacement is the combination of the prestressed and dynamic displacements. Therefore, displacement magnitude is the displacement caused by the dynamic analysis.
- For a random response analysis, the Power Spectral Density (PSD) of the
response
${S}_{xo}(f)$
, is related to the power spectral density of
the source,
${S}_{a}(f)$
, by:
(11) Where, ${H}_{xa}(f)$ is the frequency response function.$${S}_{xo}(f)={\left|{H}_{xa}(f)\right|}^{2}{S}_{a}(f)$$For better understanding, let us take an example of base excitation with acceleration as the excitation type as an input to a random response analysis.

The base excitation with amplitude is used to define the input for frequency response analysis, so the units of the acceleration excitation type would be either of below highlighted units, m/sec

^{2}; mm/sec^{2}; cm/sec^{2}; G; in/sec^{2}Units for the PSD function depend on the boundary condition. In this example, as base excitation is given as acceleration, the unit for PSD will be (mm/s^{2})^{2}/Hz.