Next Article in Journal
Effect of the Addition of Agribusiness and Industrial Wastes as a Partial Substitution of Portland Cement for the Carbonation of Mortars
Next Article in Special Issue
Dynamic Stability of Tensegrity Structures—Part I: The Time-Independent External Load
Previous Article in Journal
Polymeric Micelles: A Promising Pathway for Dermal Drug Delivery
Previous Article in Special Issue
Influence of Vehicle Number on the Dynamic Characteristics of High-Speed Train-CRTS III Slab Track-Subgrade Coupled System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Dynamic Modeling and Analysis of an L-Shaped Multi-Beam Jointed Structure with Tip Mass

1
School of Electromechanical and Automotive Engineering, Yantai University, Yantai 264005, China
2
State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
3
School of Astronautics, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(23), 7279; https://doi.org/10.3390/ma14237279
Submission received: 26 September 2021 / Revised: 11 November 2021 / Accepted: 25 November 2021 / Published: 28 November 2021
(This article belongs to the Special Issue Nonlinear Dynamics and Applications to Structural Engineering)

Abstract

:
A dynamic model of an L-shaped multi-beam joint structure is presented to investigate the nonlinear dynamic behavior of the system. Firstly, the nonlinear partial differential equations (PDEs) of motion for the beams, the governing equations of the tip mass, and their matching conditions and boundary conditions are obtained. The natural frequencies and the global mode shapes of the linearized model of the system are determined, and the orthogonality relations of the global mode shapes are established. Then, the global mode shapes and their orthogonality relations are used to derive a set of nonlinear ordinary differential equations (ODEs) that govern the motion of the L-shaped multi-beam jointed structure. The accuracy of the model is verified by the comparison of the natural frequencies solved by the frequency equation and the ANSYS. Based on the nonlinear ODEs obtained in this model, the dynamic responses are worked out to investigate the effect of the tip mass and the joint on the nonlinear dynamic characteristic of the system. The results show that the inertia of the tip mass and the nonlinear stiffness of the joints have a great influence on the nonlinear response of the system.

1. Introduction

A multi-beam jointed structure composed of multiple beams and flexible joints is widely used as a component on spacecrafts in the field of aerospace engineering, such as in solar arrays and large antenna support structures. The nonlinearities of the beams and the joints in the structure may have a significant effect on the dynamic characteristics of the entire system. Therefore, an in-depth study of the influence of the nonlinearities on dynamics of the system is necessary for structural design and vibration control.
As a typical multi-beam structure, the research on the dynamics of L-shaped beam structures has received extensive attention in the past few decades [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15]. Because of the inertial effect of the large motion, and the geometry of motion, the L-shaped beam structure is a typical quadratic nonlinearity system. Haddow et al. [4] presented a nonlinear dynamic model of an L-shaped beam structure with two-degree-freedom and considered the effect of the axial motion of the beam. Nayfeh et al. [5,6,7] investigated the nonlinear motion of the L-shaped beam structure with lumped masses, and experimentally confirmed that the structure exhibits a chaotic response when subjected to small excitation with two-to-one internal resonance. Warminski et al. [8] formulated a detailed derivation of the nonlinear PDEs of motion for the L-shaped beam structure with different flexibilities in two orthogonal directions. Cao et al. [9] conducted theoretical and experimental studies on the periodic and chaotic motion in the L-shaped beam structure under the one-to-one internal resonance. Onozato et al. [10] investigated the chaotic responses of a post-buckled L-shaped beam with an axial constraint and examined the contribution ratio of the vibration mode to the chaotic responses by using the proper orthogonal decomposition. Taking into account the inextensionality conditions and the rotary inertia of the beam, Georgiades [11,12] derived the dynamic equations of both in-plane and out-of-plane motions for the L-shaped beam structure. Yu et al. [14,15] made a further study of the global bifurcation and the multi-pulse chaotic motion of the L-shaped beam structure by using the energy-phase method and the generalized Melnikov method.
Concerning the multi-beam structure connected with joints, there have many studies on the effect of the joint on the nonlinear dynamic behavior of the system, which mainly focuses on flexible-joint manipulators [16,17,18,19,20,21,22] and deployable structures [23,24,25,26,27,28,29]. For the flexible-joint manipulator, the joint flexibility mainly comes from the harmonic reducer. Spong et al. [16] were the first to simplify the elastic joint as a torsional spring in the dynamic modeling of the flexible-joint manipulator. By using the Euler–Lagrange principle and the assumed mode method, Subudi and Morris [17] developed a dynamic modeling approach for a manipulator with any number of flexible links and flexible joints so to derive a closed-form dynamic model suitable for vibration control. Vakil et al. [18] investigated the influence of the joint flexibility on the free vibration of the manipulator and established an upper limit for the joint stiffness so to distinguish whether the joint is flexible or rigid for the manipulator. Meng et al. [19] indicated that an appropriate stiffness of the joint can reduce the overall vibration due to the coupling between the link and joint flexibilities, which can be used as a new way to explore the vibration suppression of flexible-joint manipulators. For the manipulator joint with the planetary gear train transmission, Yang et al. [20,21] studied the effect of the time-varying stiffness, the gear tooth profile error, and the backlash clearances in the joint on the dynamics of the manipulator.
In a study of the deployable structure with joints, Moon and Li [23] experimentally studied the dynamic behavior of the pin-jointed space truss structure. The experimental results showed that very small gaps in the joints led to the chaotic motion of the structures, and that the level of chaos can be reduced by increasing the compressive load of the truss structures. Folkman et al. [24] investigated the influence of pinned joints on the damping of the truss structure with an experimental test. Bowden and Dugundji [25] proposed a geometric joint participation factor to quantify the contribution of the damping in the joints to the modal damping of the whole structure. Zhang et al. [26] presented a damping ratio formulation to investigate the effects of pressure, clearance, and dynamic parameters on the damping of joints. Guo et al. [27] systematically investigated the effects of the stiffness, the damping location, and the clearance of the joint on the dynamics of the jointed deployable structure by the finite element model and experimental method. Recently, Wei et al. [28,29] proposed a dynamic modeling approach to derive a reduced-order analytical model for a multi-beam structure with nonlinear joints. Based on this model, they studied the influence of the nonlinear stiffness and the damping of the joint on the attitude and the position of the spacecraft during maneuvering.
Moreover, the effect of the tip mass on the nonlinear dynamics of the beam structure has been investigated in Refs [30,31,32,33,34]. Esmailzadeh and Jalili [30] focused on the periodic response of the cantilever beam with the tip mass under the axial base excitation, and the results indicated that the increase of the tip mass would reduce the stable periodic region. Yabuno et al. [31] theoretically analyzed the effect of the tip mass on the bifurcation points of the frequency response curve for a parametrically excited cantilever beam. Eftekhari et al. [32] conducted a study on the primary resonance of the composite cantilever beam with the tip mass subjected to base excitation both in the flapwise and the chordwise. The results showed that the addition of the tip mass can change the dynamic stability of the system. Meesala and Hajj [34] presented a sensitivity analysis of a parametrically excited cantilever beam with the tip mass. They demonstrated that very small variations in either the stiffness or the tip mass can alter the type of bifurcation, which can potentially be applied for using parametrically excited micro-cantilever beams as sensing devices.
Although a lot of work has been done in the modeling and the analysis of the L-shaped beam structure and the multi-beam jointed structure, the research on the dynamics of the L-shaped multi-beam jointed structure with the tip mass is rarely seen. This paper is devoted to obtaining a low-dimensional nonlinear dynamics model, which is both convenient for the use of the analytical method to investigate the nonlinear dynamic behavior and suitable for real-time control. Based on the nonlinear ODEs in this model, the dynamic responses of the system can be performed to investigate the effect of the tip mass and joints on the nonlinear dynamics of the L-shaped multi-beam jointed structure, which can provide support for the analysis of the payload and the nonlinear joints of deployable structures in aerospace engineering.

2. Dynamical Model

2.1. Governing Equations of Motion

Consider the planar motion of an L-shaped multi-beam jointed structure that consists of a tip mass and two beams connected with nonlinear torsional joints, as shown in Figure 1. The horizontal beam with the length L 1 is connected to the base by the first nonlinear torsional joint S 1 . Similarly, the vertical beam with the length L 2 is connected to the horizontal beam by the second nonlinear torsional joint S 2 . A tip mass with mass moment of inertia J and mass m is fixed to the end of the vertical beam. Furthermore, the L-shaped multi-beam jointed structure is subjected to the horizontal acceleration X ¨ b ( t ) = X 0 c o s Ω t and the vertical acceleration Y ¨ b ( t ) = Y 0 cos Ω t , where X 0 and Y 0 are the peak amplitudes of the horizontal and the vertical, respectively, and Ω is the excitation frequency. Assume that the two beams have the same dimensions as the cross-section and the same material properties. The mass per unit length, elastic modulus, and area moment of inertia of the beams are denoted by ρ , E , and I , respectively. The torsional deformation of the i - t h joint S i , i = 1 , 2 is described by θ i . Let ( x 1 , y 1 ) and ( x 2 , y 2 ) be the coordinates of the horizontal beam and the vertical beam, with the origin located at joint S 1 and joint S 2 , respectively. Let ( x m , y m ) be the coordinate of the tip mass, with the origin located at the center of the tip mass. x m and y m are the translational displacements of the tip mass in the horizontal and the vertical directions, respectively. θ m is used to describe the rotational displacement of the tip mass.
The deformation of arbitrary points P 10 and P 20 on the horizontal beam and the vertical beam, respectively, are depicted in Figure 1. After the deformation, the positions of P 10 and P 20 move to positions P 1 and P 2 , respectively. The location vectors of positions P 1 and P 2 are given by
r 1 = [ u 1 ( x , t ) , v 1 ( x , t ) ] ,   r 2 = [ v 2 ( x , t ) , v 1 ( L 1 , t ) + u 2 ( x , t ) ]
where u 1 ( x , t ) and v 1 ( x , t ) are the axial and transverse displacements of the horizontal beam, respectively, v 1 ( L 1 , t ) is the transverse displacement of the horizontal beam tip, and u 2 ( x , t ) and v 2 ( x , t ) are the axial and transverse displacements of the vertical beam, respectively. The axial displacements u 1 ( x , t ) and u 2 ( x , t ) are generally much smaller than the transverse displacements v 1 ( x , t ) and v 2 ( x , t ) . Then, it is reasonable that the axial displacements u 1 ( x , t ) and u 2 ( x , t ) are neglected. Then, the location vectors of the positions P 1 and P 2 can be taken as
r 1 = [ 0 , v 1 ( x , t ) ] ,   r 2 = [ v 2 ( x , t ) , v 1 ( L 1 , t ) ]
Considering the effect of the axial force, the nonlinear transverse equation that governs the motion of the horizontal beam can be written as [34]
ρ v ¨ 1 + c v v ˙ 1 + E I v 1 + η I v ˙ 1 = E I [ v 1 ( v 1 v 1 ) ] ρ 2 [ v 1 L 1 x 0 x 2 t 2 ( v 1 2 ) d x d x ] + m t 2 v 1 0 L 1 2 t 2 ( v 1 2 ) d x + ρ X ¨ b [ v 1 + ( x L 1 ) v 1 ] N 1 v 1 + F y
where an overdot denotes partial differentiation with respect to time t , a prime denotes partial differentiation with respect to x , c v and η are the external and internal damping of the beam, respectively, N 1 = m t X ¨ b + 0 L 2 ρ v ¨ 2 d x + m x ¨ m is the axial force of the horizontal beam at the joint S 2 , F y = ρ Y ¨ b is the transverse force of the horizontal beam per unit length, and m t = ρ L 2 + m is the tip mass attached to the end of the horizontal beam along the x2 axis direction. Assume that the external and internal damping of the beam are proportional to its mass and stiffness, respectively, and they are expressed as
c v = a ρ , η = b E
where a and b are the proportionality constant.
For the tip mass and the vertical beam, the equation that governs the motion in the x2 axis direction is
( ρ L 2 + m ) y ¨ m + c r y ˙ m E I v 1 ( L 1 , t ) = ( ρ L 2 + m ) Y ¨ b
where c r is the damping coefficient of the horizontal beam and the tip mass translation in the x2 axis direction, and it is expressed as
c r = a ( ρ L 2 + m )
The nonlinear transverse equation that governs the motion of the vertical beam is
ρ v ¨ 2 + c v v ˙ 2 + E I v 2   + η I v ˙ 2   = E I [ v 2 ( v 2 v 2 ) ] ρ 2 [ v 2 L 2 x 0 x 2 t 2 ( v 2 2 ) d x d x ] + m 2 v 2 0 L 2 2 t 2 ( v 2 2 ) d x + ρ ( Y ¨ b + y ¨ m ) [ v 2 + ( x L 2 ) v 2 ] N 2 v 2 + F x
where N 2 = m ( Y ¨ b + y ¨ m ) is the external force of the vertical beam on the tip mass and F x = ρ X ¨ b is the external transverse force on the vertical beam.
As for the natural matching conditions shown in Figure 2a, and according to the shear force and the bending moment of the tip mass, the equations of motion for the tip mass rotation and translation in the horizontal direction are
J θ ¨ m + c J θ ˙ m + E I v 2 ( L 2 , t ) d 2 E I v 2 ( L 2 , t ) = 0
m x ¨ m + c m x ˙ m E I v 2 ( L 2 , t ) = m X ¨ b
where d is the length of the tip mass, c J and c m are the damping of the tip mass rotational and translation in the horizontal direction, respectively, and they are expressed as
c J = a J ,   c m = a m
The nonlinear torsional joint is described as a single-degree-of-freedom nonlinear spring massless system with a damper, as shown in Figure 3. The nonlinear transmitted torque formulation of the i - t h joint can be expressed as
M i T = μ θ ˙ i + k L θ i + k N θ i 3 ,   i = 1 , 2
where k L , k N , and μ are the linear stiffness, the cubic stiffness, and the damping of the nonlinear torsional joint, respectively.
The boundary conditions at the joint S 1 are
v 1 ( 0 , t ) = 0 ,    v 1 ( 0 , t ) = θ 1 ,    E I v 1 ( 0 , t ) = M 1 T
As shown in Figure 4a, the geometric match conditions at the joint S 2 are
v 1 ( L 1 , t ) + θ 2 = v 2 ( 0 , t ) ,    v 2 ( 0 , t ) = 0
As shown in Figure 4b, the moment matching boundary conditions at the joint S 2 are
E I v 1 ( L 1 , t ) = M 2 T = E I v 2 ( 0 , t )
The geometric match conditions of the tip mass are
v 1 ( L 1 , t ) = y m ,    v 2 ( L 2 , t ) d 2 θ m = x m ,    v 2 ( L 2 , t ) = θ m

2.2. Determination of Natural Frequencies and Global Mode Shapes

In order to solve the eigenvalue problem of the system, it is assumed that the displacements of the L-shaped multi-beam jointed structure are separable in space and time. Let
v i ( x , t ) = φ i ( x ) e j ω t , θ i = Θ i e j ω t , x m = X m e j ω t ,    y m = Y m e j ω t ,    θ m = Θ m e j ω t ,    i = 1 , 2
where ω is an unknown constant corresponding to the natural frequency of the system. Substituting the separable solutions given in Equation (16) into Equations (3), (5), (7)–(9) without damping, nonlinear terms, and external force yields
ω 2 ρ φ 1 ( x ) E I φ 1 ( x ) = 0 ,
( ρ L 2 + m ) ω 2 Y m + E I φ 1 ( L 1 ) = 0 ,
ω 2 ρ φ 2 ( x ) E I φ 2 ( x ) = 0 ,
J ω 2 Θ m E I φ 2 ( L 2 ) + d 2 E I φ 2 ( L 2 ) = 0 ,
m ω 2 X m + E I φ 2 ( L 2 ) = 0 .
The boundary conditions at the joints S 1 and match conditions at the joint S 2 are reduced to
{ φ 1 ( 0 ) = 0 ,    φ 1 ( 0 ) = Θ 1 ,    E I φ 1 ( 0 ) = k L Θ 1 , φ 2 ( 0 ) = 0 ,    φ 1 ( L 1 ) + Θ 2 = φ 2    ( 0 ) , E I φ 1 ( L 1 ) = k L Θ 2 ,    E I φ 2 ( 0 ) = k L Θ 2 .
The geometric match conditions of the tip mass are reduced to
φ 1 ( L 1 ) = Y m ,    φ 2 ( L 2 ) = Θ m ,    φ 2 ( L 2 ) d 2 Θ m = X m .
The solutions of Equations (17) and (19) can be written as
{ φ 1 ( x ) = A 1 cos ( β x ) + B 1 sin ( β x ) + C 1 cosh ( β x ) + D 1 sinh ( β x ) ,    x [ 0 , L 1 ] , φ 2 ( x ) = A 2 cos ( β x ) + B 2 sin ( β x ) + C 2 cosh ( β x ) + D 2 sinh ( β x ) ,    x [ 0 , L 2 ] .
where β = ( ρ ω 2 E I ) 1 / 4 . Let
ψ = [ A 1 B 1 C 1 D 1 Θ 1 A 2 B 2 C 2 D 2 Θ 2 X m Y m Θ m ] T .
Substituting Equation (24) into the boundary and matching conditions in (22) and (23), and the dynamic equations of the tip mass in (20) and (21), yields
H ( ω ) ψ = 0 ,
where the matrix H ( ω ) R 13 × 13 is given in Appendix A.
The natural frequencies of the L-shaped multi-beam jointed structure are denoted in ascending order by ω 1 , ω 2 , , which are the positive roots of the frequency equation det ( H ( ω ) ) = 0 . Once the natural frequency ω s is obtained, the eigenvector ψ ( s ) can be obtained by Equation (26). Then, the s - t h global mode shapes for the L-shaped multi-beam jointed structure can be determined by Equation (24).

2.3. Orthogonality of the Global Mode Shapes

The global mode shapes associated with the two distinct eigenvalues ω r and ω s are denoted by ϕ r ( x ) and ϕ s ( x ) , respectively, where
ϕ r ( x ) = [ φ 1 r Θ 1 r φ 2 r Θ 2 r X m r Y m r Θ m r ] T , r = 1 , 2 , .
By Equations (17)–(21), one has
E I φ 1 r ( x ) = ω r 2 ρ φ 1 r ( x ) ,
E I φ 2 r ( L 1 ) = ω r 2 ( ρ L 2 + m ) Y m r ,
E I φ 2 r ( x ) = ω r 2 ρ φ 2 r ( x ) ,
E I φ 2 r ( L 2 ) - d 2 E I φ 2 r ( L 2 ) = ω r 2 J Θ m r ,
- E I φ 2 r ( L 2 ) = ω r 2 m X m r .
Then, multiply Equations (28) and (30) by φ 1 s and φ 2 s , integrate the resulting equations over the domain 0 x L 1 and 0 x L 2 for the horizontal and the vertical beams, respectively, and add the resulting equations to get
0 L 1 E I φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 E I φ 2 r ( x ) φ 2 s ( x ) d x = ω r 2 ( 0 L 1 ρ φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 ρ φ 2 r ( x ) φ 2 s ( x ) d x ) .
Multiply Equations (29), (31), and (32) by Y m s , Θ m s , and X m s , respectively, and add the resulting equations to get
[ E I φ 2 r ( L 2 ) - d 2 E I φ 2 r ( L 2 ) ] Θ m s - E I φ 1 r ( L 1 ) Y m s - E I φ 2 r ( L 2 ) X m s = ω r 2 [ ( ρ L 2 + m ) Y m r Y m s + m X m r X m s + J Θ m r Θ m s ] .
Integrating Equation (33) by parts, adding Equation (34), and using the boundary and matching conditions in Equations (22) and (23) and the dynamic equations in (18), (20), and (21) yields
0 L 1 E I φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 E I φ 2 r ( x ) φ 2 s ( x ) d x = 0 L 1 E I φ 1 r ( x ) φ 1 S ( x ) d x + 0 L 2 E I φ 2 r ( x ) φ 2 s ( x ) d x + k L ( Θ 1 r Θ 1 s + Θ 2 r Θ 2 s ) ω r 2 [ ( ρ L 2 + m ) Y m r Y m s + m X m r X m s + J Θ m r Θ m s ] .
Substituting Equation (35) into the left-hand side of Equation (33) yields
0 L 1 E I φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 E I φ 2 r ( x ) φ 2 s ( x ) d x + k L ( Θ 1 r Θ 1 s + Θ 2 r Θ 2 s ) = ω r 2 [ 0 L 1 ρ φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 ρ φ 2 r ( x ) φ 2 s ( x ) d x + ( ρ L 2 + m ) Y m r Y m s + m X m r X m s + J Θ m r Θ m s ] .
Exchanging the superscripts s and r in Equation (36) yields
0 L 1 E I φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 E I φ 2 r ( x ) φ 2 s ( x ) d x + k L ( Θ 1 r Θ 1 s + Θ 2 r Θ 2 s ) = ω s 2 [ 0 L 1 ρ φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 ρ φ 2 r ( x ) φ 2 s ( x ) d x + ( ρ L 2 + m ) Y m r Y m s + m X m r X m s + J Θ m r Θ m s ] .
Subtracting Equation (37) from Equation (36) yields
( ω r 2 - ω s 2 ) [ 0 L 1 ρ φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 ρ φ 2 r ( x ) φ 2 s ( x ) d x ] + ( ω r 2 - ω s 2 ) [ ( ρ L 2 + m ) Y m r Y m s + m X m r X m s + J Θ m r Θ m s ] = 0 .
From Equation (38), the first orthogonality relation can be obtained
0 L 1 ρ φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 ρ φ 2 r ( x ) φ 2 s ( x ) d x + ( ρ L 2 + m ) Y m r Y m s + m X m r X m s + J Θ m r Θ m s = M s δ r s ,
where M s is a positive constant and δ r s is the Kronecker delta. Using Equations (36) and (39), the second orthogonality relation can be obtained
0 L 1 E I φ 1 r ( x ) φ 1 s ( x ) d x + 0 L 2 E I φ 2 r ( x ) φ 2 s ( x ) d x + k L ( Θ 1 r Θ 1 s + Θ 2 r Θ 2 s ) = K s δ r s .
where K s is a positive constant.

2.4. Dynamic Model with Multi-DOF

Applying the Galerkin procedure to Equations (3), (5), and (7)–(9), the nonlinear ODEs are derived for the system. By using the first n global mode shapes, the displacements of the L-shaped multi-beam jointed structure can be written as
v 1 ( x , t ) = j = 1 n φ 1 j ( x ) q j ( t ) ,    v 2 ( x , t ) = j = 1 n φ 2 j ( x ) q j ( t ) ,    θ 1 = j = 1 n Θ 1 j q j ( t ) , θ 2 = j = 1 n Θ 2 j q j ( t ) ,    x m = j = 1 n X m j q j ( t ) ,    y m = j = 1 n Y m j q j ( t ) ,    θ m = j = 1 n Θ m j q j ( t ) .
where q j ( t ) is the j - t h modal coordinate of the system.
Substituting Equation (41) into Equations (3), (5), (7)–(9), and (11) yields
j = 1 n ρ φ 1 j q ¨ j + j = 1 n c v φ 1 j q ˙ j + j = 1 n E I φ 1 j q j + j = 1 n η I φ 1 j q ˙ j = j = 1 n ρ X ¨ b [ φ 1 j + ( x L 1 ) φ 1 j ] q j - j = 1 n k = 1 n r = 1 n E I [ φ 1 j ( φ 1 k φ 1 r ) ¢ ] ¢ q j q k q r - j = 1 n k = 1 n r = 1 n ρ { φ 1 j q j L 1 x 0 x [ φ 1 k φ 1 r ( q ˙ k q ˙ r + q ˙ k q ¨ r ) ] d x d x } + j = 1 n k = 1 n r = 1 n ( ρ L 2 + m ) φ 1 j q j 0 L 1 [ φ 1 k φ 1 r ( q ˙ k q ˙ r + q ˙ k q ¨ r ) ] d x - j = 1 n ( ρ L 2 + m ) X ¨ b φ 1 j q j + j = 1 n k = 1 n ( 0 L 2 ρ φ 2 k d x + m X m k ) q ¨ k φ 1 j q j + ρ Y ¨ b ,
j = 1 n ( ρ L 2 + m ) Y m j q ¨ j + j = 1 n c r Y m j q ˙ j - j = 1 n E I φ 1 j q j = ( ρ L 2 + m ) Y ¨ b ,
j = 1 n ρ φ 2 j q ¨ j + j = 1 n c v φ 2 j q ˙ j + j = 1 n E I φ 2 j q j + j = 1 n η I φ 2 j q ˙ j = - j = 1 n k = 1 n r = 1 n E I [ φ 2 j ( φ 2 k φ 2 r ) ¢ ] q j q k q r + j = 1 n k = 1 n ρ Y m k [ φ 2 j + ( x - L 2 ) φ 2 j ] q j q ¨ k - j = 1 n k = 1 n r = 1 n ρ { φ 2 j q j L 2 x 0 x [ φ 2 k φ 2 r ( q ˙ k q ˙ r + q ˙ k q ¨ r ) ] d x d x } + j = 1 n k = 1 n r = 1 n m φ 2 j q j 0 L 2 [ φ 2 k φ 2 r ( q ˙ k q ˙ r + q ˙ k q ¨ r ) ] d x - j = 1 n k = 1 n m Y m k φ 2 j q j q ¨ k + j = 1 n ρ Y ¨ b [ φ 2 j + ( x - L 2 ) φ 2 j ] q j - j = 1 n m Y ¨ b φ 2 j q j + ρ X ¨ b ,
j = 1 n J Θ m j q ¨ j + j = 1 n c J Θ m j q ˙ j + j = 1 n E I φ 2 j ( L 2 ) q j - d 2 j = 1 n E I φ 2 j ( L 2 ) q j = 0 ,
j = 1 n m X m j q ¨ j + j = 1 n c m X m j q ˙ j - j = 1 n E I φ 2 j q j = m X ¨ b ,
M i T = j = 1 n μ Θ i j q ˙ j + j = 1 n k L Θ i j q j + j = 1 n k = 1 n r = 1 n k N Θ i j Θ i k Θ i r q j q k q r ,    i = 1 , 2 .
Multiplying Equations (42) and (44) by φ 1 s and φ 2 s , integrating the resulting equations over the domain 0 x L 1 and 0 x L 2 for the horizontal beam and the vertical beam, multiplying Equations (43), (45), and (46) by Y m s , Θ m s , and X m s , respectively, and adding all of the resulting equations and using the matching and boundary conditions in Equations (12)–(15) and the orthogonality relations in Equations (39) and (40), we have
q ¨ s + ω s 2 q s + 2 α s ω s q ˙ s + j = 1 n γ j s q ˙ j + j = 1 n a j s q j + j = 1 n k = 1 n b j k s q j q ¨ k + j = 1 n k = 1 n r = 1 n c j k r s q j q k q r + j = 1 n k = 1 n r = 1 n d j k r s ( q j q k q ¨ r + q j q ˙ k q ˙ r ) = f s ( t ) ,    s = 1 , 2 , , n .
where α s and γ j s are the damping of the beams and joints, respectively. The terms a j s are resulted from the axial force of the beams caused by the foundation motion, the terms b j k s are resulted from the axial forces of the beams caused by the elastic motion of the beams, the terms c j k r s represent the geometric nonlinearities of the beams and the nonlinear stiffness of the joints, and the terms d j k r s represent the inertial nonlinearities of the beams and the tip mass. The relevant terms in Equation (48) are given in Appendix B.

3. Results and Discussion

In this section, a comparison of the natural frequencies obtained from the frequency equation and the finite element software ANSYS (Version 16.0, ANSYS Inc, Canonsburg, PA, USA) is presented so to verify the accuracy of the model, and the first four global mode shapes of the L-shaped multi-beam structure are also given. Then, the dynamic responses of the system with various parameters are worked out to investigate the effect of the tip mass and the joints on the nonlinear vibration of the L-shaped multi-beam structure under base excitation. The physical parameters of the L-shaped multi-beam structure are listed in Table 1.

3.1. Mode Analysis

To verify the accuracy of the dynamic model obtained in this paper, the first four natural frequencies of the L-shaped multi-beam structure with various tip mass parameters are solved by the frequency equation det ( H ( ω ) ) = 0 and compared with those obtained from the finite element software ANSYS, as shown in Figure 5. For the finite element model, the beam, the flexible joint, and the tip mass are modeled by the BEAM element, the COMBIN 14 element, and the MASS element, respectively. Each beam is divided into 300 elements. The joint linear torsional stiffness is selected as k L = E I / L , and the moment of inertia of the tip mass increases in proportion to its mass. The pink symbol ○ represents the results from the finite element software ANSYS.
As shown in Figure 5, the results, solved by the frequency equation from m = 0   kg to m = 1   kg are in good agreement with those by the finite element software ANSYS. This demonstrates the correctness of the dynamic model in this paper. The first four natural frequencies and the global mode shapes of the L-shaped multi-beam jointed structure with m = 0.5   kg are illustrated in Figure 6.

3.2. Nonlinear Responses

As the tip mass is fixed to the end of the vertical beam, it will cause a change of axial force, geometric nonlinearity, and inertial nonlinearity, and then it will affect the nonlinear dynamic response of the system. To investigate the effect of the tip mass on the nonlinear response of the system, the frequency response curves of the tip mass with various parameters are calculated by the MATLAB (Version 2017a, MathWorks, Natick, MA, USA) ode function derived from the Runge–Kutta method, as shown in Figure 7 and Figure 8. In these figures, the roman numerals I represent the results of considering only the geometric nonlinearity of the beams, which indicates that the terms a j s , b j k s and d j k r s are missed in the case I. The roman numerals II represent the results of considering only the geometric nonlinearity of the beams, the inertial nonlinearity, and the axial force of the beams caused by the foundation motion, which indicates that the terms b j k s are missed in this case II. The roman numerals III represent the results from the Equation (48), which includes all nonlinear terms. In Figure 7 and Figure 8, the damping and the nonlinear stiffness of the joints are selected as μ = 0 , k N = 0 .
From Figure 6, it is observed that the second natural frequency is much smaller than the third natural frequency. Then, one can select only the first two modes in calculating the dynamic responses of the system when Ω varies near the first natural frequency of the L-shaped multi-beam jointed structure. The results with increasing and decreasing Ω are obtained as follows: the initial conditions for the first point are assumed to be zero, and some steady-state response data at a current Ω are used as the initial conditions for the next Ω .
The spectrum plots of the steady-state response of the tip mass in the cases of m = 0 , J = 0 and m = 1 , J = 0.02 are shown in Figure 9. From this figure, it can be observed that the terms d j k r s resulted from the inertial nonlinearity, the terms c j k r s resulted from the geometric nonlinearity of the beams, which led to the generation of the third harmonics, and the terms b j k s resulted from the axial force of the beams caused by the elastic motion, which led to the generation of the second harmonics. It is important to note that the terms b j k s led to a zero frequency in the spectrum.
To investigate the effect of the nonlinear stiffness of the joint and the nonlinearities from the beams on the nonlinear responses of the system, the frequency response curves of the system with various linear stiffnesses k L and nonlinear stiffnesses k N are given in Figure 10 and Figure 11. The tip mass is selected as m = 1 , J = 0.02 . The joint linear stiffness is selected as k L = 5 E I / L and k L = E I / L in Figure 10 and Figure 11, respectively. In these figures, the roman numerals IV represent the results of considering only the nonlinearities of the joints. The roman numerals V represent the results of considering the nonlinearities of the joints and the nonlinearities from the beams, Equation (48), and including all nonlinear terms.
To investigate the effect of the nonlinear stiffness of the joint and the nonlinearities from the beams on the nonlinear responses of the system, the frequency response curves of the system with various linear stiffnesses k L and nonlinear stiffnesses k N are given in Figure 10 and Figure 11. The tip mass is selected as m = 1 , J = 0.02 . The joint linear stiffness is selected as k L = 5 E I / L and k L = E I / L in Figure 10 and Figure 11, respectively. In these figures, the roman numerals IV represent the results of considering only the nonlinearities of the joints. The roman numerals V represent the results of considering the nonlinearities of the joints and the nonlinearities from the beams, Equation (48), and including all nonlinear terms.
In the case of m = 1 , J = 0.02 , the frequency response of the system shows soft characteristics when the joint nonlinear stiffness is selected as k N = 0 . In Figure 10 and Figure 11, it is clearly seen that the frequency response of the system shows hard characteristics due to the addition of the positive cubic nonlinear stiffness of the joint. As the positive cubic nonlinear stiffness of the joint increases, the response amplitude of the system decreases gradually, and the jump occurs at a higher frequency. By comparing the results from two cases, IV and V, it can be observed that the differences obtained from the two cases become smaller and smaller. This shows that the cubic nonlinear stiffness of the joints gradually becomes the main source of the system’s nonlinearity with the increase of the positive cubic nonlinear stiffness of the joints. By comparing Figure 10 with Figure 11, when the linear stiffness of the joint is relatively small, the cubic nonlinear stiffness of the joint has a great influence on the frequency response of the system because of its large joint torsion deformation, and the influence of other nonlinearities is small or even negligible. For instance, the very small differences in the frequency responses obtained from the two cases show that the nonlinear dynamic characteristics of the system are dominated by the nonlinear stiffness of the joint in Figure 11d.
The frequency response curves of the structure with various joint dampings μ are shown in Figure 12. From this figure, it is clearly seen that, in the case without joint damping, the jump phenomenon is very obvious. With the increase of joint damping, the vibration amplitude of the tip mass decreases, and the jump phenomenon gradually diminishes until it disappears.

4. Conclusions

In this paper, by taking into account the axial force and the geometric nonlinearities of the beams, the nonlinear stiffness of the joints, and the inertial nonlinearities, a nonlinear dynamic model of an L-shaped multi-beam jointed structure has been developed by using the dynamic modeling approach proposed in Ref [28]. A comparison of the natural frequencies solved by the frequency equation and the ANSYS has been given to verify the accuracy of the model in this paper. Based on the nonlinear ODEs obtained in this model, the dynamic responses of the L-shaped multi-beam jointed structure have been presented to investigate the effect of the tip mass and the joints on the nonlinear dynamics of the system. Two conclusions are summarized as follows:
(1) In the case without the nonlinear stiffness of the joints, as the inertia of the tip mass increases, the geometric nonlinearity of the beam causes the frequency responses to change from showing the hard characteristics to showing the soft characteristics. The large inertial of the tip mass results in the frequency responses of the system without the joint nonlinearities to be dominated by the geometric nonlinearities of the beams. In addition, the axial force of the beams caused by the elastic motion leads to a zero frequency in the spectrum.
(2) With the addition of the positive cubic nonlinear stiffness of the joint, the response amplitude of the system decreases, and the jump occurs at a higher frequency. When the torsional deformation of the joint is relatively large, the nonlinear stiffness of the joints has a great influence on the frequency response of the system, and the influence of other nonlinearities is small or even negligible.
In summary, the nonlinear dynamic model derived in this paper can guarantee accuracy while at the same time greatly reduce its degree of freedom. This makes the model not only convenient for the analytical method to better understand the nonlinear dynamic behavior of the system, but also to reduce the calculation of the dynamic analysis in structural design. More importantly, this model is very suitable for real-time control due to its fewer degrees of freedom. All of these are of great significance to nonlinear dynamic analysis and the active vibration control of the multi-beam jointed structure.

Author Contributions

Conceptualization, J.W. (Jin Wei), D.J., and D.C.; methodology, J.W. (Jin Wei) and D.C.; software, J.W. (Jinjie Wang) and M.L.; investigation, writing—original draft, J.W.; writing—review and editing, J.W. (Jin Wei), T.Y., and M.L.; visualization, M.L. and J.W. (Jinjie Wang); supervision, T.Y.; project administration, J.W. (Jin Wei), D.C., and T.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 12002298 and 11732005), the China Postdoctoral Science Foundation (Grant No. 2020M681578), the Shandong Provincial Natural Science Foundation, China (Grant No. ZR2020QA038), and the Shandong Key Research and Development Program (Important Technologies Project, Grant NO. 2019JZZY020114).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Entries of the matrix H in Equation (26)
H 1 , 1 = cos ( 0 ) , H 1 , 2 = sin ( 0 ) , H 1 , 3 = cosh ( 0 ) , H 1 , 4 = sinh ( 0 ) , H 2 , 1 = - β sin ( 0 ) , H 2 , 2 = β cos ( 0 ) , H 2 , 3 = β sinh ( 0 ) , H 2 , 4 = β cosh ( 0 ) , H 3 , 1 = - E I β 2 cos ( 0 ) , H 3 , 2 = - E I β 2 sin ( 0 ) , H 3 , 3 = E I β 2 cosh ( 0 ) , H 3 , 4 = E I β 2 sinh ( 0 ) , H 3 , 5 = - k 1 , H 4 , 6 = cos ( 0 ) , H 4 , 7 = sin ( 0 ) , H 4 , 8 = cosh ( 0 ) , H 4 , 9 = sinh ( 0 ) , H 5 , 1 = - β sin ( β L 1 ) , H 5 , 2 = β cos ( β L 1 ) , H 5 , 3 = β sinh ( β L 1 ) , H 5 , 4 = β cosh ( β L 1 ) , H 5 , 6 = β sin ( 0 ) , H 5 , 7 = - β cos ( 0 ) , H 5 , 8 = - β sinh ( 0 ) , H 5 , 9 = - β cosh ( 0 ) , H 5 , 10 = 1 , H 6 , 1 = - E I β 2 cos ( β L 1 ) , H 6 , 2 = - E I β 2 sin ( β L 1 ) , H 6 , 3 = E I β 2 cosh ( β L 1 ) , H 6 , 4 = E I β 2 sinh ( β L 1 ) , H 6 , 5 = - k 1 , H 7 , 5 = - k 1 , H 7 , 6 = - E I β 2 cos ( 0 ) , H 7 , 7 = - E I β 2 sin ( 0 ) , H 7 , 8 = E I β 2 cosh ( 0 ) , H 7 , 9 = E I β 2 sinh ( 0 ) , H 8 , 1 = E I β 3 sin ( β L 1 ) , H 8 , 2 = - E I β 3 cos ( β L 1 ) , H 8 , 3 = E I β 3 sinh ( β L 1 ) , H 8 , 4 = E I β 3 cosh ( β L 1 ) , H 8 , 12 = ( ρ L 2 + m ) ω 2 , H 9 , 1 = cos ( β L 1 ) , H 9 , 2 = sin ( β L 1 ) , H 9 , 3 = cosh ( β L 1 ) , H 9 , 4 = sinh ( β L 1 ) , H 9 , 12 = - 1 , H 10 , 6 = - β sin ( β L 2 ) , H 10 , 7 = β cos ( β L 2 ) , H 10 , 8 = β sinh ( β L 2 ) , H 10 , 9 = β cosh ( β L 2 ) , H 10 , 12 = - 1 , H 11 , 6 = cos ( β L 2 ) , H 11 , 7 = sin ( β L 2 ) , H 11 , 8 = cosh ( β L 2 ) , H 11 , 9 = sinh ( β L 2 ) , H 11 , 10 = - 1 , H 11 , 13 = - d / 2 , H 12 , 6 = E I β 3 sin ( β L 2 ) , H 12 , 7 = - E I β 3 cos ( β L 2 ) , H 12 , 8 = E I β 3 sinh ( β L 2 ) , H 12 , 9 = E I β 3 cosh ( β L 2 ) , H 12 , 10 = m ω 2 , H 13 , 13 = J ω 2 , H 13 , 6 = E I β 2 cos ( β L 2 ) + d E I β 3 sin ( β L 2 ) / 2 , H 13 , 7 = E I β 2 sin ( β L 2 ) - d E I β 3 cos ( β L 2 ) / 2 , H 13 , 8 = - E I β 2 cosh ( β L 2 ) + d E I β 3 sinh ( β L 2 ) / 2 , H 13 , 9 = - E I β 2 sinh ( β L 2 ) + d E I β 3 cosh ( β L 2 ) / 2 ,
where H i , j represents the i - t h row and j - t h columns of the matrix H , and all the other elements of the matrix H are zero.

Appendix B

The relevant terms in Equation (48)
M s = 0 L 1 ρ φ 1 s 2 d x + 0 L 2 ρ φ 2 s 2 d x + ( ρ L 2 + m ) Y m s 2 + m X m s 2 + J Θ m s 2 , K s = 0 L 1 E I φ 1 s 2 d x + 0 L 2 E I φ 2 s 2 d x + k L ( Θ 1 s 2 + Θ 2 s 2 ) , ω s 2 = K s M s ,    α s = 1 2 ( a ω s + b ω s ) ,    γ j s = μ M s ( Θ 1 j Θ 1 s + Θ 2 j Θ 2 s ) ,    a j s = 1 M s { - X ¨ b 0 L 1 φ 1 s [ ρ ( φ 1 j + ( x - L 1 ) φ 1 j ) - ( ρ L 2 + m ) φ 1 j ] d x - Y ¨ b 0 L 2 φ 2 s [ ρ ( φ 2 j + ( x - L 2 ) φ 2 j ) - m φ 2 j ] d x } , b j k s = 1 M s { 0 L 1 φ 1 s φ 1 j [ 0 L 2 ρ φ 2 k d x + m X m k ] d x - 0 L 2 φ 2 s Y m k { ρ [ φ 2 j + ( x - L 2 ) φ 2 j ] - m φ 2 j } d x } , c j k r s = 1 M s { E I 0 L 1 φ 1 s [ φ 1 j ( φ 1 k φ 1 r ) ¢ ] d x + E I 0 L 2 φ 2 s [ φ 2 j ( φ 2 k φ 2 r ) ¢ ] d x + k N ( Θ 1 j Θ 1 k Θ 1 r Θ 1 s + Θ 2 j Θ 2 k Θ 2 r Θ 2 s ) } , d j k r s = 1 M s { 0 L 1 φ 1 s [ ρ ( φ 1 j L 1 x 0 x φ 1 k φ 1 r d x d x ) - ( ρ L 2 + m ) φ 1 j 0 L 1 φ 1 k φ 1 r d x ] d x + 0 L 2 φ 2 s [ ρ ( φ 2 j L 2 x 0 x φ 2 k φ 2 r d x d x ) - m φ 2 j 0 L 2 φ 2 k φ 2 r d x ] d x } , f s ( t ) = 1 M s { Y ¨ b [ 0 L 1 ρ φ 1 s d x + ( ρ L 2 + m ) Y m s ] + X ¨ b [ 0 L 2 ρ φ 2 s d x + m X m s ] } .

References

  1. Cao, Y.; Cao, D.; He, G.; Ge, X.; Hao, Y. Vibration analysis and distributed piezoelectric energy harvester design for the L-shaped beam. Eur. J. Mech. A-Solid 2021, 87, 104214. [Google Scholar] [CrossRef]
  2. Li, H.; Sun, H.; Song, B.; Zhang, D.; Shang, X.; Liu, D. Nonlinear dynamic response of an L-shaped beam-mass piezoelectric energy harvester. J. Sound Vibr. 2021, 499, 116004. [Google Scholar] [CrossRef]
  3. Sharifi-Moghaddam, S.; Tikani, R.; Ziaei-Rad, S. Energy harvesting from nonlinear vibrations of an L-shaped beam using piezoelectric patches. J. Braz. Soc. Mech. Sci. 2021, 43, 1–15. [Google Scholar] [CrossRef]
  4. Haddow, A.G.; Barr, A.D.S.; Mook, D.T. Theoretical and experimental study of modal interaction in a two-degree-of-freedom structure. J. Sound Vibr. 1984, 97, 451–473. [Google Scholar] [CrossRef]
  5. Nayfeh, A.H.; Balachandran, B.; Colbert, M.A.; Nayfeh, M.A. An Experimental Investigation of Complicated Responses of a Two-Degree-of-Freedom Structure. J. Appl. Mech.-Trans. ASME 1989, 56, 960–967. [Google Scholar] [CrossRef]
  6. Balachandran, B.; Nayfeh, A.H. Nonlinear motions of beam-mass structure. Nonlinear Dyn. 1990, 1, 39–61. [Google Scholar] [CrossRef]
  7. Balachandran, B.; Nayfeh, A.H. Observations of modal interactions in resonantly forced beam-mass structures. Nonlinear Dyn. 1991, 2, 77–117. [Google Scholar] [CrossRef]
  8. Warminski, J.; Cartmell, M.P.; Bochenski, M.; Ivanov, I. Analytical and experimental investigations of an autoparametric beam structure. J. Sound Vibr. 2008, 315, 486–508. [Google Scholar] [CrossRef]
  9. Cao, D.; Zhang, W.; Yao, M. Analytical and experimental studies on nonlinear characteristics of an L-shape beam structure. Acta Mech. Sin. 2010, 26, 967–976. [Google Scholar] [CrossRef]
  10. Onozato, N.; Nagai, K.; Maruyama, S.; Yamaguchi, T. Chaotic vibrations of a post-buckled L-shaped beam with an axial constraint. Nonlinear Dyn. 2012, 67, 2363–2379. [Google Scholar] [CrossRef]
  11. Georgiades, F.; Warminski, J.; Cartmell, M.P. Towards linear modal analysis for an L-shaped beam: Equations of motion. Mech. Res. Commun. 2013, 47, 50–60. [Google Scholar] [CrossRef] [Green Version]
  12. Georgiades, F. Nonlinear equations of motion of L-shaped beam structures. Eur. J. Mech. A-Solids 2017, 65, 91–122. [Google Scholar] [CrossRef]
  13. Wei, J.; Cao, D.; Yang, Y.; Huang, W. Nonlinear Dynamical Modeling and Vibration Responses of an L-Shaped Beam-Mass Structure. J. Appl. Nonlinear Dyn. 2017, 1, 91–104. [Google Scholar] [CrossRef]
  14. Yu, T.; Zhang, W.; Yang, X. Global bifurcations and chaotic motions of a flexible multi-beam structure. Int. J. Non-Linear Mech. 2017, 95, 264–271. [Google Scholar] [CrossRef]
  15. Yu, T.; Zhang, W.; Yang, X. Global dynamics of an autoparametric beam structure. Nonlinear Dyn. 2017, 88, 1329–1343. [Google Scholar] [CrossRef]
  16. Spong, M.W. Modeling and Control of Elastic Joint Robots. J. Dyn. Syst. Meas. Control-Trans. ASME 1987, 109, 310–319. [Google Scholar] [CrossRef]
  17. Subudhi, B.; Morris, A.S. Dynamic modelling, simulation and control of a manipulator with flexible links and joints. Robot. Auton. Syst. 2002, 41, 257–270. [Google Scholar] [CrossRef]
  18. Vakil, M.; Fotouhi, R.; Nikiforuk, P.N.; Heidari, F. A study of the free vibration of flexible-link flexible-joint manipulators. Proc. Inst. Mech. Eng. Part J-J. Eng. Tribol. 2011, 225, 1361–1371. [Google Scholar] [CrossRef]
  19. Meng, D.; She, Y.; Xu, W.; Lu, W.; Liang, B. Dynamic modeling and vibration characteristics analysis of flexible-link and flexible-joint space manipulator. Multibody Syst. Dyn. 2018, 43, 321–347. [Google Scholar] [CrossRef]
  20. Yang, T.; Yan, S.; Ma, W.; Han, Z. Joint dynamic analysis of space manipulator with planetary gear train transmission. Robotica 2016, 34, 1042–1058. [Google Scholar] [CrossRef]
  21. Yang, T.; Yan, S.; Han, Z. Nonlinear model of space manipulator joint considering time-variant stiffness and backlash. J. Sound Vibr. 2015, 341, 246–259. [Google Scholar] [CrossRef]
  22. Wei, J.; Cao, D.; Lun, L.; Huang, W. Global mode method for dynamic modeling of a flexible-link flexible-joint manipulator with tip mass. Appl. Math. Model. 2017, 48, 787–805. [Google Scholar] [CrossRef]
  23. Moon, F.C.; Li, G.X. Experimental study of chaotic vibrations in a pin-jointed space truss structure. AIAA J. 1990, 28, 915–921. [Google Scholar] [CrossRef]
  24. Folkman, S.L.; Rowsell, E.A.; Ferney, G.D. Influence of pinned joints on damping and dynamic behavior of a truss. J. Guid. Control Dyn. 1995, 18, 1398–1403. [Google Scholar] [CrossRef]
  25. Bowden, M.; Dugundji, J. Joint damping and nonlinearity in dynamics of space structures. AIAA J. 1990, 28, 740–749. [Google Scholar] [CrossRef]
  26. Zhang, J.; Guo, H.; Liu, R.; Wu, J.; Kou, Z.; Deng, Z. Damping formulations for jointed deployable space structures. Nonlinear Dyn. 2015, 81, 1969–1980. [Google Scholar] [CrossRef]
  27. Guo, H.; Zhang, J.; Liu, R.; Deng, Z. Effects of Joint on Dynamics of Space Deployable Structure. Chin. J. Mech. Eng. 2013, 26, 861–872. [Google Scholar] [CrossRef]
  28. Wei, J.; Cao, D.; Huang, H.; Wang, L.; Huang, W. Dynamics of a multi-beam structure connected with nonlinear joints: Modelling and simulation. Arch. Appl. Mech. 2018, 88, 1059–1074. [Google Scholar] [CrossRef]
  29. Wei, J.; Cao, D.; Huang, H. Nonlinear vibration phenomenon of maneuvering spacecraft with flexible jointed appendages. Nonlinear Dyn. 2018, 94, 2813–2877. [Google Scholar] [CrossRef]
  30. Esmailzadeh, A.E.; Jalili, A.N. Parametric response of cantilever Timoshenko beams with tip mass under harmonic support motion. Int. J. Non-Linear Mech. 1998, 33, 765–781. [Google Scholar] [CrossRef]
  31. Yabuno, H.; Ide, Y.; Aoshima, N. Nonlinear Analysis of a Parametrically Excited Cantilever Beam. (Effect of the Tip Mass on Stationary Response). JSME Int. J. Ser. C-Mech. Syst. Mach. Elem. Manuf. 1998, 41, 555–562. [Google Scholar] [CrossRef] [Green Version]
  32. Eftekhari, M.; Mahzoon, M.; Ziaei-Rad, S. Effect of added tip mass on the nonlinear flapwise and chordwise vibration of cantilever composite beam under base excitation. Int. J. Struct. Stab. Dyn. 2012, 12, 285–310. [Google Scholar] [CrossRef]
  33. Cheng, G.; Mei, C.; Lee, R. Large Amplitude Vibration of a Cantilever Beam with Tip Mass Under Random Base Excitation. Adv. Struct. Eng. 2001, 4, 203–210. [Google Scholar] [CrossRef]
  34. Meesala, V.C.; Hajj, M.R. Parameter sensitivity of cantilever beam with tip mass to parametric excitation. Nonlinear Dyn. 2019, 95, 3375–3384. [Google Scholar] [CrossRef]
Figure 1. Schematic of the L-shaped multi-beam jointed structure with tip mass.
Figure 1. Schematic of the L-shaped multi-beam jointed structure with tip mass.
Materials 14 07279 g001
Figure 2. Schematic of the (a) natural and (b) geometrical matching conditions on the tip mass: (a) Q is the shear force of the vertical beam on the tip mass, Q = E I v 2   ( L 2 , t ) , M is the bending moment of the vertical beam on the tip mass, M = E I v 2   ( L 2 , t ) , and F m is the external force on the tip mass; (b) v L = v 2 ( L 2 , t ) , v 2 = v 2 ( L 2 , t ) .
Figure 2. Schematic of the (a) natural and (b) geometrical matching conditions on the tip mass: (a) Q is the shear force of the vertical beam on the tip mass, Q = E I v 2   ( L 2 , t ) , M is the bending moment of the vertical beam on the tip mass, M = E I v 2   ( L 2 , t ) , and F m is the external force on the tip mass; (b) v L = v 2 ( L 2 , t ) , v 2 = v 2 ( L 2 , t ) .
Materials 14 07279 g002
Figure 3. Schematic of the joint model.
Figure 3. Schematic of the joint model.
Materials 14 07279 g003
Figure 4. Schematic of the (a) geometric and (b) force matching conditions at the joint S 2 ; v 1 = v 1 ( L 1 , t ) , v 2 = v 2 ( 0 , t ) , M h and M v are the bending moments acting on the joint S 2 , M h = E I v 1 ( L 1 , t ) ,    M v = E I v 2 ( 0 , t ) .
Figure 4. Schematic of the (a) geometric and (b) force matching conditions at the joint S 2 ; v 1 = v 1 ( L 1 , t ) , v 2 = v 2 ( 0 , t ) , M h and M v are the bending moments acting on the joint S 2 , M h = E I v 1 ( L 1 , t ) ,    M v = E I v 2 ( 0 , t ) .
Materials 14 07279 g004
Figure 5. The first four natural frequencies versus m : (ad) 1st to 4th natural frequencies.
Figure 5. The first four natural frequencies versus m : (ad) 1st to 4th natural frequencies.
Materials 14 07279 g005
Figure 6. The 1st to 4th global mode shapes of the structure from ANSYS (ad) and the presented model (eh).
Figure 6. The 1st to 4th global mode shapes of the structure from ANSYS (ad) and the presented model (eh).
Materials 14 07279 g006
Figure 7. The frequency response curves of horizontal displacement of the tip mass: (a) m = 0 , J = 0 ; (b) m = 0.2 , J = 0.004 ; (c) m = 0.5 , J = 0.01 ; (d) m = 1 , J = 0.02 .
Figure 7. The frequency response curves of horizontal displacement of the tip mass: (a) m = 0 , J = 0 ; (b) m = 0.2 , J = 0.004 ; (c) m = 0.5 , J = 0.01 ; (d) m = 1 , J = 0.02 .
Materials 14 07279 g007
Figure 8. The frequency response curves of vertical displacement of the tip mass: (a) m = 0 , J = 0 ; (b) m = 0.2 , J = 0.004 ; (c) m = 0.5 , J = 0.01 ; (d) m = 1 , J = 0.02 .
Figure 8. The frequency response curves of vertical displacement of the tip mass: (a) m = 0 , J = 0 ; (b) m = 0.2 , J = 0.004 ; (c) m = 0.5 , J = 0.01 ; (d) m = 1 , J = 0.02 .
Materials 14 07279 g008
Figure 9. The spectrum plots of the steady-state responses of the tip mass: the horizontal displacement (ac) and the vertical displacement (df) with m = 0 , J = 0 at Ω = 8.25    Hz ; The horizontal displacement (gi) and the vertical displacement (jl) with m = 1 , J = 0.02 at Ω = 8.25 Hz ; (a,d,g,j) the results from the case I; (b,e,h,k) the results from the case II; (c,f,I,l) the results from the case III.
Figure 9. The spectrum plots of the steady-state responses of the tip mass: the horizontal displacement (ac) and the vertical displacement (df) with m = 0 , J = 0 at Ω = 8.25    Hz ; The horizontal displacement (gi) and the vertical displacement (jl) with m = 1 , J = 0.02 at Ω = 8.25 Hz ; (a,d,g,j) the results from the case I; (b,e,h,k) the results from the case II; (c,f,I,l) the results from the case III.
Materials 14 07279 g009
Figure 10. The frequency response curves of the structure with k L = 5 E I / L . (a) the displacement of first joint; (b) the displacement of second joint; (c) the horizontal displacement of tip mass; (d) the vertical displacement of tip mass.
Figure 10. The frequency response curves of the structure with k L = 5 E I / L . (a) the displacement of first joint; (b) the displacement of second joint; (c) the horizontal displacement of tip mass; (d) the vertical displacement of tip mass.
Materials 14 07279 g010
Figure 11. The frequency response curves of the structure with k L = E I / L . (a) the displacement of first joint; (b) the displacement of second joint; (c) the horizontal displacement of tip mass; (d) the vertical displacement of tip mass.
Figure 11. The frequency response curves of the structure with k L = E I / L . (a) the displacement of first joint; (b) the displacement of second joint; (c) the horizontal displacement of tip mass; (d) the vertical displacement of tip mass.
Materials 14 07279 g011
Figure 12. The frequency response curves of the structure with various joint damping μ . (a) the displacement of first joint; (b) the displacement of second joint; (c) the horizontal displacement of tip mass; (d) the vertical displacement of tip mass.
Figure 12. The frequency response curves of the structure with various joint damping μ . (a) the displacement of first joint; (b) the displacement of second joint; (c) the horizontal displacement of tip mass; (d) the vertical displacement of tip mass.
Materials 14 07279 g012
Table 1. Values of the parameters of the L-shaped multi-beam jointed structure.
Table 1. Values of the parameters of the L-shaped multi-beam jointed structure.
ParameterValue
Beam   density   ρ   ( kg / m 3 ) 7800
Beam   elastic   modulus   E   ( GPa ) 200
Beam   width (m) 0.02
Beam   thickness (m) 0.005
Beam   length   L 1 = L 2 = L (m) 0.3
Mass   moment   inertia   of   tip   mas   J   ( kg m 2 ) 0.01
Joint   linear   stiffness   k L ( N m / rad ) 138.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wei, J.; Yu, T.; Jin, D.; Liu, M.; Cao, D.; Wang, J. Nonlinear Dynamic Modeling and Analysis of an L-Shaped Multi-Beam Jointed Structure with Tip Mass. Materials 2021, 14, 7279. https://doi.org/10.3390/ma14237279

AMA Style

Wei J, Yu T, Jin D, Liu M, Cao D, Wang J. Nonlinear Dynamic Modeling and Analysis of an L-Shaped Multi-Beam Jointed Structure with Tip Mass. Materials. 2021; 14(23):7279. https://doi.org/10.3390/ma14237279

Chicago/Turabian Style

Wei, Jin, Tao Yu, Dongping Jin, Mei Liu, Dengqing Cao, and Jinjie Wang. 2021. "Nonlinear Dynamic Modeling and Analysis of an L-Shaped Multi-Beam Jointed Structure with Tip Mass" Materials 14, no. 23: 7279. https://doi.org/10.3390/ma14237279

APA Style

Wei, J., Yu, T., Jin, D., Liu, M., Cao, D., & Wang, J. (2021). Nonlinear Dynamic Modeling and Analysis of an L-Shaped Multi-Beam Jointed Structure with Tip Mass. Materials, 14(23), 7279. https://doi.org/10.3390/ma14237279

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop