SAU
Fen Bıhnıleri En:;tıHlsli lkq�ı�ı 7.Cilt, 2.Sayı(
renını u/ 2001)Dlit IJI lllt·r ele 1 1 d 1
·trilm
rernin IDctltnmt.si
.Ural 1.Öı 0 \ V.Uçar
Faruk UR�L\.1,, Murat
ÖZSOY,
\'ahdct
lJ(.
R
Özet
-Günlük
hay
at
ta çok geniş kullanıln1a
;dan
ın asa
hip bir makine eJenıanı olarak,
güç ilctiıninde kullanılan diş
l
iç
arklar
ın
muhtelif
şartlarda
davranışları
bir çokaraştırn1aya
l{onu oJınuştur.Bu
ç
al
ışm
ada
silindirik düzç
arkla
rd
a diş profilin
de
ki
gerilmeler,
klasikyöntenılere de
değ
inil
ere
k Sonlu E.lemanlar MetoduiJe
çözümlemelery
ap
an
ANSYS
paket programı kullanılarak
incelendi. Modelolarak
kavranta açısı
a =20°,
di
şsayısız= 22,
dişd
ibi kavisirk
=0.9
ının, diş deri nliği t =3.75
nım, diş çıkıntJsı s =3
mmve modülü3
ol andüz
dişliçark
analiz edilmiştir.Analiz
sonucunda ol
u
şan
maksimun1 gerilıne dişlimalzenıesinin
akma sınırının al
tındak
al
mışt
ır.Anahtar kelinıele1·-
Dişli
çarklar, Sonl
u ElemanlarMetodu, Gerilme Analizi
Abstract
-In
dailytiınes, gears have
usi
ngareas as
m
a
chine elements, and their behaviors in var
ious
conditions become sub
j
e
ct
s to lots of researches.In
this research,
the stressat the
profile of tooth inc
y
lin
drica
l
spur gears ha
vebeen
examined by us
in
gANSYS
computer
code \Vith touching cl
assica
l
metbods. A s
pur gearwhich has
coupling angle a =20°,
n
um
b
er
oft
eet
b
z =22,
t ooth
botto
m bow rk =0.9
mm, tooth depth
t
=3.75
mm, tooth protr
usion
s=3
mm as the
model.At
the end
of theanalysis
maximum
stress islower
than the material yi
eld
s
tress.Keywords-
Gears, Fi nite Elements Method,Stres
Analyses
I.
GİRİŞ
İ
n
sanoğl
uyerleşik
düze nden bu
yana
tabiatla
mücadelesinde ilk önce
kendigüc
ününyanı
sı
ra
som·a da
hayvan gücünden
f
a
yda
l
a
mm şt
ır
.Daha
soma doğadan
y
ani rüzgar ve sudan yararlanmıştır.Fakat bunların
kullanıınında ortaya çıkan
zorluklard
iğer güçi
l
etim
sisteınlerini n
kul
la
n
ımın
ı zorunlu hale getirmiştir.Bu
zorluklar
k
ar
şı
sı
n
d
a insa
n
oğlu, belki deke
ş
fe
dilmiş eneski makine
elemanlarmdan
ola
n
,dişlileri
kull
ana
r
ak a şn
ıay
a
çal
ışmı
ş
tır[ 1].
F. Ura] ;
SAÜ,
Bilgiişlem
Daire Başkanlığı, SakaryaM.
Özsoy; SAÜ
Mül�endilik Fakültesi CAD CAM CAE MerkeziV. Uçar;
SAO
MülıendısJjk Fakültesi tv1akine MühendisliğiBölümü,Sakarya
Dişli,
biı· nııldendiğcrınc.:-
lıurrkel ve
g
O\
.. iletimini şekilbağı ilc sağluy�1n
rnakıra·
cleınaıııdır
.
�1illcrin konumunagöre
dişli,·:.ırkl�u·
vedişli
çark ıneknnizmaları
sın1ilandı rılır.
Günüınii7de
di�) iler, h,ı' �ı
c ı1
ı k veu7.ııy
sanayiinden tanmnıakinalaruıa kadar ht.r
alanda
yaygın olarakkulianılmaktadıı.
()/ellıkle ağır şartlarda, ya da yüksek hızda çalışannıcıkıııalaıda k
ullunıl
an dişliler kıitik biriş
l
ev iy
eri
n
egct
innekted
ir Jer.f !u
nden
le
bir dişli çjftininperformansının
dogru
olarak heliılc·ınıcsi ha
y
atı bir önem kazaıın-ıaktadır.J
)eneyse I
ve anal ı tik yöntenıJerle eldeedilecek sonuçlann sınırlı ol
n
ıa
s1 ve ol
a
sı
her di
şli çi
ftin
e uygulannıasınınınün1kiin olnıarnası
dişli
çı
ftinin dinamikkarakteristiğinin
ı node
ll
cnnıe
si
nin önenun i aıt1rmaktadır.II.
ANSYSSONLll IDLEMANI�.L\R
ANAtİZ
]•ROGRAMI
Bilim adamlan alışılını ş
an ah
t
i
knıetotlarJa
çözün1ü çok zorhatta
imkan
sızfiziksel
problemlerle sık sıkkarşılaşırlar.
Boyle bir problcn1le karşılaşan çözümleyici doğal ola
ra
k � 'sa
yısa
l
"adı verilen
çözüıne ba
şvura
caktır.Başka
metotlarlaÇÖ/.ülenıeycn
prob:cınJerin çözüı
niindekullanılabilen
çok bayıd
a
sayısal yolvardır.
S
onlu eleman met
od
ubunlardan
birid
i
r.
Sonlu
eleman
n1etoduyeni
bir çö
züm yönte1ni olup kendis
in
i
diğerlerine üs
tün kılan seçkin öze1liklere sahiptir. Bu meto
ttacismin "sonlu,
b
oyutt
a çok sayıda "eJenıana''ayrıldığı
düşünült..Metodun adı
daburadan
geln1ekteclir.
Cisim uzayda n(=
1,2,3) boyuta
sahipse, n
-bo
y
utl
u son lu elemanlar
sisternine
ayrıllr.
Sonlu elemanlar yö
n
temi, boyutları ve şe
ki
lle
rin
inesnekli ği nedeniyle, çok
ba
ğl
a
n
tılı
böl
gel
er
veya köşelerjolan bölge
l
e
r zorluk çekilmeksizinincelenebilir.
Seb
epsonuç bağlantıianna
ait
problemler gl
obal
direngenlik
matrisi ile
birbirineb
a
ğ
lan
a
ngenelleştirilmiş
"kuvve
tler" ve"yer değişt
iımeler'' cinsinden
fom1üleedilebilir.
Sonluel
e
manl
a
rmetodunun bu özelliği problemin
anlaşılınasınıhem mümkün
kılar
he nı debasitleştitir.
Günümüzde piyasada
bir
çok sonlu
ele
man yöntenli
il
eSAU Fen Bilimleri Enstitüsü Dergisi
7.Cilt, 2.Sayı (Temmuz 2003)
1970
yılından bu yana kullanılan yaygın bir analizprogramı du [ 2].
UI.DİŞLİ
Ç
ARKLARD
A
GERİLME ANALİZİ
Dişli çark boyutlarımn tayjninde gerekli ömür süre
si
nceherhangi bir hasarın meydana gelnıcsin
i
önlemek içinortaya çıkan diş zorlamaJannın ınünı.kün olduğu kadar
hassas
bir şekilde bilinmes
i gereklidir.Eş çalışan dişlilerde en büyük geriln1elcr diş tabanında
olur ve kın lmalar bu bölgede meydana gelir.
Y
orulmaözelliği taşıyan bu kırıln1alar, çekıne geriln1elerinin
bulunduğu yerde bir çatlak i1e başlar. Zanıanla ilerler ve
yükü taşıyamayacak hale gelince diş dibinde kopma olur.
Diş yüzeylerinin taşıma kabiliyeti� müsaade edilen diş yüzey basıncı ile belirlenir. Taşınabilecek basıncın aşılnıası halinde özellikle sert1eştirilmenuş yüzeylerde ve tekil kavrama noktaları civarında, malzeme
parçacıklarının kopması sonucu ufak çukurcuklar şeklinde hasar meydana gelir.
İlk
defa 1892 deI
.. cwis diş kökü muka vernetine ilişkinbasitleştirici kabuller ya
p
tı. Bu araştırmada Lewis,dişlinin yüklenmesine ilişkin en kötü konumun tüm
yükün kavrama halindek
i
tek dişlitarafından
taşınmasıhali olan ilk kavrama noktası olduğunu ve diş
in
kirişo]arak değerlendirilebileceğini varsaydı. Lewis yükün
radya1 bileşenini
ihına
l etmiştir[3].Şekil
ı. Dişleri n mukavemet hesabıFot.oelastisite yöntemleriyle elde edilen sonuçlara göre eş çalışan dişlilerde en büyük gerilmeler diş tabanında meydana gelir ve kırılmalar bu bölgede olur. YükJemniş
diş
k
ök
ünde maksinıum geri]menin tespititcınas
anındadiş üzerine gelen yük ün yönü, miktarı ve dişljnin şekliyle
101
Düz Dişlilerde Dişdibi Gedln1clerinin
incelenmesi
••
F.Ural, M.Ozsoy, V.Uçar
değişkenlik göstermektedir. Temasın ilk noktasında,
F
n diş yükü,F
r , radyal ve Fı teğ
etsel yük bileşenlerineayrılmaktadır[ 4].
(1)
F:-c
==Fnc Sin
ab(2)
bu de
ğ
erler
den yola çıkarak sırasıyla eğilme ve basmagerilmeler i cr
e
'(jb
_
F,,c Sin
a 6bsq
(3)
(4)
şeklinde hesaplanır. Burada hq; e
ğ
ilme kolunu, sq; kınlma kes
itini
ve b; diş genişliğini sembolize etmektedir.Bu verilerden yola çıkarak toplam gerilme, Gtop aşağld.aki
gibi hesaplanır. Toplam gerilme nıiktan parçamn
n1ukavemeti
için, enınıyct gerilmesinden, O"em küçük olmalıdır.(jtop == (Ye -ab
�
(jem
(5)
I\'. ANSYS
PROG
RAM
ININ
UYGULAN:\IIASI
Dişli çark profilini tayin eden ampirik bağıntılardan sonlu
elemanlar n1odelinin tesisinde referans noktaları
oluşturan bir prograın
g
eli
ştirild
i.B
.. ı noktalar ANSYSor
t
arnında değerlendirilerek n1odel oluşturuldu.Bir dişli çarkta, diş kuvvetinin diş ucuna tam değerd
e
uygularınıası en kötü haldir. Gerçekte kavrama oranına bağlı olarakdişli
çiftleri arasında paylaşılır. Dişler yi11 paylaşımını sağlayacak hassasiyette imal edilmedikleri sürece dişkuvvetinin
diş ucuna tam değerde etki ettiği kabul edilebilir. Bu nedenle diş dibinde azamig
erilmeyi sağlayan yUk tatbik noktası diş ucudur[5].lV.
1 Model Oluşturma
Bu çalışmada kullanıJan dişli ilk olarak bir bilgisayar destekli tasarım programı olan
Pro/ENGINEER
da ınodel lendi.Daha
sonra model iges transfer dosyafornıatında
ANSYS
yazıluruna aktarıldı. Kullanılan dişliSAU Fen Billınieri Enstitüsü
Dergisi
7.Cilt,
2.Sayı (Tenımuz 2003)
Malzeme
SembolO
KopO
Akma ElastisiteM
o d
Poisson Oranı : SenıentasyonÇeliği
: 1 8CrNi8
: ( 1200-1450)
x 106N/n12
: 800x 106
N/m2
O
06 NI
2:21 xl
,m
: 0.3
AKSYS yazıluruna
aktarılan
n1odclin
elen1anlarabölünmesi
i
ç
i
n3
boyutlu.6
yüzlü, 8
düğüm
noktası
veher
düğüm noktasıiçin
3 serbestlikderecesine
sah
ip
olan(
3D SOLID186
)
katıeleman kullanıldı.
l X
®�
/< l / .. _,ı 1\
®
M
1
-\
y ,1-
,tl--_
-N - ..._
Şekil
2. Solid 186 Katı Elemancv
t:7
,..r o/
®
- ---·KSonlu elenıanlar
n1odeline uygulanan sınır şartları
kinematik
ves
tatikolarak ikiye ayrılır.
Kinematiksınır
şartlan
düğüm nokt
a
la
rına tatbik edilen hareketkısıtlamalan
dır.
Dişlimodelinde, Şekil
3.'
de göıiildüğügibi
X,Y ,Z
eksenlerinde sabitlenmiştir�---_J
Şekil
3.Dişli
çarkaait
sonluelemanlar ınodeli
l>ö�: Ui�lilcrdc Uişdibl G rilmeltırinin
lncele
n
m
·
F.tJral, M.Özsoy,
v
.
uç:
()luştundan
ınodel,I
8.051
adet
vlcn1an
ve
26.368
adet
diiği.inı noktasından olu�nıaktadır.i
kidişlinin biıhirini
enÜLzla
yüzeytemasıyla kavradığı
durun1u teınsil
elınesi
arnacıyla,şekil
4.'de
1
ııumara ilegö�lterılen
alan
ü?er
ıııe 400 N nun2 hasınç uygulanrnıştır.
---- -----
·----�eki!
4 lJi�h çarka dlt yukknıe durumuV.
SONUÇLAR ve YORUM
\1• 1
X Ekseni Boyunca
Oluşan
Yer Değişimi
NOül\1. S OLUT Tt�•H STEF=l SUS al 'I'H-E,.l t;;< ( �VG 1 CMX =.02623!> SHN .,_ .025945 şı� �.192G-04 ---- -: . . .... .
• C259otı>
-
.
oıoı&.s-
. OHJ79- 0230H .017271 .01.1488 . uJ8�96 .002812
-. oo!ı70t
.
792E-o•Şekil
5. X Ekseni Boyunca OluşanYerDegiştirnıt.
Uygulanan basınç ( 400
N/ınm2
)
sonrasında Xekseni
boyuncaoluşan n1aksimum
yer değiştirme X eksen
i
negöre
ters(-)yönde 0.025946 n1n1
diJ·.SAU
Fen BilimleriEnstitüsü Dergisi
7.Cilt,
2.Sayı (Temmuz 2003)
V. 2 Y
Ekseni BoyuncaOluşan Yer
DeğişimiprtP=l
S1J6 :11
T!ME=:ı
UY IAVG) ['�t>: ; . 02 523 6 !!Wl -=-.0)8Z14 $Mli =-.OJ64J6 -.000214 •.004965 - .COb511 ı ' •.• J. ııı:cıı� -.0011" ,0015'lı OM"'Jl -. 003J1\ . 003l!ı1Şekil 6. Y Ekseni 13oyunca Oluşan Yer
Değiştiıme
.0064)6
Uygulanan
b
as
ın
ç
( 400
N/
nırn2
)sonrasında
Y
ek
se
ni
boyunca oluşan maksin1um yer değiştirme
Y
e
ksenine ters (-) yö
ndeO.
008214
mmdir.
V.
3
ZEkseni BoyuncaO
l
uş
anYer
DeğişiınİSTEr-
1
!:U5 •1
1 TIME:•!
l.l2
ı .• t..VG ı 00 -. .G:!ö2�6 SHII -.OQZC31ı
S"<io • • Cı02U:J3 , --··r·•·lıı· ... '•ı
· Oil20ll . L ! ,. .<"ı . . . � • 001110 • ���e-Ol . GIO!·Ol .001581Şekil
7. Z Ekseni Boyunca Oluşan YerDeğiştirme
. 00/U)J
Uygulanan
basınç ( 400
N/nmı2
) sonrasında Y
ek
seni
boyuncaoluşan
maksiınum yer deği
ş
tiı
me
0.002033
nun dir.V. 4 Eşdeğer Gerilme
DağılımıSTtl'- 1 !LI!s 1 'I'JY.Dol 'CQ\i (AV\i) n� • .o�..:�Jt\ ıs�: ,. .ııı:� :11 '!W. •tC6.74l l � a:•s= .Oll�>Oq 1J4.8ll 269.t:69 4011,1198 •,J<�.326 61.4/\> .t02.255 337.063 411.9];: ().)(.,7.,1
Şekil 8. Eşdeğer Gerilme Dağılımı
103
Düz Dişlilerde Dişdibi Gerihnelerinin
incelenmesi
F.Ural, M.
Ö
ı,soy, V.UçarUygulanan
b
a
sınç
( 400
N/ının2
)
soruasında oluşanmaksimum
g
eril
me değeri
606.741
N/mm2
dir.( Ş
ek
il 8.) Yapılançözümlemeler som·asında
maksimum yer
d
eğ
iş
tiını
en
i
n
X
eksenindeve
nm1'nin
binde 2.5' umertebesinde
olduğu görüldü.Model
üzerinde ise oluşan
maksimum geril
me
nıalzemenin akn1a sınırı
o
lan800 N/mm2'
nin altınd
akal dı. Fakat burada dikkat ed
ilmesige
rek
e
n konu dişliçarkın emniyetli çalışı11ası
için
b
ir en
ıniyet katsayısı
belirleınektir.
Emniyet katsayısıS=
1.25 seçildi. Buna
göre emniyet]
ig
eri
lıne sının 800/1 .25= 640
NI mm2
d
ir
.
Çözümleme
sonucunda oluşan606.741 N/nnn2'
nin640
N/mm2'
dendaha
küç
ü
k olduğu içinbu
dişliemniyet
liç
al
ış
ır.VI.
SONUÇYa
p
ı
lan
çalışmada
diş
li
çarklardadiş
kökleri üzerindekigerilmeleri hesaplamak
amacıyla ANSYS
sonluel
e
man
la
rpaket
progranukullanılarak
g
eri
ln1e analiziyapıldı.
A
s
al
geriln1eJer vasıtasıyla çeşitli kriterieregöre
ge
ri
ln1e ve
yer değiştirmeler sapt
andı
ve d
işli
nino
lmasıg
erekend
eğe
rl
er için
dek
al
dı
ğıb
el
irl
endi.
ANSYS
paket progranu kullanılarak pratik ve hızlıçö
z
ümle
rneler eldeed il diğinden,
ç
alışma
nın
geniş letil n1esine
ve değişikçap, diş
sayısı, n1odül veeğrilik yarıçapına sahip dişlileri n
ele
al
ınara
kkarşıla�tırmalı
bir
çalışma yapılm
asına
karar verildi.KAYNAKLAR
[1]
AYDENİZ,
A.
İ. ;
"D
işl
i
ÇarkMekanizmalarında
Taksimat
Ilatasımn
Statik v
eDinamik
z_,rlanmayaEtki leri", Yüksek Lisans
Tezi, İ.T.Ü.
FenBil .
Enst.,İstanbul,
2001
[2] ÖZSOY,
M.
;"TVS 2000
Yolcu Vagonlannn�Bilgisayar
Ya
rd
ın
ny
la Yapısal Analizin, Yüksek Lisans Tezi, Sakarya Üniversitesi FenBil. Enst.,
Sakarya, 1998.[3]
ÖZER,
1 I. ;"Düz D
işl
i Çarkların Sonlu ElemaıılaiYöntemiyle
Gerilıne
Analizi",Yüksek
LisansTezi,
S
aka
ry
a
Üniversitesi Fen Bil.Enst., S
ak
a
ry
a, 1996.[
4
] AKKURT, M., "Makinee
lemanları", Birsen
Yayın
e
vi
, 3. Cilt, İstanbul, 1982.[5] COLBOtJRNE,