
Inviscid Bump in a Channel
Inviscid Supersonic Wedge
Inviscid ONERA M6
Laminar Flat Plate
Laminar Cylinder
Turbulent Flat Plate
Transitional Flat Plate
Transitional Flat Plate for T3A and T3A
Turbulent ONERA M6
Unsteady NACA0012
Epistemic Uncertainty Quantification of RANS predictions of NACA 0012 airfoil
Nonideal compressible flow in a supersonic nozzle

Unconstrained shape design of a transonic inviscid airfoil at a cte. AoA
Constrained shape design of a transonic turbulent airfoil at a cte. C_{L}
Constrained shape design of a transonic inviscid wing at a cte. C_{L}
Shape Design With Multiple Objectives and Penalty Functions
Unsteady Shape Optimization NACA0012
Unconstrained shape design of a two way mixing channel
Linear Dynamics
Written by  for Version  Revised by  Revision date  Revised version 

@rsanfer  7.0.2  @rsanfer  20200128  7.0.2 
Solver: 

Uses: 

Prerequisites: 
Linear Elasticity 
Complexity: 
Basic 
Goals
Upon completion of the tutorial on Linear Elasticity, this document will guide you throught the following capabilities of SU2:
 Setting up a dynamic structural problem with small deformations
 Timedependent boundary conditions
 Newmark integration method
The problem that we will be solving consists of a vertical, slender cantilever, clamped in its base, and subject to a horizontal, timedependent load \(P\) on its left boundary. This is shown in Fig. 1.
Resources
For this tutorial, it is assumed that you have first completed the guide on Linear Elasticity guide. You can find the resources for this tutorial in the same structural_mechanics/cantilever folder in the Tutorials repository. You can reuse the mesh file mesh_cantilever.su2 but need a new config file, config_dynamic.cfg.
Background
This tutorial covers the solution of structural dynamics problems. The equation of the problem is now
\(\mathbf{M}\mathbf{\ddot{u}}+\mathbf{K}\mathbf{u}=\mathbf{F}\),
where \(\mathbf{K}\) and \(\mathbf{K}\) are respectively the stiffness and mass matrices of the cantilever, \(\mathbf{u}\) is the vector of displacements of the structural nodes, \(\mathbf{\ddot{u}}\) is the vector of accelerations of the structural nodes, and \(\mathbf{F}\) is the vector of applied forces, and the structural damping effects have been neglected.
This tutorial will be solved the for the cantilever with the following conditions:
 Youngâ€™s modulus: \(\mathrm{E = 5.0}\) \(\mathrm{GPa}\)
 Poisson ratio: \(\nu\) \(\mathrm{= 0.35}\)
 Solid density: \(\rho\) \(\mathrm{= 1000}\) \(\mathrm{kg/m^3}\)
 Height: \(\mathrm{H = 10}\) \(\mathrm{mm}\)
 Thickness: \(\mathrm{t = 0.5}\) \(\mathrm{mm}\)
Mesh Description
The cantilever is discretized using 1000 4node quad elements with linear interpolation. The boundaries are defined as follows:
Configuration File Options
We start from the configuration options used on the tutorial on [Linear Elasticity]. However, for this case set the analysis to be dynamic using
TIME_DOMAIN = YES
TIME_STEP = 0.01
TIME_ITER = 250
where the time step \(\Delta t\) has been set to 10 ms, and the total time of the simulation is \(T = 2.5\) s (250 time steps). We will be using the Newmark\(^1\) method for time integration, which relies on the approximations
\[\mathbf{u}^{t} = \mathbf{u}^{t1} + \Delta t \mathbf{\dot{u}}^{t1} + [(\frac{1}{2}  \beta)\mathbf{\ddot{u}}^{t1} + \beta \mathbf{\ddot{u}}^{t}] \Delta t^2\]and
\[\mathbf{\dot{u}}^{t} = \mathbf{\dot{u}}^{t1} + [(1  \gamma)\mathbf{\ddot{u}}^{t1} + \gamma \mathbf{\ddot{u}}^{t}] \Delta t,\]to advance the problem in time. Further details about the implementation can be found on this reference by Sanchez et al\(^2\), available for download at Spiral.
We will be introducing a 2% numerical damping, to damp out higher frequencies in the problem. This is done by setting
TIME_DISCRE_FEA = NEWMARK_IMPLICIT
NEWMARK_BETA = 0.2601
NEWMARK_GAMMA = 0.52
We maintain the structural properties from the previous example. In dynamic problems, inertial effects are incorporated throught the mass matrix \(\mathbf{M}\). Therefore, it is necessary to define the density of the material, which is done through the following option,
MATERIAL_DENSITY = 1000
Finally, we define a sinusoidal load, \(\mathrm{P}(t) = A \sin(2\pi f t+\phi)\mathrm{P}\), being \(f\) the frequency of the sine wave in Hz, \(A\) the amplitude and \(\phi\) the phase angle in rad. This introduces some dynamic effects into the problem. We do so using
SINE_LOAD = YES
SINE_LOAD_COEFF = (1.5, 2.0, 4.7124)
where SINE_LOAD_COEFF=
(\(A\), \(f\), \(\phi\)). We have set a phase angle \(\phi = 3\pi/2\) and an amplitude \(A = 1.5\), which leads to a load \(P(t) = 1.5P = 1.5 \textrm{kPa}\) at time \(t=0\). The frequency is \(f=2\) Hz, which leads to 5 full cycles of oscillation in our time span \(T = 5\) s.
Running SU2
This is a very small example that we can run in serial. To run this test case, follow these steps at a terminal command line:
 Move to the directory containing the config file (config_dynamic.cfg) and the mesh file (mesh_cantilever.su2). Make sure that the SU2 tools were compiled, installed, and that their install location was added to your path.

Run the executable by entering
$ SU2_CFD config_dynamic.cfg
at the command line.
 SU2 will print the residual of the linear solver for each time step, and the simulation will finish after 250 time steps.
 Files containing the results will be written upon exiting SU2. The flow solution can be visualized in ParaView (.vtk) or Tecplot (.dat for ASCII).
For this particular problem, the solver should have run very quickly. Please note, a large number of files (500) will be generated upon running SU2 for this case, 2 files corresponding for each time step.
The screen output of the code will go over your settings, before initializing the solver. Once done so, you will obtain an output to screen such as
++
 Time_Iter rms[DispX] rms[DispY] VonMises
++
 249 11.145342 11.309566 1.7175e+06
 Solver Exit 
Maximum number of time iterations reached (TIME_ITER = 250).

++
 File Writing Summary  Filename 
++
SU2 restart restart_linear_dynamic_00249.dat 
Paraview binary linear_dynamic_00249.vtk 
++
where rms[DispX]
and rms[DispY]
correspond to the residual of the linear solver in the x and y components of the displacements, and VonMises
corresponds to the maximum Von Mises Stress in the structural domain.
Results
The output files append the time step as linear_dynamic__*****.vtk. They contain the displacement field, the nodal tensions \(\sigma_{xx}\), \(\sigma_{yy}\) and \(\sigma_{xy}\) (as Sxx, Syy and Sxy), and the Von Misses stress. In order to visualize the deformation of the cantilever on Paraview, one can use the filter Warp By Vector applied on the displacement field.
The solution of the problem is shown next, where the horizontal displacement at the tip is plotted on the right part of the figure, and the deformed configuration on the left.
References
\(^1\) Newmark, N.M. (1959), A method of computation for structural dynamics, J Eng Mech Div, 85(3):6794
Attribution
If you are using this content for your research, please kindly cite the following reference (reference \(^2\) in the text above) in your derived works:
Sanchez, R. et al. (2016), Towards a FluidStructure Interaction solver for problems with large deformations within the opensource SU2 suite, 57th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, 48 January 2016, San Diego, California, USA. DOI: 10.2514/6.20160205

This work is licensed under a Creative Commons Attribution 4.0 International License
 Previous
 Next