Research article

Effect of Caspase Inhibitor Ac-DEVD-CHO on Apoptosis of Vascular Smooth Muscle Cells Induced by Artesunate

  • Received: 17 February 2014 Accepted: 12 May 2014 Published: 20 May 2014
  • Numerous studies have shown that the proliferation and apoptosis of vascular smooth muscle cells play a key role in restenosis. Artesunate is a triterpenoid with a peroxide structure and its antimalarial, antitumor, and antiangiogenetic activities can inhibit the proliferation and apoptosis of multifarious cells. Apoptosis is caused by the activation of a series of intracellular proteolytic enzymes, among which caspase-dependent apoptosis was the earliest to be recognized. The purpose of this article is to study the effects of caspase-3 inhibitor Ac-DEVD-CHO on proliferation and apoptosis of vascular smooth muscle cells induced by Artesunate and to explore the mechanism of Artesunate-induced apoptosis of vascular smooth muscle cells. By using the method based on methyl thiazolyl tetrazolium to observe the effects of Artesunate on the growth and proliferation of vascular smooth muscle cells; observing the change in cell shape before and after Artesunate administration by transmission electron microscopy; detecting the changes in cell cycle and apoptosis rates before and after drug administration by flow cytometry; detecting the activity of caspase-3 in the caspase apoptosis pathway by the Western Blot method, we found that Artesunate inhibits the growth and proliferation of vascular smooth muscle cells in a dose- and time-dependent manner within the concentration range of 7.5–120 μg/mL, and the inhibition rate of Artesunate can be as high as 89.49 % at a concentration of 120 μg/mL after acting for 72 hours; vascular smooth muscle cells show a typical apoptosis peak due to the effects of higher concentration of Artesunate. Compared with the control group, the higher-concentration group shows major variability, Ac-DEVD-CHO, however, can significantly decrease this induction; it has been detected by Western Blot that Artesunate can induce caspase-3 activity dramatically in vascular smooth muscle cells, but this activation may be remarkably inhibited by Ac-DEVD-CHO.

    Citation: Jingwen Zhang, Lu Wang, Huolin Chen, Tieying Yin, Yanqun Teng, Kang Zhang, Donghong Yu, Guixue Wang. Effect of Caspase Inhibitor Ac-DEVD-CHO on Apoptosis of Vascular Smooth Muscle Cells Induced by Artesunate[J]. AIMS Bioengineering, 2014, 1(1): 13-24. doi: 10.3934/bioeng.2014.1.13

    Related Papers:

    [1] Attaullah, Shahzaib Ashraf, Noor Rehman, Asghar Khan, Muhammad Naeem, Choonkil Park . Improved VIKOR methodology based on $ q $-rung orthopair hesitant fuzzy rough aggregation information: application in multi expert decision making. AIMS Mathematics, 2022, 7(5): 9524-9548. doi: 10.3934/math.2022530
    [2] Attaullah, Shahzaib Ashraf, Noor Rehman, Asghar Khan, Choonkil Park . A decision making algorithm for wind power plant based on q-rung orthopair hesitant fuzzy rough aggregation information and TOPSIS. AIMS Mathematics, 2022, 7(4): 5241-5274. doi: 10.3934/math.2022292
    [3] Jia-Bao Liu, Rashad Ismail, Muhammad Kamran, Esmail Hassan Abdullatif Al-Sabri, Shahzaib Ashraf, Ismail Naci Cangul . An optimization strategy with SV-neutrosophic quaternion information and probabilistic hesitant fuzzy rough Einstein aggregation operator. AIMS Mathematics, 2023, 8(9): 20612-20653. doi: 10.3934/math.20231051
    [4] Muhammad Saqlain, Xiao Long Xin, Rana Muhammad Zulqarnain, Imran Siddique, Sameh Askar, Ahmad M. Alshamrani . Energy supplier selection using Einstein aggregation operators in an interval-valued q-rung orthopair fuzzy hypersoft structure. AIMS Mathematics, 2024, 9(11): 31317-31365. doi: 10.3934/math.20241510
    [5] Adel Fahad Alrasheedi, Jungeun Kim, Rukhsana Kausar . q-Rung orthopair fuzzy information aggregation and their application towards material selection. AIMS Mathematics, 2023, 8(8): 18780-18808. doi: 10.3934/math.2023956
    [6] Sumbal Ali, Asad Ali, Ahmad Bin Azim, Ahmad ALoqaily, Nabil Mlaiki . Averaging aggregation operators under the environment of q-rung orthopair picture fuzzy soft sets and their applications in MADM problems. AIMS Mathematics, 2023, 8(4): 9027-9053. doi: 10.3934/math.2023452
    [7] Yasir Yasin, Muhammad Riaz, Kholood Alsager . Synergy of machine learning and the Einstein Choquet integral with LOPCOW and fuzzy measures for sustainable solid waste management. AIMS Mathematics, 2025, 10(1): 460-498. doi: 10.3934/math.2025022
    [8] Sumbal Ali, Asad Ali, Ahmad Bin Azim, Abdul Samad Khan, Fuad A. Awwad, Emad A. A. Ismail . TOPSIS method based on q-rung orthopair picture fuzzy soft environment and its application in the context of green supply chain management. AIMS Mathematics, 2024, 9(6): 15149-15171. doi: 10.3934/math.2024735
    [9] Ahmad Bin Azim, Ahmad ALoqaily, Asad Ali, Sumbal Ali, Nabil Mlaiki, Fawad Hussain . q-Spherical fuzzy rough sets and their usage in multi-attribute decision-making problems. AIMS Mathematics, 2023, 8(4): 8210-8248. doi: 10.3934/math.2023415
    [10] Ghous Ali, Kholood Alsager, Asad Ali . Novel linguistic $ q $-rung orthopair fuzzy Aczel-Alsina aggregation operators for group decision-making with applications. AIMS Mathematics, 2024, 9(11): 32328-32365. doi: 10.3934/math.20241551
  • Numerous studies have shown that the proliferation and apoptosis of vascular smooth muscle cells play a key role in restenosis. Artesunate is a triterpenoid with a peroxide structure and its antimalarial, antitumor, and antiangiogenetic activities can inhibit the proliferation and apoptosis of multifarious cells. Apoptosis is caused by the activation of a series of intracellular proteolytic enzymes, among which caspase-dependent apoptosis was the earliest to be recognized. The purpose of this article is to study the effects of caspase-3 inhibitor Ac-DEVD-CHO on proliferation and apoptosis of vascular smooth muscle cells induced by Artesunate and to explore the mechanism of Artesunate-induced apoptosis of vascular smooth muscle cells. By using the method based on methyl thiazolyl tetrazolium to observe the effects of Artesunate on the growth and proliferation of vascular smooth muscle cells; observing the change in cell shape before and after Artesunate administration by transmission electron microscopy; detecting the changes in cell cycle and apoptosis rates before and after drug administration by flow cytometry; detecting the activity of caspase-3 in the caspase apoptosis pathway by the Western Blot method, we found that Artesunate inhibits the growth and proliferation of vascular smooth muscle cells in a dose- and time-dependent manner within the concentration range of 7.5–120 μg/mL, and the inhibition rate of Artesunate can be as high as 89.49 % at a concentration of 120 μg/mL after acting for 72 hours; vascular smooth muscle cells show a typical apoptosis peak due to the effects of higher concentration of Artesunate. Compared with the control group, the higher-concentration group shows major variability, Ac-DEVD-CHO, however, can significantly decrease this induction; it has been detected by Western Blot that Artesunate can induce caspase-3 activity dramatically in vascular smooth muscle cells, but this activation may be remarkably inhibited by Ac-DEVD-CHO.


    Generally, electromagnetic transducers are employed to gauge thickness, calculate an object's rotational speed, and detect defects in materials from steel to various other alloys. Managing the nonlinear interaction between vibrating modes is crucial for advancing nano-mechanical or micro-electromechanical devices. Interconnected oscillators serve as primary models aimed at the performance over various technical, chemical, biological, and physical schemes. In nonlinear electromechanical oscillator organizations, the mechanical part functions as a sensor and is magnetically linked to an electrical part, representing a signal of the observed vibration. The space in the permanent magnet facilitates interaction between the mechanical and electrical elements. Chaos management, vacillations, and the stability of a nonlinear electromechanical structure are examined in [1]. Also, an electromechanical gyrostat system that underwent external perturbation for its chaotic behavior, synchronization, and chaotic characteristics (adaptive control, delayed feedback control) observed in [2]. Yamapi and Bowong [3] utilized a sliding mode organizer to manage the electrostatic transducer classification while exploring the dynamic and chaotic behaviors of a self-sustaining electromechanical system both with and without discontinuity. Siewe et al. [4] employed an electromechanical oscillator method to capture the vertical motion of the earth during an earthquake. The application of slight amplitude damping for managing chaos in the system was also investigated. They discovered that the damping coefficient estimation affects both the chaotic and periodic orbits. The interactions and synchronization of two systems studied in [5,6]: a magnetically linked electrical Rayleigh-Duffing oscillator with linear mechanical oscillators, along with an interconnected self-sustaining electromechanical system featuring various functions. They employed the harmonic balance and averaging methods to identify the amplitudes of the oscillatory states. The dynamics, global bifurcations, and chaotic behavior of a self-sustaining nonlinear electromechanical system exhibiting nonlinear dynamics investigated in [7,8]. Moreover, the impact of elevated nonlinearity values on the behavior and synchronization of interconnected electromechanical systems is examined in [9]. Furthermore, the chaotic dynamics and nonlinear oscillations in an electromechanical seismograph system featuring stiffness that varies over time examined in [10]. Also, different types of active controllers are experimented to minimize system oscillations and determined that negative velocity feedback is the most effective active control for the behavior of the system. The influence of noise factors, coupling coefficients, and restraining quantities on the response of an electromechanical seismograph examined in [11].

    A time-varying stiffness nonlinear electromechanical seismograph system's behavior, stability, approximate solutions, and dynamic feedback control were all investigated by Amer [12]. He also contrasted the perturbation solution with the numerical solution. Amer et al. [13] examined the performance of a twin-tail aircraft structure with both cubic and quadratic nonlinearities via an active control procedure. Sayed et al. [14,15,16,17] examined the non-linear dynamic properties of the rectangular plate within combined excitations. Also, they examined three cases of internal and primary resonances (1:2, 1:1, and 1:1:3) and matched the analytical and numerical solutions of the system. The efficiency of various control procedures in lowering the notable vibrations of a beam was examined by Hamed and Amer [18]. The oscillations and stability of the MEMS gyroscope scheme with distinct parametric forces were investigated in [19]. The frequency response equations for the concurrent resonance situation have been generated using the averaging technique. Dynamical systems motivated by parametric and external influences thoroughly examined in the works of [20,21]. The air gap of a permanent magnet, which only generates a uniform radial magnetic field, is assumed to reflect the interaction between mechanical and electrical components in most electromechanical system studies. To the best of our knowledge, the assumption of deterministic electromechanical system dynamics has been the sole method used to investigate the effects of this particular magnetic connection. However, when random disturbances are taken into account, these systems can exhibit remarkable characteristics. Similar mechanical systems have previously been shown in the literature to display a variety of behaviours and intricate chaotic dynamics when subjected to noise [22,23]. One kind of instrument that uses electromagnetic transduction to gather vibration energy is an electromechanical seismograph.

    Researchers have recently focused on the behavior of electromechanical schemes and their impact on the energy harvesting challenge. The effectiveness of vibration energy harvesting devices influenced by random ambient excitations was analyzed by Martens et al. [24] through the solution of relevant Fokker-Planck equations. Borowiec et al. [25] demonstrated that the noise element of the force influences the system's stability in their investigation of how random excitation affects the performance of an energy harvester. A stochastic averaging method in [26] was introduced to evaluate the mean square electric voltage of a nonlinear energy harvesting system. Li et al. [27] explored a piezoelectric energy harvester featuring tri-stable potential wells persuaded by external magnetic fields. They claimed that the system can be appropriately designed to enhance the frequency bandwidth for a specific deterministic or stochastic input and achieve a high harvesting efficiency at coherence resonance. The multiple scales method (MMS) [28,29,30,31,32,33] and averaging techniques [34,35] are two perturbation methods commonly employed to analyze the efficiency of parametrically excited models. These methods have proven to remain successful in forecasting the behavior of such systems, especially within the frequency range close to the significant parametric resonance [36]. Traditional MMS are shown to effectively estimate responses in very basic scenarios, such as confined frequency ranges around the central parametric resonance, moderate excitations, and minimal system features [37,38,39,40,41,42,43]. Moreover, it has been shown that several control mechanism approaches can reduce the detrimental vibrations produced by different nonlinear systems in [44,45,46,47,48,49,50,51,52].

    This study investigates and manages the nonlinear dynamics and vibration reduction of a model of a nonlinear electromagnetic transducer subjected to parametric and harmonic excitations. The response and stability of the solutions during the most unfavorable resonance cases have analyzed using the perturbation technique [53]. The dynamic response of the sandwiched functionally-graded piezoelectric semiconductor (FGPS) plate with the consideration of the initial electron density is investigated, and the natural frequencies and multi-field coupling are obtained in [54]. Based on the nonlocal piezoelectric semiconductor theory, Fang et al. [55] investigated the transient response of a piezoelectric semiconductor (PS) fiber, and analyzed the bending vibration, electric potential, and concentration of electrons along the nano-fiber with different nonlocal effects. Also, Liu et al. [56] proposed an active disturbance rejection control (ADRC) scheme for the electromagnetic docking of spacecraft in the presence of time-varying delay, fault signals, external disturbances, and elliptical eccentricity. Formerly, Lyu et al. [57] investigated an integrated predictor observer feedback control strategy for the vibration suppression of large-scale spacecraft affected by unbounded input time-delay effects.

    The main contributions of this article are summarized as follows:

    1) The study introduces a novel control mechanism Nonlinear Proportional-Derivative Velocity Feedback (NPDVF) Controller, which integrates nonlinear first- and second-order filters to improve vibration control in systems involving both mechanical and electrical components. This controller is designed to stabilize and suppress vibrations caused by complex nonlinearities and mixed forces within simultaneous resonance case ($ {\varOmega }_{2}\cong {\omega }_{1}, {\varOmega }_{4}\cong {\omega }_{2} $).

    2) The research utilizes the perturbation method to derive approximate solutions for the system equations up to the second order, allowing for a better understanding of the system's behavior under nonlinear conditions. This technique aids in analyzing the stability and response characteristics of the system.

    3) A comprehensive mathematical model is presented, which describes the coupled nonlinear ordinary differential equations governing the interaction between the mechanical and electrical components of the system. The model includes cubic and quadratic nonlinearities, which represent the complex behaviors observed in real-world systems.

    4) The proposed NPDVF controller is shown to significantly reduce harmful vibrations, mitigate unstable motion, and stabilize nonlinear bifurcations in the system. Numerical simulations demonstrate its ability to handle a variety of excitation frequencies and forcing magnitudes, eliminating self-excited vibrations and improving structural stability.

    5) Extensive numerical simulations conducted using MAPLE and MATLAB validate the performance of the controller. The simulations confirm that the NPDVF controller effectively stabilizes the system under mixed force conditions, and the results are consistent with the perturbation analysis, providing strong support for the controller's effectiveness.

    6) The examination of stability and the impact of various framework coefficients were assessed together theoretically and numerically.

    7) The various nonlinear controllers that influence the system are compared through numerical methods. The primary finding from the numerical result indicates that the new controller NPDVF, is most effective in diminishing and eradicating the model's excessive oscillations.

    A seismograph records ground movement during an earthquake. It falls into the category of electromechanical systems known as electromagnetic transducers. The system being examined is the most basic seismograph as it solely captures the vertical component of ground motion and only accommodates low-frequency movements. It, however, exhibits intricate dynamics. This apparatus is illustrated in Figure 1. In the current study, the electrical component of the seismograph includes a linear inductor $ L $, a linear capacitor $ C $, and a linear resistor $ R $, with their voltage and charge conforming to the following equations.

    $ {U}_{C} = \frac{q}{{C}_{0}} \text{ and } {U}_{R} = R\dot{q} , $ (2.1)

    where $ q $ stands the instantaneous electrical charge, $ \dot{q} $ remains its time derivative ($ \dot{q} = \frac{dq}{d\tau } = i $, where $ i $ is the current).

    Figure 1.  Sketch of an electromechanical seismograph model with associated electric circuit [11].

    The mechanical component consists of a suspended mass; its movement is influenced by the intrinsic forces of the mass-spring system and the natural forces affecting the arrangement. Let $ F\left(\tau \right) $ be the force varying with time that influences the frame because of the ground movement. We assume that damping forces (friction, air resistance, etc.) exist, and the corresponding spring exhibits nonlinear with linear stiffness characterized by $ {k}_{0} $; parameters $ {k}_{1} $ and $ {k}_{2} $ describing the nonlinearity of the stiffness based on its type. The mechanical and electrical components engage via the air-gap of a permanent magnet that generates a radial magnetic field $ \overrightarrow{B} $. The connection between the magnetic field $ B $ and the position of the coil $ y $ is considered as:

    $ B = {B}_{0}\left(1-\left(\frac{y+{y}_{0}}{{y}_{max}}({)}^{2}\right)\left(\right)\right) , $ (2.2)

    where $ {B}_{0} $ is the greatest intensity that the field $ B $ attains occurs, $ {y}_{0} $ is the armature's initial position, $ y $ remains its oscillation amplitude, and $ {y}_{max} $ represents the maximum amplitude. Subsequently, the movement of the mechanical component must consider the connection between the Laplace force and the current indicated by:

    $ {F}_{C} = {\alpha }_{0}\left(1-\left(\frac{y+{y}_{0}}{{y}_{max}}({)}^{2}\right)(\dot{)}\right) . $ (2.3)

    While in the electrical part, we must include the Lenz electromotive voltage $ {E}_{bemf} $ as :

    $ {E}_{bemf} = {k}_{0}\left(1-\left(\frac{y+{y}_{0}}{{y}_{max}}({)}^{2}\right)(\dot{)}\right). $ (2.4)

    The entire mathematical model that signifies the physical model in Figure 1 can be derived using Newton's second law and Kirchhoff's laws. It is regulated by the subsequent nonlinear differential equations:

    $ m\ddot{y}+{\mu }_{0}\dot{y}+{k}_{0}y+{k}_{1}{y}^{2}+{k}_{2}{y}^{3} = {F}_{C}+F\left(\tau \right), $ (2.5a)
    $ L\ddot{q}+R\dot{q}+\frac{1}{C}q+{E}_{bemf} = 0, $ (2.5b)

    where $ y $ is the relative displacement of the mass $ m $ with inertial forces $ m\ddot{y} $ and damping forces $ {\mu }_{0}\dot{y}, $ and $ {k}_{0}, {k}_{1}, {k}_{2} $ remain linear and nonlinear stiffness of the electromechanical oscillator system. The coupling between the above equations is assured by nonlinear radial magnetic field. The external ground motion is expected to be stochastic or periodic $ \left(F\left(\tau \right) = {F}_{0}+{F}_{1}\mathit{cos}({\varOmega }_{1}\tau)+{F}_{2}y\mathit{cos}({\varOmega }_{2}\tau)\right) $ where $ {F}_{0} $ is the critical amplitude, $ {F}_{1}\mathit{cos}({\varOmega }_{1}\tau) $ is the external force, with amplitude $ {F}_{1} $ and $ {\varOmega }_{1} $ is the excitation frequency, $ {F}_{2}y\mathit{cos}({\varOmega }_{2}\tau) $ is the parametric force, with amplitude $ {F}_{2} $ and $ {\varOmega }_{2} $ is the excitation frequency. We put Eq (2.5) into dimensionless form by setting: $ x = y/l $, $ z = q/{Q}_{0} $ where $ {Q}_{0} $ is the reference charge and $ l $ is the reference length. Let's set $ {\omega }_{e} = \sqrt{\frac{1}{LC}} $, $ {\omega }_{m} = \sqrt{\frac{{k}_{0}}{m}} $, by the time transformation $ t = {\omega }_{e}\tau $.

    Dimensionless variables introduced by scaling the system using characteristic quantities. The main goal is to eliminate the physical units and simplify the equations. Then we convert (2.5) into the following

    $ ω1=ωmωe,ω2=ωeωm,μm=μ0mωe,μe=RLωe,λ1=k1lmω2e,λ2=k2l2mω2e,f0=F0,f1=F1mlω2e,f2=F2mlω2e,γ0=α0Q0mlωe(y20y2max1),γ1=2α0y0Q0mωey2max,γ2=α0lQ0mωey2max,β0=k0lLQ0ωe(1y20y2max),β1=2k0y0l2LQ0ωey2max,β2=k0l3LQ0ωey2max.
    $

    We attain the resulting dimensionless equations for the present electromechanical seismograph as:

    $¨x+μm˙x+ω21x+λ1x2+λ2x3+(γ0+γ1x+γ2x2)˙z=f0+f1cos(Ω1t)+f2xcos(Ω2t),
    $
    (2.6a)
    $ \ddot{z}+{\mu }_{e}\dot{z}+{\omega }_{2}^{2}z+({\beta }_{0}+{\beta }_{1}x+{\beta }_{2}{x}^{2})\dot{x} = 0 . $ (2.6b)

    Assume

    $ {\mu }_{m}\to \varepsilon {\mu }_{m},{\lambda }_{1}\to \varepsilon {\lambda }_{1},{\lambda }_{2}\to \varepsilon {\lambda }_{2},{\gamma }_{0}\to \varepsilon {\gamma }_{0},{\gamma }_{1}\to \varepsilon {\gamma }_{1},{\gamma }_{2}\to \varepsilon {\gamma }_{2},{f}_{0}\to \varepsilon {f}_{0},{f}_{1}\to $
    $ \varepsilon {f}_{1},{f}_{2}\to \varepsilon {f}_{2},{\mu }_{e}\to \varepsilon {\mu }_{e},{\beta }_{0}\to \varepsilon {\beta }_{0},{\beta }_{1}\to \varepsilon {\beta }_{1},{\beta }_{2}\to \varepsilon {\beta }_{2},{f}_{3}\to \varepsilon {f}_{3},{f}_{4}\to \varepsilon {f}_{4} . $

    Where, $ \varepsilon $ is often used as a small perturbation parameter in the analysis of nonlinear systems. The idea is to treat $ \varepsilon $ as a small parameter and perform an asymptotic expansion or perturbation analysis. It ensures stability, smooths out nonlinearities, helps with perturbation analysis, and is crucial for ensuring that numerical methods work effectively.

    So, Eq (2.6) will be written as

    $¨x+εμm˙x+ω21x+ελ1x2+ελ2x3+ε(γ0+γ1x+γ2x2)˙z=ε(f0+f1cos(Ω1t)+f2xcos(Ω2t)),
    $
    (2.7a)
    $ \ddot{z}+\varepsilon {\mu }_{e}\dot{z}+{\omega }_{2}^{2}z+\varepsilon ({\beta }_{0}+{\beta }_{1}x+{\beta }_{2}{x}^{2})\dot{x} = 0 . $ (2.7b)

    The first oscillator $ x $ (mechanical part) is a forced Duffing oscillator related with nonlinear coupling term, and the second one $ z $ (electrical part) is a linear damped oscillator with nonlinear coupling term. $ \dot{x} $, $ \dot{z} $, $ \ddot{x} $ and $ \ddot{z} $ are the first and second derivative with respect to time $ t $, $ {\mu }_{m} $ and $ {\mu }_{e} $ are linear damping coefficients, $ {\lambda }_{1}, {\lambda }_{2} $ are non-linear parameters, $ \varepsilon $ is a small perturbation parameter where $ 0 < \varepsilon < < 1 $, $ {f}_{0}, {f}_{1}, {f}_{2}, {f}_{3}, {f}_{4} $, are the excitation forces amplitudes ($ {f}_{0}, {f}_{1}, {f}_{2} $ are mixed mechanical excitations and $ {f}_{3}, {f}_{4} $ are mixed electrical excitations), $ {\omega }_{1}, {\omega }_{2} $ stand the natural frequencies and $ {\varOmega }_{1}, {\varOmega }_{2}, {\varOmega }_{3}, {\varOmega }_{4} $ remain the excitation frequencies, $ {\gamma }_{j} $ and $ {\beta }_{j} $ $ (j = \mathrm{0, 1}, 2) $ are the coupling terms.

    The modified and investigated system [11] after adding mixed excitations with new controller technique is in the following form [11]:

    $¨x+εμm˙x+ω21x+ελ1x2+ελ2x3+ε(γ0+γ1x+γ2x2)˙z=ε(f0+f1cos(Ω1t)+f2xcos(Ω2t))+F1c(t),
    $
    (2.8a)
    $ \ddot{z}+\varepsilon {\mu }_{e}\dot{z}+{\omega }_{2}^{2}z+\varepsilon ({\beta }_{0}+{\beta }_{1}x+{\beta }_{2}{x}^{2})\dot{x} = \varepsilon \left({f}_{3}\mathit{cos}({\varOmega }_{3}t\right)+{f}_{4}z\mathit{cos}({\varOmega }_{4}t\left)\right)+{F}_{2c}\left(t\right). $ (2.8b)

    $ {F}_{1c}\left(t\right), {F}_{2c}\left(t\right) $ are the control inputs to reduce the vibrations that happen at principal simultaneous resonance case ($ {\varOmega }_{2}\cong {\omega }_{1}, {\varOmega }_{4}\cong {\omega }_{2} $) and this optimal controller is called Nonlinear Proportional-Derivative (NPD) within Negative Cubic Velocity Feedback Controller (NCVFC) to be a resonant innovative controller (NPDVF) as in the form:

    $ {F}_{1c}\left(t\right) = -\varepsilon ({p}_{1}x+\dot{{d}_{1}x}+{\alpha }_{1}{x}^{3}+{\alpha }_{2}{x}^{2}\dot{x}+{\alpha }_{3}x{\dot{x}}^{2}+{\dot{{G}_{1}x}}^{3}), $ (2.9a)
    $ {F}_{c}\left(t\right) = -\varepsilon \left({p}_{2}z+\dot{{d}_{2}z}+{\alpha }_{4}{z}^{3}+{\alpha }_{5}{z}^{2}\dot{z}+{\alpha }_{6}z{\dot{z}}^{2}+{\dot{{G}_{2}z}}^{3}\right), $ (2.9b)

    where $ -({p}_{1}x+\dot{{d}_{1}x}), -({p}_{2}z+\dot{{d}_{2}z}) $ are the linear control forces, and $ -({\alpha }_{1}{x}^{3}+{\alpha }_{2}{x}^{2}\dot{x}+{\alpha }_{3}x{\dot{x}}^{2}) $, $ -({\alpha }_{4}{z}^{3}+{\alpha }_{5}{z}^{2}\dot{z}+{\alpha }_{6}z{\dot{z}}^{2}) $ are the non-linear control forces, and $ {G}_{1}, {G}_{2} $ stand the gain amounts. Figure 2 depicts the block chart of the model within the novel controller.

    Figure 2.  Block graph of the controlled electromechanical seismograph model with NPDGF controller.

    The perturbation technique [53] is investigated here to get the following approximate solutions as follows:

    $ x(t,\varepsilon ) = {x}_{0}({T}_{0},{T}_{1})+\varepsilon {x}_{1}({T}_{0},{T}_{1})+O\left({\varepsilon }^{2}\right), $ (3.1a)
    $ z(t,\varepsilon ) = {z}_{0}({T}_{0},{T}_{1})+\varepsilon {z}_{1}({T}_{0},{T}_{1})+O\left({\varepsilon }^{2}\right), $ (3.1b)

    Where, $ {T}_{0} = t $ and $ {T}_{1} = \varepsilon t $ represent time scales.

    The derivatives are presented as:

    $ \frac{d}{dt} = {D}_{0}+\varepsilon {D}_{1}+... $ (3.2a)
    $ \frac{{d}^{2}}{d{t}^{2}} = {D}_{0}^{2}+2\varepsilon {D}_{0}{D}_{1}+... $ (3.2b)

    Substituting Eq (3.1) into Eq (2.8) within Eqs (2.9, 3.2) and simply solving the expressions of the order $ O\left({\varepsilon }^{0}\right) $ and $ O\left({\varepsilon }^{1}\right) $ to be in the following relation:

    $D20x0+ω21x0+εω21x1+ε(D20x1+2D0D1x0)+εμm(D0x0)+ελ1x20+ελ2x30+εγ0(D0z0)+εγ1x0(D0z0)++εγ2x20(D0z0)εf0εf1cos(Ω1t)εf2x0cos(Ω2t)+εp1x0+εd1D0x0+εα1x30+εα2(x20(D0x0))+εα3x0(D0x0)2+εG1(D0x0)3=0,
    $
    (3.3a)
    $D20z0+ω22z0+εω22z1+ε(D20z1+2D0D1z0)+εμe(D0z0)+εβ0(D0x0)+εβ1z0(D0x0)+εβ2z20(D0x0)εf3cos(Ω3t)εf4z0cos(Ω4t)+εp2z0+εd2D0z0+εα4z30+εα5(z20(D0z0))+εα6z0(D0z0)2+εG2(D0z0)3=0.
    $
    (3.3b)

    Connecting terms of the same order of $ \varepsilon $ in Eq (3.3) to get:

    $ O\left({\varepsilon }^{0}\right):({D}_{0}^{2}+{\omega }_{1}^{2}){x}_{0} = 0, $ (3.4a)
    $ ({D}_{0}^{2}+{\omega }_{2}^{2}){z}_{0} = 0 , $ (3.4b)
    $O(ε1):(D20+ω21)x1=2D0D1x0μmD0x0λ1x20λ2x30γ0D0z0γ1x0(D0z0)γ2x20(D0z0)+f0+f1cos(Ω1t)+f2x0cos(Ω2t)p1x0d1(D0x0)α1x30α2x20(D0x0)α3x0(D0x0)2G1(D0x0)3,
    $
    (3.5a)
    $(D20+ω22)z1=2D0D1z0μeD0z0β0D0x0β1z0(D0x0)β2z20(D0x0)+f3cos(Ω3t)+f4z0cos(Ω4t)p2z0d2(D0z0)α4z30α5z20(D0z0)α6z0(D0z0)2G2(D0z0)3.
    $
    (3.5b)

    The solution for Eq (3.4) is

    $ {x}_{0} = {A}_{1}\left({T}_{1}\right){e}^{i{\omega }_{1}{T}_{0}}+{\stackrel{̄}{A}}_{1}\left({T}_{1}\right){e}^{-i{\omega }_{1}{T}_{0}} , $ (3.6a)
    $ {z}_{0} = {B}_{1}\left({T}_{1}\right){e}^{i{\omega }_{2}{T}_{0}}+{\stackrel{̄}{B}}_{1}\left({T}_{1}\right){e}^{-i{\omega }_{2}{T}_{0}} . $ (3.6b)

    Substituting Eq (3.6) into Eq (3.5) yields

    $ ({D}_{0}^{2}+{\omega }_{1}^{2}){x}_{1} = \left(2iω1D1A1iω1μmA13λ2A21ˉA1+12f2ˉA1ei(Ω22ω1)T0p1A1iω1d1A13α1A21ˉA12iω1α2A21ˉA1α3ω21A21ˉA13iG1ω31A21ˉA1
    \right){e}^{i{\omega }_{1}{T}_{0}}+NST+cc, $
    (3.7a)
    $ ({D}_{0}^{2}+{\omega }_{2}^{2}){z}_{1} = \left(2iω2D1B1iω2μeB1+12f4ˉB1ei(Ω4ω2)T0p2B1iω2d2B13α4B21ˉB12iω2α5B21ˉB1α6ω22B21ˉB13iG2ω32B21ˉB1
    \right){e}^{i{\omega }_{2}{T}_{0}}+NST+cc, $
    (3.7b)

    where $ cc $ represents the complex conjugate of the preceding parameters, $ NST $ signifies non-secular terms. By means of primary parametric resonance case $ {\varOmega }_{2}\cong {\omega }_{1}, {\varOmega }_{4}\cong {\omega }_{2} $, the frequency detuning parameters $ {\sigma }_{1}, {\sigma }_{2} $ for the parametric excitation frequency considered by

    $ {\varOmega }_{2}\cong {\omega }_{1}+\varepsilon {\sigma }_{1},{\varOmega }_{4}\cong {\omega }_{2}+\varepsilon {\sigma }_{2} . $ (3.8)

    In Eq (3.7), the solvability conditions produce when terms of secular terms are cancelled via Eq (3.8). Bearing in mind Eq (3.8), removing the secular terms in Eq (3.7) we get:

    $2iω1D1A1iω1μmA13λ2A21ˉA1+12f2ˉA1eiσ1T1p1A1iω1d1A13α1A21ˉA1iω1α2A21ˉA1α3ω21A21ˉA13iG1ω31A21ˉA1=0,
    $
    (3.9a)
    $2iω2D1B1iω2μeB1+12f4ˉB1eiσ2T0p2B1iω2d2B13α4B21ˉB1iω2α5B21ˉB1α6ω22B21ˉB13iG2ω32B21ˉB1=0.
    $
    (3.9b)

    To distinct the averaging conditions that manage the elements of Eq (3.9), let definite $ {A}_{1}, {B}_{1} $ and $ {\stackrel{̄}{A}}_{1}, {\stackrel{̄}{B}}_{1} $ are expressed in the next polar expressions

    $ {A}_{1} = \frac{1}{2}{a}_{1}\left({T}_{1}\right){e}^{i{\theta }_{1}\left({T}_{1}\right)} , {\stackrel{̄}{A}}_{1} = \frac{1}{2}{a}_{1}\left({T}_{1}\right){e}^{-i{\theta }_{1}\left({T}_{1}\right)} , $ (3.10a)
    $ {B}_{1} = \frac{1}{2}{a}_{2}\left({T}_{1}\right){e}^{i{\theta }_{2}\left({T}_{1}\right)} , {\stackrel{̄}{B}}_{1} = \frac{1}{2}{a}_{2}\left({T}_{1}\right){e}^{-i{\theta }_{2}\left({T}_{1}\right)} ,$ (3.10b)

    where $ {a}_{1}, {a}_{2} $ and $ {\theta }_{1}, {\theta }_{2} $ stand the steady-state amplitudes and phases, respectively. Substitute Eq (3.10) in Eq (3.9), we acquire:

    $iω1a'1+12ω1a1(σ1γ'1)12iω1μma138λ2a31+14a1f2(cosγ1+isinγ1)12p1a112iω1d1a138α1a3118iω1α2a3118α3ω21a3138iG1ω31a31=0,
    $
    (3.11a)
    $iω2a'2+12ω2a2(σ2γ'2)12iω2μea1+14a2f4(cosγ2+isinγ2)12p2a212iω2d2a238α4a3218iω2α5a3218α6ω22a3238iG2ω32a32=0,
    $
    (3.11b)

    where $ {\gamma }_{1} = {\sigma }_{1}{T}_{1}-2{\theta }_{1} $, $ {\gamma }_{2} = {\sigma }_{2}{T}_{1}-2{\theta }_{2} $, then separating real and imaginary elements:

    $ {a}_{1}^{\text{'}} = \frac{1}{2}\left[-{d}_{1}-{\mu }_{m}+\frac{1}{2{\omega }_{1}}{f}_{2}\mathit{sin}\;{\gamma }_{1}\right]{a}_{1}-\frac{1}{8}\left[{\alpha }_{2}+3{G}_{1}{\omega }_{1}^{2}\right]{a}_{1}^{3},$ (3.12a)
    $ {a}_{1}{\gamma }_{1}^{\text{'}} = \left[{\sigma }_{1}+\frac{1}{2{\omega }_{1}}{f}_{2}\mathit{cos}\;{\gamma }_{1}-\frac{{p}_{1}}{{\omega }_{1}}\right]{a}_{1}-\frac{1}{4}\left[\frac{3}{{\omega }_{1}}{\lambda }_{2}+\frac{3}{{\omega }_{1}}{\alpha }_{1}+{\alpha }_{3}{\omega }_{1}\right]{a}_{1}^{3} ,$ (3.12b)
    $ {a}_{2}^{\text{'}} = \frac{1}{2}\left[-{d}_{2}-{\mu }_{e}+\frac{1}{2{\omega }_{2}}{f}_{4}\mathit{sin}\;{\gamma }_{2}\right]{a}_{2}-\frac{1}{8}\left[{\alpha }_{5}+3{G}_{2}{\omega }_{2}^{2}\right]{a}_{2}^{3} , $ (3.13a)
    $ {a}_{2}{\gamma }_{2}^{\text{'}} = \left[{\sigma }_{2}+\frac{1}{2{\omega }_{2}}{f}_{4}\mathit{cos}\;{\gamma }_{2}-\frac{{p}_{2}}{{\omega }_{2}}\right]{a}_{2}-\frac{1}{4}\left[\frac{3}{{\omega }_{2}}{\alpha }_{4}+{\alpha }_{6}{\omega }_{2}\right]{a}_{2}^{3}.$ (3.13b)

    For steady-state responses ($ {a}_{1}^{\text{'}} = {a}_{2}^{\text{'}} = {\gamma }_{1}^{\text{'}} = {\gamma }_{2}^{\text{'}} = 0 $), the periodic solution resultant to Eqs (3.12) and (3.13) are given as:

    $ \left[-{d}_{1}-{\mu }_{m}+\frac{1}{2{\omega }_{1}}{f}_{2}\mathit{sin}\;{\gamma }_{1}\right]{a}_{1}-\frac{1}{4}\left[{\alpha }_{2}+3{G}_{1}{\omega }_{1}^{2}\right]{a}_{1}^{3} = 0,$ (3.14a)
    $ \left[{\sigma }_{1}+\frac{1}{2{\omega }_{1}}{f}_{2}\mathit{cos}\;{\gamma }_{1}-\frac{{p}_{1}}{{\omega }_{1}}\right]{a}_{1}-\frac{1}{4}\left[\frac{3}{{\omega }_{1}}{\lambda }_{2}+\frac{3}{{\omega }_{1}}{\alpha }_{1}+{\alpha }_{3}{\omega }_{1}\right]{a}_{1}^{3} = 0 ,$ (3.14b)
    $ \left[-{d}_{2}-{\mu }_{e}+\frac{1}{2{\omega }_{2}}{f}_{4}\mathit{sin}\;{\gamma }_{2}\right]{a}_{2}-\frac{1}{4}\left[{\alpha }_{5}+3{G}_{2}{\omega }_{2}^{2}\right]{a}_{2}^{3} = 0, $ (3.15a)
    $ \left[{\sigma }_{2}+\frac{1}{2{\omega }_{2}}{f}_{4}\mathit{cos}\;{\gamma }_{2}-\frac{{p}_{2}}{{\omega }_{2}}\right]{a}_{2}-\frac{1}{4}\left[\frac{3}{{\omega }_{2}}{\alpha }_{4}+{\alpha }_{6}{\omega }_{2}\right]{a}_{2}^{3} = 0. $ (3.15b)

    The Newton-Raphson method and MATLAB software are used to determine the steady-state answers from the algebraic equations. The Lyapunov first approach is used to determine the right-hand side eigenvalues of the Jacobian matrix at Eqs. (3.15) to calculate the stability of the steady-state shell system.

    $ \left[a'1γ'1a'2γ'2
    \right] = \left[R11R12R13R14R21R22R23R24R31R32R33R34R41R42R43R44
    \right]\left[a1γ1a2γ2
    \right], $
    (3.16)

    where

    $ R11=a1a1=12[d1μm+12ω1f2sinγ1]38[α2+3G1ω21]a21,
    $
    $ {R}_{12} = \frac{\partial {a}_{1}^{\text{'}}}{\partial {\gamma }_{1}} = \frac{1}{4{\omega }_{1}}{a}_{1}{f}_{2}\mathit{cos}\;{\gamma }_{1}, $
    $ {R}_{13} = \frac{\partial {a}_{1}^{\text{'}}}{\partial {a}_{2}} = 0, $
    $ {R}_{14} = \frac{\partial {a}_{1}^{\text{'}}}{\partial {\gamma }_{2}} = 0 , $
    $ {R}_{21} = \frac{\partial {\gamma }_{1}^{\text{'}}}{\partial {a}_{1}} = \frac{1}{{a}_{1}}\left[{\sigma }_{1}+\frac{1}{2{\omega }_{1}}{f}_{2}\mathit{cos}\;{\gamma }_{1}-\frac{{p}_{1}}{{\omega }_{1}}\right]-\frac{3}{4}\left[\frac{3}{{\omega }_{1}}{\lambda }_{2}+\frac{3}{{\omega }_{1}}{\alpha }_{1}+{\alpha }_{3}{\omega }_{1}\right]{a}_{1} , $
    $ {R}_{22} = \frac{\partial {\gamma }_{1}^{\text{'}}}{\partial {\gamma }_{1}} = -\frac{1}{2{\omega }_{1}}{f}_{2}\mathit{sin}\;{\gamma }_{1}, $
    $ {R}_{23} = \frac{\partial {\gamma }_{1}^{\text{'}}}{\partial {a}_{2}} = 0, $
    $ {R}_{24} = \frac{\partial {\gamma }_{1}^{\text{'}}}{\partial {\gamma }_{2}} = 0 , $
    $ {R_{31}} = \frac{{\partial {{a'}_2}}}{{\partial {a_1}}} = 0,\,{R_{32}} = \frac{{\partial {{a'}_2}}}{{\partial {\gamma _1}}} = 0,\, $
    $ {R}_{13} = \frac{\partial {a}_{2}^{\text{'}}}{\partial {a}_{2}} = \frac{1}{2}\left[-{d}_{2}-{\mu }_{e}+\frac{1}{2{\omega }_{2}}{f}_{4}\mathit{sin}\;{\gamma }_{2}\right]-\frac{3}{8}\left[{\alpha }_{5}+3{G}_{2}{\omega }_{2}^{2}\right]{a}_{2}^{2}, $
    $ {R}_{14} = \frac{\partial {a}_{2}^{\text{'}}}{\partial {\gamma }_{2}} = \frac{1}{4{\omega }_{2}}{a}_{2}{f}_{4}\mathit{cos}\;{\gamma }_{2} , $
    $ {R}_{41} = \frac{\partial {\gamma }_{2}^{\text{'}}}{\partial {a}_{1}} = 0, $
    $ {R}_{42} = \frac{\partial {\gamma }_{2}^{\text{'}}}{\partial {\gamma }_{1}} = 0, $
    $ {R}_{43} = \frac{\partial {\gamma }_{2}^{\text{'}}}{\partial {a}_{2}} = \frac{1}{{a}_{2}}\left[{\sigma }_{2}+\frac{1}{2{\omega }_{2}}{f}_{4}\mathit{cos}\;{\gamma }_{2}-\frac{{p}_{2}}{{\omega }_{2}}\right]-\frac{3}{4}\left[\frac{3}{{\omega }_{2}}{\alpha }_{4}+{\alpha }_{6}{\omega }_{2}\right]{a}_{2}, $
    $ {R}_{44} = \frac{\partial {\gamma }_{2}^{\text{'}}}{\partial {\gamma }_{2}} = -\frac{1}{2{\omega }_{2}}{f}_{4}\mathit{sin}\;{\gamma }_{2} . $

    To determine the controlled model's stable zones, calculate the following determinant in the previous matrix.

    $ \left|R11λR1200R21R22λ0000R33λR3400R43R44λ
    \right| = 0 . $
    (3.17)

    Then,

    $ {\lambda }^{4}+{r}_{1}{\lambda }^{3}+{r}_{2}{\lambda }^{2}+{r}_{3}\lambda +{r}_{4} = 0. $ (3.18)

    Where, $ \lambda $ designates the Jacobian matrix's eigenvalue,

    $ {r}_{1} = -{R}_{11}-{R}_{22}-{R}_{33}-{R}_{44} , $
    $ {r_2} = {R_{11}}{R_{22}} + {R_{11}}{R_{33}} + {R_{11}}{R_{44}} + {R_{22}}{R_{33}} + {R_{22}}{R_{44}} + {R_{33}}{R_{44}} - {R_{12}}{R_{21}} - {R_{34}}{R_{43}}, $
    $ {r}_{3} = {R}_{12}{R}_{21}{R}_{33}+{R}_{12}{R}_{21}{R}_{44}+{R}_{11}{R}_{34}{R}_{43}+{R}_{22}{R}_{34}{R}_{43} $
    $ -{R}_{11}{R}_{22}{R}_{33}-{R}_{11}{R}_{22}{R}_{44}-{R}_{11}{R}_{33}{R}_{44}-{R}_{22}{R}_{33}{R}_{44}, $
    $ {r}_{4} = {R}_{11}{R}_{22}{R}_{33}{R}_{44}+{R}_{12}{R}_{21}{R}_{34}{R}_{43}-{R}_{11}{R}_{22}{R}_{34}{R}_{43}-{R}_{12}{R}_{21}{R}_{33}{R}_{44} $

    are the coefficient of Eq (3.18).

    Routh-Hurwitz criterion is used to examine the stability of equilibrium solutions via discovering sufficient and necessary requirements to be stable are:

    $ {r}_{1} > 0,{r}_{1}{r}_{2}-{r}_{3} > 0,{r}_{3}({r}_{1}{r}_{2}-{r}_{3})-{r}_{1}^{2}{r}_{4} > 0,{r}_{4} > 0. $ (3.19)

    Furthermore, to conclude the nature of the roots of Eq (3.18), we define the following discriminants $ \varDelta $, $ {\varDelta }_{0} $, $ {\varDelta }_{1} $, $ {\varDelta }_{2} $, and $ {\varDelta }_{3} $ [58,59] as:

    $Δ=256r34192r1r3r24128r22r24+144r2r23r427r43+144r21r2r246r21r23r480r1r22r3r4+18r1r2r33+16r42r44r32r2327r41r24+18r31r2r3r44r31r334r21r32r4+r21r22r23Δ0=8r33r21,Δ1=r31+8r234r1r2,Δ2=r223r1r3+12r4,Δ3=64r416r22+16r21r216r1r33r41.
    $

    By joining the conditions Eq (3.19) with each one of the following possible cases, one can establish the sort of roots of Eq (3.18) as:

    a) If $ \varDelta < 0 $ then the equation has two different real and two complex conjugate roots.

    b) If $ \varDelta > 0 $, then either the four roots are all complex conjugate or real according to the subsequent two cases:

    ⅰ. If $ {\varDelta }_{0} < 0 $ and $ {\varDelta }_{3} < 0 $ then all four roots are different and real.

    ⅱ. If $ {\varDelta }_{0} > 0 $ or if $ {\varDelta }_{3} > 0 $ then the roots are two pairs of complex conjugates.

    c) If $ \varDelta = 0 $ then the system has multiple roots according to the following four cases:

    ⅰ. If $ {\varDelta }_{0} < 0 $, $ {\varDelta }_{3} < 0 $ and $ {\varDelta }_{2}\ne 0 $, there are two real simple roots and a real double roots

    ⅱ. If $ {\varDelta }_{3} > 0 $ or ($ {\varDelta }_{0} > 0 $ and ($ {\varDelta }_{3}\ne 0 $ or $ {\varDelta }_{1}\ne 0 $)), there are two complex conjugate roots and two real equal roots.

    ⅲ. If $ {\varDelta }_{2} = 0 $ and $ {\varDelta }_{3}\ne 0 $, there are three real equal roots and one real different root.

    ⅳ. If $ {\varDelta }_{3} = 0 $, at that time:

    1- If $ {\varDelta }_{0} < 0 $, there are two real double roots.

    2- If $ {\varDelta }_{0} > 0 $ and $ {\varDelta }_{1} = 0 $, the roots be present two double complex conjugate roots.

    3- If $ {\varDelta }_{2} = 0 $, all four roots are equal to $ -\frac{{r}_{1}}{4} $.

    The mechanism of the NPDVF controller at the measured simultaneous resonance can be clarified with the aid of Eqs (3.12) and (3.13). It is clear from Eqs (3.12) and (3.13) that the addition of the NPDVF controller to the considered framework has modified the linear damping terms $ {\mu }_{m} $ and $ {\mu }_{e} $ to the controlled terms $ {\mu }_{m-control} $ and $ {\mu }_{e-control} $. Moreover, the detuning parameter $ {\sigma }_{1} $ and $ {\sigma }_{2} $ is modified to $ {\sigma }_{1-control} $ and $ {\sigma }_{2-control} $, where $ {\mu }_{m-control} $, $ {\mu }_{e-control} $, $ {\sigma }_{1-control} $, and $ {\sigma }_{2-control} $ are given as follows:

    $ {\mu }_{m-control} = -\frac{1}{2}\left({\mu }_{m}+{d}_{1}\right) , $ (4.1)
    $ {\sigma }_{1-control} = \left({\sigma }_{1}-\frac{{p}_{1}}{2{\omega }_{1}}\right), $ (4.2)
    $ {\mu }_{e-control} = -\frac{1}{2}\left({\mu }_{e}+{d}_{2}\right), $ (4.3)
    $ {\sigma }_{2-control} = \left({\sigma }_{2}-\frac{{p}_{2}}{2{\omega }_{2}}\right). $ (4.4)

    It is clear from Eqs (4.1)–(4.4) that $ {\mu }_{m-control} $, $ {\mu }_{e-control} $, $ {\sigma }_{1-control} $, and $ {\sigma }_{2-control} $ are periodic functions on the controller, where the controlled system has the equivalent linear damping coefficients $ {\mu }_{m-control} = -\frac{1}{2}\left({\mu }_{m}+{d}_{1}\right) $, $ {\mu }_{e-control} = -\frac{1}{2}\left({\mu }_{e}+{d}_{2}\right) $ and the detuning parameters $ {\sigma }_{1-control} = \left({\sigma }_{1}-\frac{{p}_{1}}{2{\omega }_{1}}\right) = {\varOmega }_{2}-\left({\omega }_{1}+\frac{{p}_{1}}{2{\omega }_{1}}\right) $, $ {\sigma }_{2-control} = \left({\sigma }_{2}-\frac{{p}_{2}}{2{\omega }_{2}}\right) = {\varOmega }_{4}-\left({\omega }_{2}+\frac{{p}_{2}}{2{\omega }_{2}}\right) $. This means that the linear control forces $ {p}_{1} $ and $ {p}_{2} $ are responsible for changing the system natural frequencies $ {\omega }_{1} $ and $ {\omega }_{2} $, while the velocity gain $ {d}_{1} $ and $ {d}_{2} $ is responsible for modifying the system linear damping coefficients $ {\mu }_{m} $ and $ {\mu }_{e} $. Consequently, to develop the vibration suppression efficiency of the measured scheme, the linear control forces $ {p}_{1}, {p}_{2} $ and $ {d}_{1}, {d}_{2} $ should be designated in a way that maximizes the objective function $ {\mu }_{m-control} $, $ {\mu }_{e-control} $, $ {\sigma }_{1-control} $ and $ {\sigma }_{2-control} $. By comparing the obtained results in Figures 8(b) and 9(a) with the objective function given by Eqs (4.1)–(4.4), we can notice that the best vibration suppression condition has occurred at the maximum values of function $ {\mu }_{m-control} $, $ {\mu }_{e-control} $, $ {\sigma }_{1-control} $ and $ {\sigma }_{2-control} $ as summarized in Table 1.

    Table 1.  Optimum control parameters.
    Figure $ {p}_{1} $ $ {p}_{2} $ $ {d}_{1} $ $ {d}_{2} $ $ {\omega }_{1} $ $ {\omega }_{2} $ $ {\mu }_{m-control}=\frac{1}{2}\left({\mu }_{m}+{d}_{1}\right) $ $ {\sigma }_{1-control}=\left({\sigma }_{1}-\frac{{p}_{1}}{2{\omega }_{1}}\right) $ $ {\mu }_{e-control}=\frac{1}{2}\left({\mu }_{e}+{d}_{2}\right) $ $ {\sigma }_{2-control}=\left({\sigma }_{2}-\frac{{p}_{2}}{2{\omega }_{2}}\right) $ Max
    $ {\mu }_{m-control} $
    Max
    $ {\mu }_{e-control} $
    Max
    $ {\sigma }_{1-control} $
    Max
    $ {\sigma }_{2-control} $
    16 -0.5 0 0 0 5.05 4.05 0.5$ \left({\mu }_{m}+0\right) $ $ {\sigma }_{1}+\frac{0.5}{2\left(5.05\right)} $ 0.5$ \left({\mu }_{e}+0\right) $ $ {\sigma }_{2}+0 $ 0.0155 0.2 1.05 1.0
    24 0 -0.5 0 0 5.05 4.05 0.5$ \left({\mu }_{m}+0\right) $ $ {\sigma }_{1}+0 $ 0.5$ \left({\mu }_{e}+0\right) $ $ {\sigma }_{2}+\frac{0.5}{2\left(4.05\right)} $ 0.0155 0.2 1.0 1.06
    17 0 0 2.5 0 5.05 4.05 0.5$ \left({\mu }_{m}+2.5\right) $ $ {\sigma }_{1}+0 $ 0.5$ \left({\mu }_{e}+0\right) $ $ {\sigma }_{2}+0 $ 1.2655 0.2 1.0 1.0
    25 0 0 0 0.5 5.05 4.05 0.5$ \left({\mu }_{m}+0\right) $ $ {\sigma }_{1}+0 $ 0.5$ \left({\mu }_{e}+0.5\right) $ $ {\sigma }_{2}+0 $ 0.0155 0.45 1.0 1.0

     | Show Table
    DownLoad: CSV

    The table 1 provides a list or summary of the important variables that were considered in the design process. The choice of each parameter was guided by the system's functional requirements and manufacturing constraints to avoid unstable regions of the system. This approach avoids pushing the system into regions that might otherwise lead to increased costs, manufacturing difficulties, or operational instability. The parameters listed in Tab. 1 are critical to the overall design of the system. Each of these variables was carefully chosen to balance the system's functional requirements with the practical limitations imposed by manufacturing capabilities.

    To improve the system's response, Eqs (4.1)–(4.4) introduce more sophisticated feedback mechanisms, specifically designed to manage these nonlinearities and provide better stabilization and vibration suppression. The enhanced NPDVF controller improves the system's ability to avoid resonance conditions or mitigate their effects. The improvement of the NPDCVF controller represented by Eqs (4.1)–(4.4) is essential for addressing the nonlinear dynamics, feedback instability, and resonance issues in systems with both mechanical and electrical components. The combination of these equations allows the NPDVF controller to anticipate future disturbances, adjust for system nonlinearities, and optimize control actions in real time.

    The nonlinear dynamical structure was just demonstrated by Eq (3.10) through (3.14) above. The MATLAB®18 computer programmer was then used to numerically simulate three distinct control strategies (PD, NCVF, and NPDVF) to identify which controller would minimise the destructive vibrations caused during work of the model. Figures 36 are plotted to appear the time history of the worst resonance case before and after using various controllers via the following values of the parameters:

    $ {\mu }_{m} = 0.031,{\omega }_{1} = 5.05,{\lambda }_{1} = 0.00315,{\lambda }_{2} = 0.025, $
    $ {\gamma }_{0} = -0.2,{\gamma }_{1} = -0.0015,{\gamma }_{2} = 0.035,{f}_{0} = 1.2,{f}_{1} = 3.5, $
    $ {f}_{2} = 4.7,{\varOmega }_{1} = 3.75,{\varOmega }_{2} = {\omega }_{1},{\mu }_{e} = 0.4,{\omega }_{2} = 4.05, $
    $ {\beta }_{0} = 0.0072,{\beta }_{1} = 0.005,{\beta }_{2} = 0.0015,{f}_{3} = 1.6,{f}_{4} = 3.1, $
    $ {\varOmega }_{3} = 2.45,{\varOmega }_{4} = {\omega }_{2},{p}_{1} = 2.3,{d}_{1} = 2.5, $
    $ {\alpha }_{1} = 0.09,{\alpha }_{2} = 0.04,{\alpha }_{3} = 0.4,{G}_{1} = 0.5,{p}_{2} = 1.9,{d}_{2} = 0.5, $
    $ {\alpha }_{4} = 0.05,{\alpha }_{5} = 0.03,{\alpha }_{6} = 0.2,{G}_{2} = 0.4,\varepsilon = 0.5 , $

    with the initial conditions $ x\left(0\right) = 0.5, \dot{x}\left(0\right) = 0, z\left(0\right) = 0.5, \dot{z}\left(0\right) = 0. $

    Figure 3.  The measured simultaneous resonance instance ($ {\varOmega }_{2}\cong {\omega }_{1}, {\varOmega }_{4}\cong {\omega }_{2} $) without controller.
    Figure 4.  The measured simultaneous resonance instance ($ {\varOmega }_{2}\cong {\omega }_{1}, {\varOmega }_{4}\cong {\omega }_{2} $) with PD controller.
    Figure 5.  The measured simultaneous resonance instance ($ {\varOmega }_{2}\cong {\omega }_{1}, {\varOmega }_{4}\cong {\omega }_{2} $) with NCVF controller.
    Figure 6.  The measured simultaneous resonance instance ($ {\varOmega }_{2}\cong {\omega }_{1}, {\varOmega }_{4}\cong {\omega }_{2} $) with the novel resonant controller NPDVF.

    Figure 3 presents the basic steady-state amplitudes $ x\left(t\right) $ and $ z\left(t\right) $ without controller at the simultaneous resonance situation $ {\varOmega }_{2}\cong {\omega }_{1}, {\varOmega }_{4}\cong {\omega }_{2} $ as the worst resonance case of the system. Figures 46 show the results of adding the measured controllers (PD-NCVF-NPDCVF), enabling you to indicate the best one for reducing the high vibration amplitudes. After checking, we found that the best one controller is NPDVF as it reduces the vibrations at the measured resonance in a short time. The new controller NPDCVFC is the most efficient, according to these diagrams (3–6). The controller for the observed model's active vibration is developed in this part (refer to Figure 6). The NPD and NCVFC algorithms are used for active vibration control, as mentioned earlier. With this controller rule, the depreciation rate's are enhanced for this structural model. It is obvious that the infinite norm of vibration amplitudes might be reduced with this straightforward method. The outcomes demonstrate how well the optimization strategy reduced vibrations and how quickly suitably positioned actuators and sensors were able to reduce vibrations in the model. Thus, we will select and statistically evaluate the controlled model in order to study it and examine the impact of various controlled parameters on the system.

    This part examined the effects of different parameters used the regulated model in Eqs (3.14-3.15), numerically illustrated the stable and unstable areas. The beneficial case is examined to acquire a significant quantity of parameter effects. As illustrated in Figures 726, all curves display only stable segments without any instability areas when NPDVF controller is applied to the system. This provides an additional justification for adding this supplementary controller to the system, which is a beneficial outcome for any vibrating system. Solid curves indicate stable reactions. The solid line () illustrates the stable regions. Figures 7 and 19 describe the performance of the amplitude-frequency response steady-state response curves $ {a}_{1} $ versus $ {\sigma }_{1} $ and $ {a}_{2} $ versus $ {\sigma }_{2} $, respectively, with nonstable regions, which shows only stable zones. As the linear damping coefficients $ {\mu }_{m}, {\mu }_{e} $ decreased the amplitudes $ {a}_{1}, {a}_{2} $ are increased and the stability region is increased as appeared in Figures 8, 20. By way of the excitation frequencies $ {\omega }_{1}, {\omega }_{2} $ are increased the steady-state amplitudes $ {a}_{1}, {a}_{2} $ are increased and all regions are stable as exposed in Figures 9, 10, 21. Moreover, the diagrams of the nonlinear control coefficient $ {\alpha }_{2} $ is slightly monotonic decreasing with no instability regions as depicted in Figure 11. Moreover, when the amplitudes $ {a}_{1}, {a}_{2} $ are decreased the nonlinear control parameter $ {\alpha }_{3} $ is increased as illustrated in Figures 12, 22. Also, the nonlinear coefficient $ {\alpha }_{4} $ is monotonic decreasing in which when it decreased, the amplitude $ {a}_{1} $ increased as shown in Figure 13. Also, as the parametric force coefficients $ {f}_{2}, {f}_{4} $ are increased, then the steady-state amplitudes $ {a}_{1}, {a}_{2} $ are increased with increasing in stability regions as displayed in Figures 14, 15 and 23. Besides, the curve of the linear control force coefficients $ {p}_{1}, {p}_{2} $ are shifted to right (S.R) when the values of parameters $ {p}_{1}, {p}_{2} $ are increased with no instability regions as graphed in Figures 16 and 24. On other hand, when the linear control force coefficients $ {d}_{1} $ is increased, the steady-state amplitudes $ {a}_{1}, {a}_{2} $ are increased and all zones are stable as labelled in Figures 17 and 25. In the last, as the gain coefficients $ {G}_{1}, {G}_{2} $ increased the steady-state amplitudes $ {a}_{1}, {a}_{2} $ are small with stability regions which are presented in Figures 18 and 26.

    Figure 7.  Response curve via $ {a}_{1} $ versus $ {\sigma }_{1} $ of the controlled structure.
    Figure 8.  Impact of damping coefficient $ {\mu }_{m} $ (M.D).
    Figure 9.  Effect of excitation frequency $ {\omega }_{1} $ (M.I).
    Figure 10.  Influence of excitation frequency $ {\omega }_{2} $ (S.M.I) slightly monotonic Increasing.
    Figure 11.  Influence of nonlinear control parameter $ {\alpha }_{2} $ (S.M.D) slightly monotonic decreasing.
    Figure 12.  Influence of nonlinear control parameter $ {\alpha }_{3} $ (M.D).
    Figure 13.  Influence of nonlinear control parameter $ {\alpha }_{4} $(M.D).
    Figure 14.  Influence of parametric excitation force parameter $ {f}_{2} $ (M.I).
    Figure 15.  Influence of parametric excitation force parameter $ {f}_{4} $ (M.D).
    Figure 16.  Influence of the linear control force coefficient $ {p}_{1} $ at increasing amplitude moved to right, at decreasing amplitude moved to left.
    Figure 17.  Influence of the linear control force coefficient $ {d}_{1} $(M.I).
    Figure 18.  Influence of nonlinear control gain $ {G}_{1} $ (M.D).
    Figure 19.  Influence response curves $ {a}_{2} $ versus $ {\sigma }_{2} $ of the controlled system.
    Figure 20.  Impact of the damping coefficient $ {\mu }_{e} $ (M.D).
    Figure 21.  Effect of the excitation frequency $ {\omega }_{2} $ (M.I).
    Figure 22.  Influence of nonlinear control parameter $ {\alpha }_{3} $ (M.D).
    Figure 23.  Influence of the parametric excitation force parameter $ {f}_{2} $ (M.I).
    Figure 24.  Influence of the linear control force factor $ {p}_{2} $, at increasing amplitude it moved to right, but at decreasing amplitude it shifted to left.
    Figure 25.  Influence of the linear control force coefficient $ {d}_{2} $ (M.I).
    Figure 26.  Influence of the nonlinear control gain $ {G}_{2} $ (M.D).

    As listed in the figures, in which, M.I refers to monotonic increasing and M.D refers to monotonic decreasing

    Therefore, Table 2 provides a detailed account of the effects of various parameters of the controlled system with various effects on the system to appear the stable regions.

    Table 2.  Impact of various coefficients of the controlled system.
    Parameter Symbol Status Amplitude Effect Stability Figure No.
    Linear damping $ {\mu }_{m} $ Increase $ {a}_{1} $ Decrease Stable regions occur 8
    Excitation frequency $ {\omega }_{1} $ Increase $ {a}_{1} $ Increase Stable regions occur 9
    Excitation frequency $ {\omega }_{2} $ Increase $ {a}_{1} $ Slightly Increase Stable regions occur 10
    Nonlinear control $ {\alpha }_{2} $ Increase $ {a}_{1} $ Slightly Decrease Stable regions occur 11
    Nonlinear control $ {\alpha }_{3} $ Increase $ {a}_{1} $ Decrease Stable regions occur 12
    Nonlinear control $ {\alpha }_{4} $ Increase $ {a}_{1} $ Decrease Stable regions occur 13
    Parametric excitation force $ {f}_{2} $ Increase $ {a}_{1} $ Increase Stable regions occur 14
    parametric excitation force $ {f}_{4} $ Increase $ {a}_{1} $ Decrease Stable regions occur 15
    Linear control force $ {p}_{1} $ $ \frac{\rm Increase}{\rm Decrease}$ $ {a}_{1} $ $\frac{\text{Shifted to right}}{\text{Shifted to left}} $ Stable regions occur 16
    Linear control force $ {d}_{1} $ Increase $ {a}_{1} $ Increase Stable regions occur 17
    Gain $ {G}_{1} $ Increase $ {a}_{1} $ Decrease Stable regions occur 18
    Linear damping $ {\mu }_{e} $ Increase $ {a}_{2} $ Decrease Stable regions occur 20
    Excitation frequency $ {\omega }_{2} $ Increase $ {a}_{2} $ Increase Stable regions occur 21
    Nonlinear control $ {\alpha }_{3} $ Increase $ {a}_{2} $ Decrease Stable regions occur 22
    Parametric excitation force $ {f}_{2} $ Increase $ {a}_{2} $ Increase Stable regions occur 23
    Linear control force $ {p}_{2} $ $ \frac{\text{Increase}}{\text{Decrease}}$ $ {a}_{2} $ $ \frac{\text{Shifted to right}}{\text{Shifted to left}}$ Stable regions occur 24
    Linear control force $ {d}_{2} $ Increase $ {a}_{2} $ Increase Stable regions occur 25
    Gain $ {G}_{2} $ Increase $ {a}_{2} $ Decrease Stable regions occur 26

     | Show Table
    DownLoad: CSV

    An electromechanical seismograph model associated electric circuit model system akin to Eq (2.7) was examined in [9,10,11]. However, they applied the multiple time scale process within the principal parametric resonance item to study the behaviour of the structure under mixed excitations, eliminating the need for a controller. The development of the model given in [11] is examined in the current work. I add a variety of control strategies to the vibrating structure system's modified system in order to determine which one reduces the framework structure's risk of vibration. Additionally, the upgraded system's new controller NPDVF is examined in this article. The results of this research show that the novel controller less than the other controllers decreases the high vibrational amplitude of the model exposed to parametric stimulation inside the simultaneous resonance, as shown in subsection 5.1. The perturbation approach is used to aid in the acquisition of analytical solutions. Plotting of the frequency response graphs occurs at different framework parameter levels. We end with a numerical validation of the obtained results. The comparisons show that the current approach produces findings that are remarkably similar to those found in [11] and that the discrepancies are less than 1%.

    A comparison between nonlinear controllers with the proposed NPDVF controller is presented in a summary Table 3 that compares the performance of the different controllers based on the criteria of robustness, computational cost, and control signal amplitude. The table should provide a clear visual comparison that illustrates the advantages of the proposed NPDVF controller.

    Table 3.  Comparison NPDVF controller with other controllers of nonlinear systems.
    Control Method Robustness Computational Cost Control Signal Amplitude
    NPDVF Controller Excellent (handles nonlinearities and bifurcations) Low (efficient and fast) Small (efficient suppression)
    Linear Proportional-Derivative (PD)
    Controller
    Moderate (struggles with nonlinearities) Very Low (simple) Large (less adaptive)
    Sliding Mode Controller (SMC) High (robust but chattering) Moderate (requires switching) Large (due to chattering)
    Fuzzy Logic Controller Moderate (requires rule tuning) Moderate (rule evaluation) Moderate (tuned control)
    Model Predictive Controller (MPC) High (optimal performance but computationally expensive) Very High (optimization) Moderate to Large (depends on prediction)
    Backstepping Controller High (effective for nonlinear systems) Moderate (recursive equations) Moderate to Large (strong control action)

     | Show Table
    DownLoad: CSV

    The time response of control inputs is shown in Figures 2732 where the control input is also used to eliminate the effect of disturbance mainly at the stable stage. The control inputs is not only regulating the system's behavior but is also designed to mitigate the effect of external disturbances, particularly during the stable stage of the system's operation. These figures will show how the control input varies over time as the system responds to a disturbance. Each figure shows how the control input changes over time and responds to disturbances, providing a clear visualization of the controller's effectiveness in rejecting disturbances. This suggests that the system can quickly react to external disturbances, with the linear control force $ {p}_{1} = 2.3 $ stabilizing at a specific value after some initial disturbance. The system reaches stability within a short time, indicating a fast disturbance rejection process, and the disturbance effect is effectively mitigated as in Figure 27. Also, Figure 28 shows that the system achieves a disturbance-free state after a brief period. The disturbance input $ {d}_{1} = 2.5 $ ​has been neutralized quickly, indicating that the controller effectively nullifies the impact of the disturbance on the system in a short time. The nonlinear control gain $ {G}_{1} = 0.5 $ is an important factor in reducing disturbances. When this gain is set to 0.5, it appears that the system is better able to suppress disturbances. This indicates that the nonlinear controller is tuned to counteract the disturbances effectively, ensuring that the system can reach a steady state despite the presence of external influences as in Figure 29. Additionally, the control parameters $ {p}_{2} $​ and $ {d}_{2} $ appear to be carefully selected to ensure the system reacts appropriately to disturbances while maintaining stability. The values of $ {p}_{2} = 1.9 $ and $ {d}_{2} = 0.5 $ are likely optimized for quick disturbance rejection and minimal control effort. These parameters ensure that the system reaches the desired state quickly, further reducing the influence of external disturbances as shown in Figures 30 and 31. Finally, When $ {G}_{1} = 0.4 $ is applied, the system is better able to reduce or even eliminate the impact of external disturbances. The system's response, particularly in terms of the control input, should show a clear improvement after this nonlinear gain is activated, resulting in a smoother and more stable control input over time as illustrated in Figure 32. By effectively control input parameters such as $ {p}_{1}, {p}_{2}, {d}_{1}, {d}_{2}, {G}_{1, } $ and $ {G}_{2} $​, the system is able to reject disturbances in a timely manner and ensure robust performance.

    Figure 27.  Time response of linear control force $ {p}_{1} $ ( 2.3, 3.5, 5.5).
    Figure 28.  Time response of linear control force $ {d}_{1} $ ( 2.5, 4.5, 6.5).
    Figure 29.  Time response of nonlinear control gain $ {G}_{1} $ ( 0.5, 1.5, 3.5).
    Figure 30.  Time response of linear control force $ {p}_{2} $ ( 1.9, 3.2, 4.4).
    Figure 31.  Time response of linear control force $ {d}_{2} $ ( 0.5, 1.5, 3.5).
    Figure 32.  Time response of nonlinear control gain $ {G}_{2} $ ( 0.4, 3.3, 1.5).

    This study introduces a novel nonlinear proportional-derivative cubic velocity feedback (NPDVF) Controller that effectively manages vibrations in systems with coupled mechanical and electrical components. The controller's performance has been validated through numerical simulations and perturbation methods, demonstrating its ability to reduce vibrations, stabilize nonlinear motions, and suppress unwanted oscillations in the presence of mixed forces. A number of controller design approaches (PD control, NCVF control, and NPDVF as an innovative control method) were assessed in order to ascertain which one best minimizes high amplitude vibrations during the simultaneous resonance case $ {\varOmega }_{2}\cong {\omega }_{1}, {\varOmega }_{4}\cong {\omega }_{2} $. The solution of the studied controlled model can be approximated using the perturbation approach. The key advantages of the NPDVF controller, including its robustness to parameter uncertainties and external disturbances, low computational cost, and efficient control signal amplitude, make it a promising solution for nonlinear vibration control systems. A comprehensive comparison with traditional control strategies, such as PD control, sliding mode control (SMC), and model predictive control (MPC), highlights the superior performance of the NPDVF controller in several aspects. Specifically, the NPDVF controller outperforms others in terms of computational efficiency and robustness, while also maintaining a smaller control signal amplitude, which reduces actuator wear and improves energy efficiency. These results suggest that the NPDVF controller is an effective and practical solution for real-time vibration control in complex systems. This study advances our understanding of the control dynamics in nonlinear models with combined excitations via perturbation technique for controlling chaotic behavior in such models. Future work could investigate the impact of time delay on the effectiveness of the NPDVF controller. Specifically, it would be beneficial to analyze how time delay influences the system's stability, vibration reduction, and bifurcation control.

    The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.

    The author extends his appreciation to Prince Sattam bin Abdulaziz University for funding this research work through the project number (PSAU/ 2024/01/29249).

    The author declares that he has no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    [1] Clarke MC, Figg N, Maguire JJ, et al. (2006) Apoptosis of vascular smooth muscle cells induces features of plaque vulnerability in atherosclerosis. Nat Med 12: 1075-1080. doi: 10.1038/nm1459
    [2] Scott S, O'Sullivan M, Hafizi S, et al. (2002) Human vascular smooth muscle cells from restenosis or in-stent stenosis sites demonstrate enhanced responses to p53: implications for brachytherapy and drug treatment for restenosis. Circ Res 90: 398-404. doi: 10.1161/hh0402.105900
    [3] Erl W (2005) Statin-induced vascular smooth muscle cell apoptosis: a possible role in the prevention of restenosis? Curr Drug Targets Cardiovasc Haematol Disord 5: 135-144. doi: 10.2174/1568006043586134
    [4] Liu X (2006) Recent development of Artesunate. Chinese Journal of New Drugs 15:1918-1923.
    [5] Chen HH, Zhou HJ, Wu GD, et al. (2004) Inhibitory effects of Artesunate on angiogenesis and on expressions of vascular endothelial growth factor and VEGF receptor KDR/flk-1. Pharmacology 71: 1-9. doi: 10.1159/000076256
    [6] Berger TG, Dieckmann D, Efferth T, et al. (2005) Artesunate in the treatment of metastatic uveal melanoma——first experiences. Oncol Rep 14: 1599-1603.
    [7] Wu GD, Zhou HJ, Wu XH, et al. (2004) Apoptosis of human umbilical vein endothelial cells induced by artesunate. Vascul Pharmacol 41: 205-212. doi: 10.1016/j.vph.2004.11.001
    [8] He RR, Zhou HJ (2008) Progress in research on the anti-tumor effect of Artesunate. Chin J Integr Med 14: 312-316. doi: 10.1007/s11655-008-0312-0
    [9] Li S, Xue F, Cheng Z, et al. (2009) Effect of Artesunate on inhibiting proliferation and inducing apoptosis of SP2/0 myeloma cells through affecting NFkappaB p65. Int J Hematol 90: 513-521. doi: 10.1007/s12185-009-0409-z
    [10] Zheng JS, Wang MH, Huang M, et al. (2008) Artesunate suppresses human endometrial carcinoma RL95-2 cell proliferation by inducing cell apoptosis. Nan Fang Yi Ke Da Xue Xue Bao 28: 2221-2223.
    [11] Hou J, Wang D, Zhang R, et al. (2008) Experimental therapy of hepatoma with artemisinin and its derivatives: in vitro and in vivo activity, chemosensitization, and mechanisms of action. Clin Cancer Res 14: 5519-5530. doi: 10.1158/1078-0432.CCR-08-0197
    [12] Sieber S, Gdynia G, Roth W, et al. (2009) Combination treatment of malignant B cells using the anti-CD20 antibody rituximab and the anti-malarial Artesunate. Int J Oncol 35: 149-158.
    [13] Xiao SH, Mei JY, Jiao PY (2011) Effect of mefloquine administered orally at single, multiple, or combined with artemether, Artesunate, or praziquantel in treatment of mice infected with Schistosoma japonicum. Parasitol Res 108: 399-406. doi: 10.1007/s00436-010-2080-y
    [14] Keiser J, Xiao S, Smith TA, et al. (2009) Combination Chemotherapy against Clonorchis sinensis: Experiments with Artemether, Artesunate, OZ78, Praziquantel, and Tribendimidine in a Rat Model. Antimicrob Agents Chemother 53: 3770-3776. doi: 10.1128/AAC.00452-09
    [15] Rinner B, Siegl V, Pürstner P, et al. (2004) Activity of novel plant extracts against medullary thyroid carcinoma cells. Anticancer Res 24: 495-500.
    [16] Zhou Z, Feng Y (2005) Artesunate reduces proliferation, interferes DNA replication and cell cycle and enhances apoptosis in vascular smooth muscle cells. J Huazhong Univ Sci Technolog Med Sci 25: 135-136. doi: 10.1007/BF02873558
    [17] Liao H (2006) Effects of Artesunate on the the proliferation and apoptosis of rat aortic vascular smooth muscle cells in vitro. Huazhong University of Science and Technology [Ph.D Dissertation].
    [18] Pan W, da Graca LS, Shao Y, et al. (2009) PHAPI/pp32 suppresses tumorigenesis by stimulating apoptosis. J Biol Chem 284: 6946-6954. doi: 10.1074/jbc.M805801200
    [19] LuYY, Chen TS, Qu JL, et al. (2009) Dihydroartemisinin (DHA) induces caspase-3-dependent apoptosis in human lung adenocarcinoma ASTC-a-1 cells. J Biomed Sci 16:16. doi: 10.1186/1423-0127-16-16
    [20] Bratton SB, MacFarlane M, Cain K, et al. (2000) Protein complexes activate distinct caspase cascades in death receptor and stress-induced apoptosis. Exp Cell Res 256: 27-33. doi: 10.1006/excr.2000.4835
    [21] Cho SG, Choi EJ (2002) Apoptotic signaling pathways: caspases and stress-activated protein kinases. J Biochem Mol Biol 35: 24-27. doi: 10.5483/BMBRep.2002.35.1.024
    [22] Enari M, Sakahira H, Yokoyama H, et al. (1998) A caspase-activated DNase that degrades DNA during apoptosis, and its inhibitor ICAD. Nature 391: 43-50. doi: 10.1038/34112
    [23] Korn C, Scholz SR, Gimadutdinow O, et al. (2002) Involvement of conserved histidine, lysine and tyrosine residues in the mechanism of DNA cleavage by the caspase-3 activated DNase CAD. Nucleic Acids Res 30: 1325-1332. doi: 10.1093/nar/30.6.1325
  • This article has been cited by:

    1. Cengiz Kahraman, Proportional picture fuzzy sets and their AHP extension: Application to waste disposal site selection, 2024, 238, 09574174, 122354, 10.1016/j.eswa.2023.122354
    2. Omar Barukab, Asghar Khan, Sher Afzal Khan, Fermatean fuzzy Linguistic term set based on linguistic scale function with Dombi aggregation operator and their application to multi criteria group decision –making problem, 2024, 10, 24058440, e36563, 10.1016/j.heliyon.2024.e36563
    3. Yushuo Cao, Xuzhong Wu, Ling Ding, Weizhong Wang, A Hybrid Framework for Selecting Food Waste Treatment Techniques Using q-Rung Orthopair Fuzzy CRITIC-Generalized TODIM Method, 2024, 26, 1562-2479, 1916, 10.1007/s40815-024-01714-2
    4. Ertugrul Ayyildiz, Melike Erdogan, Literature analysis of the location selection studies related to the waste facilities within MCDM approaches, 2024, 1614-7499, 10.1007/s11356-024-34370-y
    5. Mery Ellen Brandt de Oliveira, Francisco Rodrigues Lima Junior, José Marcelo Almeida Prado Cestari , Hesitant Fuzzy Vikor e suas extensões: uma revisão sistemática de literatura, 2024, 15, 2178-9010, e4287, 10.7769/gesec.v15i10.4287
    6. Zihang Jia, Junsheng Qiao, Minghao Chen, Multi-level discrimination index for intuitionistic fuzzy coverings and its applications in feature selection, 2025, 263, 09574174, 125735, 10.1016/j.eswa.2024.125735
  • Reader Comments
  • © 2014 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5071) PDF downloads(838) Cited by(1)

Figures and Tables

Figures(8)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog