Loading [MathJax]/jax/output/SVG/jax.js
Case report Special Issues

Multimodal physiotherapy treatment based on a biobehavioral approach in a patient with chronic low back pain: A case report

  • Low back pain and its chronification are among the most common causes of disabilities worldwide. This is why multiple interventions and treatment approaches have been investigated around this pathology in recent years. Currently, the evidence is increasingly leaning towards treating chronic low back pain with a treatment based on the patient's centered biopsychosocial model. This case report presents the evaluation and treatment of a 43-year-old woman with chronic low back pain using a multimodal physiotherapy based on a biobehavioral approach. In the physiotherapy intervention, the patient's pain characteristics and somatosensory, motor-functional, and affective-cognitive states were first evaluated. Subsequently, a multimodal treatment including therapeutic exercise, pain neuroscience education, and orthopedic manual and physical therapies was applied for a total of 14 sessions over a period of 9 weeks. Finally, a post-intervention evaluation and at the one-month follow-up were carried out, in which it was observed that the treatment had significantly improved the patient's symptomatology and their general state. This case suggests that a multimodal physiotherapy treatment based on a biobehavioral approach is an effective option to improve pain symptoms and somatosensory, motor-functional, and affective-cognitive aspects in the reported patient; therefore, a treatment of these characteristics may be an option for patients with chronic low back pain.

    Citation: Carlos Forner-Álvarez, Ferran Cuenca-Martínez, Rafael Moreno-Gómez-Toledano, Celia Vidal-Quevedo, Mónica Grande-Alonso. Multimodal physiotherapy treatment based on a biobehavioral approach in a patient with chronic low back pain: A case report[J]. AIMS Medical Science, 2024, 11(2): 77-89. doi: 10.3934/medsci.2024007

    Related Papers:

    [1] Kenneth H. Karlsen, Süleyman Ulusoy . On a hyperbolic Keller-Segel system with degenerate nonlinear fractional diffusion. Networks and Heterogeneous Media, 2016, 11(1): 181-201. doi: 10.3934/nhm.2016.11.181
    [2] José Antonio Carrillo, Yingping Peng, Aneta Wróblewska-Kamińska . Relative entropy method for the relaxation limit of hydrodynamic models. Networks and Heterogeneous Media, 2020, 15(3): 369-387. doi: 10.3934/nhm.2020023
    [3] Piotr Biler, Grzegorz Karch, Jacek Zienkiewicz . Morrey spaces norms and criteria for blowup in chemotaxis models. Networks and Heterogeneous Media, 2016, 11(2): 239-250. doi: 10.3934/nhm.2016.11.239
    [4] Panpan Xu, Yongbin Ge, Lin Zhang . High-order finite difference approximation of the Keller-Segel model with additional self- and cross-diffusion terms and a logistic source. Networks and Heterogeneous Media, 2023, 18(4): 1471-1492. doi: 10.3934/nhm.2023065
    [5] Steinar Evje, Kenneth H. Karlsen . Hyperbolic-elliptic models for well-reservoir flow. Networks and Heterogeneous Media, 2006, 1(4): 639-673. doi: 10.3934/nhm.2006.1.639
    [6] Raul Borsche, Axel Klar, T. N. Ha Pham . Nonlinear flux-limited models for chemotaxis on networks. Networks and Heterogeneous Media, 2017, 12(3): 381-401. doi: 10.3934/nhm.2017017
    [7] Avner Friedman . PDE problems arising in mathematical biology. Networks and Heterogeneous Media, 2012, 7(4): 691-703. doi: 10.3934/nhm.2012.7.691
    [8] Hirotada Honda . Global-in-time solution and stability of Kuramoto-Sakaguchi equation under non-local Coupling. Networks and Heterogeneous Media, 2017, 12(1): 25-57. doi: 10.3934/nhm.2017002
    [9] Caihong Gu, Yanbin Tang . Global solution to the Cauchy problem of fractional drift diffusion system with power-law nonlinearity. Networks and Heterogeneous Media, 2023, 18(1): 109-139. doi: 10.3934/nhm.2023005
    [10] Stefan Berres, Ricardo Ruiz-Baier, Hartmut Schwandt, Elmer M. Tory . An adaptive finite-volume method for a model of two-phase pedestrian flow. Networks and Heterogeneous Media, 2011, 6(3): 401-423. doi: 10.3934/nhm.2011.6.401
  • Low back pain and its chronification are among the most common causes of disabilities worldwide. This is why multiple interventions and treatment approaches have been investigated around this pathology in recent years. Currently, the evidence is increasingly leaning towards treating chronic low back pain with a treatment based on the patient's centered biopsychosocial model. This case report presents the evaluation and treatment of a 43-year-old woman with chronic low back pain using a multimodal physiotherapy based on a biobehavioral approach. In the physiotherapy intervention, the patient's pain characteristics and somatosensory, motor-functional, and affective-cognitive states were first evaluated. Subsequently, a multimodal treatment including therapeutic exercise, pain neuroscience education, and orthopedic manual and physical therapies was applied for a total of 14 sessions over a period of 9 weeks. Finally, a post-intervention evaluation and at the one-month follow-up were carried out, in which it was observed that the treatment had significantly improved the patient's symptomatology and their general state. This case suggests that a multimodal physiotherapy treatment based on a biobehavioral approach is an effective option to improve pain symptoms and somatosensory, motor-functional, and affective-cognitive aspects in the reported patient; therefore, a treatment of these characteristics may be an option for patients with chronic low back pain.



    Predator-prey model is one of basic interspecies relations for ecological and social system [1]. The more complex biochemical network structure and food chain are based on the predator-prey model [2]. The study of Lotka and Volterra [3,4] has opened the way to study the dynamics of the predator-prey systems. After that, Gause and Smaragdova also proposed a well-known Gause-type predator-prey model. Kolmogorov first focused on the qualitative analysis of this Gause-type predator-prey model in 1972. Freedman [5] introduced the generalized autonomous Gause model, which comes from accounting for periodic changes of the environment. Gause-type predator-prey models have been widely applied to describe some population models [6,7,8,9]. For example, Hasik [6] considered the generalized Gause-type predator-prey model

    {x=xg(x)yp(x),y=y[q(x)γ], (1.1)

    here g(x) represents the increase in prey density. When the natural environment is relatively bad, the mortality rate of the population is higher than its birth rate, so the g(x) here can get a negative value. p(x) represents the amount of prey consumed by a single predator per unit time. q(x)γ represents the growth rate of the predator, and the same as g(x), q(x)γ can also be taken to a negative value. Ding et al. [7] considered the periodic Gause-type predator-prey system with delay

    {x(t)=x(t)f(t,x(tτ(t)))g(t,x(t))y(tσ1(t)),y(t)=y(t)[d(t)+h(t,x(tσ2(t)))], (1.2)

    where x(0),y(0)0 are the prey and the predator and obtained the positive periodic solution of this system (1.2) by using the continuation theorem.

    For the past few years, more and more researchers are interested in the dynamic behavior of predator-prey systems with Allee effect. The Allee effect describes that the low population is affected by the positive relationship between population growth rate and density, which increases the likelihood of their extinction. Terry [10] considered predator-prey systems with Allee effect and described how to extend the traditional definition of effective components and population Allee effect for a single species model to predators in the predator-prey model. Cui et al. [11] focused on the dynamic behavior and steady-state solutions of a class of reaction-diffusion predator-prey systems with strong Allee effect. Cai et al. [12] explored the complex dynamic behavior for a Leslie-Gower predation model with additive Allee effect on prey. Without considering the influence of Allee effect on prey, the model has a unique global asymptotically stable equilibrium point. However, considering the influence of Allee reaction on prey, the model has no definite positive equilibrium point [12]. Baisad and Moonchai [13] were interested in a Gause-type predator-prey model with Holling type-Ⅲ functional response and Allee effect on prey as follows

    {dxdt=r(1xK)(xm)xsx2x2+a2y,dydt=(px2x2+a2c)y. (1.3)

    Using the linearization method, they gave the local stability of three equilibrium types and also carried out a numerical simulation of the results. Guan and Chen [14] studied the dynamical analysis of a two species amensalism model with Beddington-DeAngelis functional response and Allee effect.

    The study of the dynamics of a harvested population is a topic studied in mathematical bioeconomics [15], inside a larger chapter dealing with optimal management of renewable resources. The exploitation of biological resources and the harvesting of interacting species is applied in fisheries, forestry and fauna management [15,16,17]. Etoua and Rousseau [16] studied a generalized Gause model with both prey harvesting and Holling response function of type Ⅲ:

    {dxdt=rx(1xK)mx2yax2+bx+1h1,dydt=y(d+cmx2ax2+bx+1),x0, y0, (1.4)

    where the eight parameters: r,k,m,a,c,d,h are strictly positive and b0. Through the following linear transformation and time scaling

    (X,Y,T)=(1kx,1cky,cmk2t).

    Laurin and Rousseau[17] transformed the model (1.4) into the simplified system with the number of parameters reduced to five

    {˙x=ρx(1x)yx2αx2+βx+1λ,˙y=y(δ+x2αx2+βx+1),x0, y0, (1.5)

    with parameters

    (ρ,α,β,δ,λ)=(rcmk2,ak2,bk,dcmk2,hcmk3).

    And the Hopf bifurcation was studied in [16,17]. Du et al. [18] considered a general predator-prey model with prey and predator harvesting and proved that the predator-prey system has at least four positive periodic solutions. In addition, some other predator-prey models have been studied widely [19,20,21,22].

    In this paper, we consider a generalized Gause-type predator-prey model with Allee effect, Holling type Ⅲ functional response and harvesting terms

    {dxdt=r(t)(1x(t)K)(x(t)m(t))x(t)s(t)x2x2(t)+a2(t)y(t)H1(t),dydt=(p(t)x2(t)x2(t)+a2(t)b(t)y(t)c(t))y(t)H2(t),x(0)>0, y(0)>0,  t[0,T], (1.6)

    where x=x(t) and y=y(t) represent the population sizes of prey and predator at time t, respectively. The size can represent numbers of individuals or density in the unit space of the population. To ensure biological significance, the parameter of K is positive, and a, b, c, H1, H2, m, p, r, s are positive Tperiodic functions. The meaning of the parameters in system (1.6) is given as follows:

     a is the amount of prey at which predation rate is maximal.

     b is the predator population decays in the competition among the predators.

     c is the natural per capita death rate of the predator.

     K is the environmental capacity of the prey.

     m is the minimum viable population.

     p is the conversion efficiency of reduction rate of the predator.

     r is the growth rate of the prey.

     s is the maximum per capita consumption rate.

    In system (1.6), the Allee effect is defined by the term r(t)(1x(t)K)(x(t)m(t))x(t) and the Holling type-Ⅲ functional response is represented by the term s(t)x2x2(t)+a2(t). This Holling type-Ⅲ functional response describes a behavior in which the number of prey consumed per predator initially increases quickly as the density of prey grows and levels off with further increase in prey density [13]. H1(t) and H2(t) describe the harvesting rate of prey and predators. We consider four important assumptions as regards the interactions between prey and predator:

    the prey population is affected by the Allee effect,

    the functional response is Holling type-Ⅲ,

    the influence of artificial harvest is considered on predator and prey, and

    the predator population decays in the competition among the predators is investigated.

    In this paper, we establish some conditions to ensure that system (1.6) has at least two positive periodic solutions. We outline the format for the rest of this paper. In Sect. 2, we describe several technical lemmas. In Sect. 3, using a systematic qualitative analysis and employing the Mawhin coincidence degree theory, we obtain that system (1.6) has at least two positive Tperiodic solutions of system (1.6).

    In this section, we will give relevant definitions of the Mawhin coincidence degree theory [23] and several technical lemmas.

    Let both X and Y be Banach spaces, L: DomLXY be a linear map and N: X×[0,1]Y be a continuous map. If ImLY is closed and dimKerL=codimImL<+, then we call the operator L is a Fredholm operator of index zero. If L is a Fredholm operator with index zero and there exists continuous projections P: XX and Q: YY such that ImP=KerL and ImL=KerQ=Im(IQ), then L|DomLKerP:(IP)XImL has an inverse function, and we set it as KP. Assume that Ω×[0,1]X is an open set. If QN(¯Ω×[0,1]) is bounded and KP(IQ)N(¯Ω×[0,1])X is relatively compact, then we say that N(¯Ω×[0,1]) is Lcompact.

    Next, we will give the Mawhin coincidence degree theorem.

    Lemma 2.1. ([23,24]) Let X and Y be two Banach spaces, L: DomLXY be a Fredholm operator with index zero, ΩY be an open bounded set, and N:¯Ω×[0,1]X be Lcompact on ¯Ω×[0,1]. If all the following conditions hold

    [C1] LxλNx, for xΩDomL, λ(0,1);

    [C2] QNx0, for every xΩKerL;

    [C3] Brouwer degree deg{JQN,ΩKerL,0}0, where J: ImQKerL is an isomorphism.

    Then the equation Lx=Nx has at least one solution on ¯ΩDomL.

    Lemma 2.2. ([19]) Let x>0, y>0, z>0 and x>2yz. For functions f(x,y,z)=x+x24yz2z and g(x,y,z)=xx24yz2z, the following assertions hold:

    (ⅰ) f(x,y,z) and g(x,y,z) are monotonically increasing and monotonically decreasing with respect to the variable x(0,);

    (ⅱ) f(x,y,z) and g(x,y,z) are monotonically decreasing and monotonically increasing with respect to the variable y(0,);

    (ⅲ) f(x,y,z) and g(x,y,z) are monotonically decreasing and monotonically increasing with respect to the variable z(0,).

    Throughout this paper, we denote by C[0,T] the space of all bounded continuous functions f: RR, and denote by C+ the set of all functions fC and f0. For the convenience of statement, we use the notations as follows

    ¯f=1TT0f(t)dt, fL=mint[0,T]f(t), fM=maxt[0,T]f(t).

    In this section, we will establish the existence results of at least two positive periodic solutions for the system (1.6).

    Theorem 2.1. Assume the following conditions hold:

    (H1) (pMu+2u2+(aL)2cL)2>4bLHL2;

    (H2) (pLu2u+2+(aM)2cM)2>4bMHM2;

    (H3) the algebra equation system

    {ˉr(1euK)(euˉm)ˉseueve2u+(ˉa)2¯H1eu=0,ˉpe2ue2u+(ˉa)2ˉbevˉc¯H2ev=0,

    has finite real-valued solutions (uk,vk), k=1, 2,..., n, satisfying

    (uk,vk)detG(uk, vk)nk=1uknk=1vk0,

    where

    G(uk,vk)=(ˉreu(ˉmK+1euK)ˉseueve2u+(ˉa)2+2ˉse3uev[e2u+(ˉa)2]2+¯H1e2uˉseueve2u+(ˉa)22ˉpe2ue2u+(ˉa)2+2ˉpe4u[e2u+(ˉa)2]2ˉbev+¯H2e2v).

    Then system (1.6) has at least two positive Tperiodic solutions.

    Proof of Theorem 2.1.. Suppose (x(t),y(t))R2 is an arbitrary positive of system (1.6). Let x=eu(t), y=ev(t), it follows from system (1.6) we can obtain

    {˙u(t)=r(t)(1eu(t)K)(eu(t)m(t))s(t)eu(t)ev(t)e2u(t)+a2(t)H1(t)eu(t),˙v(t)=p(t)e2u(t)e2u(t)+a2(t)b(t)ev(t)c(t)H2(t)ev(t), (3.1)

    where ˙=ddt.

    Let

    X=Y={z(t)=(u(t), v(t))TC(R, R2):z(t+T)z(t)},

    be equipped with the norm

    z(t)=(u(t), v(t))T=maxt[0,T]|u(t)|+maxt[0,T]|v(t)|,

    where X and Y are Banach spaces, T is the transpose.

    Taking zX and then we will define operators of L, P and Q as follows.

    Firstly, let

    L: DomLXY, Lz=dzdt.

    It is clear that

    KerL={zdomL: z=c, cR2},

    that is dimKerL=dimR2=2. Next we calculate ImL. Let

    dzdt=y(t), y(t)Y.

    Integrating both sides of this equation, we have

    T0dzdtdt=T0y(t)dt,

    thus

    T0y(t)dt=z(T)z(0)=0.

    From

    X=Y={z(t)=(u(t), v(t))TC(R, R2):z(t+T)z(t)},

    we can obtain y(t)=z(t), that is

    ImL={zY: T0z(t)dt=0}

    is closed in Y. Obviously, ImLR2={0}.

    Considering P, Q are both continuous projections satisfying

    ImP=KerL, ImL=KerQ=Im(IQ).

    Let

    P: XKerL,

    then, we get P(z) is a constant. Here, we denote it by

    P(z)=1TT0z(t)dt.

    Secondly, let

    Q: YYImL, β=T0z(t)dt and Q(z)=αβ,

    then, we have

    Q(Q(z))=Q(αβ)=αT0αβdt=α2βT0dt=α2βT=Q(z)=αβ,

    i.e.

    α=1T.

    Hence,

    Q(z)=1TT0z(t)dt.

    Thirdly, for zY, z1(t)=z(t)Q(z), we're going to verify z1(t)ImL, i.e. T0z1(t)dt=0. Here

    T0z1(t)dt=T0z(t)dtT0Q(z)dt=βαT0T0z(t)dtdt=βαβT=ββ=0,

    that is

    z1(t)ImL.

    Moreover, we can obtain

    Y=ImLR2, codimImL=dimR2=2.

    i.e.,

    dimKerL=codimImL.

    So L is a Fredholm operator with index zero, which implies L has a unique inverse. We define by KP: ImLKerPDomL the inverse of L.

    By simply calculating, we have

    KP(z)=t0z(w)dw1TT0t0z(w)dwdt.

    Define N: XY by the form

    Nz=(Δ1(z(t), t)Δ2(z(t), t)),

    where

    Δ1(z(t), t)=r(t)(1eu(t)K)(eu(t)m(t))s(t)eu(t)ev(t)e2u(t)+a2(t)H1(t)eu(t),Δ2(z(t), t)=p(t)e2u(t)e2u(t)+a2(t)b(t)ev(t)c(t)H2(t)ev(t).

    Thus

    QNz=(1TT0Δ1(z(t), t)dt1TT0Δ2(z(t), t)dt),

    and

    KP(IQ)Nz=(t0Δ1(z(w), w)dwt0Δ2(z(w), w)dw)(1TT0t0Δ1(z(w), w)dwdt1TT0t0Δ2(z(w), w)dwdt) (1Tt0T0Δ1(z(w), w)dwdw1Tt0T0Δ2(z(w), w)dwdw)+(1T2T0t0T0Δ1(z(w), w)dwdwdt1T2T0t0T0Δ2(z(w), w)dwdwdt)=(t0Δ1(z(w), w)dw1TT0t0Δ1(z(w), w)dwdt(tT12)T0Δ1(z(w), w)dwt0Δ2(z(w), w)dw1TT0t0Δ2(z(w), w)dwdt(tT12)T0Δ2(z(w), w)dw). (3.2)

    Let ΩX be bounded. For zΩ, we have that zM1, |u(t)|M1 and |v(t)|M1.

    Next, we see that QN(¯Ω) is bounded.

    |1TT0Δ1(z(t), t)dt||1TT0r(t)(1eu(t)K)(eu(t)m(t))dt| +|1TT0s(t)eu(t)ev(t)e2u(t)+a2(t)dt|+|1TT0H1(t)eu(t)dt||1TT0[r(t)eu(t)r(t)m(t)r(t)e2u(t)K+r(t)m(t)eu(t)K]dt| +|1TT0s(t)eu(t)ev(t)a2(t)dt|+|1TT0H1(t)eu(t)dt|ˉreM1+¯rm+ˉre2M1K+¯rmeM1K+¯(sa2)e2M1+¯H1eM1,

    and

    |1TT0Δ2(z(t), t)dt||1TT0p(t)e2u(t)e2u(t)+a2(t)dt|+|1TT0b(t)ev(t)dt| +|1TT0c(t)dt|+|1TT0H2(t)ev(t)dt||1TT0p(t)e2u(t)+p(t)a2(t)e2u(t)+a2(t)dt|+ˉbeM1+ˉc+¯H2eM1=|1TT0p(t)dt|+ˉbeM1+ˉc+¯H2eM1=ˉp+ˉbeM1+ˉc+¯H2eM1.

    It is immediate that QN and KP(IQ)N are continuous.

    Consider a sequence of function {z}Ω. We have the following inequality for the first function of KP(IQ)Nz.

     KP(IQ)NΔ1(z(t1), t1)KP(IQ)NΔ1(z(t2), t2)=t1t2[r(w)(1eu(w)K)(eu(w)m(w))s(w)eu(w)ev(w)e2u(w)+a2(w)H1(w)eu(w)]dw (t1t2T)T0[r(w)(1eu(w)K)(eu(w)m(w))s(w)eu(w)ev(w)e2u(w)+a2(w)H1(w)eu(w)]dw(t1t2)[rMeM1rLmLrLe2M1K+rMmMeM1K] (t1t2)[ˉreM1¯rmˉre2M1K+¯rmeM1K].

    For another function, we have similar inequalities as follows

     KP(IQ)NΔ2(z(t1), t1)KP(IQ)NΔ2(z(t2), t2)=t1t2[p(w)e2u(w)e2u(w)+a2(w)b(w)ev(w)c(w)H2(w)ev(w)]dw (t1t2T)T0[p(w)e2u(w)e2u(w)+a2(w)b(w)ev(w)c(w)H2(w)ev(w)]dw(t1t2)[pMe2M1a2LbLeM1cLH2LeM1] +(t1t2)[ˉbeM1+ˉc+¯H2eM1].

    Hence the sequence {KP(IQ)Nz} is equicontinuous. Using the periodicity of the functions, we know that the sequence {KP(IQ)Nz} is uniformly bounded.

    An application of Ascoli-Arzela's theorem shows that {KP(IQ)N(¯Ω)} is compact for any bounded set ΩX. Since QN(¯Ω) is bounded, we conclude that N is Lcompact on Ω for any bounded set ΩX.

    Then, considering the operator equation Lx=λNx, as follows

    {˙u(t)=λΔ1(z(t), t),˙v(t)=λΔ2(z(t), t), (3.3)

    where λ(0,1). Let

    u(ξ1)=maxt[0,T]u(t),  u(η1)=mint[0,T]u(t),
    v(ξ2)=maxt[0,T]v(t),  v(η2)=mint[0,T]v(t).

    Through simple analysis, we have

    ˙u(ξ1)=˙u(η1)=0, ˙v(ξ2)=˙v(η2)=0.

    From (3.3), we can find that

    r(ξ1)(1eu(ξ1)K)(eu(ξ1)m(ξ1))s(ξ1)eu(ξ1)ev(ξ1)e2u(ξ1)+a2(ξ1)H1(ξ1)eu(ξ1)=0, (3.4)
    p(ξ2)e2u(ξ2)e2u(ξ2)+a2(ξ2)b(ξ2)ev(ξ2)c(ξ2)H2(ξ2)ev(ξ2)=0, (3.5)

    and

    r(η1)(1eu(η1)K)(eu(η1)m(η1))s(η1)eu(η1)ev(η1)e2u(η1)+a2(η1)H1(η1)eu(η1)=0 (3.6)
    p(η2)e2u(η2)e2u(η2)+a2(η2)b(η2)ev(η2)c(η2)H2(η2)ev(η2)=0, (3.7)

    In view of (3.4), we have

    r(ξ1)(1eu(ξ1)K)(eu(ξ1)m(ξ1))=s(ξ1)eu(ξ1)ev(ξ1)e2u(ξ1)+a2(ξ1)+H1(ξ1)eu(ξ1)>0,

    then, we get

    rLKe2u(ξ1)(rM+rMmMK)eu(ξ1)+rLmL<0,

    which implies that

    eu=:rM+rMmMK(rM+rMmMK)24(rL)2mLK2rLK<eu(ξ1)<rM+rMmMK+(rM+rMmMK)24(rL)2mLK2rLK:=eu+.

    Similarly, we can discuss the range of eu(η1) from (3.6)

    r(η1)(1eu(η1)K)(eu(η1)m(η1))=s(η1)eu(η1)ev(η1)e2u(η1)+a2(η1)+H1(η1)eu(η1)>0,

    A direct calculation gives

    rLKe2u(η1)(rM+rMmMK)eu(η1)+rLmL<0,

    so we can obtain

    eu=:rM+rMmMK(rM+rMmMK)24(rL)2mLK2rLK<eu(η1)<rM+rMmMK+(rM+rMmMK)24(rL)2mLK2rLK:=eu+.

    From (3.5), we have

    b(ξ2)e2v(ξ2)(p(ξ2)e2u(ξ2)e2u(ξ2)+a2(ξ2)c(ξ2))ev(ξ2)+H2(ξ2)=0,

    and

    bLe2v(ξ2)(pMe2u+e2u+(aL)2cL)ev(ξ2)+H2L<0,

    then, we get

    ev=:pMe2u+e2u+(aL)2cL(pMe2u+e2u+(aL)2cL)24bLH2L2bL<ev(ξ2)<pMe2u+e2u+(aL)2cL+(pMe2u+e2u+(aL)2cL)24bLH2L2bL:=ev+.

    By (3.7), we obtain

    b(η2)e2v(η2)(p(η2)e2u(η2)e2u(η2)+a2(η2)c(η2))ev(η2)+H2(η2)=0,

    and

    bLe2v(η2)(pMe2u+e2u+(aL)2cL)ev(η2)+H2L<0,

    this implies that

    ev=:pMu2+e2u+(aL)2cL(pMu2+e2u+(aL)2cL)24bLH2L2bL<ev(η2)<pMu2+e2u+(aL)2cL+(pMe2u+e2u+(aL)2cL)24bLH2L2bL:=ev+.

    And then, in view of (3.5) and (3.7) we have

    bMe2v(ξ2)(pLe2ue2u++(aM)2cM)ev(ξ2)+H2M>0,

    that is

    ev(ξ2)>pLe2ue2u++(aM)2cM+(pLe2ue2u++(aM)2cM)24bMH2M2bM:=el+

    or

    ev(ξ2)<pLe2ue2u++(aM)2cM(pLe2ue2u++(aM)2cM)24bMH2M2bM:=el.

    From (3.7), we obtain

    bMe2v(η2)(pLe2ue2u++(aM)2cM)ev(η2)+H2M>0,

    i.e.

    ev(ξ2)>l+ or ev(ξ2)<el.

    In view of Lemma 2.2, we can find that v<l<l+<v+. Thus, we have

    {v(ξ2)>l+ or v(ξ2)<l,v<v(ξ2)<v+,

    and

    {v(η2)>l+ or v(η2)<l,v<v(η2)<v+,

    that is

    v(ξ2)(v, l)(l+, v+), v(η2)(v, l)(l+, v+).

    Similarly, we get

    u(ξ1)(u, u+), u(η1)(u, u+).

    It is easy to know that u±,v±,l± are independent of λ. Consider the following two sets

    Ω1={z=(u,v)TX|u<u< u+, v<v< l},
    Ω2={z=(u,v)TX|u<u< u+, l+<v< v+}.

    Obviously, ΩiX and ¯Ω1¯Ω2=. So Ωi's (i=1, 2) satisfy the condition [C1] of Lemma 2.1.

    Next, we show that QNz0, for zΩiKerL=ΩiR2 (i=1,2).

    If it is not true, then there exists (u,v)TΩi, such that

    {T0Δ1(z(t), t)dt=0,T0Δ2(z(t), t)dt=0.

    By virtue of the mean value theorem, there exists two points tj[0,T] (j=1, 2) satisfying

    {Δ1(z(t1), t1)=0,Δ2(z(t2), t2)=0.

    So, we obtain

    u(u, u+),

    and

    v(v, l)(l+, v+),

    which contradicts (u,v)TΩiR2. So the condition [C2] in Lemma 2.1 holds.

    Then, we check the condition [C3] in Lemma 2.1. Define the homomorphism J: ImQKerL, Jzz. From [H3], we have

    deg{JQN,ΩiKerL,0}=zkQN1(0)sgnJQN(zk)=zkQN1(0)detG(uk,vk)nk=1uknk=1vk0.

    This implies that the condition [C3] in Lemma 2.1 holds too. Note that, Ω1 and Ω2 satisfies all conditions of Lemma 2.1. Therefore, system (1.6) has at least two Tperiodic solutions.

    Here, we would like to give two remarks.

    Remark 2.1. If we take H1(t)=0,H2(t)=0 and b(t)=0, i.e., system (1.6) without considering the harvesting terms of prey and predator, as well as the predator competition, we find that the system (1.3) in [13] is the system (1.6).

    Remark 2.2. In [16] and [17], the authors only considered the Gause model (1.4) and (1.5) with prey harvesting h1, respectively, but they don't investigate the influence of the predator harvesting. In fact, the influence of the predator harvesting is very important in biological populations and bioeconomics, especially in fisheries management etc.

    In this paper, we are concerned with a Gause-type predator-prey model with Allee effect, Holling type Ⅲ functional response and the artificial harvesting terms, which are very important in biological populations and bioeconomics. Four important assumptions as regards the interactions between prey and predator is considered. By applying the Mawhin coincidence degree theory, we obtain the existence of multiple positive periodic solutions for the predator-prey model.

    We express our sincere thanks to the anonymous reviewers for their valuable comments and suggestions. This work is supported by the Natural Science Foundation of China (Grant No.11771185). The first author supported by Postgraduate Research and Practice Innovation Program of Jiangsu Province (Grant Nos. KYCX18-2091, KYCX20-2082 and KYCX20-2206).

    The authors declare no conflict of interest in this paper.



    Author contributions



    Ferran Cuenca-Martínez, Mónica Grande-Alonso: conceptualization, supervision; all authors: methodology, validation, data curation, writing—original draft preparation, writing—review and editing, visualization; Mónica Grande-Alonso: investigation, project administration. All authors have read and agreed to the published version of the manuscript.

    Conflict of interest



    The authors declare no conflict of interest.

    [1] Hoy D, Bain C, Williams G, et al. (2012) A systematic review of the global prevalence of low back pain. Arthritis Rheum-Us 64: 2028-2037. https://doi.org/10.1002/art.34347
    [2] Anderssen G, Frymoyer J (1997) The epidemiology of spinal disorders. The Adult Spine: Principles and Practice, Edited by: Frymoyer JW . https://doi.org/10.12691/ajmsm-6-1-3
    [3] Nachemson A (1999) Back pain: delimiting the problem in the next millennium. Int J Law Psychiat 22: 470-490. https://doi.org/10.1016/s0160-2527(99)00022-9
    [4] Bongers P, De Winter C, Kompier M, et al. (1993) Psychosocial factors at work and musculoskeletal disease. Scand J Work Env Hea 19: 297-312. https://doi.org/10.5271/sjweh.1470
    [5] Hall H, McIntosh G (2008) Low back pain (chronic). BMJ Clin Evid .
    [6] Pincus T, Burton A, Vogel S, et al. (2002) A systematic review of psychological factors as predictors of chronicity/disability in prospective cohorts of low back pain. Spine 27: 109-120. https://doi.org/10.1097/00007632-200203010-00017
    [7] Rossignol M, Rozenberg S, Leclerc A (2009) Epidemiology of low back pain: What's new?. Joint Bone Spine 76: 608-613. https://doi.org/10.1016/j.jbspin.2009.07.003
    [8] Bialosky J, Bishop M, Pric D, et al. (2009) The mechanisms of manual therapy in the treatment of musculoskeletal pain: A comprehensive model. Manual Ther 14: 531-538. https://doi.org/10.1016/j.math.2008.09.001
    [9] Teodorczyk-Injeyan J, Injeyan H, Ruegg R (2006) Spinal manipulative therapy reduces inflammatory cytokines but not substance P production in normal subjects. J Manip Physiol Ther 29: 14-21. https://doi.org/10.1016/j.jmpt.2005.10.002
    [10] Kurt V, Aras O, Buker N (2020) Comparison of conservative treatment with and without neural mobilization for patients with low back pain: A prospective, randomized clinical trial. J Back Musculoskelet 33: 969-975. https://doi.org/10.3233/BMR-181241
    [11] Peacock M, Douglas S, Nair P (2023) Neural mobilization in low back and radicular pain: A systematic review. J Man Manip Ther 31: 4-12. https://doi.org/10.1080/10669817.2022.2065599
    [12] Dworkin B, Elbert T, Rau H, et al. (1994) Central effects of baroreceptor activation in humans: attenuation of skeletal reflexes and pain perception. P Natl Acad Sci Usa 91: 6329-6333. https://doi.org/10.1073/pnas.91.14.6329
    [13] Koltyn K (2000) Analgesia following exercise: A review. Sports Med 29: 85-98. https://doi.org/10.2165/00007256-200029020-00002
    [14] Fletcher C, Bradnam L, Barr C (2016) The relationship between knowledge of pain neurophysiology and fear avoidance in people with chronic pain: A point in time, observational study. Physiother Theor Pr 32: 271-276. https://doi.org/10.3109/09593985.2015.1138010
    [15] López-de-Uralde-Villanueva I, Beltran-Alacreu H, Fernández-Carnero J, et al. (2020) Pain management using a multimodal physiotherapy program including a biobehavioral approach for chronic nonspecific neck pain: a randomized controlled trial. Physiother Theor Pr 36: 45-62. https://doi.org/10.1080/09593985.2018.1480678
    [16] Nolan M (1985) Quantitative measure of cutaneous sensation. Two-point discrimination values for the face and trunk. Phys Ther 65: 181-185. https://doi.org/10.1093/ptj/65.2.181
    [17] Starkweather A, Heineman A, Storey S, et al. (2016) Methods to measure peripheral and central sensitization using quantitative sensory testing: A focus on individuals with low back pain. Appl Nurs Res 29: 237-241. https://doi.org/10.1016/j.apnr.2015.03.013
    [18] Maier C, Baron R, Tölle T, et al. (2010) Quantitative sensory testing in the German Research Network on Neuropathic Pain (DFNS): Somatosensory abnormalities in 1236 patients with different neuropathic pain syndromes. Pain 150: 439-450. https://doi.org/10.1016/j.pain.2010.05.002
    [19] Kinser A, Sands W, Stone M (2009) Reliability and validity of a pressure algometer. J Strength Cond Res 23: 312-314. https://doi.org/10.1519/jsc.0b013e31818f051c
    [20] Pourahmadi MR, Ebrahimi Takamjani I, Sarrafzadeh J, et al. (2017) Reliability and concurrent validity of a new iPhone® goniometric application for measuring active wrist range of motion: a cross-sectional study in asymptomatic subjects. J Anat 230: 484-495. https://doi.org/10.1111/joa.12568
    [21] Chattanooga GroupStabilizer TM pressure bio-feedback operating instructions (2002).
    [22] Roberts H, Denison H, Martin H, et al. (2011) A review of the measurement of grip strength in clinical and epidemiological studies: Towards a standardised approach. Age Ageing 40: 423-429. https://doi.org/10.1093/ageing/afr051
    [23] Duncan P, Weiner D, Chandler J, et al. (1990) Functional reach: A new clinical measure of balance. JOGE 45: M192-M197. https://doi.org/10.1093/geronj/45.6.M192
    [24] Newton RA (2001) Validity of the multi-directional reach test: a practical measure for limits of stability in older adults. J Gerontol Med Sci 56: M248-M252. https://doi.org/10.1093/gerona/56.4.m248
    [25] Shaffer S, Teyhen D, Lorenson C, et al. (2013) Y-balance test: A reliability study involving multiple raters. Mil Med 178: 1264-1270. https://doi.org/10.7205/MILMED-D-13-00222
    [26] Moreau CE, Green BN, Johnson CD, et al. (2001) Isometric back extension endurance tests: A review of the literature. J Manip Physiol Ther 24: 110-122. https://doi.org/10.1067/mmt.2001.112563
    [27] Alcántara-Bumbiedro S, Flórez-García MT, Echávarri-Pérez C, et al. (2006) Oswestry low back pain disability questionnaire. Rehabilitacion 40: 150-158. https://doi.org/10.1016/S0048-7120(06)74881-2
    [28] Gómez-Pérez L, López-Martínez AE, Ruiz-Párraga GT, et al. (2011) Psychometric properties of the spanish version of the Tampa Scale for Kinesiophobia (TSK). Journal of Pain 12: 425-435. https://doi.org/10.1016/j.jpain.2010.08.004
    [29] García Campayo J, Rodero B, Alda M, et al. (2008) Validación de la versión española de la escala de la catastrofización ante el dolor (pain catastrophizing scale) en la fibromyalgia, (Spanish) [Validation of the Spanish version of the pain catastrophizing scale in fibromyalgia]. Med Clin-Barcelona 131: 487-492. https://doi.org/10.1157/13127277
    [30] George S, Valencia C, Beneciuk J (2010) A Psychometric investigation of fear-avoidance model measures in patients with chronic low back pain. J Orthop Sport Phys 40: 197-205. https://doi.org/10.2519/jospt.2010.3298
    [31] Vilagut G, Valderas JM, Ferrer M, et al. (2008) Interpretación de los cuestionarios de salud SF-36 y SF-12 en España: componentes físico y menta, (Spanish) [Interpretation of SF-36 and SF-12 questionnaires in Spain: physical and mental components]. Med Clin 130: 726-735. https://doi.org/10.1157/13121076
    [32] Quintana J, Padierna A, Esteban C, et al. (2003) Evaluation of the psychometric characteristics of the Spanish version of the Hospital Anxiety and Depression Scale. Acta Psychiat Scand 107: 216-221. https://doi.org/10.1034/j.1600-0447.2003.00062.x
    [33] Martín-Aragón M, Pastor MA, Rodríguez-Marín J, et al. (1999) Percepción de autoeficacia en dolor crónico. Adaptación y validación de la chronic pain selfefficacy scale, (Spanish) [Perception of self-efficacy in chronic pain. Adaptation and validation of the chronic pain selfefficacy scale]. J Health Psychol 11: 51-75. https://doi.org/10.21134/pssa.v11i1.799
    [34] Grande-Alonso M, Suso-Martí L, Cuenca-Martínez F, et al. (2019) Physiotherapy based on a biobehavioral approach with or without orthopedic manual physical therapy in the treatment of nonspecific chronic low back pain: A randomized controlled trial. Pain Med 20: 2571-2587. https://doi.org/10.1093/pm/pnz093
    [35] Thomas M, Thomson OP, Kolubinski DC, et al. (2023) The attitudes and beliefs about manual therapy held by patients experiencing low back pain: a scoping review. Musculoskel Sci Prac 102752. https://doi.org/10.1016/j.msksp.2023.102752
    [36] Evers A, Colloca L, Blease C, et al. (2018) Implications of placebo and nocebo effects for clinical practice: expert consensus. Psychother Psychosom 87: 204-210. https://doi.org/10.1159/000490354
    [37] Ferrer-Peña R, Muñoz-Garcıa D, Calvo-Lobo C, et al. (2019) Pain expansion and severity reflect central sensitization in primary care patients with greater trochanteric pain syndrome. Pain Med 20: 961-970. https://doi.org/10.1093/pm/pny199
    [38] Lluch Girbés E, Dueñas L, Barbero M, et al. (2016) Expanded distribution of pain as a sign of central sensitization in individuals with symptomatic knee osteoarthritis. Phys Ther 96: 1196-1207. https://doi.org/10.2522/ptj.20150492
    [39] Searle A, Spink M, Ho A, et al. (2015) Exercise interventions for the treatment of chronic low back pain: A systematic review and meta-analysis of randomised controlled trials. Clin Rehabil 29: 1155-1167. https://doi.org/10.1177/0269215515570379
    [40] Niederer D, Mueller J (2020) Sustainability effects of motor control stabilisation exercises on pain and function in chronic nonspecific low back pain patients: A systematic review with meta-analysis and meta-regression. Plos One 15: e0227423. https://doi.org/10.1371/journal.pone.0227423
    [41] Geneen L, Moore R, Clarke C, et al. (2017) Physical activity and exercise for chronic pain in adults: An overview of Cochrane Reviews. Cochrane Db Syst Rev 4. https://doi.org/10.1002/14651858.CD011279.pub3
    [42] Alcon C, Wang-Price S (2022) Non-invasive brain stimulation and pain neuroscience education in the cognitive-affective treatment of chronic low back pain: Evidence and future directions. Front Pain Res 3: 959609. https://doi.org/10.3389/fpain.2022.959609
    [43] Lepri B, Romani D, Storari L, et al. (2023) Effectiveness of pain neuroscience education in patients with chronic musculoskeletal pain and central sensitization: A systematic review. Int J Env Res Pub Health 20: 4098. https://doi.org/10.3390/ijerph20054098
    [44] Kim K, An J, Kim JO, et al. (2022) Effects of pain neuroscience education combined with lumbar stabilization exercise on strength and pain in patients with chronic low back pain: randomized controlled trial. J Pers Med 12: 303. https://doi.org/10.3390/jpm12020303
  • Reader Comments
  • © 2024 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(1657) PDF downloads(57) Cited by(1)

Figures and Tables

Figures(3)  /  Tables(5)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog