Large Deformation Beam

Solved with COMSOL Multiphysics 4.
L arg e D e f o r m a ti on Beam
Model Definition
In this example you study the deflection of a cantilever beam undergoing very large
deflections. The model is called Straight Cantilever GNL Benchmark and is
described in detail in section 5.2 of NAFEMS Background to Finite Element Analysis
of Geometric Non-linearity Benchmarks (Ref. 1). A schematic description of the
beam and its characteristics is shown in Figure 1.
y
3.2 m
Cross section
x
0.1 m
0.1 m
Figure 1: Cantilever beam geometry.

GEOMETRY
The length of the beam is 3.2 m.

The cross section is a square with side lengths 0.1 m.
MATERIAL
The beam is linear elastic with E = 2.11011 N/m2 and = 0.

CONSTRAINTS AND LOADS
The left end is fixed. This boundary condition is compatible with beam theory
assumptions only in the case when = 0.
The right end is subjected to distributed loads with the resultants
Fx = 3.844106 N and Fy = 3.844103 N.
Results and Discussion

Due to the large compressive axial load and the slender geometry, this is a buckling
problem. If you are to study the buckling and post-buckling behavior of a symmetric
problem, it is necessary to perturb the symmetry somewhat. Here the small transversal
1 |
LARGE DEFORMATION BEAM
Solved with COMSOL Multiphysics 4.4
load serves this purpose. An alternative approach would be to introduce an initial

imperfection in the geometry.
Figure 2 below shows the final state with the 1:1 displacement scaling.
Figure 2: The effective von Mises stress of the deformed beam.

The horizontal and vertical displacements of the tip versus the compressive load
normalized by its maximum value are shown in Figure 3.
2 |
Figure 3: Horizontal (solid) and vertical (dashed) tip displacements versus normalized
compressive load.
Table 1 contains a summary of some significant results. Because the reference values
are given as graphs, an estimate of the error caused by reading this graph is added:
TABLE 1: COMPARISON BETWEEN MODEL RESULTS AND REFERENCE VALUES.
QUANTITY
COMSOL MULTIPHYSICS
REFERENCE
Maximum vertical displacement at the tip
-2.58
-2.58 0.02
Final vertical displacement at the tip
-1.34
-1.36 0.02
Final horizontal displacement at the tip
-5.07
-5.04 0.04
The results are in excellent agreement, especially considering the coarse mesh used.
The plot of the axial deflection reveals that an instability occurs at a parameter value
close to 0.1, corresponding to the compressive load 3.84105 N. It is often seen in
practice that the critical load of an imperfect structure is significantly lower than that
of the ideal structure.
This problem (without the small transverse load) is usually referred to as the Euler-1
case. The theoretical critical load is
3 |
11
EI 2.1 10 ( 0.1 12 )
5
P c = -----------= -----------------------------------------------------------------= 4.22 10 N
2
2
4L
4 3.2
(1)
Reference
1. A.A. Becker, Background to Finite Element Analysis of Geometric Non-linearity
Benchmarks, NAFEMS, Ref: -R0065, Glasgow, 1999.
Model Library path: Structural_Mechanics_Module/Verification_Models/

large_deformation_beam
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click the Model Wizard button.

MODEL WIZARD
1 In the Model Wizard window, click the 2D button.

2 In the Select physics tree, select Structural Mechanics>Solid Mechanics (solid).
3 Click the Add button.
4 Click the Study button.
5 In the tree, select Preset Studies>Stationary.
6 Click the Done button.
GEOMETRY 1
Rectangle 1
1 In the Model Builder window, under Component 1 right-click Geometry 1 and choose
Rectangle.
2 In the Rectangle settings window, locate the Size section.
3 In the Width edit field, type 3.2.
4 In the Height edit field, type 0.1.
4 |
Form Union
In the Model Builder window, under Component 1>Geometry 1 right-click Form Union
and choose Build Selected.
GLOBAL DEFINITIONS
Define parameters for the compressive and transverse load components as well as a
parameter that you will use to gradually turn up the compressive load.
Parameters
1 On the Home toolbar, click Parameters.
2 In the Parameters settings window, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
F_Lx
-3.844[MN]
-3.844E6 N
Maximum compressive load
F_Ly
1e-3*F_Lx
-3844 N
Transverse load
NCL
Normalized compressive load
By restricting the range for the parameter NCL to [0, 1], it serves as a compressive
load normalized by the maximum compressive load.
MATERIALS
Material 1
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose New Material.

2 Select Domain 1 only.
3 In the Material settings window, locate the Material Contents section.
Property
Name
Value
Unit
Property group
Young's modulus
2.1e11
Pa
Basic
Poisson's ratio
nu
Basic
Density
rho
7850
kg/m
Basic
SOLID MECHANICS
1 In the Model Builder window, under Component 1 click Solid Mechanics.

2 In the Solid Mechanics settings window, locate the 2D Approximation section.
5 |
3 From the 2D approximation list, choose Plane stress.

4 Locate the Thickness section. In the d edit field, type 0.1.
Linear Elastic Material 1

1 In the Model Builder window, expand the Solid Mechanics node, then click Linear
Elastic Material 1.
2 Select Domain 1 only.
Fixed Constraint 1
1 On the Physics toolbar, click Boundaries and choose Fixed Constraint.
2 Select Boundary 1 only.
Boundary Load 1
1 On the Physics toolbar, click Boundaries and choose Boundary Load.
2 Select Boundary 4 only.
3 In the Boundary Load settings window, locate the Force section.
4 From the Load type list, choose Total force.
5 Specify the Ftot vector as
NCL*F_Lx
F_Ly
MESH 1
In the Model Builder window, under Component 1 right-click Mesh 1 and choose
Build All.
STUDY 1
Step 1: Stationary
1 In the Model Builder window, expand the Study 1 node, then click Step 1: Stationary.
2 In the Stationary settings window, locate the Study Settings section.
3 Select the Include geometric nonlinearity check box.
Set up a continuation sweep for the NCL parameter.

4 Click to expand the Study extensions section.
5 Select the Auxiliary sweep check box.
6 Click Add.
6 |

Auxiliary parameter
Parameter value list
NCL
range(0,0.01,1)
Solver 1
1 On the Study toolbar, click Show Default Solver.
2 In the Model Builder window, expand the Study 1>Solver Configurations node.
3 In the Model Builder window, expand the Solver 1 node, then click Stationary Solver 1.
4 In the Stationary Solver settings window, locate the General section.
5 In the Relative tolerance edit field, type 1e-6.
6 In the Model Builder window, under
Study 1>Solver Configurations>Solver 1>Stationary Solver 1 click Fully Coupled 1.
7 In the Fully Coupled settings window, click to expand the Method and termination
section.
8 From the Nonlinear method list, choose Constant (Newton).
For problems with a strong influence of geometric nonlinarity, it is often more

efficient to not allow damping in the solver.
9 On the Home toolbar, click Compute.
RESULTS
Stress (solid)
1 In the Surface settings window, locate the Expression section.
2 From the Unit list, choose MPa.
3 On the Stress (solid) toolbar, click Plot.
4 Click the Zoom Extents button on the Graphics toolbar.
Compare the resulting plot with the one shown in Figure 2.

Next, add a data set to use for plotting the displacement components at the tip and for
displaying their values as functions of the normalized compressive load.
Data Sets
1 On the Results toolbar, click Cut Point 2D.
2 In the Cut Point 2D settings window, locate the Point Data section.
3 In the X edit field, type 3.2.
4 In the Y edit field, type 0.05.
7 |
5 Click the Plot button.

6 Click the Zoom Extents button on the Graphics toolbar.
Derived Values
1 On the Results toolbar, click Point Evaluation.
2 In the Point Evaluation settings window, locate the Data section.
3 From the Data set list, choose Cut Point 2D 1.
4 Click Replace Expression in the upper-right corner of the Expression section. From
the menu, choose Solid Mechanics>Displacement>Displacement field

(Material)>Displacement field, X component (u).
5 Click the Evaluate button.
As stated in Table 1, the final horizontal displacement value is roughly -5.07 m.

6 Click Replace Expression in the upper-right corner of the Expression section. From
the menu, choose Solid Mechanics>Displacement>Displacement field

(Material)>Displacement field, Y component (v).
7 Click the Evaluate button.
By scrolling in the Table window table you can verify that the maximum vertical
displacement is 2.58 m (downward) and occurs at a normalized compressive load of
0.19.
8 |
1D Plot Group 2
1 On the Home toolbar, click Add Plot Group and choose 1D Plot Group.
2 In the 1D Plot Group settings window, locate the Data section.
3 From the Data set list, choose Cut Point 2D 1.
4 On the 1D Plot Group 2 toolbar, click Point Graph.
5 In the Point Graph settings window, click Replace Expression in the upper-right
corner of the y-axis data section. From the menu, choose Solid
Mechanics>Displacement>Displacement field (Material)>Displacement field, X
component (u).
6 On the 1D Plot Group 2 toolbar, click Point Graph.
7 In the Point Graph settings window, click Replace Expression in the upper-right
corner of the y-axis data section. From the menu, choose Solid
Mechanics>Displacement>Displacement field (Material)>Displacement field, Y
component (v).
8 Click to expand the Coloring and style section. Find the Line style subsection. From
the Line list, choose Dashed.

9 In the Model Builder window, click 1D Plot Group 2.
10 In the 1D Plot Group settings window, click to expand the Title section.
11 From the Title type list, choose Manual.
12 In the Title text area, type Tip displacement components (m) vs. normalized
compressive load.
13 Locate the Plot Settings section. Select the y-axis label check box.
14 In the associated edit field, type Tip displacement.
15 On the 1D Plot Group 2 toolbar, click Plot.
9 |
10 |

Large Deformation Beam

Uploaded by

Copyright:

Available Formats

Large Deformation Beam

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Large Deformation Beam

Uploaded by

Copyright:

Available Formats

Solved with COMSOL Multiphysics 4.

Figure 1: Cantilever beam geometry.

The length of the beam is 3.2 m.

The beam is linear elastic with E = 2.11011 N/m2 and = 0.

Results and Discussion

LARGE DEFORMATION BEAM

Solved with COMSOL Multiphysics 4.4

load serves this purpose. An alternative approach would be to introduce an initial

Figure 2: The effective von Mises stress of the deformed beam.

LARGE DEFORMATION BEAM

Solved with COMSOL Multiphysics 4.4

Maximum vertical displacement at the tip

Final vertical displacement at the tip

Final horizontal displacement at the tip

LARGE DEFORMATION BEAM

Solved with COMSOL Multiphysics 4.4

Model Library path: Structural_Mechanics_Module/Verification_Models/

1 In the New window, click the Model Wizard button.

1 In the Model Wizard window, click the 2D button.

LARGE DEFORMATION BEAM

Solved with COMSOL Multiphysics 4.4

Maximum compressive load

Normalized compressive load

choose New Material.

1 In the Model Builder window, under Component 1 click Solid Mechanics.

LARGE DEFORMATION BEAM

Solved with COMSOL Multiphysics 4.4

3 From the 2D approximation list, choose Plane stress.

Linear Elastic Material 1

Set up a continuation sweep for the NCL parameter.

LARGE DEFORMATION BEAM

Solved with COMSOL Multiphysics 4.4

7 In the table, enter the following settings:

Parameter value list

For problems with a strong influence of geometric nonlinarity, it is often more

Compare the resulting plot with the one shown in Figure 2.

LARGE DEFORMATION BEAM

Solved with COMSOL Multiphysics 4.4

5 Click the Plot button.

the menu, choose Solid Mechanics>Displacement>Displacement field

As stated in Table 1, the final horizontal displacement value is roughly -5.07 m.

the menu, choose Solid Mechanics>Displacement>Displacement field

LARGE DEFORMATION BEAM

Solved with COMSOL Multiphysics 4.4

the Line list, choose Dashed.

LARGE DEFORMATION BEAM

Solved with COMSOL Multiphysics 4.4

LARGE DEFORMATION BEAM

You might also like