BAÜ Fen Bil. Ens. Derg. (1999).1 (1) FİBERLERLE TAKVİYELİ VİSKOELASTİK BİR MALZEMEYE
HOMOJEN DEFORMASYONUN UYGULANMASI Nurşen ÖNTÜRK
Trakya Üniversitesi, Çorlu Mühendislik Fakültesi, Makina Müh. Böl. Çorlu - Tekirdağ
ÖZET
Bu çalışmada, fiberlerle takviye edilmiş viskoelastik bir kompozit ortama ait gerilme bünye denklemine homojen deformasyon uygulanmıştır. Bunun için önce deformasyon geometrisi incelenmiş ve gerekli kinematik bağıntılar elde edilmiştir. Kompozit ortamda ve A
B
vektörleriyle gösterilen doğrusal fiberler kullanılmıştır.Daha sonra homojen deformasyonda kullanılacak olan büyüklükler ve kompozit ortamdaki kısıtlamaları veren ifadeler elde edilmiştir.
Bundan sonra ortaya çıkan kinematik kısıtlamalar altında homojen deformasyona maruz kalan bir kompozit ortamın gerilme analizi yapılmıştır. Elde edilen gerilme bileşenleri daha sonra sıkışmazlık ve uzamazlık kısıtlamalarından dolayı ortamda oluşan
reaksiyon gerilmelerinin bulunmasında kullanılmıştır.
p T T
,
a,
b Anahtar Kelimeler: Fiber vektörü, kompozit ortam, viskoelastik, gerilme, homojen deformasyon.A HOMOGENEOUS DEFORMATION IMPOSED ON A VISCOELASTIC CONTINUUM REINFORCED BY FIBRES ABSTRACT
In this study; the homogeneous deformations were applied to the stress constitutive equations for the viskoelastic composite medium. First, the deformation geometry was studied and the relevant kinematic relations were obtained. In the composite medium, straight fibres denoted by A and B were used.
Under the defined restrictions the stress analysis was made for the composite medium subjected to the homogeneous deformations. For the stress analysis purposes the medium was assumed to be homogeneous.
The stress components of this composite medium were used for obtaining the arbitrary values , the reaction stresses, representing incompressibility and inextensibilities, respectively.
p T T
,
a,
bKey words: Fibre vector, compozite continuum, viscoelastic, stress, homogeneous deformation.
1. GİRİŞ
Kompozit malzeme denildiği zaman üstün özelliklere sahip bir malzeme oluşturmak üzere mikroskobik düzeyde iki veya daha fazla malzemenin birleştirilmiş hali anlaşılır. Genel olarak kompozit malzemeler:
- Fiber takviyeli kompozitler, - Tabakalı kompozitler, - Parçacıklı kompozitler
olmak üzere üç türlü oluşturulurlar. Burada fiber takviyeli viskoelastik bir kompozit ortam ele alınmıştır. Böyle bir ortama ait gerilme bünye denklemi, daha önceki bir çalışmada araştırmacı tarafından elde edilmiştir [1]. Bu çalışmada; Spencer’in elastik bir kompozit ortamın homojen deformasyonları için incelediği örnek problem [2], araştırmacı tarafından viskoelastik bir kompozit ortama uyarlanmıştır.
2. DEFORMASYON GEOMETRİSİ VE KİNEMATİK BAĞLAR
Kelvin-Voigt türünde viskoelastik, fiberli bir malzeme için gerilme bünye denklemi:
t
k l= −
p
δ
k l+
T a a
a k l+
T b b
b k l+
τ
k l
=
σ
k l+
τ
k l=
σ
k l+
Eτ
k l+
Dτ
k l (2.1) şeklinde elde edilmiştir[1]. BuradaA
veB
fiber ailelerinin uzamadığı ve ortamın sıkıştıramadığı kabul edilmiştir. (2.1)denkleminde, E
τ
k l ile gösterilen elastik gerilme aşağıda görüldüğü gibi nonlineer bir şekilde elde edilmiştir [1].I
I
∂
∂
=
Σ
+
1Σ
bb c
k k lra c b
Φ
klc a a
kδ
=
+
+ 1 kl c a11 k b k 12 2λ
2
2
δ
kl+
12+
γ
X3 X[
Φ in ,Φ 0]
[
φ in ,φ 0]
X1 =(
)
E klτ
I
c
klI
c c
kr rlI
a c a a c a
k lr r l kr r∂
∂
∂
∂
∂
∂
−−
− −+
−+
−2
2
2
1 2 1 2 1 1 7 1 1Σ
Σ
(
)
(
)
+
2
2 −+
−+
2
8 1 1 9λ ∂
∂
lr r l kr rλ ∂
b∂
k l lI
b
b c b
I
a b
a b
Σ
Σ
Φ
cos
(
)
+
λ ∂
−+
−+
−+
−+
∂
λ ∂
∂
b r l kr r l lr r l kr r b k lI
a c b b c a b c a
I
b b
Σ
Σ
10 1 1 1 1 2 62
2
cos
D kl (2.2)τ
ile verilen disipatif gerilme ise, deformasyon hızları tansörüne göre lineer olarak aşağıdaki gibi elde edilmiştir [1].d
+
k[
D klτ
λ
v 1 l+
d b b d
1 k l+
21λ
bcos
2
Φ(
a b
k l+
b a
k l)
]
d
mm[
+2 + + 12 2 + 2 12µ
v kld cδ
a clλ
b k lb b fλ
bcos Φ(
a bk l +b ak l)
]
a d
i ija
j[
+ d1 2b kl +c a al +d b k lb b + f b 22 2 21 3 2λ δ
λ
λ
cos Φ(
a b b a
k l+
k l)
]
b d b
i ij j[
+
d
12λ
bcos
Φ+2
f
12λ
bcos
2
Φa a
k l+2
f
21λ
3bcos
2
Φb bk l Φ(
2
22
22 2λ β
bcos
)(
a bk l +b ak l)
]
a d bi ij j ( 2.3) Bu çalışmadaki kompozit ortamda ve A B vektör alanlarıyla gösterilen doğrusal fiberler Şekil 1 ‘de görüldüğü gibi dağılmıştır. Fiberlerin her biri deformasyondan önce X1 ekseni ileΦ
açısı yapmaktadır. ekseni düzlemine dik seçilmiştir. Ayrıca ve eksenleri de ’e dik olup farklı ailelere ait doğrusal fiberlerin açıortaylarını teşkil etmektedir. Buna göre deformasyondan önceki ve sonraki fiber dağılımı aşağıdaki denklemlerde gösterilmiştir [2].X X1 2 X1
X2 3
AT = cos ,s , BT =
[
cos , sin ,Φ− Φ 0]
,aT = cos ,s , bT =
[
cos , sin ,φ − φ 0]
(2.4)X2
Öte yandan homojen deformasyon, zamanı göstermek üzere, t
( )
şeklindedir[2]. Görülüyor ki λ λ λ1, 2,
F C c, , − 1
3
0
0
belli ise deformasyon bellidir. Burada kullanılacak olan ve gibi bazı büyüklükler: j
Deformasyon gradyanı : (2.6)
[ ]
F = xk K = , λ λ λ 1 2 3 0 0 0 0 0Green deformasyon tansörü : C = FTF = (2.7) λ λ λ 1 2 2 2 3 2 0 0 0 0 0
Cauchy deformasyon tansörü : c−1 =C (2.8) Sıkışmazlık şartı : j = det F = λ λ λ1 2 3 =1 (2.9) şeklinde, uzamazlık şartları da ;
A fiber ailesi için: C AKL K LA =A CAT =λ +λ (2.10)
1 2 2 2 2 2 1 cos Φ sin Φ= =
B fiber ailesi için: C B BKL K L=B CBT =λ +λ (2.11)
1
2 2
2
2 2 1
cos Φ sin Φ
şeklinde hesaplanır. Kompozit ortamdaki bu kısıtlamalardan dolayı λi
‘lerden sadece bir tanesi keyfi olarak verilebilir. Ortam, kinematik olarak buna müsaittir. Homojen deformasyondan dolayı Φ açısı , ortamın tüm noktalarında sabit olup Şekil 1’de görüldüğü gibi noktadan noktaya değişmemektedir.
Şekil 1. Ortamdaki fiber dağılımı.
Deformasyon homojen olduğu için deformasyondan önce doğrusal olan fiberler , deformasyondan sonra doğrusal kalırlar[3].
a= FA , (2.13) cos sin cos sin φ φ λ λ λ 0 0 0 0 0 0 0 0 1 2 3 = Φ Φ b= FB , (2.14) cos sin cos sin φ φ λ λ λ − = − 0 0 0 0 0 0 0 0 1 2 3 Φ Φ
(2.13) ve (2.14) denklemlerinden deformasyondan sonraki fiber açıları için
a1 =cosφ λ= 1cosΦ , b1 =cosφ λ= 1cosΦ a2 = sinφ = λ2 sinΦ , b2 = −sinφ = −λ2 sinΦ
a3 = , b0 3 =0 (2.15) değerleri elde edilir. Deformasyondan sonraki fiberlerin açısını gösteren
φ
değerleri de, homojen deformasyondan dolayı noktadan noktaya değişmemektedir [3].
φ
= sabit ,(
φ
≠ Φ)
(2.16) (2.9),(2.10),(2.11) ve (2.15) ifadelerinden deformasyon geometrisi ve deformasyon kinematiğine ait şu önemli noktalar tespit edilir [2].i) Eğer
λ
1≥0 ve λ2≥0 ise;λ
3 ≥ sin Φ2 ‘dir. Bu dadoğrultusundaki uzamanın sınırsız olduğunu,
X3
( )
λ
3 mini =sin2Φ olması nedeniyle doğrultusundaki kısalmanın da sonlu ve sınırlı olduğunu gösterir.Şekil 2. Fiberlerin birbirine dik olma durumu.
ii)
λ
3=
( )
λ
3 mini=
sin
2
Φ
değerini aldığı zaman fiberler birbirine dik olurlar (Şekil 2).iii) Fiberler deformasyondan önce birbirine dikse 2Φ = 90° , λ3 ≥1 ve
( )
λ
3 min i =1 elde edilir. Bu da malzemenindoğrultusunda uzayabileceğini ancak kısalamayacağını gösterir. Burada elde edilen (2.15) eşitliklerinden görüldüğü gibi
X3
a
=
b
=
A
=
B
= 1
olduğundan açıların değişmesi ancak birim fiber vektörlerinin dönmesi ile mümkün olur. Diğer taraftan verilen deformasyonda hız alanları da aşağıda görüldüğü gibi elde edilmiştir.( )
v x t X X x v 1 1 1 1 1 1 1 1 1 = = ⋅ = = ⋅∂
∂
λ
λ
λ
x (2.17)(
v x t X X x v 2 2 2 2 2 2 2 2 2 =∂
= = =∂
λ
λ
λ
x)
(2.18)(
v x t X X x v 3 3 3 3 3 3 3 3 3 =∂
= = =∂
λ
λ
λ
x)
(2.19)Elde edilen hız alanlarına bağlı olarak hız gradyanı ve deformasyon hızları tansörünün matrisi de:
[ ]
[
( ) L = vk l = v k l = , , λ λ λ λ λ λ 1 1 2 2 3 3 0 0 0 0 0 0]
(2.20)şeklinde elde edilmiştir. L matrisi simetrik olduğundan, d deformasyon hızları matrisine eşittir. Bu matris hesaplandıktan sonra sıkışmazlık ve uzamazlık şartlarının L matrisine getirdiği kısıtlamalar :
∇.v = trd =
λ
+ + =λ
λ
λ
λ
λ
1 1 2 2 3 3 0 (2.21) a T L a = a d aT =0 (2.22) b L bT = b d bT = 0 (2.23)şeklinde elde edilmiştir. Dikkat edilecek olursa (2.22) ve (2.23) denklemleri aynı sonucu vermektedir.
λ
λ
φ
λ
λ
φ
1 1 2 2 2 2 0 c o s + sin = (2.24) (2.24) denkleminden görüldüğü gibi v1 1 1 1 0 , = > λ λ ise, v2 2 2 2 0 , = < λ λ x2‘dır. Bu da x1 doğrultusundaki bir uzamaya daima doğrultusundaki bir kısalmanın eşlik ettiğini gösterir. Yukarıdaki (2.24) denkleminden
λ
λ
1λ
λ
φ
1 2 2 2 = − tg (2.25) ve (2.21) denkleminden deλ
λ
λ
λ
λ
λ
3 3 1 1 2 2 = − − (2.26) eşitlikleri elde edilir. (2.25) , (2.26) ‘da yerine konursa(
)
v3 3 3 tg tg 3 2 2 2 2 2 2 2 2 1 , = = − = −λ
λ
λ
λ
φ
λ
λ
φ
λ
λ
(2.27)bulunur. (2.25) ve (2.27) denklemleri deformasyonun herhangi bir hali için yorumlanırsa şu kinematik bağıntılar geçerli olur:
λ
λ
λ
λ
1 1 2 2 0 > , <0 ise : λ λ φ φ λ λ φ φ 3 3 2 3 3 2 0 1 4 5 0 1 4 5 > < < < > > , , , , t g t g o o (2.28)Bu da deformasyondan sonraki fiber açısına bağlı olarak doğrultusunda genleşme veya daralmanın meydana gelebileceğini göstermektedir.
x3
3. GERİLME DAĞILIMI
Ortaya çıkan belli kinematik kısıtlamalar altında homojen deformasyona maruz kalan Şekil 1 ’deki gibi bir kompozit ortamın gerilme analizi için bir basitleştirici/öğretici durum göz önüne alınmıştır. Cismin homojen olduğu ve aynı zamanda X3 = sabit düzlemine göre mekanik özellikler bakımından da ayna simetrisine (reflection) sahip olduğu varsayılmıştır. Gerilme dağılımı, bu düzlemlere göre ayna simetrisi altında invaryant kalacağından, koordinat transformasyon matrisi olmak üzere,
Q x1/ = x , , (3.1) 1 2 x3
τ
= x2/ = x x 3 / = − Q= , Q Q (3.2) − 1 0 0 0 1 0 0 0 1 Q T = = −1τ
/ = Q Qτ
T (3.3)ifadeleri elde edilir. (3.3) ifadesinin açık şekli yazıldığında
τ 3 1 = τ 3 2 = 0 (3.4)
olduğu görülür. Bu da bağımsız gerilme sayısını altıdan dörde indirir. Homojen deformasyonlar doğru çizgileri doğru çizgilere , elipsleri elipslere, elipsoidleri elipsoidlere, vs taşır [3]. Dolayısıyla
τ
kl,λ
,Φ ,φ
ve vektörleri konuma bağlı olmazlar. (2.1) denklemi ile verilen toplam gerilme ifadesi matris notasyonu ile
t
= −
p I
+
T a a
a T+
T b b
b T+
τ
(3.5) şeklinde yazılır. Buradaτ
‘nunτ
=
Eτ
+
Dτ
olduğu hatırlanır ve de değerleri yerine konursa gerilme bileşenleri verilen koordinat sisteminde aşağıdaki gibi elde edilir.a a
T,
b b
Tt
11= − +
p
(
T
a+
T
b)
cos
2φ τ
+
11 (3.6)t
22= − +
p
(
T
a+
T
b)
sin
2φ τ
+
22 (3.7)t
33= − +
p
τ
33 (3.8)t
12=
(
T
a−
T
b)
sin cos
φ
φ τ
+
12 (3.9)t
23=
τ
23=
0
(3.10)t
3 1=
τ
3 1=
0
(3.11) Diğer taraftan (2.1) denklemiyle verilen toplam gerilme ifadesi denge denklemini(
t
kl k,= 0
)
sağlamak zorundadır.
t
k l k,= −
p
,l+
T a a
a l, k l+
T b b
b l k, l= 0
(3.12) Burada , cisimsel ve atalet kuvvetleri ihmal edilmektedir. (3.12) ‘den görüldüğü gibi için seçilen keyfi sabit değerler için denge denklemi sağlanır.p
,
T
a,
T
bp
=
k
1 ,T
a=
k
2 ,T
b=
k
3 (3.13) Öte yandan,t r
τ
=
τ
k k= 0
(sıkışmazlık şartı) (3.14)a
Ta
a a
k l k lτ
=
τ
=
0
=
(A fiber ailesi için uzamazlık şartı) (3.15) (B fiber ailesi için uzamazlık şartı) (3.16)
b
Tb
b b
k l k l
τ
=
τ
0
kısıtlamalarına, (2.4)3 ve (2.4)4 eşitlikleri yerleştirilirse
τ
12=
0
(3.18)τ
12 cos2φ τ
+ 22 sin2φ
=0 (3.19)φ
denklemleri bulunur. Bu denklemlerden de;
(3.20)
τ
11= − tg
τ
22 2τ
33τ
22(
φ
)
21
=
tg
−
(3.21) eşitlikleri bulunur. (3.18), (3.20) ve (3.21) eşitlikleri (3.6) , (3.8) ve (3.9) ’da yerine yazılırsa gerilme bileşenleri aşağıdaki gibi elde edilir.t
11p
(
T
aT
b)
tg
2 22 2= − +
+
cos
φ τ
−
φ
(3.22)t
2 2p
(
T
aT
b)
2 2 2= − +
+
sin
φ
+
τ
(3.23)t
33p
22(
tg
)
21
= − +
τ
φ
−
(3.24)t
12=
(
T
a−
T
b)
sin cos
φ
φ
(3.25)Burada elde edilen gerilme bileşenlerinden sadece bir tanesi bağımsız olup diğer üç gerilme keyfidir. Bu keyfilik
keyfi reaksiyon gerilmelerinin bulunması için kullanılabilir. Bu amaçla keyfi olan üç gerilme
p T
,
a,
T
bt
22=
t
33=
t
12=
0
şeklinde seçilip (Bu, doğrultusunda tek eksenli gerilme hali demektir.) , (3.22), (3.23), (3.24), (3.25) gerilme bileşenlerinde yerine konursa X1p
=
τ
33=
τ
22(
tg
φ
−
)
21
(3.26)T
a=
T
b=
tg
−
−
2 2 2 222
φ
φ
φ
τ
cos
sin
,(
2
φ
≠
0 90
)
o,
o (3.27) reaksiyon gerilmeleri bulunur. Bulunan değerleri (3.6) ‘da yerine konursap T T
,
a,
bt
11 22(
tg
)
22
1
=
τ
−
φ
−
cot g
2φ
(3.28)gerilmesi elde edilir. O halde verilen şekilde fiberlerle takviye edilmiş olan ve X1 doğrultusunda tek eksenli bir gerilmeye maruz kompozit
ortamdaki gerilme dağılımı
t
11 22(
tg
(3.29)2
2
1
=
τ
−
φ
−
cot
g
2φ
)
t
22=
t
33=
t
12=
t
23=
t
31=
0
,(
φ
≠0
0 ,90
0)
(3.30) şeklinde elde edilir. Bu denklemlerde gözükenτ
22değeri de
τ
2 2=
Eτ
2 2+
Dτ
2 2 (3.31) formunda olduğundan, (1.2) ve (1.3) ile verilen bünye denklemlerinden bulunur.Burada ele alınan problem Spencer’in Şekil 1’de gösterilen problemiyle aynıdır. Bu problemde Spencer, elastik bir kompozit malzemeyi düşünmüş ve toplam gerilmenin bünye denklemiyle elde edilecek terimleri üzerinde herhangi bir işlem yapmamıştır. Gerilmenin bu kısmını reaksiyon gerilmeleriyle birlikte toplam gerilme ifadesinin sağ tarafında bünye denklemleriyle bulunmak üzere muhafaza etmiştir. Dolayısıyla bu safhaya kadar Spencer ve araştırmacı tarafından yapılan işlemler, malzemenin bünyesine bağlı olmayan işlemlerdir. Araştırmacı tarafından Spencer’in problemi fiber kinematiği de işlenerek viskoelastik bir kompozit malzeme için çözülmüş ve gerilme dağılımı bulunmuştur. Bundan sonra (1.3) denkleminde sıkışmazlık ve uzamazlık şartları
λ
b=1
,d
mm=0
,a da
T= 0
,b db
T= 0
(3.32) denklemlerinde yerine konup, (1.2) ile toplandığında bu problem için kullanılacak toplam gerilme ifadesi bulunur[ 4-6 ].(
)
(
)
∂
∂
∂
∂
∂
∂
=2 − +2 − − − +2 − + − 1 1 2 1 1 2 7 1 1 Σ Σ Σ I c I I c c I a a c c a a T Tτ
(
)
(
)
+2 − + − + 2 8 1 1 9∂
∂
∂
∂
Σ Σ Φ I b b c c b b I a b b a T T c o s T + T(
)
+
∂
−+
−+
−+
−+
∂
µ
Σ
Φ
I
ab c
c ba
ba c
c ab
d
T T T T v 10 1 1 1 12
2
cos
[
+ 2 d1 2 c o s2Φ I + f1 2 c o s2Φ a a T + f2 1 c o s2Φ b bT(
)(
)
]
+β
12 +γ
22 cos22Φ abT +baT a dbT (3.33) (3.33) denklemindekiµ
v,
d
12,
f
12,
f
21,
β
12veγ
22katsayıları yalnız, θ ve X ’e bağlıdır. Gene bu ifade , (2.4)
cos2Φ 3 , (2.4)4 ,(2.8) ve
(2.18) eşitliklerinin kullanılmasıyla λ1 , λ2 , λ3 cinsinden bulunur ve
buradan da τ22 gerilme bileşeni aşağıdaki gibi elde edilir.
(
)
[
]
τ
∂
∂
λ
∂
∂
λ
λ λ
∂
∂
λ
φ
∂
∂
λ
φ
22 1 2 2 2 3 2 1 2 2 2 7 2 2 2 8 2 2 2 2 2 4 4 = Σ + Σ + + Σ + Σ I I I sin I sin −2 2 − 4 2 2 9 2 1 0 2 2 2 2 2 +∂
∂
φ
∂
∂
λ
φ
µ
λ
λ
Σ Φ Σ Φ I c o s sin I c o s sin v[
+ 12(
12 21)
2 2d cos Φ+ f + f cos Φsin φ
2 2
(
)
]
−2 12 + 22 22 2 2 1 1
2
β γ cos Φ sin φ λλ cos φ (3.34) Böylece X1 doğrultusunda tek eksenli gerilmeye maruz fiberli
viskoelastik kompozit ortama ait gerilme dağılımı (3.29) ile (3.34) denklemlerinin birleştirilmesiyle aşağıdaki gibi elde edilir.
(
)
[
(
)
]
t tg g I I I 11 2 2 1 2 2 2 3 2 1 2 2 2 7 2 2 2 1 4 4 8 = − − + + + φ φ ∂ ∂ λ ∂ ∂ λ λ λ ∂ ∂ λ φ cot Σ Σ Σ sin +8 −4 2 −8 2 +4 8 2 2 2 9 2 10 2 2 2 2 2∂
∂
λ
φ
∂
∂
φ
φ
∂
∂
λ
φ
µ
λ
λ
Σ Σ Σ ΦI sin I cos sin I cos sin v
(
)
[
8 1 2 2 1 2 12 12 21 2λ
λ
cosφ
d cos Φ+ f + f cos Φsinφ
(
)
]
SEMBOLLER:
A a
,
:A- fiber ailesinin deformasyondan önceki ve sonraki fiber dağılımıB b, :B- fiber ailesinin deformasyondan önceki ve sonraki fiber dağılımı
C ,c ; C ,c :Green ve Cauchy deformasyon tansörleri ve matrisleri
C−1,c−1;C−1,c−1 :Piola ve Finger deformasyon tansörleri ve
matrisleri
d
:Deformasyon hızları tansörüE , E :Lagrange (maddesel) strain tansörü ve matrisi
e
,
e
:Euler (uzaysal) strain tansörü ve matrisiF
:Deformasyon gradyanı matrisiI i, :Maddesel (Lagrange) ve uzaysal (Euler) koordinatlarda birim vektörler
j = d e t F :Deformasyon jokobiyeni
L
:Hız gradyanı matrisiP
:A- fiber dağılımı tansörü p :Hidrostatik basınçS
:B- fiber dağılımı tansörüT T
a,
b :Fiberlerin uzamazlığından kaynaklanan reaksiyon gerilmeleriT
, ;
t
T
,
t
:Maddesel ve uzaysal koordinatlarda gerilme tansörleri ve matrisleriE
T , t
E :Maddesel ve uzaysal koordinatlarda elastikgerilme tansörleri
D
T
,
Dt
:Maddesel ve uzaysal koordinatlardadisipatif gerilme tansörleri
Q
:Transformasyon matrisiδ
KL ,δ
kl :Maddesel ve uzaysal koordinatlarda kronecker deltasıλ λ
a,
b : a ve b fiber ailelerinin uzama oranlarıλ λ λ
1,
2,
3 :Asal doğrultudaki asal germelerρ ρ0, :Deformasyondan önceki ve sonraki kütle yoğunluğu
Σ
:Gerilme potansiyeliσ
kl :Toplam reaksiyon gerilmesiτ
:Gerilme-bünye dağılımı matrisiΦ
:Deformasyondan önce fiberlerin X1 ekseni ileyaptığı açı
φ
:Deformasyondan sonra fiberlerin x1 ekseni ileyaptığı açı
XK,xk (K ,k :1,2,3) :Maddesel ve uzaysal koordinatlar KAYNAKLAR
[1] Öntürk,N.; “İki Fiber Ailesi İle Takviyeli Viskoelastik Kompozit Ortamlarda Bünye Denklemlerinin Modellenmesi”, Doktora Tezi, Gazi Üniversitesi, Mühendislik Mimarlık Fakültesi, Makina Mühendisliği Bölümü,
[2] Spencer, A.J.M.,”Deformation Of Fibre-Reinforced Materials”, Clarendon Press, Oxford,(1972).
[3] Spencer,A.J.M.,”Continuum Mechanics, London, (1980).
[4] Erıngen,A.C.,”Nonlinear Theory Of Continuous Media “, Mcgraw-Hill Book Co.,(1962).
[5] Şuhubı,E.S.,”Sürekli Ortamlar Mekaniği”, İstanbul Teknik Üniv. Rektörlüğü , (1993).
[6] Erıngen,A.C.,”Mechanics Of Continua”,Robert E. Krieger Pub. Co.,Huntington, New York,(1980).