1.6.13 Finite sliding between concentric cylinders—axisymmetric and CAXA models

Product: ABAQUS/Standard  

Element tested

ISL21A   

Features tested

*CONTACT PAIR

SLAVE, MASTER

*SLIDE LINE

*ASYMMETRIC-AXISYMMETRIC

Problem description

This example illustrates the use of ABAQUS slide line elements and contact surface definitions in an axisymmetric structure that may undergo nonlinear, nonaxisymmetric deformation. This contact problem involves the relative motion of two outer cylinders with respect to one another and with respect to an inner, constrained cylinder. The axisymmetric model is shown in Figure 1.6.13–1, where the three cylinders are identified: the inner cylinder defined by the points , the middle cylinder defined by points , and the outer cylinder defined by points . Two slide lines are used in this model: one along the outer edge of the inner cylinder, from node H through node O, and a second along the outer edge of the middle cylinder, from node L through node D. Axisymmetric contact elements for finite sliding (slide line elements) defined along edge of the middle cylinder are associated with the first slide line. Axisymmetric slide line elements defined along edge of the outer cylinder are associated with the second slide line.

The structure is subjected to localized pressurization to initiate contact between the surfaces in the three bodies, and then the two outer cylinders are forced to slide down the cylinder. These loading conditions are defined in two separate steps (pressurization followed by sliding). An additional perturbation step is created to test the *LOAD CASE option.

In the axisymmetric model the inner cylinder is restrained from motion in the z-direction along lines and . In addition, node B is restrained from radial motion. In the first step a pressure of 207 MPa (30 × 103 lb/in2) is applied to edge of the outer cylinder, while nodes L and J are restrained vertically. During the second step the pressure is maintained, and node L is displaced in the negative z-direction by 127 mm (5.0 in), while node J is displaced in the same direction by 114.3 mm (4.5 in).

In the CGAX4 model the same steps and boundary conditions that were applied in the CAX4 model are used. An additional third step is added in which the outermost cylinder is twisted by 0.1 radians about the z-axis while the innermost cylinder is prevented from twisting.

The nonaxisymmetric model is made up of CAXA elements and additional slide line elements at various locations in the -direction. The *ASYMMETRIC-AXISYMMETRIC suboption of the *INTERFACE option is used to define the area of integration for the slide line elements. The ANGLE parameter of the *INTERFACE option is used to define the angular position (measured in degrees) of the slide line elements.

In the CAXA model the boundary conditions that were applied in the axisymmetric model are kept and are extended in the -direction. The loading conditions are the same as the axisymmetric model. Any axisymmetric or nonaxisymmetric loading can be applied to the CAXA model after the second step.

Material:

Solid:

Young's modulus207 GPa (30 × 106 lb/in2)
Poisson's ratio0.3

Coefficients of friction:

Inside edge of middle cylinder0.2
Outer cylinder0.6

Results and discussion

The results from the axisymmetric and three-dimensional models match.

In Step 3 of the CGAX4 analysis (in which the outermost cylinder is twisted by 0.1 radians) the middle cylinder rotates with the outermost cylinder without slipping. Relative slip is prevented by the friction that develops between the slave/master contact pair DSURF/CSURF. However, the middle cylinder does slip with respect to the innermost cylinder. At the end of Step 2 the computed maximum torque that can be transmitted by the slave/master surface contact pair BSURF/ASURF about the z-axis is computed to be 0.2*CTRQ=41000 lb-in, where 0.2 is the friction coefficient for the contact pair. The actual moment due to frictional shear stresses transmitted by the slave/master surface contact pair BSURF/ASURF about the z-axis in Step 3 is 40600 lb-in, which is within 1% of the predicted value from Step 2.

Input files

eia2sssa.inp

Axisymmetric model with CAX4 elements using the contact surface approach, perturbation step with *LOAD CASE.

eia2sssg.inp

Axisymmetric model with CGAX4 elements using the contact surface approach.

eia2ssca.inp

Axisymmetric model with ISL21A and CAX4 elements, perturbation step with *LOAD CASE.

eia2sscn.inp

Nonaxisymmetric model with ISL21A and CAXA41 elements, perturbation step with *LOAD CASE.

Figure

Figure 1.6.13–1 Cylinder sliding model (schematic).