Professional Documents
Culture Documents
Pressed Elastoplastic Pipe
Pressed Elastoplastic Pipe
This model is licensed under the COMSOL Software License Agreement 5.3.
All trademarks are the property of their respective owners. See www.comsol.com/trademarks.
Introduction
In offshore applications, it is sometimes necessary to quickly seal a pipe as part of the
prevention of a blowout. This example shows a simulation, in which a circular pipe is
squeezed between two flat stiff indenters.
The tutorial serves as an example of an analysis with very large plastic strains and contact.
Model Definition
The pipe has an external radius, R0, of 200 mm and a wall thickness of 25 mm. The
material is a stainless steel with the following properties
TABLE 1: MATERIAL DATA
PROPERTY VALUE
The hardening curve is available as a text file containing pairs of data (plastic strain, stress)
which can be imported as a function. The function is shown in Figure 1below. This
nonlinear hardening function σh(εpe) is defined in the Materials node by an interpolation
function. The stress is measured as true (Cauchy) stress, and the strain is measured as true
(logarithmic) strain. The data can thus be used directly as input to COMSOL when large
strain plasticity is used. Note that if the strain is above the ultimate strain, then the curve
is implicitly flat, so that deformation continues under constant stress.
The pipe is compressed between two flat indenters that can be considered as rigid. The
original position of each indenter is 0.1 mm from the outer pipe wall. During the
Due to the symmetries, only one quarter of the geometry needs to be modeled. The
problem is considered as 2D with the plane strain assumption.
The distribution of the effective plastic strains is shown in Figure 4 and Figure 5. As can
be seen, the peak value (0.95) is far above the ultimate strain (0.52). All values above the
ultimate strain are however on the inside of the pipe, where the strain state is mainly in
compression. At the outer edge of the pipe, the plastic strain approximately reaches the
ultimate strain. There is thus a certain risk that cracks could start forming. The values of
ultimate stress and strain are related to the specimen used for the testing (usually a
cylindrical bar) and cannot directly be transferred to general multiaxial conditions.
Figure 5: Effective plastic strain at maximum compression, detail. The black contour
indicates the ultimate strain (0.52).
Figure 6: Deformed shape and residual stresses at the end of the process.
The load used to compress the pipe is computed from the reaction force in the indenter,
and it is shown in Figure 7.
Modeling Instructions
From the File menu, choose New.
NEW
In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D.
2 In the Select Physics tree, select Structural Mechanics>Solid Mechanics (solid).
3 Click Add.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Interpolation 1 (int1)
1 On 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 compr.
4 In the table, enter the following settings:
t f(t)
0 0
1 0.15
2 0
GEOMETRY 1
Circle 1 (c1)
1 On the Geometry toolbar, click Primitives and choose Circle.
Circle 2 (c2)
1 On the Geometry toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, locate the Size and Shape section.
3 In the Radius text field, type Ri.
4 In the Sector angle text field, type 90.
Difference 1 (dif1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object c1 only.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Active toggle button.
5 Select the object c2 only.
6 Click Build All Objects.
Rectangle 1 (r1)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type 0.05.
4 In the Height text field, type 1.3*Ro.
5 Locate the Position section. In the x text field, type Ro+0.0001.
Plasticity 1
1 On the Physics toolbar, click Attributes and choose Plasticity.
2 Select Domain 1 only.
3 In the Settings window for Plasticity, locate the Plasticity Model section.
4 From the Plasticity model list, choose Large plastic strains.
5 Find the Isotropic hardening model subsection. From the list, choose User defined.
MATERIALS
In the Model Builder window, under Component 1 (comp1) right-click Materials and choose
Blank Material.
Material 1 (mat1)
1 In the Settings window for Material, locate the Material Contents section.
Interpolation 1 (int1)
1 Right-click Component 1 (comp1)>Materials>Material 1 (mat1)>Elastoplastic material
model (ElastoplasticModel) and choose Functions>Interpolation.
2 In the Settings window for Interpolation, locate the Definition section.
3 Click Load from File.
4 Browse to the application’s Application Libraries folder and double-click the file
compressed_elastoplastic_pipe_stress_strain.txt.
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1)>Materials click Material 1
(mat1).
2 In the Settings window for Material, locate the Material Contents section.
Contact 1
1 On the Physics toolbar, in the Boundary section, click Pairs and choose Contact.
2 In the Settings window for Contact, locate the Pair Selection section.
3 In the Pairs list, select Contact Pair 1 (p1).
4 Locate the Contact Pressure Method section. From the list, choose Penalty.
Symmetry 1
1 On the Physics toolbar, click Boundaries and choose Symmetry.
2 Select Boundaries 1 and 2 only.
Prescribed Displacement 1
1 On the Physics toolbar, click Domains and choose Prescribed Displacement.
2 Select Domain 2 only.
3 In the Settings window for Prescribed Displacement, locate the Prescribed Displacement
section.
4 Select the Prescribed in x direction check box.
5 Select the Prescribed in y direction check box.
6 In the u 0x text field, type -compr(para).
MESH 1
1 In the Model Builder window, under Component 1 (comp1) click Mesh 1.
2 In the Settings window for Mesh, locate the Mesh Settings section.
3 From the Sequence type list, choose User-controlled mesh.
Free Triangular 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 right-click Free
Triangular 1 and choose Delete.
2 Click Yes to confirm.
Mapped 1
In the Model Builder window, under Component 1 (comp1)>Mesh 1 right-click Mapped 1 and
choose Distribution.
Distribution 1
1 In the Settings window for Distribution, locate the Distribution section.
2 In the Number of elements text field, type 1.
3 Select Boundaries 5 and 6 only.
One element is enough on the indenter, since the whole domain is under displacement
control.
Mapped 1
Right-click Mapped 1 and choose Distribution.
Distribution 2
1 In the Settings window for Distribution, locate the Distribution section.
2 From the Distribution properties list, choose Predefined distribution type.
3 Select Boundaries 1 and 2 only.
4 In the Number of elements text field, type 8.
5 In the Element ratio text field, type 2.
6 Select the Symmetric distribution check box.
Mapped 1
Right-click Mapped 1 and choose Distribution.
Distribution 3
1 Select Boundary 3 only.
2 In the Settings window for Distribution, locate the Distribution section.
3 In the Number of elements text field, type 80.
4 Click Build All.
STUDY 1
Step 1: Stationary
Set up an auxiliary continuation sweep for the parameter para.
1 In the Settings window for Stationary, click to expand the Study extensions section.
5 Click to expand the Results while solving section. Locate the Results While Solving section.
Select the Plot check box.
6 On the Home toolbar, click Compute.
RESULTS
Stress (solid)
Mirror the solution twice to get a full view of the pipe.
Mirror 2D 1
On the Results toolbar, click More Data Sets and choose Mirror 2D.
Mirror 2D 2
1 On the Results toolbar, click More Data Sets and choose Mirror 2D.
2 In the Settings window for Mirror 2D, locate the Data section.
3 From the Data set list, choose Mirror 2D 1.
4 Locate the Axis Data section. In row Point 2, set x to 1 and y to 0.
Stress (solid)
1 In the Model Builder window, under Results click Stress (solid).
2 In the Settings window for 2D Plot Group, locate the Data section.
3 From the Data set list, choose Mirror 2D 2.
4 On the Stress (solid) toolbar, click Plot.
5 Click the Zoom Extents button on the Graphics toolbar.
Plot the stresses at the maximum compression. This gives Figure 3.
6 From the Parameter value (para) list, choose 1.
7 On the Stress (solid) toolbar, click Plot.
Create an animation of the compression process.
8 Click Animation and choose Player.
Surface Integration 1
1 On the Results toolbar, click More Derived Values and choose Integration>Surface
Integration.
2 Select Domain 2 only.
3 In the Settings window for Surface Integration, locate the Expressions section.
4 In the table, enter the following settings:
Expression Unit
-solid.RFx*2/1E6 N
5 Click Evaluate.
TABLE
1 Go to the Table window.
2 Click Table Graph in the window toolbar.
RESULTS
1D Plot Group 2
1 In the Model Builder window, under Results click 1D Plot Group 2.
2 In the Settings window for 1D Plot Group, type Compression Force in the Label text
field.
3 Locate the Plot Settings section. Select the y-axis label check box.
4 In the associated text field, type Force/unit length [MN/m].
5 On the Compression Force toolbar, click Plot.
2D Plot Group 3
On the Home toolbar, click Add Plot Group and choose 2D Plot Group.
Surface 1
1 In the Model Builder window, right-click 2D Plot Group 3 and choose Surface.
2 In the Settings window for Surface, click Replace Expression in the upper-right corner of
the Expression section. From the menu, choose Component 1>Solid Mechanics>Strain
(Gauss points)>solid.epeGp - Effective plastic strain.
Set the plot scale limit to the ultimate strain.
Deformation 1
1 Right-click Results>2D Plot Group 3>Surface 1 and choose Deformation.
2 In the Settings window for Deformation, locate the Scale section.
3 Select the Scale factor check box.
4 In the associated text field, type 1.
Add a contour showing where the ultimate strain is exceeded.
2D Plot Group 3
In the Model Builder window, under Results right-click 2D Plot Group 3 and choose Contour.
Contour 1
1 In the Settings window for Contour, click Replace Expression in the upper-right corner
of the Expression section. From the menu, choose Component 1>Solid Mechanics>Strain
(Gauss points)>solid.epeGp - Effective plastic strain.
2 Locate the Levels section. From the Entry method list, choose Levels.
3 In the Levels text field, type 0.52.
4 Locate the Coloring and Style section. From the Coloring list, choose Uniform.
5 From the Color list, choose Black.
6 Clear the Color legend check box.
Deformation 1
1 Right-click Results>2D Plot Group 3>Contour 1 and choose Deformation.
2 In the Settings window for Deformation, locate the Scale section.
3 Select the Scale factor check box.
4 In the associated text field, type 1.
5 On the 2D Plot Group 3 toolbar, click Plot.
6 Click the Zoom Extents button on the Graphics toolbar.
2D Plot Group 3
1 In the Model Builder window, under Results click 2D Plot Group 3.
2 In the Settings window for 2D Plot Group, type Plastic Strain in the Label text field.