Distributed traction and edge loads

This problem contains basic test cases for one or more Abaqus elements and features.

The following topics are discussed:

ProductsAbaqus/StandardAbaqus/Explicit

Features tested

This section provides basic verification tests for the traction load labels TRVEC and TRSHR and the edge load labels EDLD, EDNOR, EDSHR, and EDTRA using the distributed element-based and surface loads.

Distributed shear and general traction loads

Elements tested

  • CPS3
  • CPE3
  • CPS4
  • CPE4
  • CPS6
  • CPE6
  • CPS6M
  • CPE6M
  • CPS8
  • CPE8
  • CPEG3
  • CPEG4
  • CPEG6
  • CPEG6M
  • CPEG8
  • CAX3
  • CAX4
  • CAX6
  • CAX6M
  • CAX8
  • CGAX3
  • CGAX4
  • CGAX6
  • CGAX6M
  • CGAX8
  • C3D4
  • C3D5
  • C3D8R
  • C3D6
  • C3D10
  • C3D10M
  • C3D15
  • C3D20
  • C3D27
  • CSS8
  • CCL9
  • CCL12
  • CCL18
  • CCL24
  • S3R
  • STRI3
  • S4R
  • S4R5
  • STRI65
  • S8R
  • S8R5
  • S9R5
  • SC6R
  • SC8R
  • SAX1
  • SAX2
  • RAX2
  • M3D3
  • M3D4
  • M3D6
  • M3D8
  • M3D9
  • MAX1
  • MAX2
  • MGAX1
  • MGAX2
  • MCL6
  • MCL9
  • SFMCL6
  • SFMCL9
  • SFM3D3
  • SFM3D4
  • SFM3D6
  • SFM3D8
  • SFMAX1
  • SFMAX2
  • SFMGAX1
  • SFMGAX2
  • R2D2
  • R3D3
  • R3D4
  • RAX2

Problem description

The analyses in this section test the traction load labels TRVEC and TRSHR using distributed element-based and surface loads. One-element and two-element tests are performed to verify the loading options on all the faces of supported elements. In both Abaqus/Standard and Abaqus/Explicit tests, the elements are held fixed by kinematic coupling constraints as each face of each element is loaded with a combination of distributed general tractions and shear tractions. The resultant forces at the kinematic reference nodes are output to verify that distributed loads are properly applied to each element.

Results and discussion

The results for each combination indicate that the loads are applied correctly.

Input files

Abaqus/Standard input files

tracload2d.inp

Traction loading of two-dimensional elements.

tracloadcpeg.inp

Traction loading of generalized plane strain elements.

tracloadcax.inp

Traction loading of axisymmetric elements.

tracloadcgax.inp

Traction loading of axisymmetric elements with twist.

tracload3d.inp

Traction loading of three-dimensional elements.

tracloadcss8.inp

Traction loading of continuum solid shell elements.

tracloadccl.inp

Traction loading of cylindrical elements.

tracloadshl.inp

Traction loading of shell elements.

tracloadsc.inp

Traction loading of continuum shell elements.

tracloadrsax.inp

Traction loading of axisymmetric shell elements and axisymmetric rigid link elements.

tracloadm3d.inp

Traction loading of three-dimensional membrane and surface elements.

tracloadmax.inp

Traction loading of axisymmetric membrane elements.

tracloadmgax.inp

Traction loading of axisymmetric membrane elements with twist.

tracloadmcl.inp

Traction loading of cylindrical membrane elements.

tracloadsfmax.inp

Traction loading of axisymmetric surface elements.

tracloadsfmgax.inp

Traction loading of axisymmetric surface elements with twist.

tracloadr2d.inp

Traction loading of two-dimensional rigid elements.

tracloadr3d.inp

Traction loading of three-dimensional rigid elements.

Abaqus/Explicit input files

tracload2d_xpl.inp

Traction loading of two-dimensional elements.

tracloadcax_xpl.inp

Traction loading of axisymmetric elements.

tracload3d_xpl.inp

Traction loading of three-dimensional elements.

tracloadshl_xpl.inp

Traction loading of shell, membrane, and surface elements.

tracloadsc_xpl.inp

Traction loading of continuum shell elements.

tracloadrsax_xpl.inp

Traction loading of axisymmetric shell elements and axisymmetric rigid link elements.

tracloadr2d2_xpl.inp

Traction loading of two-dimensional rigid elements.

Distributed edge loads

Elements tested

  • S3R
  • STRI3
  • S4R
  • S4R5
  • STRI65
  • S8R
  • S8R5
  • S9R5

Problem description

The analyses in this section test the edge load labels EDLD, EDNOR, EDSHR, and EDTRA using distributed element-based and surface loads. One-element and two-element tests are performed to verify the loading options on all the edges of supported shell elements. In both Abaqus/Standard and Abaqus/Explicit tests, the elements are held fixed by kinematic coupling constraints as each edge of each element is loaded with a combination of distributed edge loads. The resultant forces at the kinematic reference nodes are output to verify that distributed loads are properly applied to each element.

Results and discussion

The results for each combination indicate that the loads are applied correctly.

Input files

Abaqus/Standard input file

tracloadedge.inp

Edge loading of shell elements.

Abaqus/Explicit input file

tracloadedge_xpl.inp

Edge loading of shell elements.

Distributed shear and general traction loads in geometrically nonlinear analyses

Elements tested

  • CPS3
  • CPE3
  • CPS4
  • CPE4
  • CPS6
  • CPE6
  • CPS6M
  • CPE6M
  • CPS8
  • CPE8
  • C3D4
  • C3D5
  • C3D8R
  • C3D6
  • C3D10
  • C3D10M
  • C3D15
  • C3D20
  • CCL9
  • CCL12
  • CCL18
  • CCL24
  • S3R
  • STRI3
  • S4R
  • S4R5
  • STRI65
  • S8R
  • S8R5
  • S9R5
  • SC6R
  • SC8R
  • SAX1
  • SAX2

Problem description

The analyses in this section test the traction load labels TRVEC and TRSHR using distributed element-based and surface loads in geometrically nonlinear analyses. Tests include models under large rigid body rotations and large deformations. In the tests where elements undergo large rigid body rotations, one facet is coupled to a kinematic coupling reference node. A traction load is applied to another face. This load is kept constant as the elements are rotated by the kinematic coupling reference node. The reaction forces at the kinematic reference node are used to verify that the loads are properly applied and rotated with the element. Different combinations of the follower and non-follower surface loads and constant resultants are also used. Some of the models in the tests have cylindrical geometry. General traction or shear loadings are applied on the cylindrical surface by defining a local cylindrical coordinate system.

Results and discussion

The results for each combination indicate that the loads are applied correctly.

Input files

Abaqus/Standard input files

traclarge_rotation_2d.inp

Traction loading of two-dimensional elements.

traclarge_rotation_3d.inp

Traction loading of three-dimensional elements.

traclarge_rotation_3d_usub.inp

User-defined traction loading of three-dimensional elements.

traclarge_rotation_3d_usub.f

User subroutine used in traclarge_rotation_3d_usub.inp.

traclarge_rotation_shl.inp

Traction loading of three-dimensional shell elements.

traclarge_rotation_m3d.inp

Traction loading of three-dimensional membrane elements.

tracnlgeom_ccl9.inp

Traction loading of 9-node cylindrical element CCL9.

tracnlgeom_ccl12.inp

Traction loading of 12-node cylindrical element CCL12.

tracnlgeom_ccl12_usub.inp

User-defined traction loading of 12-node cylindrical element CCL12.

tracnlgeom_ccl12_usub.f

User subroutine used in tracnlgeom_ccl12_usub.inp.

tracnlgeom_ccl18.inp

Traction loading of 18-node cylindrical element CCL18.

tracnlgeom_ccl24.inp

Traction loading of 24-node cylindrical element CCL24.

tracnlgeom_sax.inp

Traction loading of axisymmetric shell element.

trac_cylori.inp

Traction loading of a three-dimensional cylinder.

Abaqus/Explicit input files

traclarge_rotation_2d_xpl.inp

Traction loading of two-dimensional elements.

trac_cylori_xpl.inp

Traction loading of a three-dimensional cylinder.

Distributed edge loads in a geometrically nonlinear analysis

Elements tested

  • S3R
  • STRI3
  • S4R
  • S4R5
  • STRI65
  • S8R
  • S8R5
  • S9R5

Problem description

The analyses in this section test the edge load labels EDLD, EDNOR, EDSHR, and EDTRA using distributed element-based and surface loads in geometrically nonlinear analyses. One facet is coupled to a kinematic coupling reference node. A traction load is applied to another face. This load is kept constant as the elements are rotated by the kinematic coupling reference node. The reaction forces at the kinematic reference node are used to verify that the loads are properly applied and rotated with the element. Different combinations of the follower and non-follower surface loads and constant resultants are also used.

Results and discussion

The results for each combination indicate that the loads are applied correctly.

Input files

Abaqus/Standard input files

tracedgelarge_rotation.inp

Edge loading of shell elements.

tracnlgeom_edge_usub.inp

User-defined edge loading of shell elements.

tracnlgeom_edge_usub.f

User subroutine used in tracnlgeom_edge_usub.inp.

Abaqus/Explicit input file

traclarge_rotation_edge_xpl.inp

Edge loading of shell elements.

Dead load analysis of a membrane structure using a constant resultant

Elements tested

  • M3D4

Problem description

This section provides basic verification of using constant resultants in a dead load analysis. The constant resultant method has certain advantages when a traction is used to model a distributed load with a known constant resultant.

If you choose not to have a constant resultant, the traction vector is integrated over the surface in the current configuration, a surface that in general deforms in a geometrically nonlinear analysis. The most common example of a traction that should be integrated over the current configuration is a live pressure load defined as t=-pn, where n is the normal in the current configuration. The total resultant due to a pressure load depends on the surface area in the current configuration. A live uniform normal surface traction integrated over the current surface is equivalent to applying a uniform pressure load. By default, the traction vector is integrated over the surface in the current configuration.

If you choose to have a constant resultant, the traction vector is integrated over the surface in the reference configuration, which is constant.

The analysis in this section consists of a unit planar membrane structure that is held fixed at the edges by a kinematic coupling constraint. The normal of the flat structure is in the (1,0,1) direction. A uniform dead traction load (of magnitude 4) is applied in the negative e3-direction. This could be considered a simple model of a sloped roof with a snow load.

Let So=1 and S denote the total surface area of the plate in the reference and current configurations, respectively. With no constant resultant, the total integrated load on the plate, f, is

f=StdS=S-4e2dS=-4e3S.

In this case a uniform traction leads to a resultant load that increases as the surface area of the plate increases, which is not consistent with a fixed snow load. With the constant resultant method, the total integrated load on the plate is

f=SotdSo=So-4e3dSo=-4e3So=-4e3.

In the first step the load is applied without a constant resultant. In the second step the structure is unloaded. In the third step the load is applied with a constant resultant.

Results and discussion

The magnitude of the reaction force at the kinematic coupling reference node at the end of the first step is 4.59. A reaction force greater than 4.0 reflects the fact that the surface area of the membrane is increasing with the load. The magnitude of the reaction force at the kinematic coupling reference node at the end of the third step is 4.0 as expected.

Input files

Abaqus/Standard input files

tracresultant_m3d4.inp

Testing the CONSTANT RESULTANT parameter.

tracresultant_m3d4_usub.inp

User-defined traction loading with the CONSTANT RESULTANT parameter.

tracresultant_m3d4_usub.f

User subroutine used in tracresultant_m3d4_usub.inp.

Abaqus/Explicit input file

tracresultant_m3d4_xpl.inp

Testing the CONSTANT RESULTANT parameter.