VM257

VM257
Transient Analysis of a Swing with Two Rigid Links and Beam

Overview

Reference: O.A. Bauchau, G. Damilano, and N.J. Theron Numerical Integration of Non-Linear Elastic Multi-Body Systems, International Journal for Numerical Methods in Engineering, Vol. 38, 2727-2751 (1995).
Analysis Type(s): Transient Analysis (ANTYPE = 4)
Element Type(s):
3D Linear Finite Strain Beam (BEAM188)
Multipoint Constraint Element: Rigid Link or Rigid Beam (MPC184-Link/Beam)
3D Target Segment (TARGE170)
Input Listing: vm257.dat

Test Case

The swing shown in Figure 436: Swing Consisting of Two Rigid Links and a Beam with Mid-Span Mass consists of a long aluminum beam of rectangular cross-section (width = 1mm, depth = 5 mm) and a mid-span mass (mass = 0.5 kg). The modulus of elasticity, Poisson’s ratio and density of aluminum are shown in the table below. The mass is rigidly connected to the beam at its mid-span position, labeled C in the figure. The beam is suspended at each end by two rigid links, and is initially at rest in the position as shown in the Figure 436: Swing Consisting of Two Rigid Links and a Beam with Mid-Span Mass. The rigid links impose a kinematic constraint corresponding to fixed distance between points O1 and A, and O2 and E of 0.36 and respectively. The points B and D indicate the quarter and three quarter span points of the beam, respectively. The loading of the system consists of a triangular pulse in the direction applied at the mid-span mass. This pulse starts at time t = 0 s, reaches a peak value of 2N at t =0.128 s and goes back to zero at t = 0.256 s, as shown in Figure 437: Triangular-Pulse Loading.

Figure 436: Swing Consisting of Two Rigid Links and a Beam with Mid-Span Mass

Swing Consisting of Two Rigid Links and a Beam with Mid-Span Mass

Figure 437: Triangular-Pulse Loading

Triangular-Pulse Loading

Material PropertiesGeometric PropertiesLoading
E = 73 GN/m2
ν = 0.3
ρ = 2700 kg/m2
O1A = 0.36 m
AE = 0.72 m
E O2 =
Fx = 0 at time = 0 s
Fx = 2N at time = 0.128 s
Fx = 0 at time = 0.256 s

Analysis Assumptions and Modeling Notes

The system is modeled with four equal length BEAM188 beam elements, two rigid links and a rigid mass. The dynamic response of the system was calculated over a period of 1 s using HHT method with 30% numerical damping and auto time stepping turned on with a minimum of 1000 time steps.

The system was solved twice. In the first case, the rigid links were modeled using MPC184 rigid links and in the second case the rigid links were modeled as rigid bodies using TARGE170 elements. Similar results were obtained in both analyses.

Figure 438: Predicted Time Histories for Displacement Components of Point B shows the predicted time histories for the and direction displacement components of point B and Figure 439: Calculated Time History of Axial Force at Point B shows the calculated time histories of the axial force at point B. These figures should be compared to Figures 15 and 16 respectively in the reference.

Figure 438: Predicted Time Histories for Displacement Components of Point B

Predicted Time Histories for Displacement Components of Point B

Figure 439: Calculated Time History of Axial Force at Point B

Calculated Time History of Axial Force at Point B

Results Comparison

 Target Mechanical APDL Ratio
Results using MPC184 rigid links
TIME (sec)0.64100.64000.998
DISP-UY (m)0.28000.28071.003
DISP-UX (m)0.07500.07831.043
FORCE-FX (N)112.7000107.20690.951
Results using TARGE170 rigid links
TIME (sec)0.64100.64000.998
DISP-UY (m)0.28000.28071.003
DISP-UX (m)0.07500.07831.043
FORCE-FX (N)112.7000107.19040.951