study on a fractional model of viscoelasticity of human cranial bone
Post on 06-Jan-2016
27 Views
Preview:
DESCRIPTION
TRANSCRIPT
Study on a fractional model of viscoelasticity of human cranial bone
Jiaguo Liu, Mingyu XuSchool of Mathematics,
Shandong University, Jinan, 250100, P.R. China.
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.
3
Cranial Bones: eight bones
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
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
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
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
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
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
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
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
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
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
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
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
16
By Boltzmann superposition principle
from (13) and (16) we obtain the relaxation and creep response functions
,0
dtGtt
,0
dtJtt
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)
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)
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
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.
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.
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.
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.
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.
25
Thank you very much!
Jiaguo Liuliujiaguo@sdu.edu.cn
top related