Figure 1.
Wedge friction damper.
Citation: Moazameh Golshani, Mohsen Basiri, Mohammad Shabani, Iraj Aghaei, Majid Asadi-Shekaari. Effects of erythropoietin on bile duct ligation-induced neuro-inflammation in male rats[J]. AIMS Neuroscience, 2019, 6(2): 43-53. doi: 10.3934/Neuroscience.2019.2.43
| [1] | Dike Hu, Hua Wang, Zijie Huang, Wangqiang Xiao . Research on impact reduction of flexible boundary particle damping honeycomb plate based on discrete element multi body dynamics coupling. Electronic Research Archive, 2023, 31(10): 6303-6326. doi: 10.3934/era.2023319 |
| [2] | Rafael Ortega . Variations on Lyapunov's stability criterion and periodic prey-predator systems. Electronic Research Archive, 2021, 29(6): 3995-4008. doi: 10.3934/era.2021069 |
| [3] | Qianqian Zhang, Mingye Mu, Heyuan Ji, Qiushi Wang, Xingyu Wang . An adaptive type-2 fuzzy sliding mode tracking controller for a robotic manipulator. Electronic Research Archive, 2023, 31(7): 3791-3813. doi: 10.3934/era.2023193 |
| [4] | Fang Wang, Weiguo Li, Wendi Bao, Li Liu . Greedy randomized and maximal weighted residual Kaczmarz methods with oblique projection. Electronic Research Archive, 2022, 30(4): 1158-1186. doi: 10.3934/era.2022062 |
| [5] | Shengming Hu, Yongfei Lu, Xuanchi Liu, Cheng Huang, Zhou Wang, Lei Huang, Weihang Zhang, Xiaoyang Li . Stability prediction of circular sliding failure soil slopes based on a genetic algorithm optimization of random forest algorithm. Electronic Research Archive, 2024, 32(11): 6120-6139. doi: 10.3934/era.2024284 |
| [6] | Wanxun Jia, Ling Li, Haoyan Zhang, Gengxiang Wang, Yang Liu . A novel nonlinear viscous contact model with a Newtonian fluid-filled dashpot applied for impact behavior in particle systems. Electronic Research Archive, 2025, 33(5): 3135-3157. doi: 10.3934/era.2025137 |
| [7] | Xingting Geng, Jianwen Feng, Yi Zhao, Na Li, Jingyi Wang . Fixed-time synchronization of nonlinear coupled memristive neural networks with time delays via sliding-mode control. Electronic Research Archive, 2023, 31(6): 3291-3308. doi: 10.3934/era.2023166 |
| [8] | Yang Jin, Wanxiang Li, Hongbing Zhang . Transition characteristics of the dynamic behavior of a vehicle wheel-rail vibro-impact system. Electronic Research Archive, 2023, 31(11): 7040-7060. doi: 10.3934/era.2023357 |
| [9] | Dae-Woong Lee . Homotopy comultiplications on the localization of a wedge of spheres and Moore spaces. Electronic Research Archive, 2022, 30(6): 2033-2053. doi: 10.3934/era.2022103 |
| [10] | Yichao Shao, Hengguo Yu, Chenglei Jin, Jingzhe Fang, Min Zhao . Dynamics analysis of a predator-prey model with Allee effect and harvesting effort. Electronic Research Archive, 2024, 32(10): 5682-5716. doi: 10.3934/era.2024263 |
Contact collision is a common phenomenon in engineering practice. With the deepening of scholars' research, the frontal collision model cannot describe the oblique collision system with tilt angle more accurately. Therefore, the oblique collision system with collision angle should be studied in depth. In addition, friction is inevitable in the process of oblique collision. The vibration caused by friction affects the stability of the system, accelerates the wear failure of the parts, and has potential safety hazards. Duan and Ding [1] discussed the influence of collision parameters such as mass clearance and contact surface inclination on the vibration reduction effect of turbine shroud blade group system through phase trajectory diagram, bifurcation diagram, and Poincare section diagram. Ling et al. [2] established a mechanical model of oblique collision between freight cars and trains at railway crossings and studied the dynamic response and derailment mechanism of trains in this scenario using the multi-body dynamics method. Dong and Ding [3] established the oblique collision model and equation of the whole ring bladed disk system with a shroud and analyzed the dynamic response of the blade when the gap between the shrouds and the stiffness of the blade are misaligned. In [4,5,6,7,8], the mechanical model with oblique collision was simplified for the actual engineering model, and the dynamic behaviors such as oblique collision motion and bifurcation phenomenon of the system in engineering practice were studied. Considering the dynamic friction of LuGre, Liu et al. [9] established a friction model of blade crown collision, and studied the bifurcation characteristics of contact collision and friction motion of shrouded blades. For the single-degree-of-freedom system with improved LuGre dynamic friction, Saha et al. [10] analyzed the sticking phenomenon of the system through the phase diagram, and compared and analyzed the dynamic characteristics of the system under LuGre and improved LuGre dynamic friction. Zhang et al. [11] introduced the Dankowicz dynamic friction model into the single-degree-of-freedom vibration system, analyzed the motion state and judgment conditions, and studied the friction-induced system vibration and the dynamic behavior of the system under the influence of two parameters.
The wedge friction damper is one of the key components of the three-piece bogie [12], which has an important influence on the dynamic performance of the vehicle. Li [13] analyzed the vibration characteristics of a 100 t heavy-duty truck with a wedge friction damper and found the influence of the relative friction coefficient of the damping system on the stability of the truck. Based on vehicle system dynamics, Song et al. [14] used SIMPACK software to establish a nonlinear dynamic model of six-axle flat car and analyzed the influence and variation of different friction angles on the dynamics of three-axle bogie. Li et al. [12] used SIMPACK software to establish the mechanical model of the wedge friction damper and analyzed the spatial force and vibration reduction performance of the model. Liu [15] introduced the structure of wedge friction damper for three-piece bogie of freight train bogie, analyzed the design value of the relative friction coefficient of the product, and gave the value suggestion.
On the basis of [12,13,14,15], aiming at the wedge friction damper for freight train bogie, in order to be closer to the actual working conditions and describe the friction-induced vibration characteristics more comprehensively, the dynamic friction between the wedge and the side frame is considered, and a three-degree-of-freedom oblique collision vibration system with Dankowicz dynamic friction and gap is simplified. The motion equation of the system is established, and the motion state of the system and its judgment conditions are analyzed. Based on C language programming and variable step size fourth-order Runge-Kutta algorithm, the numerical simulation analyzes the evolution process of the system motion and the frequency range of the viscous and flutter phenomena when the external excitation frequency changes, which can provide a certain reference for the selection of system structure parameters in engineering practice.
The following Figure 1 shows a wedge friction damper. The numbers 1–9 in the figure are bolster, right wedge, left wedge, intermediate bolster spring, right bolster spring, left bolster spring, right sub-wedge spring, left sub-wedge spring, and side frame.
In order to be closer to the actual working conditions and consider the influence of vibration reduction and lubrication, the bolster 1 is simplified into a mass block $ {M}_{1} $; the wedge blocks 2 and 3 are simplified into block $ {M}_{2} $ and block $ {M}_{3} $, respectively. The damping spring 7 and 8 connected by the wedge block 2 and 3 are simplified to the nonlinear spring $ {K}_{2} $ and the nonlinear spring $ {K}_{3} $, respectively. The lubricating oil between the bolster 1 and the wedge block is simplified as a nonlinear damping sum $ {C}_{2} $ and $ {C}_{3} $, respectively. The bolster spring 4, 5, 6 is simplified into a nonlinear spring $ {K}_{1} $; the air damping of the bolster is simplified to linear damping $ {C}_{1} $; the contact surfaces between the wedge block 2, 3 and the side frame 9 are A and B surfaces, and there is sliding friction force between the contact surfaces; and a dynamic model of the wedge damper with dynamic friction and gap is obtained, as shown in Figure 2.
Between two adjacent collisions, the differential equation of system motion is:
|
$ \left\{M1¨X1+F1−F2−F3−ˉF2−ˉF3+C1˙X1=Fcos(ΩT+φ)M2¨X2+F2+ˉF2=FdM3¨X3+F3+ˉF3=Fd \right. $
|
(1) |
In Eq (1): $ {F}_{d} $ is the dynamic friction generated between the contact surface of the bolster 1 and the side frame 9 during the movement of the wedge friction damper; $ {F}_{i} $ is the nonlinear spring force of the damping spring $ {K}_{i} $, where $ E $ is a small parameter; $ {\bar{F}}_{i} $ is a nonlinear damping force [16].
| $ {F}_{i} = {K}_{i}{X}_{i}+E{K}_{i}{X}_{i}^{3} $ | (2) |
In Eq (2), $ i = 1, \mathrm{ }2, \mathrm{ }3 $
| $ {\bar{F}}_{i} = {C}_{i}[1+Hsign({\dot{X}}_{i}-{\dot{X}}_{1}\left)\right]{({\dot{X}}_{i}-{\dot{X}}_{1})}^{\frac{1}{3}} $ | (3) |
In Eq (3), $ i = 2, \mathrm{ }3 $, $ \mathrm{H} $ is the asymmetry coefficient, which represents the unequal degree of damping force between the recovery stroke and the compression stroke of the shock absorber, where
|
$ \text{sign(}{\dot{X}}_{i}-{\dot{X}}_{1}) = \left\{1,˙Xi−˙X1>0−1,˙Xi−˙X1<0 \right. $
|
(4) |
Select the initial displacement $ {\mathrm{L}}_{1} = {\mathrm{x}}_{1} $(0) as the length scale, construct the time scale $ {T}_{1} = \sqrt{\frac{{M}_{1}}{{K}_{1}}} $, and take the dimensionless parameter as:
| $ {m}_{2} = \frac{{M}_{2}}{{M}_{1}} , {m}_{3} = \frac{{M}_{3}}{{M}_{1}} , {k}_{1} = \frac{{K}_{1}{T}_{1}^{2}}{{M}_{1}} , {k}_{2} = \frac{{K}_{2}{T}_{1}^{2}}{{M}_{1}} , {k}_{3} = \frac{{K}_{3}{T}_{1}^{2}}{{M}_{1}} , {\xi }_{1} = \frac{{C}_{1}{T}_{1}}{{M}_{1}} , {\xi }_{2} = \frac{{C}_{2}{T}_{1}^{\frac{5}{3}}}{{M}_{1}{L}_{1}^{\frac{2}{3}}} , {\xi }_{3} = \frac{{\mathrm{C}}_{3}{\mathrm{T}}_{1}^{\frac{5}{3}}}{{M}_{1}{\mathrm{L}}_{1}^{\frac{2}{3}}} , \eta = \frac{{L}_{1}}{{T}_{1}}H , \\ f = \frac{F{T}_{1}^{2}}{{M}_{1}{L}_{1}} , {f}_{d} = \frac{{F}_{d}{T}_{1}^{2}}{{M}_{1}{L}_{1}} , \omega = \varOmega {T}_{1} , \varepsilon = E{L}_{1}^{2} , t = \frac{T}{{T}_{1}} , {x}_{i} = \frac{{X}_{i}}{{L}_{1}} , i = 1, 2, 3. $ |
Substituting the above dimensionless parameters into Eq (1), the dimensionless equation of the system motion can be obtained:
|
$ \left\{\begin{array}{l}{\ddot{x}}_{1} = f\mathit{cos}(wt+\varphi )-{\xi }_{1}{\dot{x}}_{1}-{k}_{1}{x}_{1}-\varepsilon {k}_{1}{x}_{1}^{3}+{k}_{2}{x}_{2}+\varepsilon {k}_{2}{x}_{2}^{3}+{k}_{3}{x}_{3}+\varepsilon {k}_{3}{x}_{3}^{3}\\ +{\xi }_{2}[1+\eta \mathit{sgn}({\dot{x}}_{2}-{\dot{x}}_{1})\left]\right({\dot{x}}_{2}-{\dot{x}}_{1}{)}^{1/3}+{\xi }_{3}[1+\eta \mathit{sgn}({\dot{x}}_{3}-{\dot{x}}_{1})\left]\right({\dot{x}}_{3}-{\dot{x}}_{1}{)}^{1/3}\\ \begin{array}{c}{\ddot{x}}_{2} = \frac{1}{{m}_{2}}[-{k}_{2}{x}_{2}-\varepsilon {k}_{2}{x}_{2}^{3}+{f}_{d}-{\xi }_{2}[1+\eta \mathit{sgn}({\dot{x}}_{2}-{\dot{x}}_{1}\left)\right]({\dot{x}}_{2}-{\dot{x}}_{1}{)}^{1/3}]\\ {\ddot{x}}_{3} = \frac{1}{{m}_{3}}[-{k}_{3}{x}_{3}-\varepsilon {k}_{3}{x}_{3}^{3}+{f}_{d}-{\xi }_{3}[1+\eta \mathit{sgn}({\dot{x}}_{3}-{\dot{x}}_{1}\left)\right]({\dot{x}}_{3}-{\dot{x}}_{1}{)}^{1/3}]\end{array} \end{array}\right. $
|
(5) |
The dimensionless expression Eq (6) of Dankowicz dynamic friction model is [17]:
| $ {f}_{d} = {f}_{r}(\mu {e}^{-y/\sigma }\frac{z}{\delta }+{\alpha }_{1}{{v}_{r}}^{2}sgn({v}_{r}\left)\right(1-\frac{y}{{y}_{\infty }}\left){e}^{-y/\sigma }\right)+{\sigma }_{1}\dot{z} $ | (6) |
In Eq (6), $ {f}_{r} $ represents the normal total load on the contact surface, $ \mu $ represents the friction coefficient, $ \sigma $ represents the standard deviation of the roughness height, $ \delta $ represents the maximum allowable asperity deformation, $ {v}_{r} = {\dot{x}}_{1}-{\dot{x}}_{i} $ represents the relative slip velocity between mass 1 and mass 2 or mass 3, $ {\mathrm{y}}_{\mathrm{\infty }} $ represents the separation distance from the radius of the asperity to the normal direction, $ {\sigma }_{1} $ is the bristle damping, $ {\alpha }_{1} $ is an ordinary parameter, and $ \dot{z} $ is the internal damping term. The evolution of the internal variable $ z $ of the dynamic friction model is dominated by Eq (7):
| $ \dot{z} = {v}_{r}(1-\frac{z}{\delta }sgn({v}_{r}\left)\right) $ | (7) |
In Eq (7), for $ {v}_{r} > $0, z tends to $ \delta $; for $ {v}_{r} < $0, z tends to $ -\delta $; therefore $ z\in \left[-\delta, \delta \right] $.
The motion of the state variable y is dominated by Eq (8):
| $ \ddot{y} = {f}_{r}({\mu e}^{-y/\sigma }-1)+{f}_{r}(\beta {{v}_{r}}^{2}\sqrt{1-\frac{y}{{y}_{\infty }}}{e}^{-y/\sigma }-\gamma \left|{v}_{r}\right|\dot{y}{e}^{-y/\sigma }) $ | (8) |
In Eq (8): $ \beta > 0 $, $ \gamma > 0 $ are free parameters. The first term on the right side of Eq (8) becomes from the asperity shape in the normal direction, and the second term is the normal component of the force generated by the collision of the asperities.
When the displacement difference $ \left|{x}_{1}-{x}_{i}\right| = d $, the system collides. Since the nonlinear system contains collision and friction, the motion state of the system will be converted between chatter, sticking and collision. In order to more fully understand the changes in the motion process of the system, according to the force and velocity changes of the mass 1, the motion process of the system is divided into the following situations.
Case 1: Mass $ {M}_{1} $ moves under the action of harmonic excitation force $ F\mathit{cos}(\varOmega T+\varphi) $, and the system force analysis of mass $ {M}_{1} $ is carried out. Let $ {f}_{1} $ represent the resultant force of mass $ {M}_{1} $, the $ {f}_{1} = fcos(\omega t+\varphi)-{\xi }_{1}{\dot{x}}_{1}-{k}_{1}{x}_{1}-\varepsilon {k}_{1}{x}_{1}^{3}+{k}_{2}{x}_{2}+\varepsilon {k}_{2}{x}_{2}^{3}+{k}_{3}{x}_{3}+\varepsilon {k}_{3}{x}_{3}^{3}$+${\xi }_{2}[1+\eta sgn\left({\dot{x}}_{2}-{\dot{x}}_{1}\right)]{\left({\dot{x}}_{2}{-\dot{x}}_{1}\right)}^{1/3}+{\xi }_{3}[1+\eta sgn\left({\dot{x}}_{3}-{\dot{x}}_{1}\right)]{\left({\dot{x}}_{3}{-\dot{x}}_{1}\right)}^{1/3} $. When the difference between the motion displacement of the mass $ {M}_{1} $ and the mass $ {M}_{2} $ or the mass $ {M}_{3} $ is equal to gap d, the system collides, due to the symmetry of the left wedge and the right wedge, the resultant force on mass $ {M}_{1} $ is the resultant force in the vertical direction as shown in the Figure 3, and the judgment expression is:
| $ \left|{x}_{1}-{x}_{i}\right| = d (i = 2, 3) $ | (9) |
According to the law of conservation of momentum and the relationship between the velocity before and after the collision, it can be obtained that:
|
$ \left\{μm˙x1++˙xi+=μm˙x1−+˙xi−˙xi+−˙x1+=−R(˙xi−−˙x1−) \right. (i = 2, 3) $
|
(10) |
Among them, R is the collision recovery coefficient.
Case 2: If the resultant force of the mass $ {M}_{1} $ is greater than zero or equal to zero and the relative velocity of mass $ {M}_{1} $ and mass $ {M}_{1} $ or $ {M}_{1} $ is not zero, the mass $ {M}_{1} $ performs accelerated slip motion, and then $ {f}_{1}\ge $0, $ \left|{\dot{x}}_{1}-{\dot{x}}_{2}\right|\ne $0.
Case 3: If the resultant force of the mass $ {M}_{1} $ is less than zero and the relative velocity of mass $ {\mathrm{M}}_{1} $ and mass $ {M}_{1} $ or $ {M}_{1} $ is not zero, the mass $ {M}_{1} $ performs decelerated slip motion, and then $ {f}_{1}\le $0, $ \left|{\dot{x}}_{1}-{\dot{x}}_{2}\right|\ne $0.
Case 4: If the resultant force of mass $ {M}_{1} $ is less than zero and the relative velocity between mass $ {M}_{1} $ and mass $ {M}_{2} $ or mass $ {M}_{3} $ is zero, that is, $ {f}_{1} < $0, $ \left|{\dot{x}}_{1}-{\dot{x}}_{2}\right| = $0, mass $ {M}_{1} $ is in a sticking state.
Based on the above four cases, the numerical simulation of the wedge friction damper oblique collision system with gap and friction is carried out to analyze its dynamic behavior. In order to study the periodic motion and bifurcation characteristics of the wedge damper system, the Poincaré section $ \sigma = \left\{\left({x}_{1}, {\dot{x}}_{1}, {x}_{i}, {\dot{x}}_{i}, \theta \right)\in {\mathrm{R}}^{4}\times \mathrm{S}, {x}_{1}-{x}_{i} = d, {\dot{x}}_{1} = {\dot{x}}_{1+}, {\dot{x}}_{i} = {\dot{x}}_{i+}\right\} $ is selected.
The Poincare mapping is established:
| $ {X}^{\left(j+1\right)} = \stackrel{~}{f}\left(v, {X}^{\left(j\right)}\right) $ | (11) |
Of which:
| $ {X}^{\left(j\right)} = {\left[{\tau }^{\left(j\right)}, {{X}_{i}}^{\left(j\right)}, {\dot{X}}_{1+}^{\left(j\right)}{\dot{X}}_{i+}^{\left(j\right)}\right]}^{T} $ |
| $ {X}^{\left(j+1\right)} = {\left[{\tau }^{\left(j+1\right)}, {{X}_{i}}^{\left(j+1\right)}, {\dot{X}}_{1+}^{(j+1)}, {\dot{X}}_{i+}^{(j+1)}\right]}^{T} $ | (12) |
In Eqs (11) and (12), $ i = 2, \mathrm{ }3, X\in {\mathrm{R}}^{4} $, $ v $ is a real parameter. The periodic motion in the following discussion is represented by n-p-q-l, where n represents the number of excitation periods, p represents the number of collisions, q represents the number of chatters, and l represents the number of stickings.
In advance, through a large number of data simulation experiments, we find the different dynamic behavior of the system under different parameters, and select the system parameters of the wedge friction damper combined with the engineering practice: $ {M}_{1} $ = 489 kg, $ {M}_{2} $ = 537.9 kg, $ F $ = 112 N, $ {K}_{1} $ = 2.4 × 104 N/mm, $ {K}_{2} $ = $ {K}_{3} $ = 2.04 × 104 N/mm, $ {C}_{1} $ = 171 N∙s/m, $ {C}_{2} $ = $ {C}_{3} $ = 643 N∙s/m. Take the parameters into Eq (1), and calculate the dimensionless parameters of the system: $ {k}_{1} $ = 0.8, $ {k}_{2} $ = 0.85, $ {k}_{3} $ = 0.85, $ {m}_{2} $ = 1.1, $ {m}_{3} $ = 1.1, $ R $ = 0.85, $ {\xi }_{1} $ = 0.05, $ {\xi }_{2} $ = 0.15, $ {\xi }_{3} $ = 0.15, $ {\mu }_{m} $ = 0.4, $ \mu $ = 0.2, $ \beta $ = 100, $ {y}_{\infty } $ = 2000, $ \delta $ = 0.0001, $ \gamma $ = 3000, $ \sigma $ = 0.01, $ {\alpha }_{1} $ = 2.2, $ f $ = 1, $ {f}_{r} $ = 0.2, $ \theta $ = 53. Set the initial value: $ {x}_{1} $ = 0, $ {x}_{2} $ = 0, $ {\dot{x}}_{1} $ = 0, $ {\dot{x}}_{2} $ = 0. Taking the dimensionless frequency as the bifurcation control parameter, the global and local bifurcation diagrams of the system motion obtained by numerical simulation are shown in Figure 4. From the bifurcation diagram in Figure 4, it can be seen that as the frequency ω increases from 1 to 8, the system motion undergoes inverse periodic doubling bifurcation, periodic doubling bifurcation, Grazing bifurcation, and Hopf bifurcation to each other between single-period, multi-period, quasi-period, and chaos. As shown in Figure 4(a), when ω∈[1, 1.454], the system is in a chaotic state, and when ω∈[1.454, 3.2], the system is in period three motion. As shown in Figure 4(b), when ω∈[3.2, 4.1], the system is in chaotic and multi-period motion. As shown in Figure 4(c), when ω∈[4.3, 4.66], the system is in period six motion and period nine motion successively, and then Hopf bifurcation occurs. When ω∈[4.66, 4.93], the system is in almost periodic motion. As shown in Figure 4(a), when [4.93, 6.33], the system transitions between periodic three motions, periodic six motions and periodic twelve motions. As shown in Figure 4(d), when [7, 8], the system transitions between multi-period motion and chaos. By analyzing the motion process of the system with frequency change, it is found that the system is a stable periodic motion in the parameter interval ω∈[1.75, 3.24], ω∈[4.029, 4.65], ω∈[5, 6.33] and ω∈[6.6, 7.225]. The elastic force, damping force, and external excitation force are coupled to each other due to collision and friction in the motion system of the wedge shock absorber, which leads to chatter. Through the phase diagram, poincare section diagram, and time history diagram of the system motion with the change of bifurcation control parameters, the chatter, collision, and sticking of the wedge damper during the train operation are studied in detail.
When the system frequency is $ \omega = 1.75 $, as shown in Figure 5(a), (b), the corresponding phase diagram and displacement-time history diagram of the system are in the 1-3-2-0 periodic motion. Under the action of an external excitation, the system has three collisions and two chatters. As the external excitation frequency of the system increases, the system undergoes a period-doubling bifurcation and transits to a short period of twelve motions. When the system frequency $ \omega = 3.24 $, as shown in phase diagram 5(c) and displacement-time history diagram 5(d), the system is in a 1-12-9-0 periodic motion. Under the action of an external excitation, the system undergoes twelve collisions and nine chatters. As the system frequency continues to increase, when $ \omega \in [3.35, \mathrm{ }3.55] $ is selected, the sticking phenomenon occurs, and the system motion state is converted between sticking and chatter. As shown in Figure 5(e)–(p), the amplitude of the system chatter first decreases from large to small, reaches the minimum value at $ \omega = 3.395 $, and then decreases from small to large. Furthermore, the duration of the sticking phenomenon first increases from less to the longest at $ \omega = 3.395 $, and then the phenomenon gradually weakens until it disappears. The excitation frequency $ \omega = 3.922 $ is selected, as shown in Figure 5(q), (r). The system is in a 1-8-6-0 periodic motion. Under the action of an external excitation, the system has eight collisions and six chatters. When the excitation frequency $ \omega = 4.029 $, it can be seen from Figure 5(s), (t) that the system has a Grazing bifurcation. Because the displacement difference between mass $ {M}_{1} $ and mass $ {M}_{2} $ and mass $ {M}_{3} $ is exactly d, the direction of the resultant force of mass $ {M}_{1} $ suddenly changes, and under the action of the resultant force, it moves away from mass $ {M}_{2} $ and mass $ {M}_{3} $. At this time, the system is in a 1-6-6-0 periodic motion. Under the action of an external excitation, the system has six collisions and six chatters.
When the excitation frequency $ \omega = 4.187 $ is selected, as shown in Figure 6(a), (b), the corresponding phase diagram and displacement-time history diagram of the system are in the 1-3-3-0 periodic motion. Under the action of an external excitation, the system has three collisions and three chatters. Then, the system undergoes periodic doubling bifurcation, as shown in Figure 6(c)–(e). It can be seen that the system is in a 1-6-4-1 periodic motion. Under the action of an external excitation, the system undergoes six collisions, four chatters and one sticking. When the excitation frequency $ \omega = 4.4258 $, as shown in Figure 6(f)–(h), the system enters the period 1-9-7-2 motion, and the system appears the Grazing collision behavior through the Grazing bifurcation. When the excitation frequency $ \omega \in [4.65, \mathrm{ }4.9] $, as shown in Figure 6(i)–(l), with the increase of $ \omega $, the inverse Hopf bifurcation occurs in the system, and the system enters the periodic 1-3-2-1 motion through phase locking. Under the action of an external excitation, the system undergoes three collisions, two chatters and one sticking.
The excitation frequency $ \omega = 6.12 $ is selected, such as Figure 7(a)–(c), which shows that the system is in period 1-12-12-1 motion. Under the action of an external excitation, the system has twelve collisions, twelve chatters and one sticking. The excitation frequency $ \omega = 6.34 $ is selected, as shown in Figure 7(d), (e), which shows that the system enters the period 1-6-6-0 motion through the grazing bifurcation, and the sticking phenomenon disappears. Select the excitation frequency $ \omega = 7.05 $. As shown in Figure 7(f), (g), the system is in periodic 1-3-3-0 motion. Under the action of an external excitation, the system undergoes three collisions and three chatters. When the excitation frequency $ \omega = 7.276 $, the system enters the period 1-5-5-0 motion. Under the action of an external excitation, the system has five collisions and five chatters. As the system frequency continues to increase, when $ \omega \in [7.34, \mathrm{ }7.578] $ is selected, the system motion state is converted between sticking and chatter. As shown in Figure 7(e)–(p), the amplitude of the system chatter first decreases from large to small, reaches the minimum value at $ \omega = 3.395 $, and then decreases from small to large. At the same time, the duration of the sticking phenomenon increases from less to more, and the duration is the longest when $ \omega = 3.395 $, and then the phenomenon gradually weakens until it disappears. When the excitation frequency $ \omega = 7.877 $ is selected, the inverse Hopf bifurcation occurs in the system, and the chaotic state enters the periodic motion state through phase locking. When the excitation frequency $ \omega = 7.905 $, the system enters a period of 1-4-4-0 motion. Under the action of an external excitation, the system undergoes four collisions and four chatters.
The change of the system motion transfer process and the parameter interval of the system stable motion is analyzed further. Through numerical simulation, it is found that the periodic bubble phenomenon of the system is very sensitive to parameter $ {\mathrm{m}}_{2} $, and the slight change of parameter $ {\mathrm{m}}_{2} $ makes the number of periodic bubbles of the system change or disappear. In order to better understand the influence of parameters on the dynamic behavior of the system, the system parameters $ {m}_{2} $ = 1.15, $ {m}_{2} $ = 1.125, $ {m}_{2} $ = 1.1, $ {m}_{2} $ = 1.05, $ {m}_{2} $ = 1, $ {m}_{2} $ = 0.95, $ {m}_{2} $ = 0.9, $ {m}_{2} $ = 0.85 are taken, and the remaining parameters remain unchanged. As the mass ratio decreases, when $ {m}_{2}\in [1.1, \mathrm{ }1.15] $, as shown in Figure 8(a)–(c), the number of periodic bubble structures in the system increases from 4 to 7 and then to 11. When $ {m}_{2}\in [1, \mathrm{ }1.05] $, as shown in Figure 8(d), (e), the periodic bubble structure disappears, and the system presents six chaotic bubble structures. When $ {m}_{2}\in [0.85, \mathrm{ }0.95] $, as shown in Figure 8(f)–(h), the chaotic bubble structure disappears and the periodic bubble structure appears. The number of periodic bubble structures is 17, 6, and 3, respectively. Moreover, with the decrease of mass ratio, the transition region between flutter and viscosity of the system gradually moves to the left, and the period of period three motion window in the low frequency band is reduced. The Hopf bifurcation phenomenon before the periodic bubble phenomenon or the chaotic bubble phenomenon gradually weakens until it disappears, and the window period of the periodic motion in the corresponding frequency band gradually increases.
In this paper, the dynamic model of a kind of wedge damper with dynamic friction is simplified, and the differential equation of system motion is established. According to the change of mass force, the motion situation of the system and the judgment conditions of each situation are discussed. The dynamic characteristics of the system are analyzed by a large number of data simulation. The simulation results show that:
Under certain parameters, with the change of excitation frequency, the system undergoes periodic doubling bifurcation, inverse periodic doubling bifurcation, Grazing bifurcation and Hopf bifurcation between single period, multi period, quasi period, and chaos. In addition, when the excitation frequency of the system is between [5.2, 6.1], there is a periodic bubble phenomenon in the system motion. When the excitation frequency of the system is between [3.35, 3.55], [4.425, 6.12], and [7.34, 7.758], there are chatter and sticking phenomena in the system motion. When the excitation frequency of the system is between [1.75, 3.24] and [3.92, 4.425], the sticking phenomenon disappears, and only the chatter phenomenon exists.
When other parameters remain unchanged, as the mass ratio decreases, the number of periodic bubble structures of the system motion increases exponentially, that is, the number of cycles of the system motion increases exponentially. Then, the periodic bubble structure disappears and the chaotic bubble structure appears. Finally, the chaotic bubble structure disappears and the system enters periodic motion.
Through the theoretical analysis and numerical simulation of the wedge friction damper oblique collision system with Dankowicz dynamic friction, it provides a certain reference basis for the selection of structural parameters and system motion process control of wedge friction damper for freight train bogies in practical engineering. Furthermore, the method of applying the bifurcation theory of mapping to study the stability, bifurcation, and chaos formation process of periodic motion of multi-degree-of-freedom impact vibration system can also be used to study the periodic motion, stability and bifurcation of other types of impact vibration systems, such as double-mass impact vibration molding machine, pile driver, vibration hammer, wheel-rail collision of high-speed train, etc., and can also be applied to the dynamic analysis of some mechanical systems with clearance and elastic constraints.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The work was supported by the Gansu Provincial Science and Technology Plan Project (21YF5WA060) and the National Natural Science Foundation of China (11302092). The author is grateful for the financial support.
The authors declare that there are no conflicts of interest.
| [1] |
Felipo V (2013) Hepatic encephalopathy: effects of liver failure on brain function. Nat Rev Neurosci 14: 851–858. doi: 10.1038/nrn3587
|
| [2] |
Olde Damink SW, Jalan R, Dejong CH (2009) Interorgan ammonia trafficking in liver disease. Metab Brain Dis 24: 169–181. doi: 10.1007/s11011-008-9122-5
|
| [3] |
Butterworth RF, Norenberg MD, Felipo V, et al. (2009) Experimental models of hepatic encephalopathy: ISHEN guidelines. Liver Int 29: 783–788. doi: 10.1111/j.1478-3231.2009.02034.x
|
| [4] |
Leke R, Oliveira DL, Forgiarini LF, et al. (2013) Impairment of short term memory in rats with hepatic encephalopathy due to bile duct ligation. Metab Brain Dis 28: 187–192. doi: 10.1007/s11011-012-9347-1
|
| [5] |
Weissenborn K, Giewekemeyer K, Heidenreich S, et al. (2005) Attention, memory, and cognitive function in hepatic encephalopathy. Metab Brain Dis 20: 359–367. doi: 10.1007/s11011-005-7919-z
|
| [6] |
Ferenci P, Lockwood A, Mullen K, et al. (2002) Hepatic encephalopathy--definition, nomenclature, diagnosis, and quantification: final report of the working party at the 11th World Congresses of Gastroenterology, Vienna, 1998. Hepatology 35: 716–721. doi: 10.1053/jhep.2002.31250
|
| [7] | Patel A, Wade JB, Thacker LR, et al. (2014) 991 Brain Reserve Modulates Health-Related Quality of Life in Patients With Cirrhosis Independent of Covert Hepatic Encephalopathy and MELD Score. Gastroenterology 146: S–932. |
| [8] |
Rose CF (2012) Ammonia-lowering strategies for the treatment of hepatic encephalopathy. Clin Pharmacol Ther 92: 321–331. doi: 10.1038/clpt.2012.112
|
| [9] |
Savlan I, Liakina V, Valantinas J (2014) Concise review of current concepts on nomenclature and pathophysiology of hepatic encephalopathy. Medicina (Kaunas) 50: 75–81. doi: 10.1016/j.medici.2014.06.008
|
| [10] |
Rama Rao KV, Jayakumar AR, Norenberg MD (2012) Glutamine in the pathogenesis of acute hepatic encephalopathy. Neurochem Int 61: 575–580. doi: 10.1016/j.neuint.2012.01.012
|
| [11] |
Rothman DL, De Feyter HM, Maciejewski PK, et al. (2012) Is there in vivo evidence for amino acid shuttles carrying ammonia from neurons to astrocytes? Neurochem Res 37: 2597–2612. doi: 10.1007/s11064-012-0898-7
|
| [12] |
Jover R, Rodrigo R, Felipo V, et al. (2006) Brain edema and inflammatory activation in bile duct ligated rats with diet-induced hyperammonemia: A model of hepatic encephalopathy in cirrhosis. Hepatology 43: 1257–1266. doi: 10.1002/hep.21180
|
| [13] |
Sheen JM, Huang LT, Hsieh CS, et al. (2010) Bile duct ligation in developing rats: temporal progression of liver, kidney, and brain damage. J Pediatr Surg 45: 1650–1658. doi: 10.1016/j.jpedsurg.2009.12.019
|
| [14] | Aldridge DR, Tranah EJ, Shawcross DL (2015) Pathogenesis of hepatic encephalopathy: role of ammonia and systemic inflammation. J Clin Exp Hepatol 5: S7–S20. |
| [15] |
Rodrigo R, Cauli O, Gomez-Pinedo U, et al. (2010) Hyperammonemia induces neuroinflammation that contributes to cognitive impairment in rats with hepatic encephalopathy. Gastroenterology 139: 675–684. doi: 10.1053/j.gastro.2010.03.040
|
| [16] | Jayakumar AR, Rama Rao KV, Norenberg MD (2015) Neuroinflammation in hepatic encephalopathy: mechanistic aspects. J Clin Exp Hepatol 5: S21–28. |
| [17] |
Butterworth RF (2012) Reprint of: Neuroinflammation in acute liver failure: Mechanisms and novel therapeutic targets. Neurochem Int 60: 715–722. doi: 10.1016/j.neuint.2012.03.014
|
| [18] |
Aghaei I, Nazeri M, Shabani M, et al. (2015) Erythropoietin ameliorates the motor and cognitive function impairments in a rat model of hepatic cirrhosis. Metab Brain Dis 30: 197–204. doi: 10.1007/s11011-014-9600-x
|
| [19] |
Aghaei I, Shabani M, Doustar N, et al. (2014) Peroxisome proliferator-activated receptor-gamma activation attenuates motor and cognition impairments induced by bile duct ligation in a rat model of hepatic cirrhosis. Pharmacol Biochem Behav 120: 133–139. doi: 10.1016/j.pbb.2014.03.002
|
| [20] |
Nairz M, Schroll A, Moschen AR, et al. (2011) Erythropoietin contrastingly affects bacterial infection and experimental colitis by inhibiting nuclear factor-kappaB-inducible immune pathways. Immunity 34: 61–74. doi: 10.1016/j.immuni.2011.01.002
|
| [21] | Bond WS, Rex TS (2014) Evidence that erythropoietin modulates neuroinflammation through differential action on neurons, astrocytes, and microglia. Front Immunol 5: 523. |
| [22] |
Aghaei I, Hajali V, Dehpour A, et al. (2016) Alterations in the intrinsic electrophysiological properties of Purkinje neurons in a rat model of hepatic encephalopathy: Relative preventing effect of PPARgamma agonist. Brain Res Bull 121: 16–25. doi: 10.1016/j.brainresbull.2015.12.002
|
| [23] |
Shabani M, Ebrahimpoor F, Firouzjaei MA, et al. (2019) Modulation of sphingosine-1-phosphate receptor by FTY720 contributes in improvement of hepatic encephalopathy induced by bile duct ligation. Brain Res Bull 146: 253–269. doi: 10.1016/j.brainresbull.2019.01.012
|
| [24] |
Tahamtan M, Aghaei I, Pooladvand V, et al. (2017) Characterization of the CA1 pyramidal neurons in rat model of hepatic cirrhosis: insights into their electrophysiological properties. Metab Brain Dis 32: 881–889. doi: 10.1007/s11011-017-9966-7
|
| [25] |
Onoda A, Takeda K, Umezawa M (2017) Dose-dependent induction of astrocyte activation and reactive astrogliosis in mouse brain following maternal exposure to carbon black nanoparticle. Part Fibre Toxicol 14: 4. doi: 10.1186/s12989-017-0184-6
|
| [26] |
Matsui H, Ohgomori T, Natori T, et al. (2013) Keratan sulfate expression in microglia is diminished in the spinal cord in experimental autoimmune neuritis. Cell Death Dis 4: e946. doi: 10.1038/cddis.2013.479
|
| [27] |
Javadi-Paydar M, Ghiassy B, Ebadian S, et al. (2013) Nitric oxide mediates the beneficial effect of chronic naltrexone on cholestasis-induced memory impairment in male rats. Behav Pharmacol 24: 195–206. doi: 10.1097/FBP.0b013e3283618a8c
|
| [28] |
Nasehi M, Piri M, Abbolhasani K, et al. (2013) Involvement of opioidergic and nitrergic systems in memory acquisition and exploratory behaviors in cholestatic mice. Behav Pharmacol 24: 180–194. doi: 10.1097/FBP.0b013e3283618aab
|
| [29] |
Butterworth RF (2011) Hepatic encephalopathy: a central neuroinflammatory disorder? Hepatology 53: 1372–1376. doi: 10.1002/hep.24228
|
| [30] |
Huang LT, Tiao MM, Tain YL, et al. (2009) Melatonin ameliorates bile duct ligation-induced systemic oxidative stress and spatial memory deficits in developing rats. Pediatr Res 65: 176–180. doi: 10.1203/PDR.0b013e31818d5bc7
|
| [31] |
Brines M, Cerami A (2005) Emerging biological roles for erythropoietin in the nervous system. Nat Rev Neurosci 6: 484–494. doi: 10.1038/nrn1687
|
| [32] |
Levitt DG, Levitt MD (2016) Human serum albumin homeostasis: a new look at the roles of synthesis, catabolism, renal and gastrointestinal excretion, and the clinical value of serum albumin measurements. Int J Gen Med 9: 229–255. doi: 10.2147/IJGM.S102819
|
| [33] |
Caillaud C, Mechta M, Ainge H, et al. (2015) Chronic erythropoietin treatment improves diet-induced glucose intolerance in rats. J Endocrinol 225: 77–88. doi: 10.1530/JOE-15-0010
|
| [34] |
Dhanda S, Sandhir R (2015) Role of dopaminergic and serotonergic neurotransmitters in behavioral alterations observed in rodent model of hepatic encephalopathy. Behav Brain Res 286: 222–235. doi: 10.1016/j.bbr.2015.01.042
|
| [35] |
Su YY, Yang GF, Lu GM, et al. (2015) PET and MR imaging of neuroinflammation in hepatic encephalopathy. Metab Brain Dis 30: 31–45. doi: 10.1007/s11011-014-9633-1
|
| [36] |
Chen JR, Wang BN, Tseng GF, et al. (2014) Morphological changes of cortical pyramidal neurons in hepatic encephalopathy. BMC Neurosci 15: 15. doi: 10.1186/1471-2202-15-15
|
| [37] |
Ponce LL, Navarro JC, Ahmed O, et al. (2013) Erythropoietin neuroprotection with traumatic brain injury. Pathophysiology 20: 31–38. doi: 10.1016/j.pathophys.2012.02.005
|
| [38] |
Wenker SD, Chamorro ME, Vittori DC, et al. (2013) Protective action of erythropoietin on neuronal damage induced by activated microglia. FEBS J 280: 1630–1642. doi: 10.1111/febs.12172
|
| [39] |
McPherson RJ, Juul SE (2008) Recent trends in erythropoietin-mediated neuroprotection. Int J Dev Neurosci 26: 103–111. doi: 10.1016/j.ijdevneu.2007.08.012
|
| [40] |
Mofidi A, Bader A, Pavlica S (2011) The use of erythropoietin and its derivatives to treat spinal cord injury. Mini Rev Med Chem 11: 763–770. doi: 10.2174/138955711796355267
|
| [41] | Hu X, Liou AK, Leak RK, et al. (2014) Neurobiology of microglial action in CNS injuries: receptor-mediated signaling mechanisms and functional roles. Prog Neurobiol 119–120: 60–84. |
| [42] |
Wright G, Swain M, Annane D, et al. (2016) Neuroinflammation in liver disease: sessional talks from ISHEN. Metab Brain Dis 31: 1339–1354. doi: 10.1007/s11011-016-9918-7
|
Figures(3) / Tables(1)
Moazameh Golshani, Mohsen Basiri, Mohammad Shabani, Iraj Aghaei, Majid Asadi-Shekaari. Effects of erythropoietin on bile duct ligation-induced neuro-inflammation in male rats[J]. AIMS Neuroscience, 2019, 6(2): 43-53. doi: 10.3934/Neuroscience.2019.2.43
DownLoad:
