Models - Cfd.sports Car Fsi
Models - Cfd.sports Car Fsi
Models - Cfd.sports Car Fsi
This model is licensed under the COMSOL Software License Agreement 6.1.
All trademarks are the property of their respective owners. See www.comsol.com/trademarks.
Introduction
While a CFD analysis can give the car designer important input for drag coefficients and
turbulence around the car panels, a fluid–structure-interaction (FSI) analysis provides
insights into how the fluid forces interact with the body, thereby generating vibrations that
result in noise inside the car and potentially cause fatigue problems.
In this model, we perform a structural analysis of the front door of a sports car traveling
at a speed of 180 km/h. The forces from a transient Large Eddy Simulation (LES) model
are transformed to the frequency domain using the Fast Fourier Transform (FFT). A
structural analysis of the door and side mirror is performed in the frequency domain to
study the Fourier components of the forcing and the resulting vibrational modes of the car
door. The door geometry is simplified compared to a real door of a sports car, but the
focus in this model is to demonstrate a method for applying the loads from a CFD-LES
analysis in the time domain to a structural-mechanics analysis in the frequency domain.
Figure 1: Car-door geometry with applied loads at 50 Hz. The highest magnitude is in the
area around the side mirror attachment.
The forces are added directly from the LES solution of the complete car. Since the car door
coordinates are the same in both components, no coordinate transformation of the forces
is needed. The forces are mapped from one component to the other using an extrusion
coupling operator. This reduces the size of the mesh and makes the simulation leaner and
more efficient. The mapped time-dependent forces are transformed to the frequency
domain using a Fast Fourier Transform (FFT).
Figure 3 shows the von Mises stress on the door surface at 50 Hz. There is a concentration
of stress at the attachment of the side mirror. This indicates that the area may need
reinforcement.
Figure 4 and Figure 5 show the displacement field on the inside door panel at 50 Hz and
the large displacement of the mirror at 70 Hz, respectively. Note that the displacement is
exaggerated for visual reasons. The real deformation is in the range of millimeters at these
frequencies.
Modeling Instructions
APPLICATION LIBRARIES
1 From the File menu, choose Application Libraries.
2 In the Application Libraries window, select CFD Module>Single-Phase Flow>sports_car in
the tree.
3 Click Open.
COMPONENT 1 (COMP1)
In the Model Builder window, expand the Component 1 (comp1) node.
7 Click to expand the Advanced section. From the Mesh search method list, choose
Closest point.
GEOMETRY 2
1 In the Settings window for Geometry, locate the Advanced section.
2 From the Geometry representation list, choose CAD kernel.
Import 1 (imp1)
1 In the Home toolbar, click Import.
2 In the Settings window for Import, locate the Import section.
3 Click Browse.
4 Browse to the model’s Application Libraries folder and double-click the file
sports_car_fsi.mphbin.
ADD MATERIAL
1 Go to the Add Material window.
2 In the tree, select Built-in>Aluminum 3003-H18.
3 Click Add to Component in the window toolbar.
4 In the tree, select Built-in>Glass (quartz).
5 Click Add to Component in the window toolbar.
6 In the Home toolbar, click Add Material to close the Add Material window.
MATERIALS
To make the model slightly more efficient, it is convenient to map the surface stresses from
the fluid flow simulation onto a set of new variables. The following steps add new variables
for the fluid forces exerted on the wall, assuming that automatic wall treatment is used in
the LES-model. Then, a weak form equation is added and solved for the extruded stresses
on the car door. The time dependent stresses are Fourier transformed so that the following
structural mechanics simulations can be performed in the frequency domain.
DEFINITIONS (COMP1)
Variables 1
1 In the Model Builder window, expand the Component 1 (comp1)>Materials node, then
click Component 1 (comp1)>Definitions>Variables 1.
2 In the Settings window for Variables, locate the Variables section.
3 In the table, enter the following settings:
COMPONENT 2 (COMP2)
1 In the Model Builder window, click Component 2 (comp2).
2 In the Home toolbar, click Windows and choose Add Physics.
ADD PHYSICS
1 Go to the Add Physics window.
2 In the tree, select Mathematics>PDE Interfaces>Lower Dimensions>
Weak Form Boundary PDE (wb).
3 Find the Physics interfaces in study subsection. In the table, clear the Solve check boxes
for Study 1, Study 2, and Study 3.
4 Click Add to Component 2 in the window toolbar.
5 In the Home toolbar, click Add Physics to close the Add Physics window.
8 Click to expand the Discretization section. From the Element order list, choose Linear.
9 Click to expand the Dependent Variables section. In the Field name text field, type
Tstress.
Tstress1
ADD STUDY
1 In the Home toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
3 Find the Physics interfaces in study subsection. In the table, clear the Solve check boxes
for LES RBVM (spf), Laplace’s Equation (lpeq), and Turbulent Flow, k- 2 (spf2).
4 Find the Studies subsection. In the Select Study tree, select General Studies>
Time Dependent.
5 Click Add Study in the window toolbar.
6 In the Home toolbar, click Add Study to close the Add Study window.
STUDY 4
7 Click to expand the Mesh Selection section. In the table, enter the following settings:
Component Mesh
Component 1 No mesh
Solution 4 (sol4)
1 In the Study toolbar, click Show Default Solver.
Modify the solver settings so that the time stepping is performed for the exact same
points in time as in the solution of the fluid dynamics simulation, and set the nonlinear
solver to take only one iteration per time step.
2 In the Model Builder window, expand the Solution 4 (sol4) node, then click Time-
Dependent Solver 1.
3 In the Settings window for Time-Dependent Solver, click to expand the Time Stepping
section.
4 From the Steps taken by solver list, choose Manual.
5 In the Time step text field, type .002.
6 Locate the General section. From the Times to store list, choose Steps taken by solver.
DEFINITIONS (COMP2)
In the Model Builder window, under Component 2 (comp2) click Definitions.
Integration 2 (intop2)
1 In the Definitions toolbar, click Nonlocal Couplings and choose Integration.
2 In the Settings window for Integration, locate the Source Selection section.
3 From the Geometric entity level list, choose Boundary.
4 Select Boundary 1 only.
Integration 3 (intop3)
1 In the Definitions toolbar, click Nonlocal Couplings and choose Integration.
2 In the Settings window for Integration, locate the Source Selection section.
3 From the Geometric entity level list, choose Boundary.
4 Select Boundaries 48–53, 55, 59–64, 67–72, 85–87, 92, 93, 97–104, and 109–111
only.
Variables 2
1 Right-click Definitions and choose Variables.
2 In the Settings window for Variables, locate the Variables section.
3 In the table, enter the following settings:
MESH 3
1 In the Model Builder window, expand the Component 2 (comp2)>Mesh 3 node, then click
Mesh 3.
2 In the Settings window for Mesh, locate the Sequence Type section.
3 From the list, choose User-controlled mesh.
Size
1 In the Model Builder window, under Component 2 (comp2)>Mesh 3 click Size.
2 In the Settings window for Size, locate the Element Size section.
3 From the Calibrate for list, choose Fluid dynamics.
4 From the Predefined list, choose Fine.
Free Triangular 1
1 In the Model Builder window, click Free Triangular 1.
2 In the Settings window for Free Triangular, locate the Boundary Selection section.
3 From the Geometric entity level list, choose Entire geometry.
4 Click Build All.
STUDY 4
In the Home toolbar, click Compute.
RESULTS
1D Plot Group 13
1 In the Home toolbar, click Add Plot Group and choose 1D Plot Group.
2 In the Settings window for 1D Plot Group, locate the Data section.
3 From the Dataset list, choose Study 4/Solution 4 (5) (sol4).
Global 1
1 Right-click 1D Plot Group 13 and choose Global.
2 In the Settings window for Global, click Add Expression in the upper-right corner of the
y-Axis Data section. From the menu, choose Component 2 (comp2)>Definitions>
Variables>F_mirror - Norm of mirror force - N.
You can now clear the solution from Studies 1-3, to make the file size smaller. The results
needed for the structural mechanics analysis have already been mapped to new variables.
STUDY 3
In the Model Builder window, expand the Study 3 node.
Solution 3 (sol3)
1 In the Model Builder window, expand the Study 3>Solver Configurations node.
2 Right-click Solution 3 (sol3) and choose Clear.
STUDY 2
In the Model Builder window, expand the Study 2 node.
Solution 2 (sol2)
1 In the Model Builder window, expand the Study 2>Solver Configurations>Solution 2 (sol2)
node.
2 Right-click Solution 2 (sol2) and choose Clear.
STUDY 1
In the Model Builder window, expand the Study 1 node.
Solution 1 (sol1)
1 In the Model Builder window, expand the Study 1>Solver Configurations node.
2 Right-click Solution 1 (sol1) and choose Clear.
Look at the plot for the extruded force components to verify that the transform is working.
RESULTS
Add a structural mechanics interface to the car door structure. The car door is modeled as
surfaces, thus, using a shell interface is convenient in this context. The shell is assumed to
be thin, relative to the size of the car door.
ADD PHYSICS
1 In the Home toolbar, click Add Physics to open the Add Physics window.
2 Go to the Add Physics window.
3 In the tree, select Structural Mechanics>Shell (shell).
4 Find the Physics interfaces in study subsection. In the table, clear the Solve check boxes
for Study 1, Study 2, Study 3, and Study 4.
5 Click Add to Component 2 in the window toolbar.
6 In the Home toolbar, click Add Physics to close the Add Physics window.
MATERIALS
SHELL (SHELL)
Damping 1
1 In the Physics toolbar, click Attributes and choose Damping.
2 In the Settings window for Damping, locate the Damping Settings section.
3 From the Damping type list, choose Isotropic loss factor.
4 From the s list, choose User defined. In the associated text field, type 0.03.
3 In the Settings window for Thickness and Offset, locate the Thickness and Offset section.
4 In the d0 text field, type 1.8[mm].
Fixed Constraint 1
1 In the Physics toolbar, click Edges and choose Fixed Constraint.
Spring Foundation 1
1 In the Physics toolbar, click Boundaries and choose Spring Foundation.
2 Select Boundaries 3, 58, and 131 only.
3 In the Settings window for Spring Foundation, locate the Spring section.
Spring Foundation 2
1 In the Physics toolbar, click Edges and choose Spring Foundation.
2 Select Edge 1 only.
3 In the Settings window for Spring Foundation, locate the Spring section.
Face Load 1
In the Physics toolbar, click Boundaries and choose Face Load.
Create a selection from the weak form physics boundaries. This saves you from recreating
the selection for the boundary load on the external door surfaces.
SHELL (SHELL)
Face Load 1
1 In the Model Builder window, under Component 2 (comp2)>Shell (shell) click Face Load 1.
2 In the Settings window for Face Load, locate the Boundary Selection section.
-withsol('sol4',comp2.Tstress1,setval(freq,freq))/T0s x
-withsol('sol4',comp2.Tstress2,setval(freq,freq))/T0s y
-withsol('sol4',comp2.Tstress3,setval(freq,freq))/T0s z
ADD STUDY
1 In the Home toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
3 Find the Physics interfaces in study subsection. In the table, clear the Solve check boxes
for LES RBVM (spf), Laplace’s Equation (lpeq), Turbulent Flow, k- 2 (spf2), and
Weak Form Boundary PDE (wb).
4 Find the Studies subsection. In the Select Study tree, select
Preset Studies for Selected Physics Interfaces>Frequency Domain, Modal.
5 Click Add Study in the window toolbar.
6 In the Home toolbar, click Add Study to close the Add Study window.
STUDY 5
Step 1: Eigenfrequency
1 In the Settings window for Eigenfrequency, locate the Study Settings section.
2 Select the Desired number of eigenfrequencies check box. In the associated text field, type
50.
Component Mesh
Component 1 No mesh
Solution 6 (sol6)
1 In the Study toolbar, click Show Default Solver.
2 In the Model Builder window, expand the Solution 6 (sol6) node.
3 In the Model Builder window, expand the Study 5>Solver Configurations>
Solution 6 (sol6)>Dependent Variables 1 node, then click
Turbulent dissipation rate (comp1.ep2).
4 In the Settings window for Field, locate the General section.
5 Clear the Store in output check box.
6 In the Model Builder window, under Study 5>Solver Configurations>Solution 6 (sol6)>
Dependent Variables 1 click Turbulent kinetic energy (comp1.k2).
7 In the Settings window for Field, locate the General section.
8 Clear the Store in output check box.
9 In the Model Builder window, under Study 5>Solver Configurations>Solution 6 (sol6)>
Dependent Variables 1 click Pressure (comp1.p).
10 In the Settings window for Field, locate the General section.
11 Clear the Store in output check box.
12 In the Model Builder window, under Study 5>Solver Configurations>Solution 6 (sol6)>
Dependent Variables 1 click Pressure (comp1.p2).
13 In the Settings window for Field, locate the General section.
14 Clear the Store in output check box.
15 In the Model Builder window, under Study 5>Solver Configurations>Solution 6 (sol6)>
Dependent Variables 1 click Velocity field (comp1.u).
16 In the Settings window for Field, locate the General section.
17 Clear the Store in output check box.
Step 1: Eigenfrequency
1 In the Model Builder window, under Study 5 click Step 1: Eigenfrequency.
2 In the Settings window for Eigenfrequency, click to expand the Mesh Selection section.
3 In the table, enter the following settings:
Component Mesh
Component 1 No mesh
Component Mesh
Component 1 No mesh
RESULTS
Stress (shell)
1 In the Stress (shell) toolbar, click Plot.
2 In the Settings window for 3D Plot Group, locate the Data section.
3 From the Parameter value (freq (Hz)) list, choose 50.
4 In the Stress (shell) toolbar, click Plot.
5 In the Home toolbar, click Add Predefined Plot.
Displacement (shell)
1 In the Settings window for 3D Plot Group, locate the Data section.
2 From the Parameter value (freq (Hz)) list, choose 50.
3 In the Displacement (shell) toolbar, click Plot.
RESULTS