
Knee medial compartment osteoarthritis is effectively treated by a medial open-wedge high tibial osteotomy (MOWHTO). The feasibility and safety of MOWHTO for mild lateral meniscal tears are unknown. This study examined the feasibility and safety of knee joint weight-bearing line ratio (WBLr) adjustment during MOWHTO with lateral meniscal injuries. We used a healthy adult male's lower extremities computed tomography scans and knee joint magnetic resonance imaging images to create a normal fine element (FE) model. Based on this model, we generated nine FE models for the MOWHTO operation (WBLr: 40–80%) and 15 models for various lateral meniscal injuries. A compressive load of 650N was applied to all cases to calculate the von Mises stress (VMS), and the intact lateral meniscus' maximal VMS at 77.5% WBLr was accepted as the corrective upper limit stress. Our experimental results show that mild lateral meniscal tears can withstand MOWHTO, while severe tears cannot. Our findings expand the use of MOWHTO and provide a theoretical direction for practical decisions in patients with lateral meniscal injuries.
Citation: Lin Chen, Mingjun Wang, Zhanyu Wu, Jinbo Sun, Jianglong Li, Chun Chen, Chuan Ye. Advancements in health informatics: finite element insights into medial open-wedge high tibial osteotomy and lateral meniscal tears[J]. Mathematical Biosciences and Engineering, 2024, 21(4): 5394-5410. doi: 10.3934/mbe.2024237
[1] | Ya-Dong Zhang, Hai-Feng Huo, Hong Xiang . Dynamics of tuberculosis with fast and slow progression and media coverage. Mathematical Biosciences and Engineering, 2019, 16(3): 1150-1170. doi: 10.3934/mbe.2019055 |
[2] | H. M. Srivastava, Khaled M. Saad, J. F. Gómez-Aguilar, Abdulrhman A. Almadiy . Some new mathematical models of the fractional-order system of human immune against IAV infection. Mathematical Biosciences and Engineering, 2020, 17(5): 4942-4969. doi: 10.3934/mbe.2020268 |
[3] | Chaoyue Wang, Zhiyao Ma, Shaocheng Tong . Adaptive fuzzy output-feedback event-triggered control for fractional-order nonlinear system. Mathematical Biosciences and Engineering, 2022, 19(12): 12334-12352. doi: 10.3934/mbe.2022575 |
[4] | Tao-Li Kang, Hai-Feng Huo, Hong Xiang . Dynamics and optimal control of tuberculosis model with the combined effects of vaccination, treatment and contaminated environments. Mathematical Biosciences and Engineering, 2024, 21(4): 5308-5334. doi: 10.3934/mbe.2024234 |
[5] | Eduardo Ibargüen-Mondragón, M. Victoria Otero-Espinar, Miller Cerón Gómez . A within-host model on the interaction dynamics between innate immune cells and Mycobacterium tuberculosis. Mathematical Biosciences and Engineering, 2025, 22(3): 511-527. doi: 10.3934/mbe.2025019 |
[6] | Jimin Yu, Jiajun Yin, Shangbo Zhou, Saiao Huang, Xianzhong Xie . An image super-resolution reconstruction model based on fractional-order anisotropic diffusion equation. Mathematical Biosciences and Engineering, 2021, 18(5): 6581-6607. doi: 10.3934/mbe.2021326 |
[7] | Cristiana J. Silva, Helmut Maurer, Delfim F. M. Torres . Optimal control of a Tuberculosis model with state and control delays. Mathematical Biosciences and Engineering, 2017, 14(1): 321-337. doi: 10.3934/mbe.2017021 |
[8] | Shuai Zhang, Yongqing Yang, Xin Sui, Yanna Zhang . Synchronization of fractional-order memristive recurrent neural networks via aperiodically intermittent control. Mathematical Biosciences and Engineering, 2022, 19(11): 11717-11734. doi: 10.3934/mbe.2022545 |
[9] | Georgi Kapitanov . A double age-structured model of the co-infection of tuberculosis and HIV. Mathematical Biosciences and Engineering, 2015, 12(1): 23-40. doi: 10.3934/mbe.2015.12.23 |
[10] | Wenjie Qin, Jiamin Zhang, Zhengjun Dong . Media impact research: a discrete SIR epidemic model with threshold switching and nonlinear infection forces. Mathematical Biosciences and Engineering, 2023, 20(10): 17783-17802. doi: 10.3934/mbe.2023790 |
Knee medial compartment osteoarthritis is effectively treated by a medial open-wedge high tibial osteotomy (MOWHTO). The feasibility and safety of MOWHTO for mild lateral meniscal tears are unknown. This study examined the feasibility and safety of knee joint weight-bearing line ratio (WBLr) adjustment during MOWHTO with lateral meniscal injuries. We used a healthy adult male's lower extremities computed tomography scans and knee joint magnetic resonance imaging images to create a normal fine element (FE) model. Based on this model, we generated nine FE models for the MOWHTO operation (WBLr: 40–80%) and 15 models for various lateral meniscal injuries. A compressive load of 650N was applied to all cases to calculate the von Mises stress (VMS), and the intact lateral meniscus' maximal VMS at 77.5% WBLr was accepted as the corrective upper limit stress. Our experimental results show that mild lateral meniscal tears can withstand MOWHTO, while severe tears cannot. Our findings expand the use of MOWHTO and provide a theoretical direction for practical decisions in patients with lateral meniscal injuries.
TB is a deadly infectious disease and one of the world's top ten causes of death [1,2,3]. Most of the infected individuals can be cured by adopting a systematical diagnostic and using appropriate antibiotics. However, the impact of TB on human health and social development still cannot be ignored due to the lack of knowledge about TB prevention and diagnosis and the low treatment coverage in some areas. According to the global tuberculosis report 2021 published by the World Health Organization [4], about 9.9 million people worldwide suffered from TB in 2020, of which more than 1.3 million died. The massive destruction caused by TB to the economy and society forces us to understand the transmission mechanism of TB deeply.
Mathematical models play an essential role in understanding the dynamics of infectious diseases, including TB [5,6,7,8,9,10]. Exploring the transmission mechanism of tuberculosis is the first step to establishing a mathematical model, so references [11,12,13] have studied the transmission dynamics of tuberculosis. In addition, Castillo-Chavez and Feng [14] formulated one-strain and two-strain TB models to analyze the effects of fundamental epidemiological factors. Subsequently, the TB compartment model was proposed by Feng et al. [15] and is considered one of the basic frameworks of the transmission mechanism of TB. The authors in [16] established the TB models with fast and slow dynamics by considering multi-level contact structure. Castillo-Chavez and Song [17] comprehensively summarized the research trends of the tuberculosis dynamic models and control problems in the past decades, and proposed some exciting and challenging questions. It should be pointed out that the research results of [15] indicate that it is not enough to reduce the basic reproduction number to less than 1 for the prevention and control of TB. After that, the researchers further discussed the influence of various factors on the transmission of TB by proposing different dynamical models [18,19,20,21,22,23,24,25]. For example, Hattaf et al. [23] established a TB disease model with exogenous reinfection, and studied the optimal control problem to reduce the infectious group by using Pontryagin's maximum principle. The multi-drug resistance, and case detection and treatment were considered in [24] and [25], respectively. Meanwhile, based on the actual infection data of TB in the United States from 1988 to 2019, the authors in [26] established an SVEITR dynamical model by considering vaccination, incomplete treatment and relapse, and analyzed the sensitivity of basic reproduction number and solved the optimal control problem issue by using Pontryagin's maximum principle.
It should be pointed out that the models in all of the literature mentioned above are established by using the conventional integer-order calculus, which is a local concept. Fractional-order calculus, as an extension of integer-order one, is a global concept that can more effectively characterize the hereditary and memory properties of many physical processes and materials [27,28]. Therefore, some scholars try to describe the spread and outbreak of infectious disease by using the fractional dynamics model [29,30,31,32,33,34]. The results show that the fractional order model can better fit the actual data. Although the fractional-order model has the above advantages, the control of the fractional TB model is still in the ascendant due to less application of fractional calculus in biology, which is one of the motivations of the current paper.
Adaptive control is one of the important control methods in engineering [35,36]. Its characteristic is that the control target can be given in advance, and then the control variable can be designed so that the state variable or the related variable of the model can track the control target with small errors, so as to realize the control of the model. Recently, some results adopt the adaptive control method to realize the control of biology [37,38,39]. In particular, the authors in [39] studied the nonlinear adaptive control of COVID-19. However, the system parameters are estimated by the adaptive laws in [37,38,39], which are heavily dependent on the initial values of adaptive laws. Therefore, how to reduce this dependency is also one of the motivations of this paper.
Motivated by the above discussions, the problem of nonlinear neural networks control is studied for the fractional-order TB model. The main contributions of this paper are summarized as follows: (ⅰ) By analyzing the transmission mechanism of TB and considering the historical dependence on the spreading of TB, a fractional-order TB model under Caputo fractional derivative is established. Compared with the model in [33], in this paper we consider that the proportion of smear-positive and smear-negative from susceptible population is not equal to that from latent population, and the patients after rehabilitation may become susceptible ones again. (ⅱ) This is a successful attempt to apply the universal approximation principle of radial basis function neural networks to biological model, and thus the controllers for fractional-order TB model are proposed and the initial values dependency of adaptive laws are also reduced. Different from [37,38,39], where the model parameters are estimated by adaptive laws, the unknown functions that include model parameters are approximated in this paper, and thus the processes of estimations to model parameters are avoided.
This paper is organized as follows. The problem to be studied is formulated and some preliminaries are presented in Section 2. The adaptive neural networks controllers are designed and the stability of controlled model is also analyzed in Section 3. Some numerical examples are given in Section 4 to illustrate the proposed controllers. Finally, a brief summary is presented in Section 5.
Definition 1. [40] The Caputo derivative of function f(t)∈Cn+1([t0,+∞),R) with the order α is given as follows
Dαf(t)=1Γ(n−α)∫tt0f(n)(τ)(t−τ)α+1−ndτ, |
where Γ(⋅) denotes Gamma function, n∈N+ and n−1<α≤n.
Definition 2. [40] For α>0, the Mittag-Leffler function with a single parameter is defined as Eα(z)=∑∞i=0ziΓ(αi+1), where z is plural. In particular, E1(z)=ez.
Lemma 1. [40] Let the function f(t)∈C1([t0,+∞),R), then Dαf2(t)≤2f(t)Dαf(t) for t≥t0 and α∈(0,1].
Lemma 2. Consider the system Dαx(t)=φ(t)−x(t) with initial value x(0)>0. If 0<α<1 and function φ(t)(t≥0) is a non-negative and bounded continuous function, then the solution x(t) of this system is positive and bounded.
Proof. This Lemma can be proven by using the similar method in [44,Lemma 5], so we omit it here.
Lemma 3. [41] For a given approximation accuracy ˜ϵ>0, any continuous function f(x) defined on a compact set can be approximated by the radial basis function neural networks W∗TΨ(x) as follows
f(x)=W∗TΨ(x)+ϵ(x),|ϵ(x)|≤˜ϵ, |
where W∗=argminW[supx|f(x)−WTΨ(x)|], x represents the input vector, W∗=(w1,⋯,wl)T∈Rl denotes the ideal weight vector, l>1 represents the nodes number, Ψ(x)=(Ψ1(x),⋯,Ψl(x))T, Ψi(x) is selected as Ψi(x)=exp[−(x−σi)T(x−σi)/κ2] for i=1,⋯,l, where σi=(σi1,⋯,σin)T and κ are the center and width of the Gaussian function, respectively.
Since it usually takes four to eight weeks from the time TB enters the body to the time it shows symptoms, and people with TB can have either smear positive or smear negative. Therefore, the people carrying mycobacterium TB is divided into latent population (E), smear positive population (I1) and smear negative population (I2). The rest of the population is classified as susceptible (S) and immune or quarantined due to cure (R). Based on the above division, literature [33] established a dynamic model of TB, discussed the global stability of disease-free equilibrium point and endemic equilibrium point, and analyzed the impact of smear-negative individuals on TB transmission. Considering that the proportion of smear-positive and smear-negative from S is not equal to that from E, and the patients after rehabilitation may become susceptible ones again, we improve the model in [33] to the following one
{˙S(t)=Π−(β1I1+β2I2)S−μS+λR,˙E(t)=(1−k)(β1I1+β2I2)S−(μ+m)E,˙I1(t)=p1k(β1I1+β2I2)S+q1mE−(γ1+δ1+μ)I1,˙I2(t)=p2k(β1I1+β2I2)S+q2mE−(γ2+δ2+μ)I2,˙R(t)=γ1I1+γ2I2−(μ+λ)R, | (2.1) |
where Π represents the input rate of the population, β1 and β2 represent the infection rates of smear positive and negative population to susceptible population respectively, μ represents the natural mortality rate, γ1 and γ2 are the cure rates of smear positive and negative population respectively, δ1 and δ2 are the disease-related death rates of smear positive and negative infected persons respectively, λ denotes the transfer rate of the immune or quarantined population that loses immunity or released from quarantine and becomes susceptible again, m represents the transfer rate from latent population to infected population, k represents the proportion of susceptible people who become infected directly after carrying TB bacteria (p1 and p2 represent the proportion of smear positive and negative respectively), q1 and q2 represent the proportion of smear positive and negative respectively in the people who enter the infected population from the latent population. All parameters in model (2.1) are positive except 0≤p1,p2,q1,q2≤1 and p1+p2=q1+q2=1. The initial values of model (2.1) are S(0)>0,E(0)>0,I1(0)>0,I2(0)>0,R(0)>0.
In practice, the infection rate of the bacteria and the transfer rate of the population among the compartments depend on the infection rate and transfer rate at each historical moment. That is, it has heritability and memory, which is completely ignored by the integer-order derivative adopted in model (2.1). It is worth noting that the fractional-order derivatives can better describe the genetic and memory characteristics of the process or materials, so the fractional-order compartmental model can be established to better describe the transmission process of mycobacterium TB bacteria in the population based on model (2.1). On the other hand, effective and feasible control measures need to be taken to control the spread of the bacteria in the population. Usually, media coverage for susceptible population to reduce the infection rate and treatment for infected population are the two most commonly used control measures, which are introduced into the fractional-order compartmental model to obtain the following fractional-order control model
{DαS(t)=Π−(β1I1+β2I2)S(1−c1u1)−μS+λR,DαE(t)=(1−k)(β1I1+β2I2)S(1−c1u1)−(μ+m)E,DαI1(t)=p1k(β1I1+β2I2)S(1−c1u1)+q1mE−(γ1+δ1+μ)I1−c2u2I11+α1I1,DαI2(t)=p2k(β1I1+β2I2)S(1−c1u1)+q2mE−(γ2+δ2+μ)I2−c3u3I21+α2I2,DαR(t)=γ1I1+γ2I2+c2u2I11+α1I1+c3u3I21+α2I2−(μ+λ)R, | (2.2) |
where α∈(0,1), c1,c2,c3∈[0,1] are the effective rate of control measures, u1,u2,u3∈[0,1] denote the intensity of control measures, α1,α2>0 are half-satiation constants. The reasons why we set u1(t),u2(t),u3(t) as control variables are summarized as follows:
(ⅰ) Media coverage for susceptible population (u1): after the media coverage to susceptible population, people can take scientific and effective protective measures to reduce the infection rate of the bacteria, which is described by (1−c1u1) in (2.2).
(ⅱ) Treatment to smear positive infected population (u2) and smear negative infected population (u3): due to the limited medical resources, the saturation treatment rate is used to characterize the treatment control measures.
Remark 1. In references [29,30,31,32], fractional-order TB models were established by considering different factors that influence the transmission of TB, where the fractional derivatives with nonsingular kernel were adopted, while the Caputo fractional derivative is used in model (2.2). Although the Caputo fractional derivative contains a singular kernel, this does not affect the main purpose of this paper: extend the adaptive control and neural networks approximation method widely used in the engineering fields to the biological field, and design control variables for the fractional-order TB model. It should be pointed out that the new and interesting definition of fractional derivative [27,28,42] proposed in recent years can effectively overcome some shortcomings of the traditional fractional derivative, but the stability theory related to these new definition is still relatively lacking, which may lead to difficulties in analyzing the stability of the controlled TB model.
The control objectives of this paper are as follows: The appropriate expressions of u1,u2,u3 in model (2.2) are designed so that the susceptible population S, smear positive population I1 and smear negative population I2 can track the control targets Sd, I1d and I2d with small errors, respectively. It is assumed that Sd,I1d,I2d and their fractional-order derivative with order α are all bounded and continuously differentiable. In order to achieve the control goal, the following lemma is given.
Lemma 4. Ω={(S,E,I1,I2,R)|S≥0,E≥0,I1≥0,I2≥0,R≥0,S+E+I1+I2+R≤Πμ} is a positive invariant set of models (2.1) and (2.2).
Proof. Lemma 4 can be proved by using similar methods of [43,Theorems 1 and 2], so omitted here.
In this section, the controllers are designed for the controlled model (2.2) through three steps, and then the stability of the control system is analyzed.
For convenience, the control errors are defined as follows
˜S=S−Sd,˜I1=I1−I1d,˜I2=I2−I2d. |
Step 1: The Lyapunov function V1 is constructed as
V1=12˜S2+12ξ−11˜θ21, | (3.1) |
where ξ1>0 is a design parameter, ˜θ1=ˆθ1−θ1, ˆθ1 is the estimate of θ1=‖W∗1‖2, W∗1 will be defined later. Based on lemma 1, the fractional-order derivative of the function V1 is
DαV1≤˜SDα˜S+ξ−11˜θ1Dα˜θ1≤˜S[Π−(β1I1+β2I2)S(1−c1u1)−μS+λR−DαSd]+ξ−11˜θ1Dα˜θ1=˜S[c1u1S(β1I1+β2I2)−DαSd+F1(X1)]+ξ−11˜θ1Dα˜θ1, |
where F1(X1)=Π−(β1I1+β2I2)S−μS+λR, X1=(S,I1,I2,R)T. From Lemma 4, we know that X1 belongs to a compact set, so continuous function F1(X1) can be approximated by radial basis function neural networks. According to Lemma 3, F1(X1) can be approximated as F1(X1)=W∗1Ψ1(X1)+ϵ1(X1), ϵ1(X1) is the approximation error and satisfies |ϵ1(X1)|<˜ϵ1. Using Young's inequality, we get
˜S(W∗1Ψ1(X1)+ϵ1(X1))≤|˜S|(‖W∗1‖‖Ψ1(X1)‖+˜ϵ1)≤12a21˜S2θ1ΨT1(X1)Ψ1(X1)+12a21+12˜S2+12˜ϵ21. |
Then, DαV1 can be written as
DαV1≤˜S[c1u1S(β1I1+β2I2)+12a21˜Sθ1ΨT1(X1)Ψ1(X1)+12˜S−DαSd]+12a21+12˜ϵ21+ξ−11˜θ1Dα˜θ1. | (3.2) |
The controller u1 and Dαˆθ1 are designed as
u1=1c1S(β1I1+β2I2)[−(ρ1+0.5)˜S−12a21˜Sˆθ1ΨT1(X1)Ψ1(X1)+DαSd], | (3.3) |
Dαˆθ1=ξ12a21˜S2ΨT1(X1)Ψ1(X1)−ˆθ1, | (3.4) |
where ρ1>0 is a design parameter. By substituting (3.3) and (3.4) into (3.2) and noting that ˜θ1=ˆθ1−θ1, we can obtain
DαV1≤−ρ1˜S2−ξ−11˜θ1ˆθ1+12a21+12˜ϵ21. |
Further, since −˜θ1ˆθ1=−˜θ21−˜θ1θ1≤−12˜θ21+12θ21, we get
DαV1≤−ρ1˜S2−12ξ−11˜θ21+Δ1, | (3.5) |
where Δ1=12ξ−11θ21+12a21+12˜ϵ21.
Step 2: The Lyapunov function V2 is constructed as
V2=12˜I21+12ξ−12˜θ22, | (3.6) |
where ξ2>0 is a design parameter, ˜θ2=ˆθ2−θ2, ˆθ2 is the estimate of θ2=‖W∗2‖2, W∗2 will be defined later. Applying Young's inequality and noticing that u1∈[0,1], we can obtain
˜I1p1k(β1I1+β2I2)S(1−u1)≤12˜I21+p21k2(β1I1+β2I2)2S2+12(1−u1)2≤12˜I21+p21k2(β1I1+β2I2)2S2+12. | (3.7) |
According to Lemma 1 and using (3.7), we can easily obtain
DαV2≤˜I1[−c1u2I11+α1I1−DαI1d+F2(X2)]+ξ−12˜θ2Dα˜θ2+12, |
where F2(X2)=12˜I1p21k2(β1I1+β2I2)2S2+q1mE−(γ1+δ1+μ)I1, X2=(S,E,I1,I2,I1d)T. The function F2(X2) can be approximated by radial basis function neural networks because it is continuous and X2 belongs to a compact set. According to Lemma 3, F2(X2) can be approximated as F2(X2)=W∗2Ψ2(X2)+ϵ2(X2), where |ϵ2(X2)|<˜ϵ2 is the approximation error. Using Young's inequality, we get
˜I1F2(X2)≤|˜I1|(‖W∗2‖‖Ψ2(X2)‖+˜ϵ2)≤12a22˜I21θ2ΨT2(X2)Ψ2(X2)+12a22+12˜I21+12˜ϵ22, |
where a2>0 is a constant. Then, we have
DαV2≤˜I1[−c1u2I21+α1I1−DαI1d+12a22˜I1θ2ΨT2(X2)Ψ2(X2)+12˜I1]+ξ−12˜θ2Dα˜θ2+12+12a22+12˜ϵ22. |
Design the controller u2 and adaptive law of ˆθ2 as
u2=1+α1I1c1I1[(ρ2+12)˜I1+12a22˜I1ˆθ2ΨT2(X2)Ψ2(X2)−DαI1d], | (3.8) |
Dαˆθ2=ξ22a22˜I21ΨT2(X2)Ψ2(X2)−ˆθ2, | (3.9) |
where ρ2>0 is a design parameter. Substituting (3.8), (3.9) into DαV2, and using ˜θ2=ˆθ2−θ2 and −˜θ2ˆθ2≤−12˜θ22+12θ22, we can get
DαV2≤−ρ2˜I21−12ξ−12˜θ22+Δ2, | (3.10) |
where Δ2=12ξ−12θ22+12+12a22+12˜ϵ22.
Step 3: Similar to step 2, construct the Lyapunov function V3 as
V3=12˜I22+12ξ−13˜θ23, | (3.11) |
where ξ3>0 is a designed parameter, ˜θ3=ˆθ3−θ3, ˆθ3 is the estimate of θ3=‖W∗3‖2, W∗3 is the weight vector. Similarly, we can design u3 and Dαˆθ3 as
u3=1+α2I2c2I2[(ρ3+12)˜I2+12a23˜I2ˆθ3ΨT3(X3)Ψ3(X3)−DαI2d], | (3.12) |
Dαˆθ3=ξ32a23˜I22ΨT3(X3)Ψ3(X3)−ˆθ3. | (3.13) |
where ρ3>0 is a design parameter, X3=(S,E,I1,I2,I2d)T. Using lemma 1, 3 and Young's inequality, we have
DαV3≤−ρ3˜I22−12ξ−13˜θ23+Δ3, | (3.14) |
where Δ3=12ξ−13θ23+12+12a23+12˜ϵ23.
Remark 2. In previous works [37,38,39], the vectors to be estimated by adaptive laws were constructed by using the model parameters, with the result that the initial values of adaptive laws have an important influence on the estimated value of the parameter vector. However, in this paper, the functions that contain model parameters are approximated by neural networks, and then the norms of weight vectors are estimated by adaptive laws, so it avoids the estimations to model parameters. Thus, the initial values dependency of adaptive laws are reduced.
Theorem 1. For control model (2.2), the designed controllers (3.3), (3.8), (3.12) and the adaptive laws (3.4), (3.9), (3.13) can ensure that S, I1 and I2 track the control target Sd, I1d and I2d with small errors, respectively. That is, the control goal can be achieved.
Proof. According to (3.1), (3.6) and (3.11), constructing Lyapunov function as V(t)=V1+V2+V3. Thus, based on (3.5), (3.10) and (3.14), we have
DαV≤−ρ1˜S2−ρ2˜I21−ρ3˜I22−123∑i=1ξ−1i˜θ2i+Δ≤−ρV+Δ, |
where ρ=min{2ρ1,2ρ2,2ρ3,1} and Δ=Δ1+Δ2+Δ3 are positive constants. Since there is 0≤Eα(−ρtα)<1 for any t≥0, according to [44,Lemma 5], it can be obtained
V(t)≤(V(0)−Δρ)Eα(−ρtα)+Δρ≤V(0)Eα(−ρtα)+Δρ. | (3.15) |
From the form of Lyapunov function V1,V2,V3, limt→∞Eα(−ρtα)=0 and (3.15), the control errors ˜S, ˜I1 and ˜I2 can enter into a small neighborhood of the origin by choosing appropriate parameters Δ and ρ. In other words, under the designed control variables, the control objective of TB model (2.2) can be achieved with small error.
Remark 3. According to Eqs (3.3), (3.8) and (3.12), the values of u1,u2,u3 will lie outside the interval [0,1] at some time, so, the value of u1,u2,u3 can be adjusted as follows
ui(t)={1,ui(t)>1,ui(t),0≤ui(t)≤1,0,ui(t)<0,i=1,2,3. |
The above adjustment is reasonable, which is mainly because the restriction on the range of control variables will reduce the control effect of the model to a certain extent, but will not affect the stability of the model, which is explained as follows by taking the "Step 1" in Subsection 3.1 as an example. In fact, if u1 is restrict into interval [0,1], then based on Lemma 4 and the assumptions for Sd in Subsection 2.2, we can conclude that the term c1u1S(β1I1+β2I2)+12˜Sθ1ΨT1(Z1)Ψ1(Z1)+12˜S−DαSd+ρ1˜S is bounded by a positive constant ζ1. Thus, (3.2) can be rewritten as
DαV1≤|˜S|ζ1+˜S[−12˜S˜θ1ΨT1(Z1)Ψ1(Z1)+12˜S˜θ1ΨT1(Z1)Ψ1(Z1)−ρ1˜S]+12a21+12˜ϵ21+ξ−11˜θ1Dα˜θ1. |
According to Lemma 2, we know that ˆθ1 is bounded and thus ˜θ1=ˆθ1−θ is bounded, which means that 12˜S˜θ1ΨT1(Z1)Ψ1(Z1) is bounded. Further, by using the boundedness of ˜S, DαV1 can be calculated as
DαV1≤−ρ1˜S2−12ξ−11˜θ21+ζ2, |
where ζ2>0 is a constant. Therefore, the stability of model (2.2) can be ensured by using the similar argument as in Subsection 3.2. Besides, this restriction also makes it possible to use the control measures in practice, because the control input outside the interval [0,1] is meaningless in practice.
Remark 4. Although there are many TB models which were built by using different fractional derivative definitions [29,30,31,32], these studies mainly focus on the existence and uniqueness of the solution, dynamic behavior, numerical methods and so on. However, the research results on the control problem of fractional TB model using adaptive control method are relatively scarce. In this paper, we successfully designed the expressions of control variables for the model by using adaptive control method and radial basis function neural networks, and analyzed the stability of the model under control variables.
This section is devoted to illustrating the theoretical results in previous section by some numerical examples.
The values of parameters in model (2.2) are chosen as: α=0.95, Π=1.4,β1=5.07394×10−3,β2=1.01479×10−3, μ=1/66,m=0.001,k=0.001, p1=0.35,p2=0.65, q1=0.4,q2=0.6, γ1=0.2038,γ2=0.4038, δ1=δ2=0.22, c1=0.7,c2=0.8,c3=0.7, α1=1,α2=1. The designed parameters in the controllers are selected as: ρ1=250,ρ2=ρ3=20, a1=a3=1,a2=2, ξ1=ξ2=ξ3=1. The initial values of model (2.2) and adaptive laws (3.4), (3.9) and (3.13) are (S(0),E(0),I1(0),I2(0),R(0))=(5000,80,40,40,20) and (ˆθ1(0),ˆθ2(0),ˆθ3(0))=(10−4,10−5,10−5), respectively. Finally, the desired values are
Iid=Ii(0)e−bit+ri∫tt−τ[pik(β1I1(w)+β2I2(w))×S(w)(1−c1u1(w))+qimE(w)]dw,(i=1,2)Sd=b3+(S(0)−b3)e−b4t, |
where τ=3,r1=0.3,r2=0.2, b1=b2=0.4,b3=1500,b4=0.3. The parameters in radial basis function neural networks are chosen as σi=(−30,−20,−10,0,10,20,30)T,κ=10.
Remark 5. Since the number of infected people at the current time depends on the cumulative number of infected people in the recent period and the initial values, the desired values I1d and I2d are selected as the above forms. Thus, r1 (or r2) represents the proportion of people who enter I1(or I2) from t−τ to t and have not recovered at t; I1(0)e−b1t(or I2(0)e−b2t) refers to the number of patients at the initial time who have not recovered at the current time t. For susceptible populations, our goal is to try to prevent them from contracting the bacteria, that is, to keep as many people as possible in S, which can be achieve by adjusting b3 and b4 in Sd. So, b3 represents the number of people who are still susceptible in the end, and b4 means the rate from S(0) to b3.
Example 1. Tracking performance of S, I1, I2 to Sd, I1d, I2d, respectively.
As presented in Section 2, the control objective of this study is to design the controllers such that S, I1, I2 can track desired values Sd, I1d, I2d with small errors, respectively. Thus, we consider the following two cases to verify the effectiveness of the controller and the tracking performance of S,I1,I2 to Sd,I1d,I2d, respectively: Case Ⅰ, τ=2,r1=0.3,r2=0.3, b1=b2=0.4,b3=1500,b4=0.3; Case Ⅱ, τ=2,r1=0.25,r2=0.2, b1=0.5,b2=0.3,b3=2000,b4=0.2. The main reason for choosing the above parameters is that we hope that most people can change their behavior and enter the compartment R through media coverage, and most people can recover after treatment. The simulation results are presented in Figure 1. From Figure 1, we see that the trajectories of S, I1, I2 can track the desired values Sd, I1d, I2d with small errors in both Case Ⅰ and Ⅱ, respectively, which indicate the effectiveness of the proposed controllers. In addition, compared with the case of without control, the number of infected people decreased significantly in both Case Ⅰ and Ⅱ.
Example 2. The effect of the fractional-order derivative order α.
The order α is an important parameter in fractional-order model (2.2), which can be adjusted based on the actual demand. In order to exhibit the effect of α in model (2.2), we choose α=1 (i.e., integer-order model), 0.9,0.8,0.7,0.6, and run the numerical simulations. The results are shown in Figure 2. It can be seen from Figure 2 that with the decrease of α, the convergence speed of the trajectory of S, I1 and I2 are decreasing, that is, the period of TB epidemic is prolonged. This shows that the fractional-order model is more flexible than the integer-order model and has the ability to meet the needs of different situations.
Example 3. The effect of different control combinations.
The purpose of setting control variables in model (2.2) is to prevent as many people as possible from contracting TB. In order to better understand the role of each control measure, we consider the following five cases: Case Ⅰ, u1,u2,u3=0; Case Ⅱ, u1,u3=0,u2≠0; Case Ⅲ, u1,u2≠0,u3=0; Case Ⅳ, u1=0,u2,u3≠0; Case Ⅴ, u1,u2,u3≠0, where ui=0 denotes ui is in absence, and ui≠0 denotes that control measure ui is adopted (this does not mean that ui has never been equal to zero).
The simulation results are presented in Figure 3. It follows from Figure 3(a) that, compared with the cases without media coverage (i.e., Cases Ⅰ, Ⅱ, Ⅳ), there are more people who can remain in compartment S when media coverage measures are taken for the susceptible people (i.e., Cases Ⅲ, Ⅳ), which shows that media coverage (u1) is an effective control measure to curb the spread of TB. Figure 3(b) shows the exposed individuals under five cases, which also indicate the significance of media coverage. From Figure 3(c), (d), we know that, when only control measures u2 or u3 (i.e., treatment to the smear positive and negative individuals) are adopted, their effect on reducing the number of TB infections is extremely limited. However, when media coverage and treatment measures are adopted at the same time, the infectious scale of TB epidemic can be greatly reduced. Therefore, we can conclude that comprehensive control measures can be more effective than single ones to prevent the large-scale spread of the bacteria when the TB outbreak occurs. Figure 3(e) presents the trajectories of R. The trajectories of control variables u1,u2 and u3 in Case Ⅴ are given in Figure 3(f).
The adaptive neural networks control problem has been investigated for the fractional-order TB model. First, based on the spreading mechanism of TB, a fractional-order compartment model of TB has been established. Then, by using the universal approximation principle of neural networks and constructing Lyapunov functions, adaptive controllers have been designed, which can ensure the control objective can be realized. Finally, some numerical simulations have been presented to illustrate the effectiveness of the proposed controllers. It should be noted that some interesting and meaningful new definitions for fractional-order derivative, including Caputo-Fabrizio, Atangana, Atangana-Baleanu, Losada-Nieto and so on [28,42], have emerged in recent years, which, compared with Caputo derivative adopted in this paper, have nonsingular kernel. Thus, we need to consider the modeling and control problem of TB under these new definition in the future. Meanwhile, some comparisons should be made to explore the advantage and disadvantage of Caputo derivative and other new definitions, which is also one of our future works.
This work was supported by the Natural Science Foundation of China (12201330).
The author declares there is no conflicts of interest.
[1] |
L. N. Reynard, M. J. Barter, Osteoarthritis year in review 2019: Genetics, genomics and epigenetics, Osteoarthritis Cartilage, 28 (2020), 275–284. http://doi.org/10.1016/j.joca.2019.11.010 doi: 10.1016/j.joca.2019.11.010
![]() |
[2] |
D. J. Hunter, S. Bierma-Zeinstra, Osteoarthritis, Lancet, 393 (2019), 1745–1759. http://doi.org/10.1016/S0140-6736(19)30417-9 doi: 10.1016/S0140-6736(19)30417-9
![]() |
[3] |
J. W. Bijlsma, F. Berenbaum, F. P. Lafeber, Osteoarthritis: an update with relevance for clinical practice, Lancet, 377 (2011), 2115–2126. http://doi.org/10.1016/S0140-6736(11)60243-2 doi: 10.1016/S0140-6736(11)60243-2
![]() |
[4] |
Q. Liu, S. Wang, J. Lin, Y. Zhang, The burden for knee osteoarthritis among Chinese elderly: estimates from a nationally representative study, Osteoarthritis Cartilage, 26 (2018), 1636–1642. http://doi.org/10.1016/j.joca.2018.07.019 doi: 10.1016/j.joca.2018.07.019
![]() |
[5] |
S. Safiri, A. Kolahi, E. Smith, C. Hill, D. Bettampadi, M. A. Mansournia, et al., Global, regional and national burden of osteoarthritis 1990–2017: A systematic analysis of the Global Burden of Disease Study 2017, Ann. Rheum. Dis., 79 (2020), 819–828. http://doi.org/10.1136/annrheumdis-2019-216515 doi: 10.1136/annrheumdis-2019-216515
![]() |
[6] |
J. C. Stoddart, O. Dandridge, A. Garner, J. Cobb, R. J. van Arkel, The compartmental distribution of knee osteoarthritis-a systematic review and meta-analysis, Osteoarthritis Cartilage, 29 (2021), 445–455. http://doi.org/10.1016/j.joca.2020.10.011 doi: 10.1016/j.joca.2020.10.011
![]() |
[7] |
C. Jacquet, F. Gulagaci, A. Schmidt, A. Pendse, S. Parratte, J. Argenson, et al., Opening wedge high tibial osteotomy allows better outcomes than unicompartmental knee arthroplasty in patients expecting to return to impact sports, Knee Surg. Sports Traumatol. Arthroscopy, 28 (2020), 3849–3857. http://doi.org/10.1007/s00167-020-05857-1 doi: 10.1007/s00167-020-05857-1
![]() |
[8] |
A. Agarwalla, D. R. Christian, J. N. Liu, G. H. Garcia, M. L. Redondo, A. B. Yanke, et al., Return to work following isolated opening wedge high tibial osteotomy, Cartilage, 12 (2021), 468–474. http://doi.org/10.1177/1947603519852417 doi: 10.1177/1947603519852417
![]() |
[9] |
M. Loia, S. Vanni, F. Rosso, D. Bonasia, M. Bruzzone, F. Dettoni, et al., High tibial osteotomy in varus knees: Indications and limits, Joints, 4 (2016), 98–110. http://doi.org/10.11138/jts/2016.4.2.098 doi: 10.11138/jts/2016.4.2.098
![]() |
[10] |
B. Zampogna, S. Vasta, R. Papalia, Patient evaluation and indications for osteotomy around the knee, Clin. Sports Med., 38 (2019), 305–315. http://doi.org/10.1016/j.csm.2019.02.011 doi: 10.1016/j.csm.2019.02.011
![]() |
[11] |
D. H. Kim, S. C. Kim, J. S. Yoon, Y. S. Lee, Are there harmful effects of preoperative mild lateral or patellofemoral degeneration on the outcomes of open wedge high tibial osteotomy for medial compartmental osteoarthritis?, Orthop. J. Sports Med., 8 (2020). http://doi.org/10.1177/2325967120927481 doi: 10.1177/2325967120927481
![]() |
[12] |
C. Jin, E. Song, A. Santoso, P. S. Ingale, I. Choi, J. Seon, Survival and risk factor analysis of medial open wedge high tibial osteotomy for unicompartment knee osteoarthritis, J. Arthroscopic Relat. Surg., 36 (2020), 535–543. http://doi.org/10.1016/j.arthro.2019.08.040 doi: 10.1016/j.arthro.2019.08.040
![]() |
[13] |
R. S. Khakha, H. R. B. A. Razak, K. Kley, R. van Heerwaarden, A. J. Wilson, Role of high tibial osteotomy in medial compartment osteoarthritis of the knee: Indications, surgical technique and outcomes, J. Clin. Orthop. Trauma, 23 (2021), 101618. http://doi.org/10.1016/j.jcot.2021.101618 doi: 10.1016/j.jcot.2021.101618
![]() |
[14] |
J. L. Martay, A. J. Palmer, N. K. Bangerter, S. Clare, A. P. Monk, C. P. Brown, et al., A preliminary modeling investigation into the safe correction zone for high tibial osteotomy, Knee, 25 (2018), 286–295. http://doi.org/10.1016/j.knee.2017.12.006 doi: 10.1016/j.knee.2017.12.006
![]() |
[15] | Z. Trad, A. Barkaoui, M. Chafra, J. M. R. Tavares, Finite element analysis of the effect of high tibial osteotomy correction angle on articular cartilage loading, in Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine, 232 (2018), 553–564. http://doi.org/10.1177/0954411918770706 |
[16] |
Y. Dong, G. Hu, Y. Dong, Y. Hu, Q. Xu, The effect of meniscal tears and resultant partial meniscectomies on the knee contact stresses: a finite element analysis, Comput. Methods Biomech. Biomed. Eng., 17 (2013), 1452–1463. http://doi.org/10.1080/10255842.2012.753063 doi: 10.1080/10255842.2012.753063
![]() |
[17] |
L. Li, L. Yang, K. Zhang, L. Zhu, X. Wang, Q. Jiang, Three-dimensional finite-element analysis of aggravating medial meniscus tears on knee osteoarthritis, J. Orthop. Transl., 20 (2020), 47–55. http://doi.org/10.1016/j.jot.2019.06.007 doi: 10.1016/j.jot.2019.06.007
![]() |
[18] |
P. Jiang, J. Cui, Z. Chen, Z. Dai, Y. Zhang, G. Yi, Biomechanical study of medial meniscus after posterior horn injury: a finite element analysis, Comput. Methods Biomech. Biomed. Eng., 23 (2020), 127–137. http://doi.org/10.1080/10255842.2019.1702167 doi: 10.1080/10255842.2019.1702167
![]() |
[19] |
B. P. Nguyen, C. Chui, S. Ong, S. Chang, An efficient compression scheme for 4-D medical images using hierarchical vector quantization and motion compensation, Comput. Biol. Med., 41 (2011), 843–856. http://doi.org/10.1016/j.compbiomed.2011.07.003 doi: 10.1016/j.compbiomed.2011.07.003
![]() |
[20] |
Q. H. Nguyen, B. P. Nguyen, M. T. Nguyen, M. C. H. Chua, T. T. T. Do, N. Nghiem, Bone age assessment and sex determination using transfer learning, Expert Syst. Appl., 200 (2022), 116926. http://doi.org/10.1016/j.eswa.2022.116926 doi: 10.1016/j.eswa.2022.116926
![]() |
[21] |
X. Chen, B. P. Nguyen, C. Chui, S. Ong, An automatic framework for multi-label brain tumor segmentation based on kernel sparse representation, Acta Polytech. Hung., 14 (2017). http://doi.org/10.12700/aph.14.1.2017.1.3 doi: 10.12700/aph.14.1.2017.1.3
![]() |
[22] |
S. Lu, B. Yang, Y. Xiao, S. Liu, M. Liu, L. Yin, et al., Iterative reconstruction of low-dose CT based on differential sparse, Biomed. Signal Process. Control, 79 (2023), 104204. http://doi.org/10.1016/j.bspc.2022.104204 doi: 10.1016/j.bspc.2022.104204
![]() |
[23] |
L. Cai, W. Tay, B. P. Nguyen, C. Chui, S. Ong, Automatic transfer function design for medical visualization using visibility distributions and projective color mapping, Comput. Med. Imaging Graphics, 37 (2013), 450–458. http://doi.org/10.1016/j.compmedimag.2013.08.008 doi: 10.1016/j.compmedimag.2013.08.008
![]() |
[24] |
Y. Xu, F. Zhang, W. Zhai, S. Cheng, J. Li, Y. Wang, Unraveling of advances in 3D-printed polymer-based bone scaffolds, Polymers, 14 (2022), 566. http://doi.org/10.3390/polym14030566 doi: 10.3390/polym14030566
![]() |
[25] |
S. Su, J. He, C. Wang, F. Gao, D. Zhong, P. Lei, A new dressing system reduces the number of dressing changes in the primary total knee arthroplasty: A randomized controlled trial, Front. Surg., 9 (2022). http://doi.org/10.3389/fsurg.2022.800850 doi: 10.3389/fsurg.2022.800850
![]() |
[26] |
Y. Miao, X. Wang, S. Wang, R. Li, Adaptive switching control based on dynamic zero-moment point for versatile hip exoskeleton under hybrid locomotion, IEEE Trans. Ind. Electron., 70 (2023), 11443–11452. http://doi.org/10.1109/TIE.2022.3229343 doi: 10.1109/TIE.2022.3229343
![]() |
[27] |
Y. Chen, J. Xiang, Z. Wang, Y. Xiao, D. Zhang, X. Chen, et al., Associations of bone mineral density with lean mass, fat mass, and dietary patterns in postmenopausal chinese women: A 2-year prospective study, PLOS ONE, 10 (2015), 0137097. http://doi.org/10.1371/journal.pone.0137097 doi: 10.1371/journal.pone.0137097
![]() |
[28] |
N. M. Luís, R. Varatojo, Radiological assessment of lower limb alignment, Efort Open Rev., 6 (2021), 487–494. http://doi.org/10.1302/2058-5241.6.210015 doi: 10.1302/2058-5241.6.210015
![]() |
[29] |
D. M. Allen, L. Li, M. D. Crema, M. D. Marra, A. Guermazi, B. T. Wyman, et al., The relationship between meniscal tears and meniscal position, Ther. Adv. Musculoskeletal Dis., 2 (2010), 315–323. http://doi.org/10.1177/1759720X10383198 doi: 10.1177/1759720X10383198
![]() |
[30] |
M. E. Wells, J. P. Scanaliato, J. C. Dunn, E. J. Garcia, Meniscal injuries: Mechanism, classification, Sports Med. Arthroscopy Rev., 29 (2021), 154–157. http://doi.org/10.1097/JSA.0000000000000311 doi: 10.1097/JSA.0000000000000311
![]() |
[31] |
P. Beaufils, R. Becker, S. Kopf, O. Matthieu, N. Pujol, The knee meniscus: management of traumatic tears and degenerative lesions, Efort Open Rev., 2 (2017), 195–203. http://doi.org/10.1302/2058-5241.2.160056 doi: 10.1302/2058-5241.2.160056
![]() |
[32] |
K. Ding, W. Yang, H. Wang, S. Zhan, P. Hu, J. Bai, et al., Finite element analysis of biomechanical effects of residual varus/valgus malunion after femoral fracture on knee joint, Int. Orthop., 45 (2021), 1827–1835. http://doi.org/10.1007/s00264-021-05039-9 doi: 10.1007/s00264-021-05039-9
![]() |
[33] |
H. R. C. Bao, D. Zhu, G. S. Gu, H. Gong, The effect of complete radial lateral meniscus posterior root tear on the knee contact mechanics: a finite element analysis, J. Orthop. Sci., 18 (2013), 256–263. http://doi.org/10.1007/s00776-012-0334-5 doi: 10.1007/s00776-012-0334-5
![]() |
[34] |
T. Kozaki, D. Fukui, E. Yamamoto, D. Nishiyama, M. Yamanaka, A. Murata, et al., Medial meniscus extrusion and varus tilt of joint line convergence angle increase stress in the medial compartment of the knee joint in the knee extension position-finite element analysis-, J. Exp. Orthop., 9 (2022). http://doi.org/10.1186/s40634-022-00490-y doi: 10.1186/s40634-022-00490-y
![]() |
[35] |
C. Pan, X. Wang, L. Ding, X. Zhu, W. Xu, L. Huang, The best position of bone grafts in the medial open-wedge high tibial osteotomy: A finite element analysis, Comput. Methods Programs Biomed., 228 (2023), 107253. http://doi.org/10.1016/j.cmpb.2022.107253 doi: 10.1016/j.cmpb.2022.107253
![]() |
[36] |
A. E. Kedgley, T. Saw, N. A. Segal, U. N. Hansen, A. M. J. Bull, S. D. Masouros, Predicting meniscal tear stability across knee-joint flexion using finite-element analysis, Knee Surg. Sports Traumatol. Arthroscopy, 27 (2019), 206–214. http://doi.org/10.1007/s00167-018-5090-4 doi: 10.1007/s00167-018-5090-4
![]() |
[37] |
Y. Wu, X. Jin, X. Zhao, Y. Wang, H. Bai, B. Lu, et al., Computer-aided design of distal femoral osteotomy for the valgus knee and effect of correction angle on joint loading by finite element analysis, Orthop. Surg., 14 (2022), 2904–2913. http://doi.org/10.1111/os.13440 doi: 10.1111/os.13440
![]() |
[38] |
B. Ollivier, P. Berger, C. Depuydt, H. Vandenneucker, Good long-term survival and patient-reported outcomes after high tibial osteotomy for medial compartment osteoarthritis, Knee Surg. Sports Traumatol. Arthroscopy, 29 (2020), 3569–3584. http://doi.org/10.1007/s00167-020-06262-4 doi: 10.1007/s00167-020-06262-4
![]() |
[39] |
M. Englund, A. Guermazi, S. L. Lohmander, The role of the meniscus in knee osteoarthritis: A cause or consequence?, Radiol. Clin., 47 (2009), 703–712. http://doi.org/10.1016/j.rcl.2009.03.003 doi: 10.1016/j.rcl.2009.03.003
![]() |
[40] |
J. B. Guimaraes, M. C. Nevitt, C. E. McCulloch, B. J. Schwaiger, A. S. Gersing, L. Facchetti, et al., Association of weight change with progression of meniscal intrasubstance degeneration over 48 months: Data from the Osteoarthritis Initiative, Eur. Radiol., 28 (2017), 953–962. http://doi.org/10.1007/s00330-017-5054-y doi: 10.1007/s00330-017-5054-y
![]() |
[41] |
M. A. Wesdorp, S. M. Eijgenraam, D. E. Meuffels, S. M. A. Bierma-Zeinstra, G. Kleinrensink, Y. M. Bastiaansen-Jenniskens, et al., Traumatic meniscal tears are associated with meniscal degeneration, Am. J. Sports Med., 48 (2020), 2345–2352. http://doi.org/10.1177/0363546520934766 doi: 10.1177/0363546520934766
![]() |
[42] |
R. J. Ward, J. B. Driban, J. W. MacKay, T. E. McAlindon, B. Lu, C. B. Eaton, et al., Meniscal degeneration is prognostic of destabilzing meniscal tear and accelerated knee osteoarthritis: Data from the osteoarthritis initiative, J. Orthop. Res., 41 (2023), 2418–2423. http://doi.org/10.1002/jor.25575 doi: 10.1002/jor.25575
![]() |
[43] |
J. Moore, L. Mychaltchouk, F. Lavoie, Applicability of a modified angular correction measurement method for open-wedge high tibial osteotomy, Knee Surg. Sports Traumatol. Arthroscopy, 25 (2016), 846–852. http://doi.org/10.1007/s00167-015-3954-4 doi: 10.1007/s00167-015-3954-4
![]() |
[44] |
K. Briem, D. K. Ramsey, W. Newcomb, K. S. Rudolph, L. Snyder-Mackler, Effects of the amount of valgus correction for medial compartment knee osteoarthritis on clinical outcome, knee kinetics and muscle co‐contraction after opening wedge high tibial osteotomy, J. Orthop. Res., 25 (2007), 311–318. http://doi.org/10.1002/jor.20326 doi: 10.1002/jor.20326
![]() |
[45] |
Y. Yin, S. Li, R. Zhang, J. Guo, Z. Hou, Y. Zhang, What is the relationship between the "fujisawa point" and postoperative knee valgus angle? A theoretical, computer-based study, Knee, 27 (2020), 183–191. http://doi.org/10.1016/j.knee.2019.10.018 doi: 10.1016/j.knee.2019.10.018
![]() |
[46] |
X. Liu, Z. Chen, Y. Gao, J. zhang, Z. Jin, High tibial osteotomy: Review of techniques and biomechanics, J. Healthcare Eng., 2019 (2019), 8363128. http://doi.org/10.1155/2019/8363128 doi: 10.1155/2019/8363128
![]() |
1. | Ariyanatchi M., Vijayalakshmi G.M., Nonlinear Robust Adaptive Sliding Mode Control Strategies Involve a Fractional Ordered Approach to Reducing Dengue Vectors, 2024, 14, 26667207, 100406, 10.1016/j.rico.2024.100406 |