Next Article in Journal
Performance Evaluation of Automobile Fuel Consumption Using a Fuzzy-Based Granular Model with Coverage and Specificity
Next Article in Special Issue
Critically-Finite Dynamics on the Icosahedron
Previous Article in Journal
Constraints on the String T-Duality Propagator from the Hydrogen Atom
Previous Article in Special Issue
An E-Sequence Approach to the 3x + 1 Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Dynamic Modelling of Two-Point and Symmetrically Supported Pipeline Brackets with Elastic-Porous Metal Rubber Damper

1
Engineering Research Center for Metal Rubber, School of Mechanical Engineering and Automation, Fuzhou University, Fuzhou 350116, China
2
Dongfeng Nissan Technical Center, Guangzhou 510800, China
*
Authors to whom correspondence should be addressed.
Symmetry 2019, 11(12), 1479; https://doi.org/10.3390/sym11121479
Submission received: 18 November 2019 / Revised: 27 November 2019 / Accepted: 1 December 2019 / Published: 4 December 2019
(This article belongs to the Special Issue Symmetry and Dynamical Systems)

Abstract

:
This paper aims to investigate the nonlinear dynamic properties of a two-point and symmetrically supported pipeline bracket system coated with the damping element using an elastic-porous metal rubber. The dynamic model of the studied two-point and symmetric pipeline system was established based on impulse response matrix for accurate and reliable description on its nonlinear behaviours, e.g., energy dissipation and loss factor. The experimental verification of the developed model was performed by means of dynamic test as well as the analyses of nonlinear damping characteristics. The experimental results show a good agreement with the prediction results obtained from the proposed dynamic model. This work provides an alternative method to investigate the dynamics of pipeline vibration system equipped with a damping structure.

1. Introduction

Flexible and vibrating structures of pipeline system are ubiquitous in extensive engineering areas and industrial applications of civil, aerospace, marine and so on. These structures can experiment harmful vibrations resulting from its interaction with the environment, as resonance and flutter oscillations, that if are not correctly controlled, can lead to undesirable noise and even the collapsing of the structure. In this scenario, the use of feasible damping materials for vibration control is an attractive option to mitigate such vibrations. Modelling and/or control of such damping vibration systems with nonlinear dynamic behaviours should be a challenge due to the complexity of the characteristic polynomial.
Many researchers have extensively attempted to establish dynamic models and analyse the vibration characteristics in terms of various pipeline systems. Li et al. [1] proposed a semi-analytical method base on the He’s variational iteration method to analyse the conservative free vibration problem of conveying fluid pipe. For almost the similar pipeline system, Zhou et al. [2] considered the cantilevered fluid-transporting pipe as a non-conservative system and developed a model control strategy involving the mass, damping, stiffness and suppression region using the idea of nonlinear energy sink. Zhang et al. [3] provided a feasible solution for the nonlinear vibration problem of a fluid-conveying pipeline system with general boundary working conditions, in which the pipeline could be divided into the straight pipe and the bent pipe. Rong et al. [4] analysed the dynamics of flexible beams with large deformations via the absolute nodal coordinate transfer matrix method. Some other scholars focused on the control strategies of vibration reduction associated with the nonlinear dynamics of fluid conveying curved pipeline systems [5,6,7,8,9]. The pipeline systems with two-point and symmetrically supported brackets are widely applied in some precise apparatus and have also made some efforts in the past decade. For example, Modarres-Sadeghi et al. [10] attempted to explore the dynamic characteristics of straight pipes which were supported at both ends. However, the previous efforts have tended to focus on the rigidly supported pipeline system or linear vibration, rather than on the elastic support with damping performance.
In the view of pipeline vibration system considering damping structure, some efforts on the theoretic model have been made in terms of individual cases. Amabili et al. [11] studied the nonlinear damping performance under large amplitude conditions and proposed a fractional linear solid model with the geometrically nonlinear viscoelasticity. Kong et al. [12] attempted to explore the effects of damping, the dynamical behaviours of coupled systems of nonlinear energy sink (NES) including the bifurcations and strongly modulated regimes. In particular, the influence of a 2-dof NES on the response regimes and performance of vibration suppression was analysed in detail. Rezaiee-Pajand and Sarafrazi [13] introduced a novel dynamic relaxation method with zero damping terms. They indicated that a proper time-step ratio is able to guarantee the convergence of the proposed model. Zhai et al. [14] reported that the active constrained layer damping treatment was adopted to investigate the vibration control of an aero pipeline system in terms of the vibration and stress distribution. Furthermore, a three-dimensional finite element model of such a pipeline with active constrained layer damping (ACLD) patches was developed. Recently, there has been growing interest in the high static and low dynamic characteristics with a geometric nonlinear damping system. For example, a novel design of the geometric nonlinear damping (GND) was presented to enhance the vibration isolation performance of the high-static—low-dynamic stiffness (HSLDS) isolator under both base and force excitations [15]. The nonlinear dynamic stability of Euler-Bernoulli beam-columns with damping structure was investigated by combining the analytical and numerical methods [16]. Wang et al. utilized Riccati discrete time transfer matrix method to analyse the dynamics of an underwater towed system and prove its capability [17]. Some researchers further applied Riccati discrete time transfer matrix method to analyse the nonlinear dynamics of pipeline system [18,19]. Their results indicate that dynamic relaxation methods can be an effective nonlinear dynamic analysis method for the truss and frame structures.
However, there is no comprehensive research so far available to well describe the nonlinear dynamic characteristics of the supported pipeline system that contains a damping structure. Elastic-porous metal rubber, also known as entangled wire mesh, exhibits typical nonlinear characteristics with a hysteresis loop of displacement-force curve under dynamic loading conditions [20,21,22,23]. It has been proved to effectively reduce the transmission of vibration in a variety of industrial structures such as coated structure [24], smart rotor support [25] and the pipeline systems [26]. Wu et al. [24] designed a new coated structure with large-size metal rubber bellows and investigated the energy dissipation characteristics. Their results indicated that the dynamic vibration characteristics coated structure with large-size metal rubber bellows has great dependence on the vibration amplitude and frequency. Ma et al. [25] designed a smart rotor support with shape memory alloy metal rubber (SMA-MR) and studied the mechanical performance of smart rotor support with SMA-MR under high temperature variable amplitude vibrations. The results revealed that the smart rotor support with SMA-MR can provide variable stiffness and damping characteristics with temperature, excitation amplitude and frequency. Kwon et al. [26] developed a ring-type SMA mesh washer isolator that can be used in both launch vibration and on-orbit micro-jitter environments without requiring a launch lock device. The manufacturing processes of metal rubber generally include several steps including encircled tight helix, coiled wire with a particular pitch, a rough blank in a crisscross pattern, press formation and post-treatment [27,28].
The objective of this paper is to analyse the dynamic properties of the two-point and symmetrically supported pipeline system coated with a damping structure (elastic-porous metal rubber) and establish the dynamic model for the two-point and symmetrically supported pipeline brackets. The nonlinear dynamic analysis of long and straight pipeline with low-order modes can be performed as well as the experimental verification of the developed dynamic model. In Section 2, the simplification of the physical model of vibration isolation system for two-point and symmetrically supported pipeline brackets is introduced. In Section 3, the dynamic model based on impulse response matrix for accurate and reliable description on its nonlinear behaviors is established. In Section 4, the experiment for two-point and symmetrically supported pipeline brackets under different excitation conditions are carried out to analyze the nonlinear characteristics of the coated pipeline system. In Section 5, the main conclusions of this paper are drawn.

2. Description of Pipeline Vibration Isolation System

For the pipeline coated system, due to the modal characteristics of the long pipeline and thin-plate installation base, the whole pipeline system can be regarded as a physical structure of active vibration isolation. As shown in Figure 1, the pipeline B that connected with the installation base A through the coated damping structure G (the support) is subjected to external excitation force fB, which can be weakened through the coated damping G. Thus, the pipeline system is able to realize the purpose of the vibration reduction.
In this paper, the damping structure G contains metal rubber, and pipeline B connected with installation base A through the two-point and symmetrically supported pipeline brackets. Based on these assumptions, the studies of the dynamics for the coated pipeline system were performed.

2.1. Model Simplifications

To install the installation base A with concrete size, simplifying the installation base A can be simplified as a regular elastic support with uniform thickness and considered as a uniform rectangular plate (a × b). The material properties of the installation base A mainly include elastic modulus EA, Poisson’s ratio μA, density ρA, etc. As a result of the restraint mode of the installation base A that is supported and fixed at both ends, it has only an elastic mode and no rigid body displacement.
The pipeline B can be simplified as a geometric structure with length c, and the external diameter d. For its material properties, the pipeline B has specific mechanical behaviours, such as elastic modulus EB, Poisson’s ratio μB, density ρB, etc. The pipeline B connected on the installation base A by the coated damping structure G in a manner using the two-point and symmetrically supported brackets which have elastic and rigid modes when it is suffered by an external excitation.
In terms of the coordinate system, the earth is chosen and set as the absolute coordinate system, and the left side of the installation base A is considered as the origin of the coordinate system.

2.2. Junction Definition

A corresponds relationship is established one by one between the connection point on the installation base A, and the connection point on the pipeline B are denoted as Ai (I = 1, 2) and Bi (I = 1, 2), respectively. It is assumed that the coated damping structure has three-dimensional (x, y, z) stiffness and damping characteristics. The isotropic mechanical behaviours are also assumed with the same stiffness and damping characteristics in the three directions.
For the installation base A, the set of connection degrees of freedom I is defined as:
{ x A 1 , y A 1 , z A 1 } ,   { x A 2 , y A 2 , z A 2 } ,   { f x A 1 , f y A 1 , f z A 1 } ,   { f x A 2 , f y A 2 , f z A 2 }
For the pipeline B, the set of connection degrees of freedom J can be denoted as follows:
{ x B 1 , y B 1 , z B 1 } ,   { x B 2 , y B 2 , z B 2 } ,   { f x B 1 , f y B 1 , f z B 1 } ,   { f x B 2 , f y B 2 , f z B 2 }
Based on the above two definitions, the set of connection degrees of freedom I is in one-to-one correspondence with that of J, and the number of degrees of freedom for both A and B is 6. { x A i , y A i , z A i } , { x B i , y B i , z B i } correspond to the absolute displacement components of the connecting points { A i , B i } respectively; { f x A i , f y A i , f z A i } , { f x B i , f y B i , f z B i } represent the components of the connecting force applied to A and B at the connecting points { A i , B i } ; { f x A i , f y A i , f z A i } , { f x B i , f y B i , f z B i } represent the component of force applied to the coated damping structure G at the connecting points { A i , B i } .
For the absolute displacement component, a relation generally exists as follows:
{ x A i , y A i , z A i } { x B i , y B i , z B i } ( i = 1 , 2 )
This can be attributed to the elastic connection between installation base A and pipeline B, which leads to the existence of relative displacement between freedom I and J. The relationship between acting force and reaction force between joints is then given by the following relation:
{ f x A i , f y A i , f z A i } = { f x B i , f y B i , f z B i } . ( i = 1 , 2 )

3. Dynamic Modelling

Based on the above model simplification and junction definition, the dynamic model can be established by analysing the two-point elastic partial impulse response function matrix in this section. Firstly, the mathematical model of elastic partial impulse response function matrix on I and J is put forward respectively. Then, response analysis of arbitrary points on A and B caused by connecting force is provided. Finally, the dynamic model of two-point and symmetrically supported pipeline brackets with elastic-porous metal rubber damper is proposed.

3.1. Elastic Partial Impulse Response Function Matrix on I

For the viscous damped nA DOF discrete vibration system with installation base A, the differential equation of motion can be written as follows:
M A x ¨ A + C A x ˙ A + K A x A = f A ( t ) .
where, M A , C A , K A R n A × n A ; nA is the set of discrete degree of freedom (InA) of the installation base A; f A ( t ) is the external excitation (including connection force { f x A i , f y A i , f z A i } ,   ( i = 1 , 2 ) , column vector on nA; x ¨ A , x ˙ A and x A correspond to the acceleration, velocity and displacement column vectors on the set of discrete degrees of freedom nA, respectively.
At present, the main work is to find the elastic partial impulse response function matrix on I. Thus, only the external excitation (the connection force) on I is considered and the f A ( t ) is given by
f A ( t ) R n A × 1 = { 0 , , f x A 1 , f y A 1 , f z A 1 , , f x A 2 , f y A 2 , f z A 2 , 0 } T .
The position n A { f x A 1 } of mark f x A 1 in f A ( t ) is denoted as n ˜ 1 A , then n A { f y A 1 } = n ˜ 1 A + 1 , n A { z y A 1 } = n ˜ 1 A + 2 . Similarly, the position n A { f x A 2 } of mark f x A 2 in f A ( t ) is expressed as, n ˜ 2 A , then n A { f y A 2 } = n ˜ 2 A + 1 , n A { z y A 2 } = n ˜ 2 A + 2 .
In the same manner, the velocity and the displacement for the column vectors of acceleration on the set of discrete degrees of freedom nA. They can be expressed as:
x A R n A × 1 = { , x A 1 , y A 1 , z A 1 , , x A 2 , y A 2 , z A 2 , } T
x ˙ A R n A × 1 = { , x ˙ A 1 , y ˙ A 1 , z ˙ A 1 , , x ˙ A 2 , y ˙ A 2 , z ˙ A 2 , } T
x ¨ A R n A × 1 = { , x ¨ A 1 , y ¨ A 1 , z ¨ A 1 , , x ¨ A 2 , y ¨ A 2 , z ¨ A 2 , } T
n A { x A i } = n A { x ˙ A i } = n A { x ¨ A i } = n ˜ i A ( i = 1 , 2 )
Considering the consistency of finite element analysis and modal experimental analysis for A, it is suggested that the two methods should have combined in an identical way. Each grid node is described by { x , y , z } coordinate (degree of freedom) and labelled continuously according to node and order of x y z . In this case, n ˜ i A can be determined uniquely. Then, block matrix is introduced as follows:
R = [ C A M A M A 0 ] , S = [ K A 0 0 M A ] , y = { x A x ˙ A } , z = { f A ( t ) 0 }
Equation (3) is algebraically treated according to Equation (7) and solved by homogeneous breeze equation. Then, its eigenvalues are obtained as well as feature vectors, marked as λ 1 , λ 2 , …, λ n A , λ 1 * , λ 2 * , …, λ n A * , and φ 1 , φ 2 , …, φ n A , φ 1 * , φ 2 * , …, φ n A * respectively.
By solving the equation and transforming some columns, the displacement expression of the system in the physical coordinate system under the zero initial condition is obtained as:
x A ( t ) = i = 1 n A { φ i φ i T a i 0 t e λ i ( t τ ) f A ( τ ) d τ + φ i * φ i H a i * 0 t e λ i * ( t τ ) f A ( τ ) d τ } .
where ai, called complex mode mass, is consistent with the expression of natural frequency. By applying Laplacian transformation on both sides of Equation (8) and using the properties of convolution Laplacian transformation, Equation (8) can be further transformed into:
L [ x A ( t ) ] = X A ( s ) = i = 1 n A { φ i φ i T a i ( s λ i ) + φ i * φ i H a i * ( s λ i * ) } L [ f A ( t ) ] .
In order to solve the transfer function matrix from I to the whole A, a series of transformations are requested for Equation (9). They mainly include expanding and extracting the elements in f A ( t ) of Equation (9) to form a new matrix f ˜ A ( t ) , and the concept of φ I i , φ I i * . Therefore,
{ φ i φ i T a i ( s λ i ) + φ i * φ i H a i * ( s λ i * ) } L [ f A ( t ) ] = { φ i φ I i T a i ( s λ i ) + φ i * φ I i H a i * ( s λ i * ) } L [ f ˜ A ( t ) ] .
where φ I i , φ I i * represent the components of φ i , φ i * on the mode vector, respectively. By substituting Equation (10) into Equation (9), one obtains:
L [ x A ( t ) ] = X A ( s ) = i = 1 n A { φ i φ I i T a i ( s λ i ) + φ i * φ I i H a i * ( s λ i * ) } L [ f ˜ A ( t ) ] .
In this case, the transfer function matrix from I to the whole A can be obtained as:
H A I ( s ) = X A ( s ) L [ f ˜ A ( t ) ] = X A ( s ) F ˜ A ( s ) = i = 1 n A { φ i φ I i T a i ( s λ i ) + φ i * φ I i H a i * ( s λ i * ) } .
Based on Equation (12), the impulse response function matrix from I to the whole A can be derived as:
h A I ( t ) = i = 1 n A φ i φ I i T a i e λ i t + φ i * φ I i H a i * e λ i * t .
By expanding, the Equation (11) can be transformed into:
L { , x A 1 , y A 1 , z A 1 , , x A 2 , y A 2 , z A 2 , } T = i = 1 n A { 1 a i ( s λ i ) { φ i φ I i 1 φ i φ I i 2 φ i φ I i 6 } + 1 a i * ( s λ i * ) { φ i * φ I i 1 * φ i * φ I i 2 * φ i * φ I i 6 * } } F ˜ A ( s ) .
Similarly, we extracted and reorganized the { n ˜ i A , n ˜ i A + 2 } ,   ( i = 1 , 2 ) element in x A ( t ) to obtain a new displacement column vector x ˜ A ( t ) on I, that is
x ˜ A ( t ) = { x A 1 , y A 1 , z A 1 , x A 2 , y A 2 , z A 2 } T .
Based on Equations (14a) and (14b), the I I impulse response function matrix h I I ( t ) can be derived as:
h I I ( t ) = L 1 [ H I I ( s ) ] = L 1 [ X ˜ A ( s ) F ˜ A ( s ) ] = L 1 [ i = 1 n A { φ I i φ I i T a i ( s λ i ) + φ I i * φ I i H a i * ( s λ i * ) } ] = i = 1 n A φ I i φ I i T a i e λ i t + φ I i * φ I i H a i * e λ i * t .
Due to the necessity of distinguishing the parameters A and B, a i a A i , λ i λ A i , a i * a A i * , λ i * λ A i * are introduced, and Equation (15a) can be rewritten as:
h I I ( t ) = i = 1 n A φ I i φ I i T a A i e λ A i t + φ I i * φ I i H a A i * e λ A i * t .

3.2. Impulse Response Function Matrix on J

The pipeline B and installation base A elastically connected each other by the coated damping structure G in a manner of two-point and symmetrically supported brackets which have elastic and rigid modes when they are suffered external excitation. Then, elastic modes are analysed firstly. The differential equation of motion for the viscous damped nB DOF discrete vibration system of B can be achieved according to the analysis of A as follows:
M B x ¨ B + C B x ˙ B + K B x B = f B ( t ) .
where, M B , C B , K B R n B × n B ; nB is the set of discrete degree of freedom (InA) of the base B; f B ( t ) is the external excitation (including connection force { f x B i , f y B i , f z B i } ,   ( i = 1 , 2 ) , column vector on nA; x ¨ B , x ˙ B and x B correspond to the acceleration, velocity and displacement column vectors on the set of discrete degrees of freedom nB, respectively.
To achieve the impulse response function matrix from J to B, only the external excitation on J is considered when f B ( t ) is set. Then,
f B ( t ) R n B × 1 = { 0 , , f x B 1 , f y B 1 , f z B 1 , , f x B 2 , f y B 2 , f z B 2 , 0 } T .
Referring to the derivation process of A, and introducing similar vector element position symbols, block matrices, and vectors as well as meshing B in the same way, the elastic partial transfer function matrix from J to B can be shown as below after a series of deductions:
H ˜ B J ( s ) = X B ( s ) F ˜ B ( s ) = i = 1 n B { φ i φ J i T a i ( s λ i ) + φ i * φ J i H a i * ( s λ i * ) } .
After simplification, the elastic partial transfer function matrix from J to B can be obtained:
h ˜ B J ( t ) = i = 1 n B φ i φ J i T a i e λ i t + φ i * φ J i H a i * e λ i * t .
Then, referring to the deduction process and considering the same parameter discrimination between A and B, a i a B i , λ i λ B i , a i * a B i * , λ i * λ B i * are introduced. Therefore, the elastic partial impulse response function matrix on J to J can be represented as:
h ˜ J J ( t ) = i = 1 n B φ J i φ J i T a B i e λ B i t + φ J i * φ J i H a B i * e λ B i * t .
The relevant symbolic parameters in Equation (18) to Equation (20) are specified in the derivation process of A.
In order to improve the reliability of the analysis about the rigid modes, we need to locate the coordinate system on B (mass mB, moment of inertia { J x , J y , J z } ), and establish the origin of this coordinate system at the centre of mass of B. As shown in Figure 2, the coordinate is parallel to the coordinate system { x , y , z } set up earlier. Then, a series of assumptions are made, including assuming L at any point on B. L { x l , y l , z l , θ x l , θ y l , θ z l } denotes six degrees of freedom of L, and L { x L , y L , z L } represents the position of L point in B coordinate system. Also, p J denotes as any degree of freedom over J; f p f B ( t ) is expressed as a force exerted on the degree of p-freedom; and P { x P , y P , z P } { B 1 , B 2 } denotes the position of point P in the coordinate system.
The transfer function of coordinate x l can be derived when f P is parallel to x c direction. The equation of motion for the centre of mass B can be written as follows:
f P ( t ) = m B x ¨ c = m B x ¨ l { 1 } .
By applying Laplacian transformation on Equation (21) under zero initial condition, Equation (21) can be reformulated as:
H l p { 1 } ( s ) = X l { 1 } ( s ) F p ( s ) = 1 m B s 2 .
where F p ( s ) , X l { 1 } ( s ) are the Laplace transformation of f p ( t ) , x l { 1 } ( t ) , respectively. H l p { 1 } ( s ) is generated by displacement x l { 1 } of L point along x c direction under the action of f p .
The equation of motion of B centre of mass around z c can be expressed as follows:
M f P = f p ( t ) z p = J z θ ¨ z l .
Then, a series of deductive transformations is applied for Equation (23), and a formula can be finally obtained by Laplace transformation. That is
H l p { 2 } ( s ) = X l { 2 } ( s ) F p ( s ) = z p y L J z s 2 .
where, F p ( s ) , X l { 2 } ( s ) denote the Laplace transformation of f p ( t ) , x l { 2 } ( t ) , respectively. H l p { 2 } ( s ) is achieved when L point rotates θ z l angle along z c direction under the action of f p .
The equation of motion of B centre of mass around y c can be written as follows:
M f P = f p ( t ) y p = J y θ ¨ y l .
Similarly, based on the above, one can obtain a formula as follows:
H l p { 3 } ( s ) = X l { 3 } ( s ) F p ( s ) = y p z L J y s 2 .
where F p ( s ) , X l { 3 } ( s ) represent the Laplace transformation of x l { 3 } ( t ) and f p ( t ) , respectively. H l p { 3 } ( s ) is generated when L point rotates θ y l angle along y c direction under the action of f p .
Thus, based on Equations (22), (24) and (26), one obtains:
H l p ( s ) = X l ( s ) F p ( s ) = { H l p { 1 } ( s ) + H l p { 2 } ( s ) + H l p { 3 } ( s ) } = { 1 m B s 2 + z p y L J z s 2 + y p z L J y s 2 } .
This means the transfer function of the physical coordinate is x l x B ( t ) when the force of f p ( t ) f B ( t ) is applied to the p-DOF, p J . Due to the multiple direction of the action of p, Equation (27) can be rewritten to make it easy to distinguish as follows:
H x l p x ( s ) = { 1 m B s 2 + z p x y L J z s 2 + y p x z L J y s 2 } .
In the Equation (28a), H x l p x ( s ) denotes the transfer function of P’s force vector f p x applied along x c in x l coordinate; z p x , y p x denote the vertical distance between f p x and z c , y c .
Similarly, the transfer functions of physical coordinates y l and z l can be obtained:
H y l p x ( s ) = { z p x x L J z s 2 } .
H z l p x ( s ) = { y p x x L J y s 2 } .
Meanwhile, on the basis of the deductions above, the transfer function of force f p in the same direction as y c and direction z c can be obtained, where f p ( t ) f B ( t ) .
The f p is parallel to y c direction:
H x l p y ( s ) = { z p y y L J z s 2 } .
H y l p y ( s ) = { 1 m B s 2 + z p y x L J z s 2 + x p y z L J x s 2 } .
The f p is parallel to z c direction:
H x l p z ( s ) = { y p z z L J y s 2 } .
H y l p z ( s ) = { x p z z L J x s 2 } .
H z l p z ( s ) = { 1 m B s 2 + x p z y L J x s 2 + y p z x L J y s 2 } .
In order to derivate conveniently, the following assumptions are proposed: (i) Assuming that B is discretized into NB elements or grid nodes, and defining nB = 3 × NB, node L { { x L , y L , z L } , { x l , y l , z l } x B ( t ) } N B , where l is related to L, {xl, yl, zl} is marked as {xLl, yLl, zLl}, L = 1, 2, …, NB. (ii) P{{xP, yP, zP} ∈ J, P = 1, 2; { x p x , x p y , x p z } , { y p x , y p y , y p z } , { z p x , z p y , z p z } } ∈ NB, p = 1, 2, where p is related to P.
Actually, the partial transfer function matrix H ¯ B J ( s ) from J to the whole B rigid body should be a nB × 6 order matrix. According to the difference of excitation and measurement points, 6 items in H ¯ B J ( s ) can be obtained respectively. Then, the partial transfer function matrix H ¯ B J ( s ) from J to the whole B rigid body can be obtained:
H ¯ B J ( s ) n B × 6 = [ H ¯ B J ( s ) | 1 , H ¯ B J ( s ) | 2 , H ¯ B J ( s ) | 3 , H ¯ B J ( s ) | 4 , H ¯ B J ( s ) | 5 , H ¯ B J ( s ) | 6 ] .
where matrix elements H ¯ B J ( s ) | i are represented as:
H ¯ B J ( s ) | 1 = { H x 1 l 1 x ( s ) , H y 1 l 1 x ( s ) , H z 1 l 1 x ( s ) , , H x N B l 1 x ( s ) , H y N B l 1 x ( s ) , H z N B l 1 x ( s ) } T H ¯ B J ( s ) | 2 = { H x 1 l 1 y ( s ) , H y 1 l 1 y ( s ) , H z 1 l 1 y ( s ) , , H x N B l 1 y ( s ) , H y N B l 1 y ( s ) , H z N B l 1 y ( s ) } T H ¯ B J ( s ) | 3 = { H x 1 l 1 z ( s ) , H y 1 l 1 z ( s ) , H z 1 l 1 z ( s ) , , H x N B l 1 z ( s ) , H y N B l 1 z ( s ) , H z N B l 1 z ( s ) } T H ¯ B J ( s ) | 4 = { H x 1 l 2 x ( s ) , H y 1 l 2 x ( s ) , H z 1 l 2 x ( s ) , , H x N B l 2 x ( s ) , H y N B l 2 x ( s ) , H z N B l 2 x ( s ) } T H ¯ B J ( s ) | 5 = { H x 1 l 2 y ( s ) , H y 1 l 2 y ( s ) , H z 1 l 2 y ( s ) , , H x N B l 2 y ( s ) , H y N B l 2 y ( s ) , H z N B l 2 y ( s ) } T H ¯ B J ( s ) | 6 = { H x 1 l 2 z ( s ) , H y 1 l 2 z ( s ) , H z 1 l 2 z ( s ) , , H x N B l 2 z ( s ) , H y N B l 2 z ( s ) , H z N B l 2 z ( s ) } T
Based on Equations (31) and (32), the impulse response function matrix h ¯ B J ( t ) can be obtained from J to the whole rigid part of B, which can be expressed as:
h ¯ B J ( t ) n B × 6 = L 1 [ H ¯ B J ( s ) ] = L 1 [ H ¯ B J ( s ) | 1 , H ¯ B J ( s ) | 2 , , H ¯ B J ( s ) | 6 ] .
By examining each column in Equation (33) separately, the impulse response function matrix h ¯ B J ( t ) from J to the whole rigid part of B can be reformulated after itemized expansion, collation and reorganization.
h ¯ B J ( t ) n B × 6 = r = 1 6 A r = r = 1 6 ϕ r ϕ ˜ r T m r t ( t 0 ) , ( r = 1 ~ 3 , m r = m B ; r = 4 , m r = J x ; r = 5 , m r = J y ; r = 6 , m r = J z ) .
in which, modes with different order can be achieved depending on the level of r. Meanwhile, the impulse response function matrix h ¯ J J ( t ) of J J rigid part can be given as:
h ¯ J J ( t ) 6 × 6 = r = 1 6 ϕ J r ϕ J r T m r t ( t 0 ) , ( r = 1 ~ 3 , m r = m B ; r = 4 , m r = J x ; r = 5 , m r = J y ; r = 6 , m r = J z ) .
Then, integrating the h ˜ J J ( t ) and h ¯ J J ( t ) given by Equations (20) and (35), the impulse response function matrix h J J ( t ) 6 × 6 of J J can be obtained as follows:
h J J ( t ) = h ˜ J J ( t ) + h ¯ J J ( t ) = i = 1 n B φ J i φ J i T a B i e λ B i t + φ J i * φ J i H a B i * e λ B i * t + r = 1 6 ϕ J r ϕ J r T m r ( t 0 ) , ( r = 1 ~ 3 , m r = m B ; r = 4 , m r = J x ; r = 5 , m r = J y ; r = 6 , m r = J z ) .

3.3. Response Analysis of Arbitrary Points on A and B Caused by Connecting Force

Due to the consistency between the parameters of A and B, we are denoting a i a A i , λ i λ A i , a i * a A i * , λ i * λ A i * to distinguish them conveniently. Thus, Equation (8) can be written as:
x A ( t ) = i = 1 n A { φ i φ i T a A i 0 t e λ A i ( t τ ) f A ( τ ) d τ + φ i * φ i H a A i * 0 t e λ A i * ( t τ ) f A ( τ ) d τ } .
For the complex conjugate pairs of φ i φ i T a A i 0 t e λ A i ( t τ ) f A ( τ ) and φ i * φ i H a A i * 0 t e λ A i * ( t τ ) f A ( τ ) , Equation (37) can be reformulated as:
x A ( t ) = 2 R e { i = 1 n A { φ i φ i T a A i 0 t e λ A i ( t τ ) f A ( τ ) d τ } } .
Expanding the sum term of Equation (38), the element f A ( t ) can be extracted and recombined and further rewritten as follows:
x A ( t ) = 2 R e { i = 1 n A { 1 a A i 0 t { j = 1 6 φ I i j φ i f ˜ A j ( τ ) } e λ A i ( t τ ) d τ } } .
where f ˜ A ( t ) is composed of { n ˜ i A , n ˜ i A + 2 } ( i = 1 ,   2 ,   3 ,   4 ) and elements extracted from f A ( t ) , f ˜ A j ( t ) represents the j-th elements in f ˜ A ( t ) . Similarly, φ I i is the composition of { n ˜ i A , n ˜ i A + 2 } ( i = 1 , 2 , 3 , 4 ) elements extracted from φ i , φ I i j ( t ) , which means the j-th element in φ I i .
Therefore, the physical coordinate response of any displacement response of a node caused by I on the k ∈ {1, nA} of A can be calculated by the following:
x A k ( t ) = 2 R e { i = 1 n A { 1 a A i 0 t { j = 1 6 φ I i j φ i k f ˜ A j ( τ ) } e λ A i ( t τ ) d τ } } .
where x A k ( t ) denotes the k-th element in x A ( t ) , and φ i k means the k-th element in φ i .
The following works are to analyse the response of arbitrary points on B caused by connecting force. Due to the convenience of distinguishing, a i a B i , λ i λ B i , a i * a B i * and λ i * λ B i * are marked. Then, rewriting Equation (19) and according to Equation (34), one can be obtained as follows:
h B J ( t ) = h ˜ B J ( t ) + h ¯ B J ( t ) = i = 1 n B φ i φ J i H a B i e λ B i t + φ i * φ J i H a B i * e λ B i * t + i = 1 6 ϕ r ϕ J r T m r = 2 R e { i = 1 n B φ i φ J i H a B i e λ B i t } + i = 1 6 ϕ r ϕ J r T m r t ( t 0 ) , ( r = 1 ~ 3 , m r = m B ; r = 4 , m r = J x ; r = 5 , m r = J y ; r = 6 , m r = J z ) .
The response vector x B ( t ) on B caused by J can be obtained by convolution h B J ( t ) f ˜ B ( t ) of matrix and vector, where f ˜ B ( t ) and f ˜ A ( t ) are identical. Additionally, the convolution h B J ( t ) f ˜ B ( t ) of matrix and vector is given as:
x B ( t ) = h B J ( t ) f ˜ B ( t ) = { 2 R e { i = 1 n B φ i φ J i T a B i e λ B i t } + r = 1 6 ϕ r ϕ J r T m r t } f ˜ B ( t ) = 0 t 2 R e { i = 1 n B φ i φ J i T a B i e λ B i ( t τ ) f ˜ B ( τ ) d τ } + 0 t r = 1 6 ϕ r ϕ J r T m r ( t τ ) f ˜ B ( τ ) d τ = 2 R e { i = 1 n B 1 a B i 0 t { φ i φ J i T f ˜ B ( τ ) } e λ B i ( t τ ) d τ } + r = 1 6 1 m r 0 t { ϕ r ϕ J r T f ˜ B ( τ ) } ( t τ ) d τ .
By simplifying Equation (42), a formula with a form convenient for calculation is obtained as:
x B ( t ) = 2 R e { i = 1 n B { 1 a B i 0 t { j = 1 6 φ i φ J i j T f ˜ B j ( t ) } e λ B i ( t τ ) d τ } } + r = 1 6 1 m r 0 t { j = 1 6 ϕ r ϕ J r j T f ˜ B j ( τ ) } ( t τ ) d τ ( t 0 ) , ( r = 1 ~ 3 , m r = m B ; r = 4 , m r = J x ; r = 5 , m r = J y ; r = 6 , m r = J z ) .
where φ J i j T denotes the j-th element in φ J i T , φ J r j T represents the j-th element in φ J r T , and f ˜ B j represents the j-th element in f ˜ B .
Based on Equation (43), the k-th physical coordinate response on any B caused by J can be calculated as:
x B k ( t ) = 2 R e { i = 1 n B { 1 a B i 0 t { j = 1 6 φ i k φ J i j T f ˜ B j ( t ) } e λ B i ( t τ ) d τ } } + r = 1 6 1 m r 0 t { j = 1 6 ϕ r k ϕ J r j T f ˜ B j ( τ ) } ( t τ ) d τ ( t 0 ) , ( r = 1 ~ 3 , m r = m B ; r = 4 , m r = J x ; r = 5 , m r = J y ; r = 6 , m r = J z ) .
where x B k ( t ) is the k-th element in x B ( t ) , φ i k is the k-th element in φ i , and ϕ r k denotes the k-th element in ϕ r .

3.4. Dynamics Model

Displacement response vectors of linear structures A and B subjected to external excitation are recorded as x A ( t ) R n A × 1 , x B ( t ) R n B × 1 , in this case, where A, B and other structures are mutually coupled. Due to A and B being connected by G, the connecting forces f I ( t ) and f J ( t ) act on the degrees of freedom I and J, which makes the displacement response become x A ( t ) R n A × 1 , x B ( t ) R n B × 1 .
x A ( t ) = x A ( t ) 0 t h A I ( t ζ ) f I ( ζ ) d ζ ,   I = A G .
x B ( t ) = x B ( t ) 0 t h B J ( t ζ ) f J ( ζ ) d ζ ,   I = B G .
where
f I ( t ) = f ˜ A ( t ) R 6 × 1 = { f x A 1 , f y A 1 , f z A 1 , f x A 2 , f y A 2 , f z A 2 } .
in which, f J ( t ) and f I ( t ) are equal in size and opposite in direction. Unit impulse response matrix h A I ( t ) and h B J ( t ) are determined by Equations (13) and (41), respectively, or measured by experiments.
The equation of motion of damped structure G can be described as:
g N ( y ( t ) , y ˙ ( t ) , t ) + f I ( t ) = 0
where g N ( y ( t ) , y ˙ ( t ) , t ) R 6 × 1 means the constitutive relation vector of G, y ( t ) and y ˙ ( t ) has the following forms:
y ( t ) R 6 × 1 = x J ( t ) x I ( t ) = { x B 1 , y B 1 , z B 1 , x B 2 , y B 2 , z B 2 } { x A 1 , y A 1 , z A 1 , x A 2 , y A 2 , z A 2 } .
y ˙ ( t ) R 6 × 1 = x ˙ J ( t ) x ˙ I ( t ) = { x ˙ B 1 , y ˙ B 1 , z ˙ B 1 , x ˙ B 2 , y ˙ B 2 , z ˙ B 2 } { x ˙ A 1 , y ˙ A 1 , z ˙ A 1 , x ˙ A 2 , y ˙ A 2 , z ˙ A 2 } .
g N ( y ( t ) , y ˙ ( t ) , t ) R 6 × 1 = g 0 { y ( t ) , y ˙ ( t ) } + z ( t ) . g 0 { y ( t ) , y ˙ ( t ) } = a 0 sgn { y ( t ) } + n = 1 n 1 a n | y ( t ) | n 1 y ( t ) + b 0 sgn { y ˙ ( t ) } + n = 1 n 2 b n | y ˙ ( t ) | n 1 y ˙ ( t ) . d z ( t ) = k s 2 [ 1 + sgn { z s | z ( t ) | } ] d y ( t ) , k s = z s y s .
where { g 0 , z ( t ) , a 0 , a n , n 1 , b 0 , b n , n 2 , k s , z s , y s } R 6 × 1 .
Equation (45a) to (49) show that when the characteristics of each part of the pipeline structure are determined and can be measured and solved, the independent unknown quantity with the smallest dimension in the pipeline structure is the physical quantity at both ends of the non-linear element, namely the relative displacement or the connecting force. Therefore, the response of the pipeline structure can be obtained by solving the independent unknown variables and the equation.
Utilizing the Fourier transform of Equations (45a), (45b) and (48a), the corresponding Fourier spectrum of the original time domain function can be calculated as:
X A ( ω ) = X A ( ω ) H A I ( ω ) F I ( ω ) .
X B ( ω ) = X B ( ω ) + H B J ( ω ) F J ( ω ) .
Y ( ω ) = X J ( ω ) X I ( ω ) .
If only the set of connection degrees of freedom for Equations (50a) and (50b) are considered, the two formulas can be rewritten:
X I ( ω ) = X I ( ω ) H I I ( ω ) F I ( ω ) .
X J ( ω ) = X J ( ω ) + H J J ( ω ) F J ( ω ) .
where H I I ( ω ) and H J J ( ω ) are determined analytically by Equations (15b) and (36) or measured experimentally. Considering the Equation (51), and F I ( ω ) = F J ( ω ) , an expression is shown as follows:
Y ( ω ) = X J ( ω ) X I ( ω ) = X J ( ω ) X I ( ω ) + H J J ( ω ) F J ( ω ) + H I I ( ω ) F I ( ω ) = X J ( ω ) X I ( ω ) + { H I I ( ω ) + H J J ( ω ) } F I ( ω )
Equation (53) can be further transformed into:
{ H I I ( ω ) + H J J ( ω ) } F I ( ω ) = Y ( ω ) + X I ( ω ) X J ( ω ) .
By further simplifying and multiplying the compound matrix { H I I ( ω ) + H J J ( ω ) } 1 on both sides of Equation (54):
F I ( ω ) = { H I I ( ω ) + H J J ( ω ) } 1 { Y ( ω ) + X I ( ω ) X J ( ω ) } .
introducing the time domain vector p ( t ) and defining that p ( t ) R 6 × 6 = F 1 = { { H I I ( ω ) + H J J ( ω ) } 1 } . Meanwhile, defining y ( t ) as the relative displacement response vector on the connecting surface caused by external excitation while linear structure A and B are not coupled with each other, so that y ( t ) R 6 × 1 = x J ( t ) x I ( t ) . Thus, the time-domain form of Equation (55) can be written as:
f I ( t ) = + p ( t ζ ) { y ( ζ ) y ( ζ ) } d ζ .
By substituting Equation (56) into Equation (47), the dynamic equations regarding the displacement response as unknown variables are obtained:
0 + p ( t ζ ) { y ( ζ ) y ( ζ ) } d ζ + g N ( y ( t ) , y ˙ ( t ) , t ) = 0 .
From Equation (54):
Y ( ω ) = { H I I ( ω ) + H J J ( ω ) } F I ( ω ) X I ( ω ) + X J ( ω ) .
The inverse Fourier transform is performed on both sides of Equation (59). Then,
y ( t ) = 0 t { h I I ( t ζ ) + h J J ( t ζ ) } f I ( ζ ) d ζ x I ( t ) + x J ( t ) = 0 t q ( t ζ ) f I ( ζ ) d ζ + y ( t ) .
where
q ( t ) = h I I ( t ) + h J J ( t ) .
The following equation can be obtained by using integral derivative method with parametric variables on both sides of Equation (59).
y ˙ ( t ) = 0 t q ˙ ( t ζ ) f I ( ζ ) d ζ + y ˙ ( t )
Then, substituting Equations (59) and (61) into Equation (47), the dynamic equations regarding connection force as an unknown variable can be obtained:
g N { 0 t q ( t ζ ) f I ( ζ ) d ζ + y ( t ) , 0 t q ˙ ( t ζ ) f I ( ζ ) d ζ + y ˙ ( t ) , t } + f I ( t ) = 0 .
By substituting Equation (47) into Equation (61) and utilizing the calculus relation, one can obtain as follows:
y ( t ) = 0 t y ˙ ( ζ ) d ζ .
y ˙ ( t ) = 0 t q ˙ ( t ζ ) g N ( y ( ζ ) , y ˙ ( ζ ) , ζ ) d ζ + y ˙ ( t ) .
As a consequence of the above, the dynamic equations on the connecting surface have been established completely. To verify the reliability of the dynamic equation, the necessity of experimental verification should be considered.

4. Experimental Verification

In this section, the reliability of the dynamic model presented in Section 3 is considered to be verified by the physical experiment. The main section context includes the introduction of the experiment equipment, dynamic testing methods and result analysis for the two-point and symmetrically pipeline system.

4.1. Case Description

To further explore the nonlinear dynamic model proposed in this paper, according to the installation method of engineering application, the experimental bench of the coated pipeline system as shown in Figure 3 was established. The experimental system mainly contained the following five parts: The excitation system, data acquisition system and measurement system; the pipeline with two-point and symmetrically supported brackets; the coated damping structure; the installation base; and the rigid connection. The excitation system consists of JZK-50 vibration exciter (maximum excitation force 500 N, amplitude ±10 mm) and E5874A power amplifier. The data acquisition system consists of VT-900X vibration controller with data acquisition and analysis software. The measurement system consists of YD-303 piezoelectric quartz force sensor and KD9004 eddy current displacement senor. The coated damping structure was composed of elastic-porous metal rubber, coated rings, and the specific coated form. The related experiment equipment and effective test method can be found in the previous works [27]. It should be pointed out that the coated damping structure with elastic-porous metal rubber has a sound effect on vibration reduction of the pipeline system [28]. In this work, the application of elastic-porous metal rubber for the damping element enwrapped on the outer wall of the pipeline was due to its excellent vibration absorption performance [29,30].

4.2. Nonlinearity and Energy Consumption

During the experimental test, the sinusoidal excitation force with various frequencies was applied to the coated pipeline system and its displacement response was measured. For the unification of excitation magnitude, an excitation condition of 40 mvpp for the pipeline system with the voltage loading amplitude was adopted. The curve of response force F and response displacement y can be plotted as Figure 4 when the excitation frequencies of 15 Hz and 20 Hz were used respectively.
It can be observed from Figure 4 that the coated pipeline system has an obvious nonlinearity. This can mainly be caused by the used elastic-porous metal rubber in the system as a nonlinear damping structure with energy dissipation characteristics. It should be mentioned that the total restoring force of elastic-porous metal rubber coated damped structure is composed of nonlinear elastic restoring force and nonlinear damping force. Nonlinear elastic restoring force denoted by an inclined curve in Figure 4 plays a restoring role in the working process of elastic-porous metal rubber, and nonlinear damping force mainly dissipates the damping energy. The measured dynamic results show a good accordance with the analytical expression of Equation (49). In Equation (49), z(t) as a nonlinear component of elastic-porous metal rubber has a hysteresis characteristic. In the view of physical mechanism based on material microstructure, z(t) is associated with the slip force and slip displacement produced by the dynamic movement between the coiled wires, which possesses both stiffness and damping.
According to the dynamic parameter identification of loss factor η, energy dissipation ∆W and maximum elastic potential energy W in [30,31], the energy dissipation characteristics of the corresponding coated pipeline system under 15 Hz and 20 Hz experimental conditions are obtained, as listed in Table 1.
Table 1 indicates that a reduction existing in the energy dissipation ∆W with the increase of frequency. This phenomenon can be explained by the displacement of the coated damping structure decreases with the excitation frequency increases under a constant excitation force. From the microscopic level, with the decrease of displacement, the slip distance between metal wires in metal rubber is reduced, and the capacity for energy dissipation of elastic-porous metal rubber becomes weaker. Meanwhile, it is known that the maximum elastic potential energy is proportional to displacement, and the maximum elastic potential energy of elastic-porous metal rubber may degrade with the reduction of displacement. The loss factor in Table 1 represents the ratio of the energy dissipation to the maximum elastic potential energy. The maximum elastic potential energy shows a faster speed of degradation than the energy dissipation, which leads to an augmentation of loss factor.
The experimental dynamic results associated with the energy dissipation characteristics would be nonlinear in terms of the change of frequency. This is attributed to the resonance frequency of the coated pipeline system. Thus, the dynamic test with the sinusoidal sweep frequency should be considered in the following work. It should be noted that one limitation of the conducted research is that it is difficult to accurately identify dynamic properties of the used coating damping structure. This is due to the damping element, namely elastic-porous metal rubber, which exhibits obvious nonlinear characteristics and has not been determined very well so far.

4.3. Dynamic Test and Its Result

The dynamic test with a broadband sweep frequency range (5–200 Hz) for the coated pipeline system was performed to improve the reliability of the experiment. Two layers of elastic-porous metal rubber were used, and the magnitude of excitation applied is set as 20 N. The experimentally measured force transfer rate-frequency curve can be plotted as Figure 5. It can be seen that the system exists three modes in the frequency range of 5–200 Hz. They are 14.610 Hz, 93.657 Hz, 186.76 Hz, respectively. The corresponding response force F and the response displacement y are not consistent under the conditions of 15 Hz and 20 Hz as mentioned in the Section 4.2. Due to the condition of 15 Hz as the resonance frequency band, the response force of the whole system is far greater than the basic force. This is reasonable for the dynamic analysis. Due to the increase of the response displacement caused by the resonance, the displacement under 15 Hz is larger than other frequencies. However, when the excitation frequency is up to 20 Hz, the non-resonant station of the coated pipeline system occurs, resulting in the degradation of the response force and the corresponding response displacement.
Figure 6 shows that the peak of the force transfer rate TAm decreases gradually and the structural loss factor η linearly increases with the increase of the magnitude of the excitation. This is due to that the vibration damping amplitude of the elastic-porous metal rubber increases with the increase of excitation magnitude. The slip amplitude between the inner metal wires increases, leading to the increase of energy dissipation. The energy dissipation of elastic-porous metal rubber is generated by the sliding dry friction between the inner hooked wires. When the excitation magnitude increases from 10 N to 80 N, the force transfer rate obviously decreases and the loss factor shows an obvious increase trend. The excitation magnitude plays a considerable role on the energy dissipation of elastic-porous metal rubber and its vibration damping system.
Integrated with the excitation magnitude of 80 N, different pre-tightening conditions (0.5 mm, 2 mm, 3 mm) were applied to the coated pipeline system. Figure 7 shows the corresponding force transfer rate curves. The peak of the force transfer rate TAm can be obtained and then the structural loss factor η is calculated based on the half-power method, as listed in Table 2.
As can be seen from Figure 7 and Table 2, the variation of pre-tightening conditions not only affects the natural frequency of the system, but also affects the energy dissipation characteristics of the coated pipeline system. The natural frequency increases with the increase of pre-tightening. The peak TAm of the force transfer rate increases, the structural loss factor η decreases, and the energy dissipation characteristics degrades. The source of this phenomenon may be that with the increase of pre-tightening, the contact pairs between metal wires are obviously dense and the slipping movements between wires become relatively difficult. When the excitation condition is always 80 N, the slip amplitude is relatively slight under the condition of large pre-tightening with the same exciting magnitude. Additionally, the slip amplitude decreases with the increase of pre-tightening.

5. Conclusions

In this paper, a dynamic analysis based on impulse response matrix to describe two-point and symmetrically supported pipeline bracket system was performed by means of mathematically modelling. The dynamic model involved of rigid pipeline, nonlinear elastic connection structure and coated damping structure is simplified and established. As for the experimental verification, the sinusoidal sweep dynamic tests were carried out to obtain the first three modes of the coated pipeline damping system and compared with the analytical results. The experimental response force-displacement curves and the energy dissipation characteristics under different excitation conditions were obtained to analyse the nonlinear characteristics of the coated pipeline system. The dynamic experimental results show a good agreement with the developed dynamic model. This work should be of great significance for the dynamic analysis of pipeline system with a damping structure. The obtained dynamic model should be further explored and validated by means of more application cases with different boundary configurations. For example, two-point and symmetrically supported pipeline bracket systems with different damping elements or different service environments.

Author Contributions

X.X. and K.X. wrote the main original draft; S.R. and A.L. performed the experiments; S.R. and K.X. analyzed the data; H.B. and K.X. developed the theoretical model; X.X. and H.B. are the project administration and review the draft.

Funding

This research was funded by the National Natural Science Foundation of China (NSFC) through project number [51705080], and the Natural Science Foundation of Fujian Province through project number [2018J01764]. And the APC was funded by Xin Xue grant number [0020-510726].

Acknowledgments

The authors gratefully acknowledge the support of the National Natural Science Foundation of China (NSFC) through Grant No. 51705080, and the Natural Science Foundation of Fujian Province (No. 2018J01764).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this article.

References

  1. Li, Y.D.; Yang, Y.R. Vibration analysis of conveying fluid pipe via He’s variational iteration method. Appl. Math. Model. 2016, 43, 409–420. [Google Scholar]
  2. Zhou, K.; Xiong, F.R.; Jiang, N.B.; Dai, H.L.; Yan, H.; Wang, L.; Ni, Q. Nonlinear vibration control of a cantilevered fluid-conveying pipe using the idea of nonlinear energy sink. Nonlinear Dyn. 2019, 95, 1435–1456. [Google Scholar] [CrossRef]
  3. Zhang, T.; Qu, Y.H.; Zhang, Y.O.; Lv, B.L. Nonlinear dynamics of straight fluid-conveying pipes with general boundary conditions and additional springs and masses. Appl. Math. Model. 2016, 40, 7880–7900. [Google Scholar] [CrossRef]
  4. Bao, R. Efficient dynamics analysis of large-deformation flexible beams by using the absolute nodal coordinate transfer matrix method. Multibody Syst. Dyn. 2014, 32, 535–549. [Google Scholar]
  5. Lin, W.; Qiao, N. Nonlinear dynamics of a fluid-conveying curved pipe subjected to motion-limiting constraints and a harmonic excitation. J. Fluids Struct. 2008, 24, 96–110. [Google Scholar] [CrossRef]
  6. Zhang, Z.W.; Wang, Y.J.; Wang, W.; Tian, Y.L. Periodic Solution of the Strongly Nonlinear Asymmetry System with the Dynamic Frequency Method. Symmetry 2019, 11, 676. [Google Scholar] [CrossRef] [Green Version]
  7. Meng, D.; Guo, H.Y.; Xu, S.P. Non-linear dynamic model of a fluid-conveying pipe undergoing overall motions. Appl. Math. Model. 2011, 35, 781–796. [Google Scholar] [CrossRef]
  8. Yu, D.L.; Wen, J.H.; Zhao, H.G. Vibration reduction by using the idea of phononic crystals in a pipe-conveying fluid. J. Sound Vib. 2008, 318, 193–205. [Google Scholar] [CrossRef]
  9. Bi, K.M.; Hao, H. Using pipe-in-pipe systems for subsea pipeline vibration control. Eng. Struct. 2016, 109, 75–84. [Google Scholar] [CrossRef] [Green Version]
  10. Modarres-Sadeghi, Y.; Paıdoussis, M.P. Nonlinear dynamics of extensible fluid-conveying pipes, supported at both ends. J. Fluids Struct. 2009, 25, 535–543. [Google Scholar] [CrossRef]
  11. Amabili, M. Nonlinear damping in large-amplitude vibrations: Modelling and experiments. Nonlinear Dyn. 2018, 93, 5–18. [Google Scholar] [CrossRef]
  12. Kong, X.R.; Li, H.Q.; Wu, C. Dynamics of 1-dof and 2-dof energy sink with geometrically nonlinear damping: Application to vibration suppression. Nonlinear Dyn. 2018, 91, 733–754. [Google Scholar] [CrossRef]
  13. Rezaiee-Pajand, M.; Sarafrazi, S.R. Nonlinear dynamic structural analysis using dynamic relaxation with zero damping. Eng. Struct. 2011, 89, 1274–1285. [Google Scholar] [CrossRef]
  14. Zhai, J.; Li, J.; Wei, D.; Gao, P.; Yan, Y.; Han, Q. Vibration control of an aero pipeline system with active constraint layer damping treatment. Appl. Sci. 2019, 9, 2094. [Google Scholar] [CrossRef] [Green Version]
  15. Dong, G.X.; Zhang, Y.H.; Luo, Y.J.; Xie, S.; Zhang, X.N. Enhanced isolation performance of a high-static-low-dynamic stiffness isolator with geometric nonlinear damping. Nonlinear Dyn. 2018, 93, 2339–2356. [Google Scholar] [CrossRef]
  16. Gao, K.; Gao, W.; Wu, D. Nonlinear dynamic stability analysis of Euler-Bernoulli beam-columns with damping effects under thermal environment. Nonlinear Dyn. 2017, 90, 2423–2444. [Google Scholar] [CrossRef]
  17. Wang, G.P.; Rong, B.; Tao, L.; Rui, X.T. Riccati discrete time transfer matrix method for dynamic modeling and simulation of an underwater towed system. J. Appl. Mech. 2012, 79, 041014. [Google Scholar] [CrossRef]
  18. Rong, B.; Lu, K.; Rui, X.T.; Li, X.J.; Tao, L.; Wang, G.P. Nonlinear dynamics analysis of pipe conveying fluid by Riccati absolute nodal coordinate transfer matrix method. Nonlinear Dyn. 2018, 92, 699–708. [Google Scholar] [CrossRef]
  19. Rezaiee-Pajand, M.; Sarafrazi, S.R.; Rezaiee, H. Efficiency of dynamic relaxation methods in nonlinear analysis of truss and frame structures. Comput. Struct. 2012, 112, 295–310. [Google Scholar] [CrossRef]
  20. Yang, P.; Bai, H.B.; Xue, X.; Xiao, K.; Zhao, X. Vibration reliability characterization and damping capability of annular periodic metal rubber in the non-molding direction. Mech. Syst. Signal Process. 2019, 132, 622–639. [Google Scholar] [CrossRef]
  21. Chandrasekhar, K.; Rongong, J.; Cross, E. Mechanical behaviour of tangled metal wire devices. Mech. Syst. Signal Process. 2019, 118, 13–29. [Google Scholar] [CrossRef] [Green Version]
  22. Zhang, D.Y.; Fabrizio, S.; Ma, Y.H. Dynamic mechanical behavior of nickel-based superalloy metal rubber. Mater. Des. 2014, 56, 69–77. [Google Scholar] [CrossRef]
  23. Wang, H.; Rongong, J.A.; Tomlinson, G.R. Nonlinear static and dynamic properties of metal rubber dampers. Int. Conf. Noise Vib. Eng. 2010, 10, 1301–1316. [Google Scholar]
  24. Wu, K.N.; Bai, H.B.; Xue, X.; Li, T.; Li, M. Energy dissipation characteristics and dynamic modeling of the coated damping structure for metal rubber of bellows. Metals 2018, 8, 562. [Google Scholar] [CrossRef] [Green Version]
  25. Ma, Y.H.; Zhang, Q.C.; Zhang, D.Y.; Scarpa, F.; Liu, B.L.; Hong, J. A novel smart rotor support with shape memory alloy metal rubber for high temperatures and variable amplitude vibrations. Smart Mater. Struct. 2014, 23, 125016. [Google Scholar] [CrossRef]
  26. Kwon, S.C.; Jo, M.S.; Oh, H.U. Experimental Validation of Fly-Wheel Passive Launch and On-Orbit Vibration Isolation System by Using a Superelastic SMA Mesh Washer Isolator. Int. J. Aerosp. Eng. 2017, 2017, 1–16. [Google Scholar] [CrossRef] [Green Version]
  27. Xiao, K.; Bai, H.B.; Xue, X.; Wu, Y.W. Damping characteristics of metal rubber in the pipeline coating system. Shock. Vib. 2018, 42, 1–11. [Google Scholar] [CrossRef] [Green Version]
  28. Ma, Y.H.; Zhang, Q.C.; Zhang, D.Y.; Scarpa, F.; Gao, D.; Hong, J. Size-dependent mechanical behavior and boundary layer effects in entangled metallic wire material systems. J. Mater. Sci. 2017, 52, 3741–3756. [Google Scholar] [CrossRef] [Green Version]
  29. Catania, G.; Strozzi, M. Damping oriented design of thin-walled mechanical components by means of multi-layer coating technology. Coatings 2018, 8, 73. [Google Scholar] [CrossRef] [Green Version]
  30. Hou, J.F.; Bai, H.B.; Li, D.W. Damping capacity measurement of elastic porous wire-mesh material in wide temperature range. J. Mater. Process. Technol. 2008, 206, 412–418. [Google Scholar] [CrossRef]
  31. Zhang, D.Y.; Fabrizio, S.; Ma, Y.H.; Boba, K.; Hong, J.; Lu, H.W. Compression mechanics of nickel-based superalloy metal rubber. Mater. Sci. Eng. A 2013, 580, 305–312. [Google Scholar] [CrossRef]
Figure 1. Physical model of vibration isolation system for the two-point and symmetrically supported pipeline bracket.
Figure 1. Physical model of vibration isolation system for the two-point and symmetrically supported pipeline bracket.
Symmetry 11 01479 g001
Figure 2. The coordinate of pipeline B.
Figure 2. The coordinate of pipeline B.
Symmetry 11 01479 g002
Figure 3. The experimental bench of coated pipeline system.
Figure 3. The experimental bench of coated pipeline system.
Symmetry 11 01479 g003
Figure 4. The response curves of displacement and force.
Figure 4. The response curves of displacement and force.
Symmetry 11 01479 g004
Figure 5. The experimental measured force transfer rate-frequency curves of coated pipeline damping system.
Figure 5. The experimental measured force transfer rate-frequency curves of coated pipeline damping system.
Symmetry 11 01479 g005
Figure 6. The evolution of force transfer rate and loss factor with excitation level.
Figure 6. The evolution of force transfer rate and loss factor with excitation level.
Symmetry 11 01479 g006
Figure 7. The force transfer rate-frequency curves.
Figure 7. The force transfer rate-frequency curves.
Symmetry 11 01479 g007
Table 1. The energy dissipation characteristics of the coated pipeline system.
Table 1. The energy dissipation characteristics of the coated pipeline system.
Frequency (Hz)Energy Dissipation ΔW (N·mm)Maximum Elastic Potential Energy W (N·mm)Loss Factor η
151.0161.1360.142
200.8220.7520.174
Table 2. Energy dissipation characteristics at different pre-tightening conditions.
Table 2. Energy dissipation characteristics at different pre-tightening conditions.
Energy Dissipation CharacteristicsPre-Tightening Conditions (mm)
0.5 23
Natural Frequency ωn (Hz)13.99614.42014.538
Peak of the Force Transfer Rate TAm17.68219.02520.182
Structural Loss Factor η0.05330.04540.0435

Share and Cite

MDPI and ACS Style

Xue, X.; Ruan, S.; Li, A.; Bai, H.; Xiao, K. Nonlinear Dynamic Modelling of Two-Point and Symmetrically Supported Pipeline Brackets with Elastic-Porous Metal Rubber Damper. Symmetry 2019, 11, 1479. https://doi.org/10.3390/sym11121479

AMA Style

Xue X, Ruan S, Li A, Bai H, Xiao K. Nonlinear Dynamic Modelling of Two-Point and Symmetrically Supported Pipeline Brackets with Elastic-Porous Metal Rubber Damper. Symmetry. 2019; 11(12):1479. https://doi.org/10.3390/sym11121479

Chicago/Turabian Style

Xue, Xin, Shixin Ruan, Angxi Li, Hongbai Bai, and Kun Xiao. 2019. "Nonlinear Dynamic Modelling of Two-Point and Symmetrically Supported Pipeline Brackets with Elastic-Porous Metal Rubber Damper" Symmetry 11, no. 12: 1479. https://doi.org/10.3390/sym11121479

APA Style

Xue, X., Ruan, S., Li, A., Bai, H., & Xiao, K. (2019). Nonlinear Dynamic Modelling of Two-Point and Symmetrically Supported Pipeline Brackets with Elastic-Porous Metal Rubber Damper. Symmetry, 11(12), 1479. https://doi.org/10.3390/sym11121479

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