Open AccessArticle

Stability Analysis of a Mathematical Model for Adolescent Idiopathic Scoliosis from the Perspective of Physical and Health Integration

Yuhua Zhang

¹ and

Haiyin Li

^2,*

Institute of Physical Education, Henan University of Economics and Law, Zhengzhou 450046, China

Department of Mathematics and Information, Henan University of Economics and Law, Zhengzhou 450046, China

Author to whom correspondence should be addressed.

Symmetry 2023, 15(8), 1609; https://doi.org/10.3390/sym15081609

Submission received: 6 July 2023 / Revised: 17 August 2023 / Accepted: 18 August 2023 / Published: 20 August 2023

(This article belongs to the Special Issue Mathematical Modeling in Biology and Life Sciences)

Download

Browse Figures

Versions Notes

Abstract

In this paper, we take physical and health integration as the entry point. Firstly, based on the transformation mechanism of adolescent idiopathic scoliosis we construct a time delay differential model. Moreover, using the theory of characteristic equation we discuss the stability of a positive equilibrium under the delays of

τ = 0

and

τ \neq 0

. Furthermore, through numerical simulation, it has been verified the delay,

τ

, exceeds a critical value, the positive equilibrium loses its stability and Hopf bifurcation occurs. Lastly, we determine that sports have a positive effect on adolescent idiopathic scoliosis, directly reducing the number of people with adolescent idiopathic scoliosis.

Keywords:

adolescent idiopathic scoliosis; physical and health integration; time delay; stability; Hopf bifurcation

1. Introduction

In the report of the 20th National Congress, President Xi Jinping pointed out the need to “Extensively carry out national fitness activities, strengthen youth sport work, promote the development of mass sport, and accelerate the construction of a sporting country”. Adolescent health has drawn top-level national attention. In 2019, an excellent rate for high-quality health standards among students aged 6–22 in China of 23.8% was achieved, and the areas with a highly excellent rate were the economically developed eastern and coastal areas [1]. The rate of abnormal spinal curvatures among middle school students is 2.8%, and this is mainly due to abnormal postural curvature [2]. In 2022, China Central Television News Channel’s Chaowen Tianxia program reported that adolescent idiopathic scoliosis is the third-largest disease among adolescents after obesity and myopia. The number of primary and secondary school students with scoliosis in China is expected to exceed 5 million, with an annual increase of 300,000 [3]. Therefore, it is urgent to enhance the physique of adolescents and improve their physical health. How to rehabilitate adolescent idiopathic spinal curvatures has become a new topic for scholars.

Adolescent idiopathic scoliosis is a form of structural scoliosis of the spine, and accounts for over 80% of all scoliosis populations [4]. Adolescent idiopathic scoliosis is a common orthopedic disease, and it affects the alignment, growth and function of the spine, and the physical and mental health of adolescents especially during the growth and development period. If there is not a timely intervention made for adolescent idiopathic scoliosis, with the increase in lateral deformity, it will be accompanied by abnormal posture, lower back pain, decreased cardiovascular function, decreased respiratory function, decreased athletic ability and other functional problems, and then create conditions such as depression, anxiety and other psychological problems [5]. So far, the specific etiology of adolescent idiopathic scoliosis is unclear, and scholars from various academic circles, such as genetics, hormones, growth abnormalities, biomechanics, and neuromuscular theories [6], have conducted extensive research on it and formed many theories but it is recognized somewhat that long durations of abnormal posture in adolescents are a cause. Moreover, Chen [7] believes that high-intensity studying and a lack of sleep are the causes of rising incidence and increased severity of scoliosis among children and adolescents in China, while a lack of participation in sports is the direct cause of an increasing number of scoliosis cases.

The International Society for Scoliosis Research used Cobb’s method to measure the angle of the lateral curvature of the spine in X-ray films. X-ray films are taken with the patient in an upright position. Cobb angles exceeding 10° define scoliosis. For the angle of mild adolescent idiopathic scoliosis, there is no exact value. Generally speaking, when the Cobb angle is between 10° and 25°, it indicates mild adolescent idiopathic scoliosis; angles of 25° to 45° indicate to moderate adolescent idiopathic scoliosis; an angle greater than 45° indicates severe adolescent idiopathic scoliosis. In medicine, individuals without scoliosis and those whose Cobb angles are less than 10° are defined as healthy individuals, while moderate adolescent idiopathic scoliosis usually requires the use of corrective braces for treatment, and severe adolescent idiopathic scoliosis requires surgery.

Most of the literature on adolescent idiopathic scoliosis is from the perspective of pathogenesis and biomechanics [8,9]. The theory of differential equation plays a very important role in not only practical applications, but also theoretical research such as that of Slyn’ko and Tunç [10] and Tunç and Tunç [11]; it is commonly used to study infectious diseases such as COVID-19 and HIV [12,13,14], as well as skin cancer, leukemia and other non-communicable diseases [15,16]. However, there are few researchers who have studied adolescent idiopathic scoliosis by constructing differential models and using stability theory. Physical and health integration is the in-depth combination of two major systems for sports and health, and is the interaction and sharing of resources such as technology, talents, and facilities in the two systems. Compared to physical and medical integration, it has the characteristics of extensive coverage, a diversity of actions and effective results. It is one of the important ways to realize the healthy China strategy. Furthermore, starting from physical and health integration, we investigate the trend of adolescent idiopathic scoliosis and the influence of sports on adolescent idiopathic scoliosis from a mathematical perspective through theoretical analysis and numerical simulation, and provide a new research method for studying adolescent idiopathic scoliosis. In this paper, we refer individuals without scoliosis and those with a Cobb angle less than of 10 degrees as susceptible individuals, those with a Cobb angle between 10 degrees and 45 degrees as affected individuals, and those who require surgical treatment as surgical individuals. Therefore, adolescents are classified into three groups: susceptible individuals, infected individuals and surgical individuals.

This paper is organized as follows: in Section 2, we construct the model of adolescent idiopathic scoliosis; in Section 3, we investigate stability including stability switches of the system (2); in Section 4, applications with numerical simulations are described to illustrate our results; in Section 5, we discuss the effect of sports on adolescent idiopathic scoliosis; in Section 6, we conclude the paper.

2. Formulation of Mathematical Model

Adolescent idiopathic scoliosis, a disease with a persistent deviation of the spine from the midline of the body, making the spine curved or “S”-shaped to the side, is not an infectious disease. Let susceptible individuals be denoted by

h (t)

and affected individuals be denoted by

c (t)

. For adolescents, from childhood to adolescence,

h (t)

increases logistically. Under acquired or congenital influence, some adolescents have scoliosis, and the transfer rate from

h (t)

c (t)

μ

. Adolescent idiopathic scoliosis is a reversible disease; most

c (t)

can recover with the help of a brace and sports. The recovery rate with a brace is

β^{'}

, and that with sports is

α

. However, a small number of patients cannot use a brace effectively and do not participate in appropriate sports, so they require surgery after the aggravation of the condition. It takes about 3 months to fully recover after surgery. After recovery, teenagers are aware of the importance of daily sitting posture and sports, both psychologically and physically, which will make them change bad habits and strengthen sports performance, and make them no longer suffer from scoliosis in the adolescent stage. Therefore, the surgical individual is not considered in the model’s construction. This means that the conversion process between susceptible individuals and affected individuals is represented by the flowchart provided in Figure 1.

Furthermore, time delay refers to the phenomenon of delay, and reflects the dependence of the system on both the current and the past state. In the models for many biological processes, time delay is usually introduced to more accurately describe the change in objective things [17,18,19,20]. When we construct a model of adolescent idiopathic scoliosis, in order to make our research more practical, we need to consider the impact of time delay. During the onset of adolescent idiopathic scoliosis, the transition period from individuals being susceptible to affected is above 3–6 months, and the recovery period for individuals from being affected to susceptible is above 3–6 months under the combined action of sports and brace use (except for those who need surgical treatment, from affected individuals to surgical individuals). Therefore, we use a time delay,

τ

, to represent the transition period and the recovery period in this paper.

Based on the above analysis and the transformation mechanism of adolescent idiopathic scoliosis, we obtain the following system,

{\begin{matrix} h^{'} (t) = r h (t) (1 - \frac{h (t)}{K}) + (α + β^{'}) c (t - τ) - μ h (t - τ), \\ c^{'} (t) = μ h (t - τ) - (α + β^{'}) c (t - τ) - δ c (t), \end{matrix}

(1)

where

μ

is the transmission rate,

α

is the recovery rate from sports,

β^{'}

is the recovery rate from brace use, and

δ

is the surgery rate.

h (t)

is assumed to grow logistically with the growth rate,

r

, and carrying capacity,

K

τ

is the time delay, and indicates the transmission period and recovery period. All parameters in

μ, α, β^{'}, δ, r, K

are positive.

The initial conditions for the delayed system above take the form of

h_{0} (θ) = ϕ_{1} (θ) \geq 0, c_{0} (θ) = ϕ_{2} (θ) \geq 0, θ \in [- τ, 0], h (0) > 0, c (0) > 0,

where

ϕ = (ϕ_{1}, ϕ_{2}) \in C ([- τ, 0], R_{+}^{2}), R_{+}^{2} = {(h, c) : h \geq 0, c \geq 0}

‖ ϕ ‖ = \max {| ϕ (θ) | : θ \in

[- τ, 0]}

, and

| ϕ |

are any norm in

R_{+}^{2}

. As usual, we use the conventional notation

h_{t} (θ) = h (t + θ)

for

θ \in [- τ, 0]

. In this paper, for the biological reasons, we only consider that and

h (θ) > 0, c (θ) > 0

are all continuously differentiable in

- τ \leq θ \leq 0

and

h (0) > 0, c (0) > 0

System (1) also can be rewritten as

{\begin{matrix} h^{'} (t) = r h (t) (1 - \frac{h (t)}{K}) + β c (t - τ) - μ h (t - τ), \\ c^{'} (t) = μ h (t - τ) - β c (t - τ) - δ c (t), \end{matrix}

(2)

where

β = α + β^{'}

is the recovery rate aided by the combined action of sports and brace use.

The steady states occur by setting the left-hand side of system (2) to zero, i.e.,

{\begin{matrix} r h (t) (1 - \frac{h (t)}{K}) + β c (t) - μ h (t) = 0, \\ μ h (t) - (β + δ) c (t) = 0 . \end{matrix}

so we can obtain the following conclusion.

Lemma 1.

We assume that the functions

c (t)

and

h (t)

are continuously differentiable. When the condition

\frac{μ δ}{r (β + δ)} < 1 (H_{0})

holds, system (2) has a unique positive equilibrium,

E^{*} (h^{*}, c^{*})

, where

h^{*} = K (1 - \frac{μ δ}{r (β + δ)}), c^{*} = \frac{μ h^{*}}{β + δ} .

Obviously, we can obtain

h^{*} < K

from the expression

h^{*}

3. Stability Analysis

In this section, we mainly provide the stability of the positive equilibrium and the stability switch of the system (2). Firstly, the characteristic equation of the linearized system (2) is given, and then the stability of the positive equilibrium and stability switch of the system (2) are determined via the sign of the characteristic root.

Let

h (t) = h^{*} + H (t)

and

c (t) = c^{*} + C (t)

. Then, the linearized system (2) is

{\begin{matrix} H^{'} (t) = h^{*} [(r - \frac{2 r h^{*}}{K}) H (t) - μ H (t - τ) + β C (t - τ)], \\ C^{'} (t) = c^{*} [μ H (t - τ) - β C (t - τ) - δ C (t)] . \end{matrix}

(3)

Then, the corresponding Jacobian matrix of system (3) is

(\begin{matrix} r - \frac{2 r h^{*}}{K} - μ e^{- λ τ} & β e^{- λ τ} \\ μ e^{- λ τ} & - β e^{- λ τ} - δ \end{matrix}),

so the characteristic equation of system (3) is

| \begin{matrix} λ - r + \frac{2 r h^{*}}{K} + μ e^{- λ τ} & - β e^{- λ τ} \\ - μ e^{- λ τ} & λ + β e^{- λ τ} + δ \end{matrix} | = 0,

that is,

λ^{2} + (β + μ) λ e^{- λ τ} + [μ δ + β (\frac{2 r h^{*}}{K} - r)] e^{- λ τ} + (δ + \frac{2 r h^{*}}{K} - r) λ + δ (\frac{2 r h^{*}}{K} - r) = 0 .

(4)

Let

a_{1} = β + μ > 0, a_{2} = μ δ + β (\frac{2 r h^{*}}{K} - r), a_{3} = δ + \frac{2 r h^{*}}{K} - r, a_{4} = δ (\frac{2 r h^{*}}{K} - r)

. Then, Equation (4) becomes

λ^{2} + a_{1} λ e^{- λ τ} + a_{2} e^{- λ τ} + a_{3} λ + a_{4} = 0 .

(5)

When

τ = 0

, Equation (4) is

λ^{2} + (a_{1} + a_{3}) λ + a_{2} + a_{4} = 0 .

(6)

Therefore, if

a_{1} + a_{3} > 0, a_{2} + a_{4} > 0

(7)

holds, all roots of Equation (6) have negative real parts.

Note 1.

h^{*} > \frac{K}{2}

, then

a_{2} > 0, a_{3} > 0, a_{4} > 0

, and condition (7) holds. Meanwhile, the condition

h^{*} > \frac{K}{2}

is equivalent to

\frac{μ δ}{r (β + δ)} < \frac{1}{2} . (H_{1})

Obviously, when

(H_{1})

holds,

(H_{0})

holds. That is, there is a positive equilibrium

E^{*}

Based on the above analysis and Routh Hurwitz criterion [21], we have the following conclusion.

Theorem 1.

If conditions

(H_{0})

and (7) hold, then the positive equilibrium,

E^{*}

, of system (2) is asymptotically stable for

τ = 0

Note 2.

At the same time, when

(H_{1})

holds, then

(H_{0})

and (7) both hold.

Next, we want to confirm whether or not the real part of the root,

λ

, for the characteristic equation, Equation (5), will increase to zero and eventually become positive as

τ

changes. When

τ \neq 0

, let

λ = i ω

(

ω \neq 0

, since

a_{2} + a_{4} > 0

); through Equation (5) we can obtain

{\begin{array}{l} a_{1} ω \sin ω τ + a_{2} \cos ω τ = ω^{2} - a_{4}, \\ a_{1} ω \cos ω τ - a_{2} \sin ω τ = - a_{3} ω . \end{array}

(8)

Thus, we have

{(a_{1} ω)}^{2} + a_{2}^{2} = {(ω^{2} - a_{4})}^{2} + {(a_{3} ω)}^{2}

, that is

ω^{4} + (a_{3}^{2} - a_{1}^{2} - 2 a_{4}) ω^{2} + a_{4}^{2} - a_{2}^{2} = 0,

(9)

and the root of Equation (9) is

ω_{\pm}^{2} = \frac{(a_{1}^{2} - a_{3}^{2} + 2 a_{4}) \pm \sqrt{{(a_{1}^{2} - a_{3}^{2} + 2 a_{4})}^{2} - 4 (a_{4}^{2} - a_{2}^{2})}}{2} .

(10)

Therefore, when the condition

a_{1}^{2} - a_{3}^{2} + 2 a_{4} < 0 a n d a_{4}^{2} - a_{2}^{2} > 0 o r {(a_{1}^{2} - a_{3}^{2} + 2 a_{4})}^{2} < 4 (a_{4}^{2} - a_{2}^{2})

(11)

holds,

ω_{\pm}^{2}

values are all not positive, that is, (9) has no positive root, which means (5) has no pure imaginary root. Since Equation (6) has negative real parts, according to Rouche’s theorem, the characteristic equation, Equation (5), also has negative real parts for the roots. We can draw the following conclusion.

Theorem 2.

If conditions (7) and (11) hold, all roots of the characteristic equation, Equation (5), have negative real parts for the arbitrary

τ > 0 .

Otherwise, when

a_{4}^{2} - a_{2}^{2} < 0 o r {\begin{matrix} a_{1}^{2} - a_{3}^{2} + 2 a_{4} > 0 a n d \\ {(a_{1}^{2} - a_{3}^{2} + 2 a_{4})}^{2} = 4 (a_{4}^{2} - a_{2}^{2}) \end{matrix}

(12)

holds, Equation (9) has a positive root,

ω_{+}^{2} .

a_{1}^{2} - a_{3}^{2} + 2 a_{4} > 0, a_{4}^{2} - a_{2}^{2} > 0 a n d {(a_{1}^{2} - a_{3}^{2} + 2 a_{4})}^{2} > 4 (a_{4}^{2} - a_{2}^{2})

(13)

holds, Equation (9) has two positive roots

ω_{\pm}^{2} .

If the conditions (12) and (13) hold, Equation (5) has pure imaginary root when

τ

reaches a certain value. The critical value

τ_{j}^{\pm}

τ

determined by (8) is

τ_{j}^{\pm} = \frac{1}{ω_{\pm}} \arccos {\frac{(a_{2} - a_{1} a_{3}) ω_{\pm}^{2} - a_{2} a_{4}}{a_{1}^{2} ω_{\pm}^{2} + a_{2}^{2}}} + \frac{2 j π}{ω_{\pm}}, j = 0, 1, 2, \dots .

(14)

From the above analysis, we can draw the following conclusion.

Lemma 2.

(1) When (7) and (12) hold and

τ = τ_{j}^{+}

, (5) has a pair of pure roots,

\pm i ω_{+} .

(2) When (7) and (13) hold and

τ = τ_{j}^{+}

(

τ = τ_{j}^{-}

), (5) has a pair of pure roots,

\pm i ω_{+}

(

\pm i ω_{-}

When

τ > τ_{j}^{+}

and

τ < τ_{j}^{+}

, we would expect that the real parts of certain characteristic roots of (5) become positive. For convenience, we denote

λ_{j}^{\pm} = α_{j}^{\pm} (τ) +

i ω_{j}^{\pm} (τ),

j = 0, 1, 2, \dots

; meanwhile, the roots of (5) satisfy

α_{j}^{\pm} (τ_{j}^{\pm}) = 0

and

ω_{j}^{\pm} (τ_{j}^{\pm}) = ω_{\pm} .

We can verify that the transversality conditions

\frac{d}{d τ} Re λ_{j}^{+} (τ_{j}^{+}) > 0

and

\frac{d}{d τ} Re λ_{j}^{-} (τ_{j}^{-}) < 0

hold. It follows that

τ_{j}^{\pm}

represents bifurcation values. Therefore, we have the following theorem on the distribution of eigenvalues.

Theorem 3.

Let

τ_{j}^{\pm}

be as defined in (14).

(1): If (7) and (12) hold, all roots of Equation (5) have negative real parts for $τ \geq 0$ .
(2): If (7) and (12) hold, all roots of Equation (5) have negative real parts for $τ \in [0, τ_{0}^{+})$ ; Equation (5) has a pair of pure roots, $\pm i ω_{+}$ , for $τ = τ_{0}^{+}$ ; Equation (5) has at least one root with a positive real part for $τ > τ_{0}^{+}$ .
(3): If (7) and (13) hold, then there exists a positive integer, $k$ , such that there are $k$ switches from stability to unstability to stability. That is, when

$τ \in [0, τ_{0}^{+}], (τ_{0}^{-}, τ_{1}^{+}), \dots, (τ_{k - 1}^{-}, τ_{k}^{+}),$

all roots of Equation (5) have negative real parts; when

$τ \in [τ_{0}^{+}, τ_{0}^{-}), [τ_{1}^{+}, τ_{1}^{-}), \dots, [τ_{k - 1}^{+}, τ_{k - 1}^{-}) a n d τ > τ_{k}^{+},$

Equation (5) has at least one root with a positive real part.

Note 3.

The case (3) of Theorem 3 shows that when the time delay,

τ

, exceeds the critical value,

τ_{j}^{+},

j = 0, 1, 2, \dots, k - 1

, the positive equilibrium,

E^{*}

, of system (2) loses stability, and system (2) exists as a Hopf bifurcation [22].

4. Numerical Simulation

In this part, we use numerical simulation to show the interesting dynamical behavior of system (2). From the following examples, it can be seen that with the change in the time delay,

τ

, the stability of the positive equilibrium of system (2) will also change, and there is a Hopf bifurcation, which is also a numerical verification of the conclusion in the second part.

Application 1.

Let

r = 0.02, K = 10^{8}, β = 0.7, δ = 0.3, μ = 0.02

; then, system (2) becomes

{\begin{matrix} h^{'} (t) = 0.02 h (t) (1 - 10^{- 8} h (t)) + 0.7 c (t) - 0.02 h (t), \\ c^{'} (t) = 0.02 h (t) - 0.7 c (t) - 0.3 c (t) . \end{matrix}

(15)

Via simple numerical computation, the positive equilibrium of system (15) is

E^{*} = (7 \times 10^{7}, 1.4 \times 10^{6})

, and

a_{1} = 0.7200, a_{2} = 0.0116, a_{3} = 0.3080,

a_{4} = 0.0024,

\frac{μ δ}{r (β + δ)} = 0.3 < \frac{1}{2} .

Therefore, the condition

(H_{1})

holds; that is,

(H_{0})

and (7) are satisfied, and then Theorem 1 holds.

From Figure 2, it can be seen that

(h (t), c (t))

tends towards

(h^{*}, c^{*}) = (7 \times 10^{7}, 1.4 \times 10^{6})

t \to + \infty

for

τ = 0

. That is, the positive equilibrium,

E^{*}

, of system (2) is asymptotically stable for

τ = 0

. Meanwhile

a_{1}^{2} - a_{3}^{2} + 2 a_{4} = 0.4283, a_{4}^{2} - a_{2}^{2} = - 1.288 \times 10^{- 4},

and

{(a_{1}^{2} - a_{3}^{2} + 2 a_{4})}^{2} - 4 (a_{4}^{2} - a_{2}^{2}) = 0.1840,

so condition (12) holds. Via (10) and (14), we obtain

τ_{0}^{+} = 3.0366 .

Application 2.

Based on Application 1, when

τ = 3

, the corresponding system is

{\begin{array}{l} h^{'} (t) = 0.02 h (t) (1 - 10^{- 8} h (t)) + 0.7 c (t - 3) - 0.02 h (t - 3), \\ c^{'} (t) = 0.02 h (t - 3) - 0.7 c (t - 3) - 0.3 c (t) . \end{array}

(16)

where

τ = 3 < τ_{0}^{+}

, the curves of

t

and

h (t)

t

and

c (t)

are shown in Figure 3.

From Figure 3, it can be seen that

(h (t), c (t))

still tends towards

(h^{*}, c^{*}) = (7 \times 10^{7}, 1.4 \times 10^{6})

t \to + \infty

for

τ = 3

, indicating that the positive equilibrium,

E^{*}

, of system (2) is still asymptotically stable.

Application 3.

Based on Application 1, when

τ = 3.1

, the corresponding system is

{\begin{array}{l} h^{'} (t) = 0.02 h (t) (1 - 10^{- 8} h (t)) + 0.7 c (t - 3.1) - 0.02 h (t - 3.1), \\ c^{'} (t) = 0.02 h (t - 3.1) - 0.7 c (t - 3.1) - 0.3 c (t) . \end{array}

(17)

When

τ = 3.1 > τ_{0}^{+}

the curves of

t

and

h (t)

, and

t

and

c (t)

are shown in Figure 4.

Moreover, the orbits of

h (t)

and

c (t)

under the same conditions are in Figure 5. Further, we extend

t

2000

and obtain the orbits of

h (t)

and

c (t)

, as shown in in Figure 6.

From Applications 1–3, it can be seen that the positive equilibrium,

E^{*}

, is stable for

τ \in [0, τ_{0}^{+}]

, and the positive equilibrium,

E^{*}

, is unstable for

τ > τ_{0}^{+}

, which verifies the conclusion of Theorem 3(3). Therefore,

τ_{0}^{+}

is a critical value, and Hopf bifurcation occurs at

τ = τ_{0}^{+}

5. Discussion

The bifurcation problem is an important problem in the study of dynamical systems and nonlinear differential equations, and its research objects are structurally unstable systems. The bifurcation phenomenon refers to the sudden changes in certain properties of the system (such as the equilibrium state, periodic phenomena, stability, etc.) when parameters change and exceed certain critical values in a parameter-dependent system. In this paper, time delay,

τ

, is not only the transition period from individuals being susceptible to affected but also the recovery period from individuals being affected to susceptible under the combined action of sports and brace use. In Applications 1–3, when

τ \in [0, τ_{0}^{+}]

the positive equilibrium,

E^{*}

, is stable, which means that the number of susceptible individuals

h (t) \to 7 \times 10^{7}

t \to + \infty

, and the number of affected individuals

c (t) \to 1.4 \times 10^{6}

t \to + \infty

. When

τ > τ_{0}^{+}

, the positive equilibrium,

E^{*}

, is unstable, which means that

h (t)

and

c (t)

do not tend toward

(7 \times 10^{7}, 1.4 \times 10^{6})

for a large

t > 0

. The transmission period and recovery period are 3–6 months, which is basically greater than 3.0366, and

E^{*}

is unstable, which shows sports have a positive effect on adolescent idiopathic scoliosis. Therefore, with the help of sports and brace use, the number of affected individuals is not fixed at a large value, but rather far away from this large fixed value.

In system (1), the recovery rate is denoted by

β = α + β^{'}

, where

α

is the recovery rate of movement,

β^{'}

is the recovery rate from brace use; since the recovery rate from brace use is relatively fixed, we mainly refer to the recovery rate from sports when we discuss recovery rate in this article. Next, in order to observe the impact of

β

on stability, the parameters

r, K, β, δ, μ

in Applications 1–3 remain unchanged, but the value of

β

is changed.

We consider the case of

β < 0.7

. Assuming

β = 0.6

, we determine that

τ_{0}^{+} = 3.7916

, so

E^{*} (6.6667 \times 10^{7}, 1.4815 \times 10^{6})

is stable for

τ \in [0, τ_{0}^{+}]

, and

E^{*} (6.6667 \times 10^{7}, 1.4815 \times 10^{6})

is unstable for

τ > τ_{0}^{+}

; with

β = 0.5

τ_{0}^{+} = 5.0772

E^{*} (6.2500 \times 10^{7}, 1.5625 \times 10^{6})

is stable for

τ \in [0, 5.0772]

, and

E^{*} (6.2500 \times 10^{7}, 1.5625 \times 10^{6})

is unstable for

τ > 5.0772

. It is easy to see that as the parameter

β

decreases and

τ_{0}^{+}

becomes larger, the chance of

τ > τ_{0}^{+}

decreases.

However, when

β > 0.7

, let

β = 0.8

and

τ_{0}^{+} = 2.5364

;

E^{*} (7.2727 \times 10^{7}, 1.3223 \times 10^{6})

is stable for

τ \in [0, 2.5364]

, and

E^{*} (7.2727 \times 10^{7}, 1.3223 \times 10^{6})

is unstable for

τ > 2.5364

; let

β = 0.85

and

τ_{0}^{+} = 2.3442

;

E^{*} (7.3913 \times 10^{7}, 1.2854 \times 10^{6})

is stable for

τ \in [0, 2.3442]

, and

E^{*} (7.3913 \times 10^{7}, 1.2854 \times 10^{6})

is unstable for

τ > 2.3442

. It is easy to see that as the parameter

β

increases and

τ_{0}^{+}

becomes smaller, the chance of

τ > τ_{0}^{+}

decreases. That is, participating in an adolescent idiopathic scoliosis intervention process through sports can improve the overall condition of adolescent idiopathic scoliosis.

Moreover, it also can be found that the positive equilibrium

E^{*} (h^{*}, c^{*})

of system (2) changes with the continuous increase in parameter

β

during the change process of parameter

β

above; the number of susceptible individuals,

h (t)

, continuously increases from

6.2500 \times 10^{7} (β = 0.5)

7.3913 \times 10^{7} (β = 0.85)

; the number of infected individuals

c (t)

continuously decreases from

1.5625 \times 10^{6} (β = 0.5)

1.2854 \times 10^{6} (β = 0.85)

. This shows that participation in sports directly affects the prevalence of adolescent idiopathic scoliosis. However, participation in sports here requires scientific sports prescriptions as a medium, such as corrective exercises, asymmetrical crawling and breathing training and other methods. In the specific intervention implementation process, it is also necessary to combine scientific training with the patient’s bending degree, the comfort of the support, and the patient’s physical fitness.

6. Conclusions

In this paper, we have investigated the stability of an adolescent idiopathic scoliosis system from the perspective of physical and health integration. We started with the generation mechanism of adolescent idiopathic scoliosis, and obtained a differential system for the transformation of adolescent idiopathic scoliosis. Afterwards, we qualitatively analyzed the stability of the positive equilibrium and the stability switches of the system with respect to the delay parameter. Moreover, it was also interesting to investigate the stability switches and Hopf bifurcation with respect to the delay parameter via numerical simulations. Lastly, we discussed the effect of sports on adolescent idiopathic scoliosis.

Author Contributions

Y.Z. conceived the idea, designed the numerical simulation, and wrote the paper; H.L. provided a stability analysis and conducted article correction. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Science and Technology Project of Henan Province, grant number 232102320306; Science and Technology Project of Henan Province, grant number 222102320192; Key Research Projects of Higher Education Institutions in Henan Province, grant number 23A110005. The APC was funded by 23A110005.

Data Availability Statement

Date will be available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

People’s Information. The 8th National Student Physical Fitness and Health Survey Results Announce That Development Indicators Such as Height and Weight of Students Continue to Improve. Available online: https://baijiahao.baidu.com/s?id=1709867909792215991&wfr=spider&for=pc (accessed on 3 September 2021).
China Daily CCTV News. The Detection Rate of Abnormal Spinal Curvature among Primary and Secondary School Students in China is 2.8%. Available online: http://t.m.china.com.cn/convert/c_zqCLyD6z.html (accessed on 14 July 2021).
CCTV-13 News Channel News Live Room. Pay Attention to Adolescent Idiopathic Scoliosis—The Number of Primary and Secondary School Students with Scoliosis in China Is Expected to Exceed 5 Million. Available online: https://tv.cctv.com/2022/06/27/VIDEYBWgcLdi5Sn88bBVBxCB220627.shtml (accessed on 27 June 2022).
Hefti, F. Pathogenesis and biomechanics of adolescent idiopathic scoliosis(AIS). J. Child. Orthop. 2013, 7, 17–24. [Google Scholar] [CrossRef] [PubMed]
Machida, M.; Weinstein, S.L.; Dubousset, J. Pathogenesis of Idiopathic Scoliosis; Springer Japan K.K.: Tokyo, Japan, 2018. [Google Scholar]
Anggriani, N.; Ndii, M.Z.; Amelia, R. etc. A mathematical COVID-19 model considering asymptomatic and symptomatic classes with waning immunity. Alex. Eng. J. 2022, 61, 113–124. [Google Scholar] [CrossRef]
Batistela, C.M.; Correa, D.P.F.; Bueno, A.M. etc. SIRSi compartmental model for COVID-19 pandemic with immunity loss. Chaos Soliton. Frac. 2021, 142, 110388. [Google Scholar] [CrossRef] [PubMed]
Yang, Y.; Huang, G.; Dong, Y.P. Stability and Hopf bifurcation of an HIV infection model with two time delays. Math. Biosci. Eng. 2023, 20, 1938–1959. [Google Scholar] [CrossRef] [PubMed]
Parvin, T.; Biswas, M.H.A.; Datta, B.K. Mathematical analysis of the transmission dynamics of skin cancer caused by UV radiation. J. Appl. Math. 2022, 2022, 5445281. [Google Scholar] [CrossRef]
Vitalii, S.; Cemil, T. Stability of abstract linear switched impulsive differential equations. Automatica 2019, 107, 433–441. [Google Scholar]
Cemil, T.; Osman, T. On the stability, integrability and boundedness analyses of systems of integro-differential equations with time-delay retardation. RACSAM 2021, 115, 115. [Google Scholar]
Besse, A.; Clapp, G.D.; Bernard, S. Stability analysis of a model of interaction between the immune system and cancer cells in chronic myelogenous leukemia. Bull. Math. Biol. 2018, 80, 1084–1110. [Google Scholar] [CrossRef] [PubMed]
Negrini, S.; Donzelli, S.; Aulisa, A.G.; Czaprowski, D.; Schreiber, S.; de Mauroy, J.C.; Diers, H.; Grivas, T.B.; Knott, P.; Kotwicki, T.; et al. 2016 SOSORT guidelines:orthopaedic and rehabilitation treatment of idiopathic scoliosis during growth. Scoliosis 2018, 13, 1–48. [Google Scholar]
Chinese Preventive Medicine Association; Spinal Disease Prevention and Control Committee; Chinese Preventive Medicine Association, Sports Exercise and Population. Implementation pathway of exercise therapy for adolescent idiopathic scoliosis guided by sports and medicine integration: Based on expert consensus. Chi. J. Prev. Med. 2023, 24, 7–14. [Google Scholar]
Peter, O.N. Adolescent Idiopathic Scoliosis; People’s Health Publishing House: Beijing, China, 2006. [Google Scholar]
Chen, Y.; Li, L.; Yang, H.; Hu, W.; Jia, F.; Zhai, F. Current status and influencing factors of scoliosis of children in Shijiazhuang. Chin. J. Sch. Health 2021, 42, 1674–1678. [Google Scholar]
Huo, H.F.; Li, W.T.; Nieto, J.J. Periodic solutions delayed predator–prey model with the Beddington–DeAngelis functional response. Chaos Soliton. Frac. 2007, 33, 505–512. [Google Scholar] [CrossRef]
Liu, S.; Beretta, E.; Breda, D. Predator–prey model of Beddington–DeAngelis type with maturation and gestation delays. Nonlinear Anal. Real World Appl. 2010, 11, 4072–4091. [Google Scholar] [CrossRef]
Sun, X.K.; Huo, H.F.; Xiang, H. Bifurcations and stability analysis in predator–prey model with a stage-structure for predator. Nonlinear Dynam. 2009, 58, 497–513. [Google Scholar] [CrossRef]
Zeng, G.; Wang, F.; Nieto, J.J. Complexity and a delayed predator–prey model with impulsive harvest and Holling type II functional response. Adv. Complex Syst. 2008, 11, 77–97. [Google Scholar] [CrossRef]
Ma, Z.E.; Zhou, Y.C. Qualitative and Stability Methods for Ordinary Differential Equation; Science Press: Beijing, China, 2013; pp. 174–198. [Google Scholar]
Hassard, B.D.; Kazarinoff, N.D.; Wan, Y.H. Theory and Applications of Hopf Bifurcation; Cambridge University Press: Cambridge, UK, 1981. [Google Scholar]

Figure 1. Transformation process between susceptible individuals

h (t)

and affected individuals

c (t)

Figure 1. Transformation process between susceptible individuals

h (t)

and affected individuals

c (t)

Figure 2. The curves of system (15), where the blue and green curves represent

h (t)

and

c (t)

, respectively.

Figure 2. The curves of system (15), where the blue and green curves represent

h (t)

and

c (t)

, respectively.

Figure 3. The curves of

t

and

h (t)

, and

t

and

c (t)

in system (16) for

τ = 3

, where the blue and green curves respresent

h (t)

and

c (t)

, respectively.

Figure 3. The curves of

t

and

h (t)

, and

t

and

c (t)

in system (16) for

τ = 3

, where the blue and green curves respresent

h (t)

and

c (t)

, respectively.

Figure 4. The curves of

t

and

h (t)

, and

t

and

c (t)

in system (17) with

τ = 3.1

and

t = 1000

, where the blue and green curves respresent

h (t)

and

c (t)

, respectively.

Figure 4. The curves of

t

and

h (t)

, and

t

and

c (t)

in system (17) with

τ = 3.1

and

t = 1000

, where the blue and green curves respresent

h (t)

and

c (t)

, respectively.

Figure 5. Orbits of

h (t)

and

c (t)

for system (17) with

τ = 3.1

and

t = 1000

Figure 5. Orbits of

h (t)

and

c (t)

for system (17) with

τ = 3.1

and

t = 1000

Figure 6. Orbits of

h (t)

and

c (t)

for system (17) with

τ = 3.1

and

t = 2000

Figure 6. Orbits of

h (t)

and

c (t)

for system (17) with

τ = 3.1

and

t = 2000

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Zhang, Y.; Li, H. Stability Analysis of a Mathematical Model for Adolescent Idiopathic Scoliosis from the Perspective of Physical and Health Integration. Symmetry 2023, 15, 1609. https://doi.org/10.3390/sym15081609

AMA Style

Zhang Y, Li H. Stability Analysis of a Mathematical Model for Adolescent Idiopathic Scoliosis from the Perspective of Physical and Health Integration. Symmetry. 2023; 15(8):1609. https://doi.org/10.3390/sym15081609

Chicago/Turabian Style

Zhang, Yuhua, and Haiyin Li. 2023. "Stability Analysis of a Mathematical Model for Adolescent Idiopathic Scoliosis from the Perspective of Physical and Health Integration" Symmetry 15, no. 8: 1609. https://doi.org/10.3390/sym15081609

APA Style

Zhang, Y., & Li, H. (2023). Stability Analysis of a Mathematical Model for Adolescent Idiopathic Scoliosis from the Perspective of Physical and Health Integration. Symmetry, 15(8), 1609. https://doi.org/10.3390/sym15081609

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

Article Menu

Stability Analysis of a Mathematical Model for Adolescent Idiopathic Scoliosis from the Perspective of Physical and Health Integration

Abstract

1. Introduction

2. Formulation of Mathematical Model

3. Stability Analysis

4. Numerical Simulation

5. Discussion

6. Conclusions

Author Contributions

Funding

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI