Models - Sme.vibrating Deep Beam
Models - Sme.vibrating Deep Beam
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
This model studies free and forced vibrations of a deep beam. With the increase of the ratio
of section area per beam length, shear deformations and rotational inertia effects can no
longer be neglected as it is done in the Euler-Bernoulli theory.
The solution for Eigenfrequency, frequency response and transient analysis are computed
using a Timoshenko beam and compared with analytical results.
Model Definition
The model studied in this example consists of a simply supported beam with a square cross
section. One end is pinned and has a constrained rotation along the beam axis. At the other
end, the displacements in the plane of the beam cross section are constrained.
For the forced vibration cases, a load in the y direction is applied all along the beam.
GEOMETRY
• Beam length, L = 10 m
• Beam cross section dimension l = 2 m
MATERIAL
• Young’s modulus, E = 200 GPa
• Poisson’s ratio, ν = 0.3
• Mass density, ρ = 8000 kg/m3
DAMPING
The Rayleigh damping coefficients are α = 5.36 s-1, β = 7.46·10−5 m/s. The damping
ratio curve is shown in Figure 1.
CONSTRAINTS
At x = 0, the displacements are u = v = w = 0 and the rotation around the x-axis is fixed at
θx = 0
At x = 10, v = w = 0
LOAD CASES
The load with a magnitude F0 = 106 N/m is applied and oriented in the local positive
y direction. The following forced vibration cases are studied:
MESH
To satisfy the benchmark specifications, five edge elements are used to mesh the geometry.
Free vibration
In Table 1 the computed results are compared with the analytical results for the free
vibration case. The agreement is good. The accuracy decreases with increasing complexity
of the mode shape because the possibility for the relatively coarse mesh to describe such a
shape is limited.
TABLE 1: COMPARISON BETWEEN ANALYTICAL AND COMPUTED NATURAL FREQUENCIES.
Reference
1. J. Maguire, D.J. Dawswell, and L. Gould, “Selected Benchmarks for Forced Vibration”,
NAFEMS R0016, 1989.
NEW
In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select Physics tree, select Structural Mechanics>Beam (beam).
3 Click Add.
4 Click Study.
5 In the Select Study tree, select General Studies>Eigenfrequency.
6 Click Done.
GLOBAL DEFINITIONS
Parameters 1
1 In the Model Builder window, under Global Definitions click Parameters 1.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
GEOMETRY 1
Polygon 1 (pol1)
1 In the Geometry toolbar, click More Primitives and choose Polygon.
2 In the Settings window for Polygon, locate the Coordinates section.
3 In the table, enter the following settings:
MATERIALS
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, locate the Material Contents section.
3 In the table, enter the following settings:
BEAM (BEAM)
1 In the Model Builder window, under Component 1 (comp1) click Beam (beam).
2 In the Settings window for Beam, locate the Beam Formulation section.
3 From the list, choose Timoshenko.
Cross-Section Data 1
1 In the Model Builder window, under Component 1 (comp1)>Beam (beam) click Cross-
Section Data 1.
2 In the Settings window for Cross-Section Data, locate the Cross-Section Definition section.
3 From the list, choose Common sections.
4 In the hy text field, type 2.
5 In the hz text field, type 2.
Section Orientation 1
1 In the Model Builder window, expand the Cross-Section Data 1 node, then click
Section Orientation 1.
2 In the Settings window for Section Orientation, locate the Section Orientation section.
3 From the Orientation method list, choose Orientation vector.
0 X
0 Y
1 Z
Prescribed Displacement/Rotation 1
1 In the Physics toolbar, click Points and choose Prescribed Displacement/Rotation.
2 Select Point 1 only.
3 In the Settings window for Prescribed Displacement/Rotation, locate the
Prescribed Displacement section.
4 From the Displacement in x direction list, choose Prescribed.
5 From the Displacement in y direction list, choose Prescribed.
6 From the Displacement in z direction list, choose Prescribed.
7 Locate the Prescribed Rotation section. From the list, choose Rotation.
8 Select the Free rotation around y direction check box.
9 Select the Free rotation around z direction check box.
Prescribed Displacement/Rotation 2
1 In the Physics toolbar, click Points and choose Prescribed Displacement/Rotation.
2 Select Point 2 only.
3 In the Settings window for Prescribed Displacement/Rotation, locate the
Prescribed Displacement section.
4 From the Displacement in y direction list, choose Prescribed.
5 From the Displacement in z direction list, choose Prescribed.
MESH 1
Edge 1
1 In the Mesh toolbar, click Boundary and choose Edge.
2 In the Settings window for Edge, locate the Edge Selection section.
3 From the Geometric entity level list, choose Entire geometry.
Distribution 1
1 Right-click Edge 1 and choose Distribution.
2 Right-click Distribution 1 and choose Build All.
Step 1: Eigenfrequency
1 In the Model Builder window, under Free Vibration click Step 1: Eigenfrequency.
2 In the Settings window for Eigenfrequency, locate the Study Settings section.
3 Select the Desired number of eigenfrequencies check box. In the associated text field, type
10.
Solution 1 (sol1)
1 In the Study toolbar, click Show Default Solver.
2 In the Model Builder window, expand the Solution 1 (sol1) node, then click
Eigenvalue Solver 1.
3 In the Settings window for Eigenvalue Solver, locate the General section.
4 In the Relative tolerance text field, type 1.0E-15.
5 Click Compute to Selected.
RESULTS
Line 1
1 In the Model Builder window, expand the Mode Shape (beam) node, then click Line 1.
2 In the Settings window for Line, locate the Coloring and Style section.
3 In the Radius scale factor text field, type 0.1.
BEAM (BEAM)
Edge Load 1
1 In the Physics toolbar, click Edges and choose Edge Load.
2 In the Settings window for Edge Load, locate the Edge Selection section.
3 From the Selection list, choose All edges.
0 X
0 Y
F0 Z
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 In the αdM text field, type 5.36.
4 In the βdK text field, type 7.46e-5.
In order to visualize the damping ratio curve, create the Damping Ratio plot through an
action button from the Damping Settings section.
5 Click Damping Ratio Preview in the upper-right corner of the Damping Settings section.
From the menu, choose Create Damping Ratio Plot.
RESULTS
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 Studies subsection. In the Select Study tree, select
Preset Studies for Selected Physics Interfaces>Frequency Domain, Modal.
4 Click Add Study in the window toolbar.
5 In the Home toolbar, click Add Study to close the Add Study window.
RESULTS
Cut Point 3D 1
1 In the Results toolbar, click Cut Point 3D.
2 In the Settings window for Cut Point 3D, locate the Data section.
3 From the Dataset list, choose Harmonic Forced Vibration/Solution 2 (sol2).
4 Locate the Point Data section. In the X text field, type L/2.
5 In the Y text field, type 0.
6 In the Z text field, type 0.
Point Graph 1
1 Right-click Harmonic Peak Displacement and choose Point Graph.
2 In the Settings window for Point Graph, locate the y-Axis Data section.
3 In the Expression text field, type beam.uAmpZ.
4 From the Unit list, choose mm.
5 In the Harmonic Peak Displacement toolbar, click Plot.
Point Graph 1
1 In the Model Builder window, expand the Harmonic Peak Stress node, then click
Point Graph 1.
2 In the Settings window for Point Graph, locate the y-Axis Data section.
3 In the Expression text field, type abs(beam.sb1).
4 From the Unit list, choose MPa.
5 In the Harmonic Peak Stress toolbar, click Plot.
8 Click Evaluate.
Expression Unit
abs(beam.sb1) MPa
DEFINITIONS
Interpolation 1 (int1)
1 In the Home toolbar, click Functions and choose Global>Interpolation.
2 In the Settings window for Interpolation, locate the Definition section.
3 In the Function name text field, type phase.
4 In the table, enter the following settings:
t f(t)
20 -pi/2
60 pi/2
5 Locate the Interpolation and Extrapolation section. From the Interpolation list, choose
Nearest neighbor.
6 Locate the Units section. In the Argument table, enter the following settings:
Argument Unit
t Hz
Function Unit
phase rad
BEAM (BEAM)
Edge Load 2
1 In the Physics toolbar, click Edges and choose Edge Load.
2 Select Edge 1 only.
3 In the Settings window for Edge Load, locate the Force section.
0 X
0 Y
F0 Z
Phase 1
1 In the Physics toolbar, click Attributes and choose Phase.
2 In the Settings window for Phase, locate the Load Phase section.
3 Specify the φ vector as
0 X
0 Y
phase(freq) 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 Studies subsection. In the Select Study tree, select General Studies>
Frequency Domain.
4 Right-click and choose Add Study.
5 In the Home toolbar, click Add Study to close the Add Study window.
STUDY 3
Solution 4 (sol4)
1 In the Study toolbar, click Show Default Solver.
2 In the Model Builder window, expand the Solution 4 (sol4) node, then click FFT Solver 1.
3 In the Settings window for FFT Solver, locate the General section.
4 From the Defined by study step list, choose User defined.
5 In the Model Builder window, click Study 3.
6 In the Settings window for Study, type Periodic Forced Vibration in the Label text
field.
7 Locate the Study Settings section. Clear the Generate default plots check box.
8 In the Study toolbar, click Compute.
RESULTS
Cut Point 3D 2
1 In the Results toolbar, click Cut Point 3D.
2 In the Settings window for Cut Point 3D, locate the Data section.
3 From the Dataset list, choose Periodic Forced Vibration/Solution 4 (sol4).
4 Locate the Point Data section. In the X text field, type L/2.
5 In the Y text field, type 0.
6 In the Z text field, type 0.
Point Graph 1
1 Right-click Periodic Peak Displacement and choose Point Graph.
Point Graph 1
1 In the Model Builder window, expand the Periodic Peak Stress node, then click
Point Graph 1.
2 In the Settings window for Point Graph, locate the y-Axis Data section.
3 In the Expression text field, type beam.sb1.
4 From the Unit list, choose MPa.
5 In the Periodic Peak Stress toolbar, click Plot.
Expression Unit
abs(w) mm
6 Locate the Data Series Operation section. From the Transformation list, choose Maximum.
7 Click Evaluate.
Expression Unit
abs(beam.sb1) MPa
BEAM (BEAM)
Edge Load 3
1 In the Physics toolbar, click Edges and choose Edge Load.
2 In the Settings window for Edge Load, locate the Edge Selection section.
3 From the Selection list, choose All edges.
4 Locate the Force section. Specify the FL vector as
0 X
0 Y
F0 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 Studies subsection. In the Select Study tree, select General Studies>
Time Dependent.
4 Click Add Study in the window toolbar.
5 In the Home toolbar, click Add Study to close the Add Study window.
STUDY 4
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 Transient Forced Vibration>
Solver Configurations>Solution 6 (sol6)>Dependent Variables 1 node, then click
Displacement field (comp1.u).
4 In the Settings window for Field, locate the Scaling section.
5 In the Scale text field, type 1e-4.
6 In the Model Builder window, under Transient Forced Vibration>Solver Configurations>
Solution 6 (sol6) click Time-Dependent Solver 1.
7 In the Settings window for Time-Dependent Solver, click to expand the Time Stepping
section.
8 From the Steps taken by solver list, choose Intermediate.
9 Click Compute.
RESULTS
Cut Point 3D 3
1 In the Results toolbar, click Cut Point 3D.
2 In the Settings window for Cut Point 3D, locate the Data section.
3 From the Dataset list, choose Transient Forced Vibration/Solution 6 (sol6).
4 Locate the Point Data section. In the X text field, type L/2.
5 In the Y text field, type 0.
6 In the Z text field, type 0.
Point Graph 1
1 Right-click Transient Peak Displacement and choose Point Graph.
2 In the Settings window for Point Graph, locate the y-Axis Data section.
3 In the Expression text field, type w.
4 From the Unit list, choose mm.
5 In the Transient Peak Displacement toolbar, click Plot.
6 Click Plot.
3 Locate the Data section. From the Dataset list, choose Transient Forced Vibration/
Solution 6 (sol6).
4 From the Time selection list, choose From list.
5 In the Times (s) list, select 0.0117.
6 Select Edge 1 only.
8 Locate the Data Series Operation section. From the Transformation list, choose Maximum.
9 Click Evaluate.
3 Locate the Expressions section. In the table, enter the following settings:
Expression Unit
abs(beam.sb1) MPa
3 Locate the Data section. From the Time selection list, choose Interpolated.
4 In the Times (s) text field, type range(1.1,1e-3,1.2).
5 Locate the Expressions section. In the table, enter the following settings:
6 Locate the Data Series Operation section. From the Transformation list, choose Average.
7 Click next to Evaluate, then choose New Table.
TABLE
1 Go to the Table window.
The vibration analysis is now finished. If you want to save the model and run the study
with the configuration set in the previous steps, you need to disable some of the nodes
FREE VIBRATION
Step 1: Eigenfrequency
1 In the Model Builder window, under Free Vibration click Step 1: Eigenfrequency.
2 In the Settings window for Eigenfrequency, locate the Physics and Variables Selection
section.
3 Select the Modify model configuration for study step check box.
4 In the tree, select Component 1 (comp1)>Beam (beam)>Linear Elastic Material 1>
Damping 1.
5 Right-click and choose Disable.