4.10.8 3DNLG-8: Collapse of a straight pipe segment under pure bending

Product: ABAQUS/Standard  

Elements tested

S3    S3R    S4    S4R    S4R5    S8R    S8R5    S9R5   

STRI3    STRI65   

Problem description

Elastic and elastic-plastic cases of a beam under pure bending are tested. Nonlinear multipoint constraints and section ovalization are also tested.

Model:

A quarter model is used because of symmetry.

Material:

Young's modulus = 1.93 × 105, Poisson's ratio = 0.26, plastic (isotropic hardening).


Equiv. StressEquiv. Plastic Strain
2720.0
3460.00473
3790.01264
4040.02836
4240.0491

Boundary conditions:

Symmetry on plane = 0 (; symmetry on plane  = 0 (; symmetry on rotated plane: i) Shell nodes must remain in rotated plane (Nonlinear MPC). ii) The -rotation of shell nodes in the rotated plane must be the same as the master node (Linear MPC). iii) The component of rotation about the local -axis must be zero (Nonlinear MPC). The master node is constrained to move on the -axis only.

Loading:

-rotation applied to pipe end ( = 0.004).

Reference solution

This is a test recommended by the National Agency for Finite Element Methods and Standards (U.K.): Test 3DNLG-8 from NAFEMS Publication R0024 “A Review of Benchmark Problems for Geometric Non-linear Behaviour of 3-D Beams and Shells (SUMMARY).”

The published results of this problem were obtained with ABAQUS. Thus, a comparison of ABAQUS and NAFEMS results is not an independent verification of ABAQUS. The NAFEMS study includes results from other sources for comparison that may provide a basis for verification of this problem.

Results and discussion

In the following table, the end moments for the elastic and elastic-plastic cases are reported at three rotation values of the master node. All the meshes have the same nodal spacing.

Element = 0.0012 = 0.0024 = 0.004
 Moment × 10–9Moment × 10–9Moment × 10–9
 ElasticPlasticElasticPlasticElasticPlastic
S3/S3R2.0560.51463.4880.54023.7040.5066
S42.0320.50883.4520.53503.6760.5042
S4R2.0200.50803.4360.53323.6640.5024
S4R52.0180.50783.4200.53343.5780.5008
S8R2.0540.51163.4480.53803.5680.5024
S8R52.0540.51183.4520.53803.5820.5030
S9R52.0540.51183.4520.53803.5820.5030
STRI32.0560.51443.4740.53903.6620.5036
STRI652.0540.51223.4520.53883.5920.5046

Response predicted by ABAQUS

Essentially identical moment-rotation curves (at the master node) are obtained for all test cases. The response predicted using S8R5 elements is shown below.

Input files

Elasticity:


n3g8xf31.inp

S3/S3R elements.

n3g8xf31.f

User subroutine MPC used in n3g8xf31.inp.

n3g8xe41.inp

S4 elements.

n3g8xe41.f

User subroutine MPC used in n3g8xe41.inp.

n3g8xf41.inp

S4R elements.

n3g8xf41.f

User subroutine MPC used in n3g8xf41.inp.

n3g8x541.inp

S4R5 elements.

n3g8x541.f

User subroutine MPC used in n3g8x541.inp.

n3g8x681.inp

S8R elements.

n3g8x681.f

User subroutine MPC used in n3g8x681.inp.

n3g8x581.inp

S8R5 elements.

n3g8x581.f

User subroutine MPC used in n3g8x581.inp.

n3g8x591.inp

S9R5 elements.

n3g8x591.f

User subroutine MPC used in n3g8x591.inp.

n3g8x631.inp

STRI3 elements.

n3g8x631.f

User subroutine MPC used in n3g8x631.inp.

n3g8x561.inp

STRI65 elements.

n3g8x561.f

User subroutine MPC used in n3g8x561.inp.

Elastoplasticity:


n3g8xf32.inp

S3/S3R elements.

n3g8xf32.f

User subroutine MPC used in n3g8xf32.inp.

n3g8xe42.inp

S4 elements.

n3g8xe42.f

User subroutine MPC used in n3g8xe42.inp.

n3g8xf42.inp

S4R elements.

n3g8xf42.f

User subroutine MPC used in n3g8xf42.inp.

n3g8x542.inp

S4R5 elements.

n3g8x542.f

User subroutine MPC used in n3g8x542.inp.

n3g8x682.inp

S8R elements.

n3g8x682.f

User subroutine MPC used in n3g8x682.inp.

n3g8x582.inp

S8R5 elements.

n3g8x582.f

User subroutine MPC used in n3g8x582.inp.

n3g8x592.inp

S9R5 elements.

n3g8x592.f

User subroutine MPC used in n3g8x592.inp.

n3g8x632.inp

STRI3 elements.

n3g8x632.f

User subroutine MPC used in n3g8x632.inp.

n3g8x562.inp

STRI65 elements.

n3g8x562.f

User subroutine MPC used in n3g8x562.inp.