1.4.17 Distributed traction and edge loads

Products: ABAQUS/Standard  ABAQUS/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 *DLOAD and *DSLOAD options.

I. 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    C3D8R    C3D6    C3D10    C3D10M    C3D15    C3D20    C3D27   

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 the *DLOAD and *DSLOAD options. 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.

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.

II. 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 the *DLOAD and *DSLOAD options. 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.

III. Distributed shear and general traction loads in geometrically nonlinear analyses

Elements tested

CPS3    CPE3    CPS4    CPE4    CPS6    CPE6    CPS6M    CPE6M    CPS8    CPE8   

C3D4    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 the *DLOAD and *DSLOAD options 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 CONSTANT RESULTANT parameters 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 with the ORIENTATION paremeter.

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.

IV. 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 the *DLOAD and *DSLOAD options 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 CONSTANT RESULTANT parameters are also used.

Results and discussion

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

Input files

ABAQUS/Standard input files

traclarge_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_xpl.inp

Edge loading of shell elements.

V. Dead load analysis of a membrane structure using the CONSTANT RESULTANT parameter

Element tested

M3D4   

Problem description

This section provides basic verification of the CONSTANT RESULTANT parameter 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 , where 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 direction. A uniform dead traction load (of magnitude 4) is applied in the negative -direction. This could be considered a simple model of a sloped roof with a snow load.

Let 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, , is

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

In the first step the load is applied with CONSTANT RESULTANT=NO. In the second step the structure is unloaded. In the third step the load is applied with CONSTANT RESULTANT=YES.

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.