Geri Dön

Elastik krank-biyel mekanizması titreşimlerin sonlu elemanlar yöntemiyle incelenmesi

Vibration analysis of cranck-connecting rod-slider mechanism by using finite element method

  1. Tez No: 46549
  2. Yazar: ERHAN DELİGÖZ
  3. Danışmanlar: DOÇ.DR. VEDAT KARADAĞ
  4. Tez Türü: Yüksek Lisans
  5. Konular: Makine Mühendisliği, Mechanical Engineering
  6. Anahtar Kelimeler: Belirtilmemiş.
  7. Yıl: 1995
  8. Dil: Türkçe
  9. Üniversite: İstanbul Teknik Üniversitesi
  10. Enstitü: Fen Bilimleri Enstitüsü
  11. Ana Bilim Dalı: Belirtilmemiş.
  12. Bilim Dalı: Belirtilmemiş.
  13. Sayfa Sayısı: 70

Özet

ÖZET Bu çalışmada amaç, elastik krank-biyel mekanizması nın titreşim karakteristiklerinin sonlu elemanlar yönte miyle incelenmesi ve sistemde yer İTİ an çeşitli,frekans par amet reler inin yine dizayn sırasında kullanılacak olan sistem parametrelerinden(literatürlerde daha önce söz konusu edilebilen çalışmalara ek olarak) piston ile si lindir arasındaki lineer yay etkileri, boyuna yöndeki yer değiştirmeler, boyuna kuvvet etkileri ve krank elas tikliği dahil analiz edilmesidir. Bu amaca bağlı kalarak, öncelikle krank biyel mekanizmasının davranışlarını yete rince ifade edecek bir matematiksel model oluşturulmuştur. Bu model oluşturulurken, krank açısal hızının sabit olduğu kabul edilmiş, piston kütlesinden dolayı krank ve biyele etkiyen normal kuvvetlerin mekanizmaya etkileri göz önüne alınmıştır. Ayrıca piston ile silindir arasındaki lineer yay etkileri de göz önüne alınmıştır. Titreşim analizi yapılırken, rijit krank elastik biyel ve elastik krank biyel olmak üzere inceleme yapılmış ve bu şekilde elde edilen titreşim karakteristikleri < titreşim öz frekans ları ve titreşim biçimleri) arasında karşılaştırma yapı larak, çeşitli parametrelerin Cpiston yay katsayısı, piston kütlesi, devir sayısı, krank açısı değişimi) etkileri incelenmiştir. Her iki karşılaştırma modeli için fortran IV programlama diliyle programlar hazırlan mış ve sonuçlar gerek tablolar, gerekse grafik şeklinde gösterilerek analiz edilmiştir. - K -

Özet (Çeviri)

VIBRATION ANALYSIS OF CRANK-CONNECTING ROD-SLIDER MECHANISM BY USING FINITE ELEMENT METHOD SUMMARY In the present work, vibration characteristics of crank -connecting rod-slider mechanism, are analysed con sidering bending and longitudinal deformations of crank and coupler, and thin beam theory, by the finite element method. First, some details about application of the fini te element method, and then, the details of the vibration analysis of the elastic crank -si i der mechanism are summa rized as follows ; General Considerations About Finite Element Method The increasing complexity of structures and sophis tication of digital computers have been instrumental in the development of new methods of analysis, particularly of the finite element method. The idea behind the finite element method is to provide a formulation which can ex ploit digital computer automation for the analysis of ir regular systems. To this end, the method regards a comp lex structure as an assemblage of finite elements, where the points are known as nodes, the entire structure is compelled to act as one entity. Although the finite element method considers conti nuous individual elements, it is in essence a discretiza tion procudure, as it expresses the displacement at any point of the continuous element in terms of finite number of displacements at the nodal points multiplied by given interpolation functions. To illustrate the idea, we refer to the one- dimensional system shown in fig. 1. The system is divided into a finite number n of elements of width 1, where n. 1 * L, and the motion of the system is defined in terms of the“nodal displacement”UiCO Ci =1, 2,..n3. The advantage of finite element method over any other method is that the equations of motion for the sys tem can be derived by first deriving the equations of mo tion for a typical finite element and then assembling the individual elements 'equations of motion. The motion at any point inside the element is obtained by means of in terpolation, where the interpolation functions are gene rally low-degree polynomials and they are same for &sr&ry element. X Fig. 1. Interpolation Functions - XIThe finite element method» as practiced today, began as a method of structural analysis, being related to di rect stifness method. This direct approach may be satis factory for static problems, but encounters difficulties in handling dynamic problems, such as in vibrations of con tinuous media. Such problems are treated better by a va riational approach. In fact, the finite element method can be regarded as a special case of Rayleigh-Rltz method, although since its inception the method has acquired a li fe of its own, going well beyond the original structural appl i cat i ons. Our purpose is to present some of basic ideas invol ved in the use of the finite element method for vibration problems rather than a detailed treatment of the subject. We shall be concerned only one-dimensional elements, al though the concepts and developments presented quite ge neral and can readily be applied to two-and three-di men- si onal el emen t s. - - Derivation Of The Element Mass» Stiffness Matrix By The Direct Approach We adopt an approach by deriving the element mass and stiffness matrix as the matrix relating the nodal for ce vector. To this end, we consider a rod in axial vib ration and derive the mass and stiffness matrix for a ty pical element such as in fig. 2. We carry out the task in two steps. In the first step we derive an expression for the axial displacement of an arbitrary point inside the element in terms of the nodal displacements and in the second step we use this expression to relate the nodal displacements to the nodal forces. Although, in vibrati ons all these quantities are functions of time, for the purpose of deriving the stiffness matrix, they can all be regarded as constant in time. i Ç(t) U, Ui“* f2ct>. fta/t) A Fig. 2. Derivation Of Axial Stiffness And Mass Matrix The axial stiffness can be assumed to be constant over the element, so that the differential equation for the axial displacement UCx) is : EA 32U ax2 = O CO < x ”Uz- Ui Ci =. ; Cz = Ui inserting the constants just obtained into eq. C2Z>, we ob- obtain the expression for the axial displacement ; U = CI --*->. Ui. + *-Uz. C5Z> We can expression kinetic and potential energy for the e- lement as ;“ _ 1 f! «U.2 _ 1 fV ÖU _: Ek ~ 2”mJC -q^-) i Ep - 2-EAU -^-) C6D Using Lagrange equation C73 in eq. C5D, C6Z> we obtain C M ] and C K 3 as follows ; d f 3L dt^ Oq ÖL 7> --g- = Ffc Ck = 1,2..nZ> ; L = E^ -E C73 [ M 3 = -T2 'l.,,,.jaT1 -1! 6^1 2J ' K J TT^-i lj C8D where C M 3 is mass matrix and L K 3 is stiffness matrix The another case of particular interest is the bar in bending vibration. We propose to use the same approach as above to derive the corresponding element shown in fig. 3 I i i.u*«> itKx,t> U{*> U Fig. 3. Derivation Of Bending Stiffness And Mass Matrix For uniform bending stifness, the differential equation for the displacement UCx.tZ) is EI ox4 CO < x four times, we have UCxD = -|-Ci.xa+-|-C2.x2+ Ca.x + C*. C1CÖ To determine these constants, we refer to fig. 3 and write UCOZ) = Ui; UC03 = Uz; UC13 = Ua; I^ID = U* C11D where Ui. Ua are nodal displacements and Ua, U4 are nodal rotations, or nodal angular displacements. Introducina eq.ClOT into eq. CUD we obtain the expression for the ben ding displacement as follows ; U = CI ~ 3C -fo2 + 2C -f Z>a)Ui +Çx -21C J^O2 + 1C J*0 9;> u2 +C3C^D2 - 2C^^)Ua +Ç-lC^-f +1C^>U* C12> In equations C6D, C125 using lagrange eq. C73, we obtain t M j and t K 1 as follows : CM3 = ml 420 156 221 54 -131 221 41 2 131 -31 2 54 131 1S6 -221 -131 -312-221 41 2 ;CK1 = The another case is axial and bending vibration of a bar. To this end, we will use a way which called super position method. At this method we consider the bar is consisting axial and bending vibrations, so we seperate these vibrations and derive mass and stiffness matrix I“, iUCt) i.lK%-t> Cx/fc)6 jr uwjiyi^ X Fig. 4. Derivation Of Axial 4 Bending Stiffness And Mass Matrix If we consider the bar as consisting only axial vibration the mass and stiffness matrix are as follows ; xiv -C140 If we consider the bar Is consisting only bending vibra tion, the mass and stiffness matrix are as follows ; C153 If we sum eq. CI 43 and C15D using superpozision method, t M 3 = C M 1+ i M ] b e £M3' ml 420 140 O O 70 O O t K 3 = t K O 156 221 O 54 -131 -31 221 41 2 O 131 2 70 O O O 54 131 140 O O 156 O -221 O -131 -31 2 O -221 41 2 ]+ C K 3 CK3 El Cl/rD O O O 121 61 C-l/rD O O -12 O 61 O 61 41 : O -61 21 : c-ı/rr o o ClXr32 O O -12 -61 O 12 -61 O 61 21 ; O -61 41 C163 C where r is rotation radius of cross section area r = I/A> - xv -Analysis Of Elastic Mechanism The finite element theory of structural analysis has been applied to modelling an elastic linkage using force method as well as displacement method. The lumped parame ter approach has also been a useful technique. The rigid-body kinematic analysis of the linkage is dependent on its input motion characteristic as well as on its geometry. Due to elastic deformations of the members, the effective link lengths as well as their resulting angu lar positions vary from those of a purely rigid linkage. If small elastic deformations are assumed» the effects of elastic deformations on the rigid-body kinematic equations are negligible. The assumption that the elastic motions of the members do not affect the rigid-body kinematics of the linkage, has been made extensively by researches. Each member of linkage is idealized by one beam element only. This, in conjuction with the assumed mode shapes, permits at most one inflection point across the length of the element and limits the number of modes can be correct ly computed. Also, one-element Idealization results in lost of higher model frequencies. While one element idea lization is used for simplicity, mul ti -el ement idealization of each member must be used for a more accurate modelling. In the force finite el ement method, the internal for ces are assumed to be unknowns and equilibrium conditions are used to generate equations. Additional equations are sometimes essential and are derived with the help of com patibility conditions. The displacement finite element met hod, on the other hand, utilizes the nodal displacements as the unknowns of the problem. The compatibility conditi ons in and among the elements are initially satisfied. The equations are expressed in terms of these nodal disp lacements utilizing the equilibrium conditions at each no dal point. The displacement finite element method formu lation is in general found to be much simpler for a larg er number of structural problems. For this reason, this latter method will be used here to develop the equations of motion. A general beam element is shown in fig. 3 in two fra mes of reference. These are fixed (O-X-Y) and the rotating (o-x-y) frames, the rotating frame of reference being such that the x-axis parallel to the beam element axis. The elastic deformations of the beam element may be complete ly described by six nodal displacements, Ui through U = S at and < U > = S r t X A y*A ©.. X”B e“ B i »1,2..6 C22Z> By combining eqns C19D-C21D > C23D Equation C23D may be rewritten as ; < U*”> = < U*> + < U*> + < U*> + < u“> + < U”> a r net C24D Where the vectors from left to right represent the absolute, rigid-body, elastic, normal, coriolis, and tan gential accelerations respectively. The last three vec tors in eq. C233 are the rigid-body and elastic motion co upling terms. If small elastic displacements are considered, and if the rigid-body velocity and acceleration generally small compared to the velocity and acceleration of the elastic no dal deformation, then in eq. C24D the product terms in vec tors , < U'*> and < \J"> are considered small compared en t to corresponding terms in < U > +< U>. The coriolis, r normal and tangentianal acceleration terms are therefore neglected. These have also been neglected by researc hers and in the present work. Equation C24D may be modi fied as ; =+ a r C2SD Vibration Analysis Of Elastic Crank-Connecting Rod Slider Mechanism By Using Finite Element Method Because the crank -connecting rod mechanism is consis ting axial and bending displacements, we can use mass and - xvin -stiffness matrix of axi al -bendi ng displacement state derivating in eq. 16. Also, crank has a rotating movement, therefore the use of rotation matrix is needed C 4 3. In the present work, the following matrices are obtained as ; ££Ü££&^^^Ş£ W777777777777777, Fig 8. Elastic Crane-Rod Mechanism CM3=i:R3CM3CR3 C26Z> where C R 3 - C M 3 = C27Z> C28D EK3=ER3EK3ER3 r CXl/rD2+12y2 X^Cl/r32-12X^ C29D £ K 3 » r El -61 £i -CXl/r32-12y' -XyCl/rD2+12XfJ -61fci X^jCl/rD -12Xy C^Jİ/rD2+12X2 61 X -XyCl/rD2+12Xfj -C^l/rD2-12X2 61 X -61 £* 61 X 41 2 61 fj -61X 21 2 - xıx -where C M ed stiffness matrix In order to add the additional analytical strains, the following discussion is necessary. The beam element transverse displacements are assumed as independent of axial displacements. In reality, however, a compressive axial deformation will tend to increase transverse defor mation of the beam, thus effectively decreasing its trans verse stiffness. Similiarly.a tensial axial deformation will have opposite effect. An approximate method to cor rect this defficiency, in the stiffness matrix E K 3, a geometric matrix E K 3, is simply added to the matrix E K above. The stiffness matrix C ment as ; K 3 is given for a beam ele- g E K 3 g C31D By using eqns C30D and C31D we can derive stiffness mat rix for connecting-rod slider mechanism as follows ; [ K 3 = - XX-X^C 1 /r Z> 2 +1 2Xfci -61 y -CXl/rZ>Z-12y2 - 1.2P12/EI 61X + 0.1P13/EI -61X - 0. 1P13/EI 21 2 - Pl*/30 X^Cl/r32-12Xy 61fc* Cyl/rD2+12X2 + 1.2P12/EI -61X - 0. 1P13/EI -61 X - O. IPlVEI 41 2 + 2PlVl5EI C32Z) C M 3 < U*> + C K ] < U> =< F> C33Z> s a a where CM] and I K 3 are system matrix of* connecting- s & rod mechanism CCM3=CM3 D a r Conclusion t In the work, the related numerical values are calcu lated and diagrams relating to the vibration frequencies and mode shapes: are constructed to demonstrate the effects of several parameters such as elasticity of crank and coup ler, mass of piston and the stiffness of the spring bet ween piston and cylinder. It has been also, shown that crank and coupler elasticity longitudinal deflection and forces are important parameters in the design of crank-con necting rod-slider mechanism. - xxx -

Benzer Tezler

  1. Dynamic analysis of planar flexible mechanisms

    Düzlemsel esnek mekanizmaların dinamik analizi

    DUYGU ÇÖMEN

    Yüksek Lisans

    İngilizce

    İngilizce

    2011

    Makine Mühendisliğiİzmir Yüksek Teknoloji Enstitüsü

    Makine Mühendisliği Ana Bilim Dalı

    DOÇ. DR. BÜLENT YARDIMOĞLU

  2. Elastik uzuvlu makinaların dinamik kararlılığı

    Dynamics stability of mechanisms with elastic members

    GÖKHAN BULUT

    Yüksek Lisans

    Türkçe

    Türkçe

    2003

    Makine Mühendisliğiİstanbul Teknik Üniversitesi

    Makine Mühendisliği Ana Bilim Dalı

    DOÇ. DR. ÖZGÜR TURHAN

  3. Soğutucularda kullanılan hermetik pistonlu kompresörler için özgün krank biyel mekanizması tasarımı ve imalatı

    Design and manufacturing of novel crank connecting rod mechanism for hermetic piston compressors used in refrigerators

    MOHAMED MUSTAFA AHMED IZZAT AL-DALAWI

    Yüksek Lisans

    Türkçe

    Türkçe

    2024

    Endüstri Ürünleri TasarımıGazi Üniversitesi

    Otomotiv Ana Bilim Dalı

    PROF. DR. MELİH OKUR

  4. Mekanizmaların sonlu elemanlar yöntemiyle dinamik analizi

    Dynamic analysis of mechanisms with finite elements method

    SABRİ İLKSÖZ

    Yüksek Lisans

    Türkçe

    Türkçe

    2008

    Makine MühendisliğiYıldız Teknik Üniversitesi

    Makine Mühendisliği Bölümü

    YRD. DOÇ. DR. CİHAN DEMİR

  5. Dynamic investigation of diesel engine gear train system under fluctuated torque excitation

    Dizel motor zamanlama dişli sistemi dinamiğinin salınımlı tork etkisi altında incelenmesi

    OSMAN SUBAŞI

    Yüksek Lisans

    İngilizce

    İngilizce

    2017

    Makine Mühendisliğiİstanbul Teknik Üniversitesi

    Makine Mühendisliği Ana Bilim Dalı

    YRD. DOÇ. DR. OSMAN AKIN KUTLAR