study on a fractional model of viscoelasticity of human cranial bone
DESCRIPTION
Study on a fractional model of viscoelasticity of human cranial bone. Jiaguo Liu, Mingyu Xu School of Mathematics, Shandong University, Jinan, 250100, P.R. China. 1. Introduction. Bone is anisotropic and viscoelastic - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/1.jpg)
Study on a fractional model of viscoelasticity of human cranial bone
Jiaguo Liu, Mingyu XuSchool of Mathematics,
Shandong University, Jinan, 250100, P.R. China.
![Page 2: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/2.jpg)
2
1. IntroductionBone is anisotropic and viscoelasticStudy on mechanical behavior of Cranial
bone is the basic work of research on craniocerebral injury.
The researches on dynamic behavior of bones are important in guiding orthopaedics diseases, cure of bone injure, substitutive materials and healing study.
![Page 3: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/3.jpg)
3
Cranial Bones: eight bones
![Page 4: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/4.jpg)
4
• Zhu et al’ study on the behavior of cranial bone by classical St.Venant model• Classical Maxwell and Zener model’s fractional order generalizations
![Page 5: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/5.jpg)
5
2. Fractional generalization of classical St.Venant model2.1 The classical (integer order) St.Venant
model is shown as follows
Its constitutive equation is
(1))()()()( 1
11 tEtEtt rr
![Page 6: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/6.jpg)
6
)(t where and denote the stress and strain, is the elastic coefficients, and is the viscosity.
(2) Obviously, (3)
)(t
,21
21
EE
EEE
,
21 EEr
2E
d
21, EE
.1 dr EE
)()()()( 1 tEtEtt rr
![Page 7: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/7.jpg)
7
2.2 Fractional generalization of St.Venant model
Riemann-Liouville fractional operators:
(4)
(5)
,0Re,)(:)(
0
1
0
dft
tfDt
t
.0,0,)(:)( 00 qnqtfDdt
dtfD nq
tn
nqt
![Page 8: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/8.jpg)
8
Let ,
(6) .
(7)
is equivalent with Eq. (1).Integrals from to give
, (8)
, (9)
where and are initial values of and respectively.
)()()]([ 1 ttt r
)()()]([ 11 tEtEt r
01
0
110 )()( ttDD trt
))(()( 011
0
110 tEtDED trt
)]([)]([ tt
0 0 )(t )(t
t0
![Page 9: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/9.jpg)
9
Substituting by in (8),
and by in (9), we obtain
results the fractional St.Venant model:
(10)
,)()(~
00 ttD q
t
q
r
,))(()(~010
tEtDE tr
)(0 tD qt
qr )(1
01 tDtr
)(10
1 tDtr )(0 tD qt
qr
.)1,0( q
~~
))(()()()( 01000 tEtDEttD tr
qt
qr
![Page 10: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/10.jpg)
10
3. Solutions of fractional St.Venant model3.1 Relaxation and creep function of fractional
St.Venant model
Laplace transform of (10) gives
. (11)
Let , where is the Heaviside unit step function, from (11) we obtain
(12)
1011
10 )(ˆ)(ˆ)(ˆ)(ˆ pEpEppEpppp r
qqr
)()( 0 tt
.1
)(ˆ11
0qq
r
r
p
pEpp
)(t
![Page 11: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/11.jpg)
11
The discrete inverse Laplace transform of (12) give the relaxation modulus of fractional St.Venant model
, (13)
where is the H-Fox function.
0 0
11
0
01 )1()1()(k k
qkqkr
kqkqkr
k pEpLtG
qk
rk k
kqk
r
k t
qkE
t
qkE
0 01 )1(
)1(
)1(
)1(
q
qrr
q
qr
tH
t
q
EtH
q
E1
,0
1,;)1
,0(
1,12,1
1,0
1,0;)1
,0(
1,12,1
1
)(1,12,1 xH
![Page 12: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/12.jpg)
12
The deduction uses the following properties of the H-Fox function:
(14)
where is also called Maitland’s generalized hypergeometric function.
(15)
,)(;),(
),(
)1();1,0(
),1(
)(!
)()(
,
,11,
0
1
1
z
Bb
Aa
Bb
AazH
nBbn
nAaz
ppqp
pppqp
nq
jjj
p
jjj
n
.),(
),(
),(
),(1 ,,
,,
ppKnmqp
ppnmqp Kb
KazH
b
azH
K
)( zqp
![Page 13: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/13.jpg)
13
In a similar way, the creep compliance of fractional St.Venant model can be obtained
(16)
where
,
11)(
1,0
1,;)1
,0(
1,12,1
1
1,0
1,0;)1
,0(
1,12,1
1
qr
q
rr
tWH
t
E
tWH
EtJ
.
1
1
E
EW
![Page 14: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/14.jpg)
14
When , (13) and (16) reduce to the relaxation modulus and creep compliance of classical (integer order) St.Venant model
(17)
(18)
1q
1,01,1;)1,0(
1,12,1
1,01,0;)1,0(
1,12,11)(
rrr
tH
tE
tHEtG
r
tEEE exp1
1,01,1;)1,0(
1
1,12,1
1
1,01,0;)1,0(
1
1,12,1
1
11)(
rrr E
EtH
t
EE
EtH
EtJ
d
tEE exp111
21
![Page 15: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/15.jpg)
15
3.2 The fractional relaxation and creep functions under quasi-static loading
The loading processes of relaxation and creep tests are, respectively,
(19)
(20)
where and are constant strain and stress rates, respectively, and and are constants.
,
0
1
1
1 tt
ttAtt
,
0
1
1
1 tt
ttBtt
A B
11
![Page 16: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/16.jpg)
16
By Boltzmann superposition principle
from (13) and (16) we obtain the relaxation and creep response functions
,0
dtGtt
,0
dtJtt
![Page 17: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/17.jpg)
17
1
1,0
1,1;)1
,0(
11,12,1
11
1,0
1,1;)1
,0(
1,12,1
1,0
1,1;)1
,0(
11,12,1
11
1,0
1,1;)1
,0(
1,12,1
1
1
1,0
1,1;)1
,0(
1,12,1
1,0
1,1;)1
,0(
1,12,1
1
0
,
tt
ttH
tt
q
ttAEtH
t
q
AEt
ttH
q
ttAEtH
q
tAE
tt
tH
t
q
AEttH
q
tAE
t
q
qrr
q
qrr
q
qr
q
qr
q
qrr
q
qr
(21)
![Page 18: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/18.jpg)
18
1
1,0
1,1;)1
,0(
11,12,1
1
1
1
1,0
1,1;)1
,0(
1,12,1
1
1,0
1,1;)1
,0(
11,12,1
1
1
1,0
1,1;)1
,0(
1,12,1
1
1
1,0
1,1;)1
,0(
1,12,1
1
1,0
1,1;)1
,0(
1,12,1
1
0
,
tt
ttWH
tt
E
ttBWtH
t
E
Bt
ttWH
E
ttBWtH
E
Bt
tt
WtH
t
E
BtWtH
E
Bt
t
qr
q
rq
r
q
r
rr
qr
q
rr
(22)
![Page 19: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/19.jpg)
19
When , the relaxation and creep
functions of classical St.Venant model are
(23)
(24)
1q
1
111
11
,exp1exp
0,exp1
ttttEEAAEt
tttEEAAEtt
rtr
rr
1
1
2
1
12
,exp1exp
0,exp1
ttttE
B
E
Bt
tttE
B
E
Bt
t
dd
d
d
d
![Page 20: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/20.jpg)
20
4. Data fitting and comparison
The relaxation and creep functions (21) and
(22) are fitted with the experimental data
from Zhu et al’s, and we take parameters A,
B, E1, E2, t1 , τr, τd the same values as Zhu
et al’s.
![Page 21: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/21.jpg)
21
(***): Relaxation experimental data from [5]; (---): the relaxation function (23) of the standard St.Venant model; (—):the relaxation function (21) of the fractional St.Venant model. Here, q=0.965, μ=0.96.
![Page 22: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/22.jpg)
22
(***): Creep experimental data from [5]; (---): the creep function (24) of the standard St.Venant model; (—):the creep function (22) of the fractional St.Venant model. Here, q=0.5, μ=0.47.
![Page 23: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/23.jpg)
23
It is shown that, the fractional St.Venant model is more efficient than the standard St.Venant model with integer order in describing the stress-strain constitutive relations for the viscoelasticity of human cranial bone.
![Page 24: Study on a fractional model of viscoelasticity of human cranial bone](https://reader035.vdocument.in/reader035/viewer/2022062723/56813bac550346895da4deec/html5/thumbnails/24.jpg)
24
5. Conclusion
The fractional St.Venant model is more efficient than the classical model in describing the stress-strain constitutive relations for the viscoelasticity of human cranial bone.
It is efficient that applying fractional calculus method to describe constitutive relations of biological viscoelastic materials.