# Instantaneous Characteristics of Nonlinear Torsion Pendulum and Its Application in Parameter Estimation of Nonlinear System

Authors: Yan Zhao; Baofeng Zhang; Fangfang Han; Huan Tian; Xiao Yu; Junchao Zhu
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

The nonlinear model of torsion pendulum is presented by considering the nonlinear damping force and nonlinear restoring force. The analytic solution of the nonlinear model is calculated to analyze the relationship between the characteristics of torsion pendulum and the nonlinear factors. The instantaneous characteristics of nonlinear torsion pendulum are analyzed by instantaneous undamped natural frequency and instantaneous damping coefficient. The instantaneous characteristics can be used for the parameter estimation of nonlinear torsion pendulum system. The nonlinear characteristics of the torsion pendulum are validated by the torsion pendulum based on the air-hovered turntable. The parameter estimation method based on the instantaneous characteristics is validated by the moment of inertia measurement system based on the torsion pendulum.


---

## Body

## 1. Introduction

The torsion pendulum system has been widely used to measure the moment of inertia of mechanical systems such as missile and engine [1–4].The torsion pendulum model mainly contains three factors: moment of inertia, damping force, and restoring force. With the increasing complexity of the torsion pendulum, the linear model has not been able to accurately describe the actual torsion pendulum. For the large measured object with complex shape, the air damping is a nonlinear force related to linear terms and square terms of speed [5]. When the torsion bar is alloy or composite material, the restoring force of the torsion bar is a nonlinear force related to linear terms and cubic terms of torsional angle [6].The linear model cannot accurately describe the complex torsion pendulum system. So it is needed to establish the torsion pendulum model based on the nonlinear dynamic system. Because the nonlinear system and the linear system have significantly different characteristics, it is necessary to analyze the influence of nonlinear factors on the characteristics of the torsion pendulum.For the quantitative analysis of the nonlinear torsion pendulum, the approximate analytic solution of the nonlinear differential equation of the torsion pendulum is calculated. The solving methods of the nonlinear differential equation mainly include small parameter method, multiple scales method, evolutionary method, and average method [7–10]. The average method regards the amplitude and phase of the solution as slowly changing with time and uses the derivative of amplitude and phase to compute the approximate analytic solution of the nonlinear differential equation. This method is of high computational accuracy and is suitable for solving nonlinear dynamical model.The logarithmic decrement and vibration period are the important parameters of the torsion pendulum. The logarithmic decrement refers to the ratio of the amplitude of adjacent periods. For the linear torsion pendulum, the logarithmic decrement and vibration period are constant. But for the nonlinear torsion pendulum, the logarithmic decrement and vibration period are time-varying which are associated with nonlinear factors.The instantaneous characteristics of the nonlinear torsion pendulum can be described by instantaneous undamped natural frequency and instantaneous damping coefficient which are computed by the instantaneous amplitude and instantaneous frequency [11, 12].To calculate the instantaneous frequency, many methods are developed such as short-time Fourier transformation [13, 14], wavelet transformation [15, 16], Wigner-Vile Distribution [17, 18], and Hilbert transform [19]. Reference [11] presented the construction method of analytic signal based on Hilbert transform, which is used for the calculation of instantaneous frequency.Considering the influence of nonlinear damping force and nonlinear restoring force, the torsion pendulum model based on the nonlinear dynamic system is established. The analytic solution of the nonlinear model is used to analyze the relationship between the characteristics of torsion pendulum and the nonlinear factors. The instantaneous characteristics of the nonlinear system is analyzed by the instantaneous undamped natural frequency and instantaneous damping coefficient. At last, the nonlinear torsion pendulum experimental system is used to verify the presented model and method.


## 2. Nonlinear Model of Torsion Pendulum

The modeling methods based on the differential equation include Newton second law, Darren Bell principle method, and Lagrange equation method. Among them, the Lagrange equation method is based on the energy principle and it only needs to analyze the main power of the system. Therefore, the modeling method based on the Lagrange equation is more suitable for the complex system.As shown in Figure1, the torsion pendulum system consists of torsion bar, turntable, and test object. The test object is placed on the turntable that is connected to the torsion bar. With impulse excitation, the test object and turntable make torsion pendulum movement together while the restoring force is generated by the torsion bar. The torsion pendulum system in Figure 1 has one freedom degree and the torsion angle θ is regarded as the generalized coordinate.Figure 1
Schematic diagram of torsion pendulum.The second Lagrange equation can be expressed as(1)ddt∂L∂θ˙-∂L∂θ=Q,where L is the Lagrange function, L=T-V; T is the kinetic energy of the system; V is the potential energy of the system; Q is the generalized force of the system; and t is the time.The kinetic energy of the system is(2)T=12Jω2=12Jθ˙2,where J is the moment of inertia with respect to the rotation axis; ω is the angular speed of the torsion pendulum; θ˙ is the derivative of the torsion angle θ.The potential energy of the system is(3)V=∫0θK1θ+K3θ3dθ=12K1θ2+14K3θ4,where K1 is the linear stiffness coefficient of the torsion bar; K3 is the nonlinear stiffness coefficient of the torsion bar; when K3≠0, the restoring force of the torsion bar is a nonlinear force.The generalized force of the system is(4)Q=-C1θ˙-C2θ˙θ˙,where C1 is the linear damping coefficient; C2 is the nonlinear damping coefficient.Substituting (2)–(4) into (1), we can get the nonlinear model of the torsion pendulum:(5)Jθ¨+C1θ˙+C2θ˙θ˙+K1θ+K3θ3=0.Each side of (5) is divided by J; the normalization model is expressed as(6)θ¨+c1θ˙+c2θ˙θ˙+k1θ+k3θ3=0,where c1=C1/J; c2=C2/J; k1=K1/J; and k3=K3/J.


## 3. Analytic Solution of Nonlinear Torsion Pendulum Model

In order to analyze the characteristics of the nonlinear torsion pendulum, the analytic solution of (5) is calculated by the average method. If the influence of nonlinear factors is weak, a small parameter will appear in the differential equation and the differential equation is typically expressed as(7)θ¨+ω02θ=εfθ,θ˙,where ε (ε≪1) is a small parameter; f(θ,θ˙) is the nonlinear function of θ.Scholar Rylov and Bogoliubov wrote the solution as [20](8)θ=atcos⁡ω0t+φt.According to (8), the amplitude a(t) and phase φ(t) are regarded as the functions of time. If the influence of the nonlinear factors is weak, the solution can be set as simple harmonic form.The first derivative of (8) with respect to time is expressed as(9)θ˙=a˙tcos⁡ω0t+φt-atω0sin⁡ω0t+φt-atφ˙tsin⁡ω0t+φt.Becausea(t) and φ(t) are the slow-variant functions of time, we can get the following equation:(10)a˙tcos⁡ω0t+φt-atφ˙tsin⁡ω0t+φt=0.The first derivative ofθ(t) is expressed as(11)θ˙=-atω0sin⁡ω0t+φt.The second derivative ofθ(t) is expressed as(12)θ¨=-atω02cos⁡ω0t+φt-a˙tω0sin⁡ω0t+φt-atω0φ˙tcos⁡ω0t+φt.Substituting (8), (11), and (12) into (7) to obtain(13)-a˙tω0sin⁡ψ-atω0φ˙tcos⁡ψ=εfacos⁡ψ,-aω0sin⁡ψ,where ψ=ω0t+φ(t).According to (10) and (13), the first-order equation system is expressed as(14)a˙t=-εω0facos⁡ψ,-aω0sin⁡ψsin⁡ψφ˙t=-εatω0facos⁡ψ,-aω0sin⁡ψcos⁡ψ.Becausea(t) and φ(t) are the slow-variant functions of time, the right side of (14) is the periodic function of ψ and they can be expanded as Fourier series. The first term of Fourier series is taken as the approximate value of a˙(t) and φ˙(t):(15)a˙t≈-ε2πω0∫02πfacos⁡ψ,-aω0sin⁡ψsin⁡ψdψ(16)φ˙t≈-ε2πaω0∫02πfacos⁡ψ,-aω0sin⁡ψcos⁡ψdψ.If the initial condition is given, the value ofa(t) and φ(t) can be obtained by the integration of (15) and (16).If (5) is expressed as the form of (7), then(17)ω02=K1J(18)εfθ,θ˙=-C1Jθ˙-C2Jθ˙θ˙-K3Jθ3.Substitute (18) into (15) to obtain(19)a˙t=-ε2πω0∫02πfacos⁡ψ,-aω0sin⁡ψsin⁡ψdψ=-ε2πJC1aπ+83C2a2ω0.Substituting (18) into (16) to obtain(20)φ˙t=-ε2πaω0∫02πfacos⁡ψ,-aω0sin⁡ψcos⁡ψdψ=-ε3K3a28ω0J.Whent=0, the initial value of a(t) is a0. By the integration of (19), we can obtain(21)at=3a0πC1e-εC1/2Jt8a0C2ω0+3πC1-8a0C2ω0e-εC1/2Jt.Whent=0, the initial value of φ(t) is φ0. By the integration of (20), we can obtain(22)φt=-ε3K3a28ω0Jt+φ0.Substitute (21) and (22) into (8) to obtain the angular displacement θ(t):(23)θt=atcos⁡ω0t+φt=3a0πC1e-εC1/2Jt8a0C2ω0+3πC1-8a0C2ω0e-εC1/2Jtcos⁡ω0-ε3K3a2t8JK1t+φ0.According to (21), it is known that the amplitude of the torsion pendulum is related to the damping coefficients C1 and C2. According to (22), the period of the torsion pendulum is related to the stiffness coefficients of the torsion bars K1 and K3.


## 4. Instantaneous Characteristics of Torsion Pendulum

### 4.1. Construction Method of Analytic Signal

If the angular displacement signal of the torsion pendulum isθ(t), its Hilbert transform θ~(t) is defined as follows:(24)θ~t=Hθt=1π∫-∞+∞θτt-τdτ.The angular displacement signalθ(t) and its Hilbert transform θ~(t) can form the analytic signal β(t):(25)βt=θt+jθ~t=AtejΦt,where instantaneous phase(26)Φt=arctan⁡θ~tθt.instantaneous angular frequency(27)ωt=dΦdt,and instantaneous envelope(28)At=θ2t+θ~2t.The first derivative of analytic signalβ(t) is as follows:(29)β˙=βA˙tAt+jωt.The second derivative of analytic signalβ(t) is as follows: (30)β¨=βA¨A-ω2+2jA˙ωA+jω˙.The calculation formula of the instantaneous angular frequency is as follows:(31)ωt=Imβ˙tβt=θtθ~˙t-θ˙tθ~tθ2t+θ~2t.The calculation formula of the first derivative of instantaneous envelope is as follows:(32)A˙t=AtReβ˙tβt=θtθ˙t+θ~tθ~˙tθ2t+θ~2t.


### 4.2. Analysis Method Based on Hilbert Transform

The nonlinear model of the torsion pendulum can be expressed as follows:(33)θ¨+∑j=1ncjθ˙j·sign⁡θ˙+∑i=1mkiθi·sign⁡θ=0,where n is the maximum order of the angular velocity and m is the maximum order of the angular displacement.For convenience, (33) can be rewritten as [11](34)θ¨+2h0tθ˙+ω02tθ=0.Ifθ(t) and h0(t) are signals with nonoverlapping spectrum, we can use Hilbert transform for both sides of (34), multiplying each side of the obtained new equation by j and adding it to the corresponding side of (34)(35)β¨+2h0tβ˙+ω02tβ=0.Substitute (29) and (30) into (35) to obtain (36)βA¨A-ω2+ω02+2h0A˙A+j2A˙Aω+ω˙+2h0ω=0.Solving (36) we can obtain the expressions of instantaneous undamped natural frequency ω0(t) and instantaneous damping coefficient h0(t) [11]: (37)ω0t=ω2-A¨A+2A˙2A2+A˙ω˙Aω1/2(38)h0t=-A˙A-ω˙2ω.The instantaneous characteristicsω0(t) and h0(t) can be used for the parameter estimation of the nonlinear torsion pendulum system. The relationship between the parameters of (6) and ω0(t) can be expressed as [21](39)ω0A=∑i=1m2πΓi/2+1Γi+1/2+1ki·Ai-1t1/2.The parametersk1 and k3 of (6) can be estimated by the curve fitting of ω0(t). The moment of inertia can be calculated by J=K1/k1, if the value of K1 is known.The relationship between the parameters of (6) and h0(t) can be expressed as [21](40)h0Aω=∑j=0n1πΓj/2+1Γj+1/2+1cj·Atωtj-1.The parametersc1 and c2 of (6) can be estimated by the curve fitting of h0(t).


## 4.1. Construction Method of Analytic Signal

If the angular displacement signal of the torsion pendulum isθ(t), its Hilbert transform θ~(t) is defined as follows:(24)θ~t=Hθt=1π∫-∞+∞θτt-τdτ.The angular displacement signalθ(t) and its Hilbert transform θ~(t) can form the analytic signal β(t):(25)βt=θt+jθ~t=AtejΦt,where instantaneous phase(26)Φt=arctan⁡θ~tθt.instantaneous angular frequency(27)ωt=dΦdt,and instantaneous envelope(28)At=θ2t+θ~2t.The first derivative of analytic signalβ(t) is as follows:(29)β˙=βA˙tAt+jωt.The second derivative of analytic signalβ(t) is as follows: (30)β¨=βA¨A-ω2+2jA˙ωA+jω˙.The calculation formula of the instantaneous angular frequency is as follows:(31)ωt=Imβ˙tβt=θtθ~˙t-θ˙tθ~tθ2t+θ~2t.The calculation formula of the first derivative of instantaneous envelope is as follows:(32)A˙t=AtReβ˙tβt=θtθ˙t+θ~tθ~˙tθ2t+θ~2t.


## 4.2. Analysis Method Based on Hilbert Transform

The nonlinear model of the torsion pendulum can be expressed as follows:(33)θ¨+∑j=1ncjθ˙j·sign⁡θ˙+∑i=1mkiθi·sign⁡θ=0,where n is the maximum order of the angular velocity and m is the maximum order of the angular displacement.For convenience, (33) can be rewritten as [11](34)θ¨+2h0tθ˙+ω02tθ=0.Ifθ(t) and h0(t) are signals with nonoverlapping spectrum, we can use Hilbert transform for both sides of (34), multiplying each side of the obtained new equation by j and adding it to the corresponding side of (34)(35)β¨+2h0tβ˙+ω02tβ=0.Substitute (29) and (30) into (35) to obtain (36)βA¨A-ω2+ω02+2h0A˙A+j2A˙Aω+ω˙+2h0ω=0.Solving (36) we can obtain the expressions of instantaneous undamped natural frequency ω0(t) and instantaneous damping coefficient h0(t) [11]: (37)ω0t=ω2-A¨A+2A˙2A2+A˙ω˙Aω1/2(38)h0t=-A˙A-ω˙2ω.The instantaneous characteristicsω0(t) and h0(t) can be used for the parameter estimation of the nonlinear torsion pendulum system. The relationship between the parameters of (6) and ω0(t) can be expressed as [21](39)ω0A=∑i=1m2πΓi/2+1Γi+1/2+1ki·Ai-1t1/2.The parametersk1 and k3 of (6) can be estimated by the curve fitting of ω0(t). The moment of inertia can be calculated by J=K1/k1, if the value of K1 is known.The relationship between the parameters of (6) and h0(t) can be expressed as [21](40)h0Aω=∑j=0n1πΓj/2+1Γj+1/2+1cj·Atωtj-1.The parametersc1 and c2 of (6) can be estimated by the curve fitting of h0(t).


## 5. Simulation Analysis

### 5.1. Influence of Nonlinear Damping Force

It is known from (21) that the amplitude of the torsion pendulum is related to the damping factor. To analyze the influence of nonlinear damping force, the logarithmic decrement of the angular displacement signals is compared under different damping conditions. By forth-order Runge-Kutta method, (5) is solved to obtain the simulation signals of angular displacement with the parameters in Table 1. The initial displacement is 0.06 rad.Table 1
Parameter value of simulation signals under different damping conditions.
J C 1 C 2 K 1 K 3 (kg⋅m2) (kg⋅m2⋅s−1) (kg⋅m2⋅rad−1) (N⋅m⋅rad−1) (N⋅m⋅rad−3) θ 1 ( t ) 1 0 0 5 0 θ 2 ( t ) 1 0.01 0 5 0 θ 3 ( t ) 1 0.01 0.1 5 0The logarithmic decrement refers to the ratio of the amplitude of adjacent periods, reflecting the degree of amplitude attenuation. The logarithmic decrement of the angular displacement simulation signals is shown in Figure2. Without the damping force, the logarithmic decrement of θ1(t) is zero, indicating that the amplitude of the angular displacement does not attenuate. With the linear damping force, the logarithmic decrement of θ2(t) is 0.014, indicating that the amplitude of the angular displacement attenuates at a constant rate. With the nonlinear damping force, the logarithmic decrement of θ3(t) decreases gradually, indicating that the degree of the amplitude attenuation is variable.Figure 2
Logarithmic decrement of simulation signals.


### 5.2. Influence of Nonlinear Restoring Force

It is known from (22) that the period of torsion pendulum is related to the stiffness coefficient of the torsion bar. To analyze the influence of the nonlinear restoring force, periods of angular displacement signals are compared under different stiffness coefficient conditions. By forth-order Runge-Kutta method, (5) is solved to obtain the simulation signals of angular displacement with parameters in Table 2. The initial displacement is 0.06 rad.Table 2
Parameter value of simulation signals under different stiffness coefficient conditions.
J C 1 C 2 K 1 K 3 (kg⋅m2) (kg⋅m2⋅s−1) (kg⋅m2⋅rad−1) (N⋅m⋅rad−1) (N⋅m⋅rad−3) θ 4 ( t ) 1 0.01 0 5 0 θ 5 ( t ) 1 0.01 0.1 5 20 θ 6 ( t ) 1 0.01 0.1 5 −20The periods of simulation signals under different stiffness coefficient conditions are shown in Figure3. When K3=0, the restoring force of the torsion pendulum is a linear force and the time interval of every period of θ4(t) is the same. When K3>0, the restoring force of the torsion pendulum is a nonlinear force and the time interval of every period of θ5(t) increases as the amplitude decreases. When K3<0, the restoring force of the torsion pendulum is a nonlinear force and the time interval of every period of θ6(t) decreases as the amplitude decreases. It can be seen that when there is a nonlinear restoring force in the torsion pendulum system, the period of the torsion pendulum will change with the amplitude.Figure 3
Periods of simulation signal.


### 5.3. Instantaneous Characteristics

The simulation parameter value of (5) is shown in Table 3. The initial angular displacement is 0.05 rad.Table 3
Simulation parameter value of torsion pendulum.
J C 1 C 2 K 1 K 3 (kg⋅m2) (kg⋅m2⋅s−1) (kg⋅m2⋅rad−1) (N⋅m⋅rad−1) (N⋅m⋅rad−3) θ 7 ( t ) 100 1 10 400 −600The simulation signal of the angular displacement is shown in Figure4.Figure 4
Simulation signal of angular displacement.According to (25), the analytic signal is formed by the angular displacement signal θ7(t) and its Hilbert transform. The three-dimensional stereogram of the analytic signal is shown in Figure 5(a). The projection of the analytic signal in each coordinate plane is shown in Figure 5(b). The angular displacement signal θ7(t) is shown in the plane consisting of real-axis and time-axis. The Hilbert transform of the angular displacement signal is shown in the plane consisting of imaginary-axis and time-axis. The instantaneous envelope obtained according to (28) is shown in Figure 4 (dashed line).Figure 5
Stereogram and projection graph of analytic signal.
(a)
Stereogram of analytic signal (b)
Projection graph of analytic signalAccording to (37), the instantaneous undamped natural frequency ω0(t) is calculated by the instantaneous envelope and instantaneous angular frequency, which is shown in Figure 6.Figure 6
Instantaneous undamped natural frequency.According to (38), the instantaneous damping coefficient h0(t) is calculated by the instantaneous envelope and instantaneous angular frequency, which is shown in Figure 7.Figure 7
Instantaneous damping coefficient.The instantaneous characteristicsω0(t) and h0(t) can be used for the parameter estimation of the nonlinear torsion pendulum system. According to (39) and (40), the normalized parameters of the torsion pendulum are estimated by the curve fitting of the curves in Figures 6 and 7. The estimation result is shown in Table 4. If the stiffness coefficient of the torsion bar K1 is known, the moment of inertia can be calculated by J=K1/k1.Table 4
Estimation result of nonlinear torsion pendulum system.
Standard value Estimated value c 1 (s−1) 0.0100 0.0100 c 2 (rad−1) 0.1000 0.1001 k 1 (rad−1⋅s−2) 4.0000 3.9999 k 3 (rad−3⋅s−2) −6.0000 −6.0043The estimation results under different SNR noise are shown in Table5.Table 5
Estimation results under different SNR noise.
SNR c 1 (s−1) c 2 (rad−1) k 1 (rad−1⋅s−2) k 3 (rad−3⋅s−2) 60 dB 0.0100 0.1000 4.0001 −6.0173 50 dB 0.0100 0.1002 3.9994 −5.9731 40 dB 0.0098 0.1032 3.9979 −5.9592


## 5.1. Influence of Nonlinear Damping Force

It is known from (21) that the amplitude of the torsion pendulum is related to the damping factor. To analyze the influence of nonlinear damping force, the logarithmic decrement of the angular displacement signals is compared under different damping conditions. By forth-order Runge-Kutta method, (5) is solved to obtain the simulation signals of angular displacement with the parameters in Table 1. The initial displacement is 0.06 rad.Table 1
Parameter value of simulation signals under different damping conditions.
J C 1 C 2 K 1 K 3 (kg⋅m2) (kg⋅m2⋅s−1) (kg⋅m2⋅rad−1) (N⋅m⋅rad−1) (N⋅m⋅rad−3) θ 1 ( t ) 1 0 0 5 0 θ 2 ( t ) 1 0.01 0 5 0 θ 3 ( t ) 1 0.01 0.1 5 0The logarithmic decrement refers to the ratio of the amplitude of adjacent periods, reflecting the degree of amplitude attenuation. The logarithmic decrement of the angular displacement simulation signals is shown in Figure2. Without the damping force, the logarithmic decrement of θ1(t) is zero, indicating that the amplitude of the angular displacement does not attenuate. With the linear damping force, the logarithmic decrement of θ2(t) is 0.014, indicating that the amplitude of the angular displacement attenuates at a constant rate. With the nonlinear damping force, the logarithmic decrement of θ3(t) decreases gradually, indicating that the degree of the amplitude attenuation is variable.Figure 2
Logarithmic decrement of simulation signals.


## 5.2. Influence of Nonlinear Restoring Force

It is known from (22) that the period of torsion pendulum is related to the stiffness coefficient of the torsion bar. To analyze the influence of the nonlinear restoring force, periods of angular displacement signals are compared under different stiffness coefficient conditions. By forth-order Runge-Kutta method, (5) is solved to obtain the simulation signals of angular displacement with parameters in Table 2. The initial displacement is 0.06 rad.Table 2
Parameter value of simulation signals under different stiffness coefficient conditions.
J C 1 C 2 K 1 K 3 (kg⋅m2) (kg⋅m2⋅s−1) (kg⋅m2⋅rad−1) (N⋅m⋅rad−1) (N⋅m⋅rad−3) θ 4 ( t ) 1 0.01 0 5 0 θ 5 ( t ) 1 0.01 0.1 5 20 θ 6 ( t ) 1 0.01 0.1 5 −20The periods of simulation signals under different stiffness coefficient conditions are shown in Figure3. When K3=0, the restoring force of the torsion pendulum is a linear force and the time interval of every period of θ4(t) is the same. When K3>0, the restoring force of the torsion pendulum is a nonlinear force and the time interval of every period of θ5(t) increases as the amplitude decreases. When K3<0, the restoring force of the torsion pendulum is a nonlinear force and the time interval of every period of θ6(t) decreases as the amplitude decreases. It can be seen that when there is a nonlinear restoring force in the torsion pendulum system, the period of the torsion pendulum will change with the amplitude.Figure 3
Periods of simulation signal.


## 5.3. Instantaneous Characteristics

The simulation parameter value of (5) is shown in Table 3. The initial angular displacement is 0.05 rad.Table 3
Simulation parameter value of torsion pendulum.
J C 1 C 2 K 1 K 3 (kg⋅m2) (kg⋅m2⋅s−1) (kg⋅m2⋅rad−1) (N⋅m⋅rad−1) (N⋅m⋅rad−3) θ 7 ( t ) 100 1 10 400 −600The simulation signal of the angular displacement is shown in Figure4.Figure 4
Simulation signal of angular displacement.According to (25), the analytic signal is formed by the angular displacement signal θ7(t) and its Hilbert transform. The three-dimensional stereogram of the analytic signal is shown in Figure 5(a). The projection of the analytic signal in each coordinate plane is shown in Figure 5(b). The angular displacement signal θ7(t) is shown in the plane consisting of real-axis and time-axis. The Hilbert transform of the angular displacement signal is shown in the plane consisting of imaginary-axis and time-axis. The instantaneous envelope obtained according to (28) is shown in Figure 4 (dashed line).Figure 5
Stereogram and projection graph of analytic signal.
(a)
Stereogram of analytic signal (b)
Projection graph of analytic signalAccording to (37), the instantaneous undamped natural frequency ω0(t) is calculated by the instantaneous envelope and instantaneous angular frequency, which is shown in Figure 6.Figure 6
Instantaneous undamped natural frequency.According to (38), the instantaneous damping coefficient h0(t) is calculated by the instantaneous envelope and instantaneous angular frequency, which is shown in Figure 7.Figure 7
Instantaneous damping coefficient.The instantaneous characteristicsω0(t) and h0(t) can be used for the parameter estimation of the nonlinear torsion pendulum system. According to (39) and (40), the normalized parameters of the torsion pendulum are estimated by the curve fitting of the curves in Figures 6 and 7. The estimation result is shown in Table 4. If the stiffness coefficient of the torsion bar K1 is known, the moment of inertia can be calculated by J=K1/k1.Table 4
Estimation result of nonlinear torsion pendulum system.
Standard value Estimated value c 1 (s−1) 0.0100 0.0100 c 2 (rad−1) 0.1000 0.1001 k 1 (rad−1⋅s−2) 4.0000 3.9999 k 3 (rad−3⋅s−2) −6.0000 −6.0043The estimation results under different SNR noise are shown in Table5.Table 5
Estimation results under different SNR noise.
SNR c 1 (s−1) c 2 (rad−1) k 1 (rad−1⋅s−2) k 3 (rad−3⋅s−2) 60 dB 0.0100 0.1000 4.0001 −6.0173 50 dB 0.0100 0.1002 3.9994 −5.9731 40 dB 0.0098 0.1032 3.9979 −5.9592


## 6. Experimental Verification

To verify the correctness of the proposed method, the moment inertia of the experimental sample was measured by the instantaneous characteristics of the torsion pendulum system based on air-hovered turntable. Before the moment of inertia measurement of the experimental sample, the parameters of the torsion pendulum system need to be calibrated, including the stiffness coefficient of the torsion bar and the moment of inertia of the turntable under the no-load condition. The angular displacement signal of the turntable is measured by an incremental encoder using a grating disk. The sampling frequency is 300 Hz.


### 6.1. Parameter Calibration of Torsion Pendulum

As shown in Figure8, the parameters of the torsion pendulum system are calibrated by the standard weight with known moment of inertia. The weight is made of iron cylinder with uniform density. Because of the uniform density, regular shape and high machining accuracy, the weight moment of inertia can be accurately calculated according to the external dimensions and quality of the weight. By the theoretical calculation, the total moment of inertia of weight A and weight B about the rotation axis is 57.3002 kg·m2.Figure 8
Calibration system of the torsion pendulum.The normalized parameters of the torsion pendulum model under the no-load condition can be estimated by the instantaneous characteristics of the torsion pendulum. If the moment of inertia of the torsion pendulum under the no-load condition isJ0, then(41)k11=K1J0k31=K3J0,where k11 and k31 are the normalized stiffness coefficients of the torsion bar under the no-load condition.As shown in Figure8, two identical standard weights are symmetrically placed on either side of the shaft. The normalized parameters of the torsion pendulum model under the standard weight condition can be estimated by the instantaneous characteristics of the torsion pendulum. If the moment of inertia of weight A and weight B is expressed as JA and JB, respectively, then(42)k12=K1J0+JA+JBk32=K3J0+JA+JB,where k12 and k32 are the normalized stiffness coefficients of the torsion bar under the weight condition.According to (41) and (42), the parameters of the torsion pendulum system can be expressed as(43)K1=k11·k12·JA+JBk11-k12K3=k31·k12·JA+JBk11-k12J0=k12·JA+JBk11-k12.The calibration results of the torsion pendulum system are shown in Table6.Table 6
Calibration results of torsion pendulum.
K 1/(N⋅m⋅rad−1) K 3/(N⋅m⋅rad−3) J 0/(kg⋅m2) ( 1 ) 467.0433 −2212.7707 94.6332 ( 2 ) 467.1702 −2211.7868 94.6616 ( 3 ) 467.2111 −2200.9289 94.6730


### 6.2. Experimental Result

As shown in Figure9, two identical samples with large outline size are symmetrically placed on the both sides of the rotation axis, which can enhance the air damping force. By the theoretical calculation, the total moment of inertia of the samples about the rotation axis is 53.9056 kg·m2. By the calibration of the standard weight, the stiffness coefficient of the torsion bar is K1 = 467.1702 N·m·rad−1, K3 = −2211.7868 N·m·rad−3, and the moment of inertia of the turntable is 94.6616 kg·m2.Figure 9
Moment of inertia measurement system.The angular displacement signal of the turntable is shown in Figure10. To analyze the influence of nonlinear damping force and nonlinear restoring force, it is necessary to calculate the period and logarithmic decrement of the angular displacement signal.Figure 10
Experimental angular signal of turntable.The logarithmic decrement of the experimental signal is shown in Figure11. It is known from Figure 11 that the logarithmic decrement of the experimental signal decreases gradually, which is consistent with the simulation result of the nonlinear model of the torsion pendulum.Figure 11
Logarithmic decrement of experimental signal.The time intervals of the experimental signal are shown in Figure12. It is known from Figure 12 that the periods of the experimental signal decrease gradually, which is consistent with the simulation result of the nonlinear model of the torsion pendulum.Figure 12
Periods of experimental signal.The instantaneous undamped natural frequencyω0(t) and instantaneous damping coefficient h0(t) of the experimental system are shown in Figures 13 and 14.Figure 13
Instantaneous undamped natural frequency.Figure 14
Instantaneous damping coefficient.The parameter estimation method based on the instantaneous characteristics can be validated by the measurement of the moment of inertia of the samples. According to (39), the measured value of k1 is 3.1368 by the curve fitting of the instantaneous undamped natural frequency in Figure 13. As the stiffness coefficient K1 is known, the measured value of the total moment of inertia of samples and turntable is 148.9321 kg·m2. As the moment of inertia of the turntable is known, the measured value of the total moment of inertia of samples is 54.2705 kg·m2.The moment of inertia of the samples is measured by 5 times. As shown in the result in Table7, the relative error of the moment of inertia of the samples is within 0.7%, which verifies the effectiveness of the parameter estimation method based on the instantaneous characteristics.Table 7
Relative error of moment of inertia.
Estimated value Relative error ( 1 ) 54.2722 kg⋅m2 0.68% ( 2 ) 54.2236 kg⋅m2 0.59% ( 3 ) 54.2021 kg⋅m2 0.55% ( 4 ) 54.2398 kg⋅m2 0.62% ( 5 ) 54.2129 kg⋅m2 0.57%The main factors affecting the measurement results include the calibration of the torsion pendulum, the theoretical value of the moment of inertia of the samples, the environmental noise, and the sampling accuracy. The method of improving measurement accuracy will be carried out in future research.


## 6.1. Parameter Calibration of Torsion Pendulum

As shown in Figure8, the parameters of the torsion pendulum system are calibrated by the standard weight with known moment of inertia. The weight is made of iron cylinder with uniform density. Because of the uniform density, regular shape and high machining accuracy, the weight moment of inertia can be accurately calculated according to the external dimensions and quality of the weight. By the theoretical calculation, the total moment of inertia of weight A and weight B about the rotation axis is 57.3002 kg·m2.Figure 8
Calibration system of the torsion pendulum.The normalized parameters of the torsion pendulum model under the no-load condition can be estimated by the instantaneous characteristics of the torsion pendulum. If the moment of inertia of the torsion pendulum under the no-load condition isJ0, then(41)k11=K1J0k31=K3J0,where k11 and k31 are the normalized stiffness coefficients of the torsion bar under the no-load condition.As shown in Figure8, two identical standard weights are symmetrically placed on either side of the shaft. The normalized parameters of the torsion pendulum model under the standard weight condition can be estimated by the instantaneous characteristics of the torsion pendulum. If the moment of inertia of weight A and weight B is expressed as JA and JB, respectively, then(42)k12=K1J0+JA+JBk32=K3J0+JA+JB,where k12 and k32 are the normalized stiffness coefficients of the torsion bar under the weight condition.According to (41) and (42), the parameters of the torsion pendulum system can be expressed as(43)K1=k11·k12·JA+JBk11-k12K3=k31·k12·JA+JBk11-k12J0=k12·JA+JBk11-k12.The calibration results of the torsion pendulum system are shown in Table6.Table 6
Calibration results of torsion pendulum.
K 1/(N⋅m⋅rad−1) K 3/(N⋅m⋅rad−3) J 0/(kg⋅m2) ( 1 ) 467.0433 −2212.7707 94.6332 ( 2 ) 467.1702 −2211.7868 94.6616 ( 3 ) 467.2111 −2200.9289 94.6730


## 6.2. Experimental Result

As shown in Figure9, two identical samples with large outline size are symmetrically placed on the both sides of the rotation axis, which can enhance the air damping force. By the theoretical calculation, the total moment of inertia of the samples about the rotation axis is 53.9056 kg·m2. By the calibration of the standard weight, the stiffness coefficient of the torsion bar is K1 = 467.1702 N·m·rad−1, K3 = −2211.7868 N·m·rad−3, and the moment of inertia of the turntable is 94.6616 kg·m2.Figure 9
Moment of inertia measurement system.The angular displacement signal of the turntable is shown in Figure10. To analyze the influence of nonlinear damping force and nonlinear restoring force, it is necessary to calculate the period and logarithmic decrement of the angular displacement signal.Figure 10
Experimental angular signal of turntable.The logarithmic decrement of the experimental signal is shown in Figure11. It is known from Figure 11 that the logarithmic decrement of the experimental signal decreases gradually, which is consistent with the simulation result of the nonlinear model of the torsion pendulum.Figure 11
Logarithmic decrement of experimental signal.The time intervals of the experimental signal are shown in Figure12. It is known from Figure 12 that the periods of the experimental signal decrease gradually, which is consistent with the simulation result of the nonlinear model of the torsion pendulum.Figure 12
Periods of experimental signal.The instantaneous undamped natural frequencyω0(t) and instantaneous damping coefficient h0(t) of the experimental system are shown in Figures 13 and 14.Figure 13
Instantaneous undamped natural frequency.Figure 14
Instantaneous damping coefficient.The parameter estimation method based on the instantaneous characteristics can be validated by the measurement of the moment of inertia of the samples. According to (39), the measured value of k1 is 3.1368 by the curve fitting of the instantaneous undamped natural frequency in Figure 13. As the stiffness coefficient K1 is known, the measured value of the total moment of inertia of samples and turntable is 148.9321 kg·m2. As the moment of inertia of the turntable is known, the measured value of the total moment of inertia of samples is 54.2705 kg·m2.The moment of inertia of the samples is measured by 5 times. As shown in the result in Table7, the relative error of the moment of inertia of the samples is within 0.7%, which verifies the effectiveness of the parameter estimation method based on the instantaneous characteristics.Table 7
Relative error of moment of inertia.
Estimated value Relative error ( 1 ) 54.2722 kg⋅m2 0.68% ( 2 ) 54.2236 kg⋅m2 0.59% ( 3 ) 54.2021 kg⋅m2 0.55% ( 4 ) 54.2398 kg⋅m2 0.62% ( 5 ) 54.2129 kg⋅m2 0.57%The main factors affecting the measurement results include the calibration of the torsion pendulum, the theoretical value of the moment of inertia of the samples, the environmental noise, and the sampling accuracy. The method of improving measurement accuracy will be carried out in future research.


## 7. Conclusion

Considering the influence of nonlinear damping force and nonlinear restoring force, the nonlinear torsion pendulum model is established. The analytic solution of the nonlinear model is solved by the average method. The simulation results show that the logarithmic decrement of the torsion pendulum varies by the nonlinear damping force and the period of the torsion pendulum varies by the nonlinear restoring force. The instantaneous characteristics can be used for the parameter estimation of the nonlinear torsion pendulum system. The torsion pendulum system based on the air-hovered turntable verifies the nonlinear characteristics of the torsion pendulum.



---

*Source: 7858403-2018-05-10.xml*