63 VI. VERIFICAREA LA REZISTENTA A ARBORELUI COTIT Verificarea presiunii specifice maxime si medii Pe baza calculelor la CCMAI au fost stabilite valorile maximă şi medie ale for ţ ei care acţ ionează înlag maneton. Prin urmare, presiunea specifica maximă pe fusul maneton este: p Mmax R Mmax 10 d M 100 × l M × 100 × := p Mmax 105.473 = daN cm 2 La rândul ei, presiunea specifică medie are valoarea: p Mmed R Mmed 10 d M 100 × l M × 100 × := p Mmed 37.044 = daN cm 2 Ambele valori nu depăşesc limitele admisibile indicate: p Mmax 100...300 := daN cm 2 p Mmed 30...100 := daN cm 2 R p.max R p.max 0 ¬ R p.max max R pmax j 1 + R pmax j , æ è ö ø ¬ R p.max max R pmax j 1 + R pmax j , æ è ö ø < if j 5 5 .. Î for R p.max := R p.max 6.047 10 4 ´ = N p pmax R p.max 10 d p 100 × l p × 100 × := p pmax 140.933 = daN cm 2 R p.med R p.med 0 ¬ R p.med max R pmed j 1 + R pmed j , æ è ö ø ¬ R p.med max R pmed j 1 + R pmed j , æ è ö ø < if j 5 5 .. Î for R p.med := R p.med 2.714 10 4 ´ = N
Proiect Vibratii Navale : Pef-448 kw , 6 cilindrii in linie
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
63
VI. VERIFICAREA LA REZISTENTA A ARBORELUI COTIT
Verificarea presiunii specifice maxime si medii
Pe baza calculelor la CCMAI au fost stabilite valorile maximă şi medie ale forţei care acţionează înlagărulmaneton. Prin urmare, presiunea specifica maximă pe fusul maneton este:
pMmax
RMmax
10
dM 100× lM× 100×:= pMmax 105.473=
daN
cm2
La rândul ei, presiunea specifică medie are valoarea:
pMmed
RMmed
10
dM 100× lM× 100×:= pMmed 37.044=
daN
cm2
Ambele valori nu depăşesc limitele admisibile indicate:
pMmax 100...300:=daN
cm2pMmed 30...100:=
daN
cm2
Rp.max Rp.max 0¬
Rp.max max Rpmaxj 1+Rpmaxj
, æè
öø
¬ Rp.max max Rpmaxj 1+Rpmaxj
, æè
öø
<if
j 5 5..Îfor
Rp.max
:=
Rp.max 6.047 104´= N
ppmax
Rp.max
10
dp 100× lp× 100×:=
ppmax 140.933=daN
cm2
Rp.med Rp.med 0¬
Rp.med max Rpmedj 1+Rpmedj
, æè
öø
¬ Rp.med max Rpmedj 1+Rpmedj
, æè
öø
<if
j 5 5..Îfor
Rp.med
:=
Rp.med 2.714 104´= N
64
daN
cm2ppmed
Rp.med
10
dp 100× lp× 100×:= ppmed 63.24=
Valoarea nu depăşeşte limitele admisibile indicate:
ppmed 20...75:=daN
cm2
Verificarea la incalzire
Pentru verificarea fusurilor la încălzire, se determină iniţial valoarea coeficientului de corecţie ξ.
ξ ξ 1.052¬ λd 0.25=if
ξ
:=
ξ 1.052=
În aceste condiţii coeficientul de uzură pentru cele două fusuri înregistrează valorile:
KM pMmed ξπ dM× n×
60×
æçè
ö÷ø
3
×:= KM 249.132=
Kp ppmedπ dp× n×
60
æçè
ö÷ø
3
×:= Kp 243.399=
Valoarile nu depăşeşc limitele admisibile indicate, pentru motoare rapide:
KM 150...250:= Kp 150...250:=
Valorile momentelor de răsucire care solicită fusurile palier
Pentru verificarea la oboseală a fusurilor palier, au fost calculeate iniţ ial momentele de răsucire ale acestoraşi au fost apoi calculate valorile extreme şi amplitudinea:
Mr.II.max Mr.II.max 0¬
Mr.II.max max Mr.IIj 1+Mr.IIj
, æè
öø
¬ Mr.II.max max Mr.IIj 1+Mr.IIj
, æè
öø
<if
j 0 71..Îfor
Mr.II.max
:=
Mr.II.min Mr.II.min Mr.II.max¬
Mr.II.min min Mr.IIj 1+Mr.IIj
, æè
öø
¬ Mr.II.min min Mr.IIj 1+Mr.IIj
, æè
öø
³if
j 0 71..Îfor
Mr.II.min
:=
AMr.II Mr.II.max Mr.II.min-:=
Mr.III.max Mr.III.max 0¬
Mr.III.max max Mr.IIIj 1+Mr.IIIj
, æè
öø
¬ Mr.III.max max Mr.IIIj 1+Mr.IIIj
, æè
öø
<if
j 0 71..Îfor
Mr.III.max
:=
Mr.III.min Mr.III.min Mr.III.max¬
Mr.III.min min Mr.IIIj 1+Mr.IIIj
, æè
öø
¬ Mr.III.min min Mr.IIIj 1+Mr.IIIj
, æè
öø
³if
j 0 71..Îfor
Mr.III.min
:=
AMr.III Mr.III.max Mr.III.min-:=
Mr.IV.max Mr.IV.max 0¬
Mr.IV.max max Mr.IVj 1+Mr.IVj
, æè
öø
¬ Mr.IV.max max Mr.IVj 1+Mr.IVj
, æè
öø
<if
j 0 71..Îfor
Mr.IV.max
:=
Mr.IV.min Mr.IV.min Mr.IV.max¬
Mr.IV.min min Mr.IVj 1+Mr.IVj
, æè
öø
¬ Mr.IV.min min Mr.IVj 1+Mr.IVj
, æè
öø
³if
j 0 71..Îfor
Mr.IV.min
:=
AMr.IV Mr.IV.max Mr.IV.min-:=
69
Mr.V.max Mr.V.max 0¬
Mr.V.max max Mr.Vj 1+Mr.Vj
, æè
öø
¬ Mr.V.max max Mr.Vj 1+Mr.Vj
, æè
öø
<if
j 0 71..Îfor
Mr.V.max
:=
Mr.V.min Mr.V.min Mr.V.max¬
Mr.V.min min Mr.Vj 1+Mr.Vj
, æè
öø
¬ Mr.V.min min Mr.Vj 1+Mr.Vj
, æè
öø
³if
j 0 71..Îfor
Mr.V.min
:=
AMr.V Mr.V.max Mr.V.min-:=
Mr.VI.max Mr.VI.max 0¬
Mr.VI.max max Mr.VIj 1+Mr.VIj
, æè
öø
¬ Mr.VI.max max Mr.VIj 1+Mr.VIj
, æè
öø
<if
j 0 71..Îfor
Mr.VI.max
:=
Mr.VI.min Mr.VI.min Mr.VI.max¬
Mr.VI.min min Mr.VIj 1+Mr.VIj
, æè
öø
¬ Mr.VI.min min Mr.VIj 1+Mr.VIj
, æè
öø
³if
j 0 71..Îfor
Mr.VI.min
:=
AMr.VI Mr.VI.max Mr.VI.min-:=
Mr.VII.max Mr.VII.max 0¬
Mr.VII.max max Mr.VIIj 1+Mr.VIIj
, æè
öø
¬ Mr.VII.max max Mr.VIIj 1+Mr.VIIj
, æè
öø
<if
j 0 71..Îfor
Mr.VII.max
:=
Mr.VII.min Mr.VII.min Mr.VII.max¬
Mr.VII.min min Mr.VIIj 1+Mr.VIIj
, æè
öø
¬ Mr.VII.min min Mr.VIIj 1+Mr.VIIj
, æè
öø
³if
j 0 71..Îfor
Mr.VII.min
:=
AMr.VII Mr.VII.max Mr.VII.min-:=
Mr.II.max 2.653 103´= Mr.II.min 2.062- 103
´= AMr.II 4.715 103´=
Mr.III.max 4.304 103´= Mr.III.min 2.341- 103
´= AMr.III 6.645 103´=
Mr.IV.max 2.892 103´= Mr.IV.min 2.3- 103
´= AMr.IV 5.192 103´=
Mr.V.max 5.492 103´= Mr.V.min 1.575- 103
´= AMr.V 7.068 103´=
70
Mr.VI.max 7.143 103´= Mr.VI.min 1.899- 103
´= AMr.VI 9.042 103´=
Mr.VII.max 5.731 103´= Mr.VII.min 3.079- 103
´= AMr.VII 8.81 103´=
Modulul de rezistenţă polar al secţiunii transversale a fusului palier are valoare:
WpPπ dp
3×
16:= WpP 1.053 10 4-
´= m3
Rezultă valorile extreme ale efortului unitar de răsucire:
τP.maxMr.VI.max
WpP10 5-
×:= τP.max 678.266=daN
cm2
τP.minMr.VI.min
WpP10 5-
×:= τP.min 180.273-=daN
cm2
Amplitudinea şi media acestor eforturi sunt:
τPvτP.max τP.min-
2:= τPv 429.269=
daN
cm2
τPmτP.max τP.min+
2:= τPm 248.997=
daN
cm2
Conform indicaţiilor din manual pentru materialul de construcţie al arborelui cotit (OLC45), se adoptărezistenţa la oboseală precum şi valorile coeficienţilor de calcul astfel:
γ 0.9:= Rβ.ε 2:=τ_1 2300:=daN
cm2
ψτ 0.09:=Rβ.ε
βkτετ
=
Cu ajutorul lor se determină coeficientul de siguranţă la oboseală pentru fusul palier:
cPτ_1
Rβ.εγ
τPv× ψτ τPm×+
:=
cP 2.356=
71
cp 2...3:=Valoarea obţinută este superioară celei minim admisibile indicată mai jos:
Verificarea la oboseala pentru fusul manetonPentru verificarea la oboseală a fusului maneton este necesară determinarea iniţială a forţei radiale Zj,normalălameneton. Aceasta se obţine prin însumarea algebrică a forţei radiale Zβ cu forţele de inerţie alefusului maneton şi ale masei bielei aferentă acestuia calculate anterior. Rezultatele vor fi prezentatecentralizat in tabelul de mai jos.Zj
ZjjZBj
Fibm+ FR+¬
j 0 72..Îfor
Zj
:=
kN
Rezultă astfel reacţiunea din reazemul stâng al fusului, rezultatele find cuprinse în tabelul de mai jos:Zsj
Zsjj12
Zjj× Fbr+ Fcg+¬
j 0 72..Îfor
Zsj
:=
kN
Totodată s-a calculat momentul de încovoiere în planul cotului, cu relaţia de mai jos, rezultatele findcuprinse centralizat în tabel.MZ
Pentru determinarea momentului de încovoiere din planul tangenţial se utilizează valorile forţei tangenţialeTj. Pe baza lao a fost determinata variaţia reacţiunii tangenţiale.
TjTjj
Tj¬
j 0 72..Îfor
Tj
:= Tsj
Tsjj
Tjj2
¬
j 0 72..Îfor
Tsj
:=
kN kN
Aceste reacţiuni determină momentul de încovoiere.
Fusul maneton fiind prevăzut cu orificiu de ungere, momentul rezultant se determină cu ajutorul relaţiei:
MoσMoσj
MZjcos deg ϕ×( )× MTj
sin deg ϕ×( )×-¬
j 0 72..Îfor
Moσ
:=ϕ ϕM17
:=
Rezultatele calculelor sunt cuprinse în tabelul de mai jos, rezultând că valoarea maximă a momentuluiîncovoietor total se înregistrează în cazul prezeţei contragreutăţii.
Conform celor prezentate se consideră un ciclu simetric de încărcare, în care Mσmin=- Mσmax. Înaceste condiţii, modulul de rezistenţă polar al secţiunii transversale a fusului manetor.
WpMπ dM
3×
16:=
WpM 1.618 10 4-´= m3
Rezultă eforturile unitare extreme.
σmaxMoσ.max
WpM10 5-
×:= σmax 46.396=daN
cm2
σminMoσ.min
WpM10 5-
×:= σmin 141.399-=daN
cm2
Amplitudinea şi media acestor eforturi sunt:
σvσmax σmin-
2:= σv 93.898=
daN
cm2
σmσmax σmin+
2:= σm 47.502-=
daN
cm2
78
Pe baza indicaţiilor de mai sus pentru fusuri din OLC45 se adoptă:
σ_1 2300:=daN
cm2
Totodată se adoptă şi parametrii:
γ 0.9= βkσ 2:= εσ 0.9:= ψσ 0.1:=
Rezultă coeficientul de siguranţă la oboseală pentru solicitarea de încovoiere a manetonului.
cσσ_1
βkσγ εσ×
σv× ψσ σm×+
:= cσ 10.128=
Pentru verificarea solicitarii de torsiune este luat în consideraţie fusul maneton al celui mai solicitat cotal arborelui cotit. Pe baza datelor calculate au fost utilizate valorile momentului de intrare şi ale forţeitangenţiale. Cu ajutorul lor poate fi calculat momentul de torsiune al fusului maneton.
Se calculează eforturile unitare corespunzătoare momentelor de torsiune extreme.
τmaxMτ6.max
WpM10 5-
×:= τmax 397.886=daN
cm2
τminMτ6.min
WpM10 5-
×:= τmin 148.234-=daN
cm2
Amplitudinea şi media acestor eforturi sunt:
τvτmax τmin-
2:= τv 273.06=
daN
cm2
τmτmax τmin+
2:= τm 124.826=
daN
cm2
Prin urmare menţinându-se valorile adoptate la verificarea fusului palier, coeficientul de siguranţă laoboseală pentru solicitarea de torsiune va fi:
cττ_1
Rβ.εγ
τv× ψτ τm×+
:= cτ 3.721=
Rezultă astfel coeficientul global de siguranţă la oboseală a fusului maneton:
81
cMcσ cτ×
cσ2 cτ
2+
:=
cM 3.493=
Pentru verificarea la rezistenţă a braţului arborelui cotit, se pleacă de la valorile maximă şi minimă aleeforturilor unitare de încovoiere şi compresiune. Pe baza valorilor extreme ale reacţiunii radiale rezultăeforturile unitare:
Zsj.max Zsj.max 0¬
Zsj.max max Zsjj 1+Zsjj
, æè
öø
¬ Zsj.max max Zsjj 1+Zsjj
, æè
öø
<if
j 0 71..Îfor
Zsj.max
:=
Zsj.min Zsj.min Zsj.max¬
Zsj.min min Zsjj 1+Zsjj
, æè
öø
¬ Zsj.min min Zsjj 1+Zsjj
, æè
öø
³if
j 0 71..Îfor
Zsj.min
:=
σbr.max Zsj.max3 lp×
br hbr2
×
1br hbr×
+æççè
ö÷÷ø
× 10 5-×:= σbr.max 207.396=
daN
cm2
σbr.min Zsj.min3 lp×
br hbr2
×
1br hbr×
+æççè
ö÷÷ø
× 10 5-×:= σbr.min 561.737-=
daN
cm2
Amplitudinea şi media eforturilor vor fi:
σbr.vσbr.max σbr.min-
2:= σbr.v 384.566=
daN
cm2
σbr.mσbr.max σbr.min+
2:= σbr.m 177.17-=
daN
cm2
Cu ajutorul acestor valori anterior adoptate (în cadrul verificării fusului maneton), rezultă coeficientul desiguranţă la oboseală pentru solicitările de încovoiere şi compresiune:
cbr.σσ_1
βkσγ εσ×
σbr.v× ψσ σbr.m×+
:= cbr.σ 2.468=
Pentru verificarea solicitării de torsiune a braţului este necesar să se determine mai întâi modulul derezistenţă al acestuia. În acest scop avem:
¬:=
82
Tj.max Tj.max 0¬
Tj.max max Tjj 1+Tjj
, æè
öø
¬ Tj.max max Tjj 1+Tjj
, æè
öø
<if
j 0 35..Îfor
Tj.max
:=
Tj.min Tj.min Tj.max¬
Tj.min min Tjj 1+Tjj
, æè
öø
¬ Tj.min min Tjj 1+Tjj
, æè
öø
³if
j 0 35..Îfor
Tj.min
:=
brhbr
5= din nomogramă
avem apoi coeficientul Saint-Venant K 0.292:=
Prin urmare modulul de rezistenţă va fi:
Wd K br× hbr2
×:= Wd 5.158 10 5-´= m3
Utilizând valorile extreme ale forţei tangenţiale rezultă valorile extreme ale efortului unitar de torsiune lacare este supus braţul arborelui cotit:
τbr.maxlp Tj.max×
4 Wd×10 5-
×:= τbr.max 63.809=daN
cm2
τbr.minlp Tj.min×
4 Wd×10 5-
×:= τbr.min 75.417-=daN
cm2
Amplitudinea şi media acestor eforturi sunt:
τbr.vτbr.max τbr.min-
2:= τbr.v 69.613=
daN
cm2
τbr.mτbr.max τbr.min+
2:= τbr.m 5.804-=
daN
cm2
Pe baza valorilor adoptate încadrul calculului fusului palier, coeficientul de siguranţă la oboseală pentrusolicitarea de torsiune va fi:
cbr.ττ_1
Rβ.εγ
τbr.v× ψτ τbr.m×+
:= cbr.τ 14.918=
În aceste condiţii coeficientul global de siguranţă la oboseală are valorile:
cbrcbr.σ cbr.τ×
cbr.σ2 cbr.τ
2+
:=cbr 2.435=
84
VII. CALCULUL PRELIMINAR AL VIBRAŢIILOR TORSIONALE
Pe baza calculelor anterioare şi a documentaţiei tehnice a motorului de referinţă, poate fi stabilit sistemulosculant echivalent din figura de mai jos.Pentru determinarea momentului de inerţie al unui cot al arborelui cotit, se calculeazăiniţial momentul deinerţie al fusului palier.
Jpπ32
dp4
× lp× ρOLC45× 1.762 10 3-´=:= Nms2
La rândul său, momentul de inerţie al fusului maneton redus la axa de rotaţie al arborelui cotit arevaloarea:
JM0 JM mM R2×+:=
JM0π32
dM2
× lM× ρ× dM 8 R2×+æ
èöø× 1.821 10 8-
´=:= Nms2
Sistemul oscilant echivalent al instalaţiei de propulsie
Momentul de inerţie al braţelor fără contragreutăţi:
Jb 1043:= Nms2
J'b 478:= Nms2
Momentul de inertie al bratului fara contragreutate, redus la axa de rotate este:
J'b0 J'b mbr rGbr 10 3-×æ
èöø
2×+ 478=:= Nms2
Pentru braţul cu contragreutate, centrul de masă al ansamblului braţ - contragreutate va fi amplasat ladistanţa următoare faţă de axa de rotaţie:
85
mrbrcgmbr rGbr× mcg rGcg×-
mbr mcg+6.215 10 3-
´=:=
Prin urmare, pentru acest tip de braţ, momentul de inerţie redus la axa de rotaţie este:
Jb0 Jb mbr mcg+( ) rbrcg2
×+ 1.043 103´=:= Nms2
Rezulta astfel momentele de inertie ale coturilor arborelui cotit:-pentru coturile 1, 3, 5 prevazute cu contragreutati:
Jcot Jp JM0+ 2 Jb0×+ 2.086 103´=:= Nms2
-pentru coturile 2, 4,6 prevazute fara contragreutati:
J'cot Jp JM0+ 2 J'b0×+ 956.002=:= Nms2
Pentru stabilirea momentelor de inertie ale discurilor (volantilor) arborelui cotit, se mai determinasi momentul de inertie al maselor in miscare, redus la axa de rotatie:
Jm0 mbm 0.5 mit×+( ) R2× 0.044=:= Nms2
În consecință, momentele de inertie ale discurilor oscilant vor fi:J1=J3=J5 si J2=J4=J6
J1 Jcot Jm0+ 2.086 103´=:= Nms2
J2 J'cot Jm0+ 956.046=:= Nms2
Pentru volantul motorului, plecand de la dimensiunile sale principale, rezulta momentul deinertie:Dv 0.976:= m
hv 0.42:= m
Jvπ32
Dv4
× hv× ρOLC45× 291.839=:= Nms2
Pentru determinarea momentului de inerție al elicei, se recurge la o metodă simplificată decalcul. Prin prelucrarea diagramelor Troost, diametrul elicei este dat de relația:
Del 15.25Pef
0.2 103×
n0.6× 1.979 103
´=:= mm
Grosimea celor patru pale ale elicei variază între 0.3 și 2.8% din diametrul Del. Substituimelicea cu un disc cu următoarea grosime:
hel 1.2 10 2-× Del× 23.75=:= mm
Având în vedere că elicea este confecționată din bronz cu o densitate mai mare, avem:
ρel 8160:=Kg
m3
86
Jelπ32
Del4
× hel× ρel× 10 15-× 291.954=:= Nms2
Arborele echivalent are diametrul egal cu cel al fusului palier, iar arborele port -elice are diametruldel=180 mm. Rezultă astfel și rigiditățile porțiunilor de arbore ale sistemului de mai sus:k12=k23=k34=k45=k56
del 0.180:= m
G 8.4 1010×:= -modulul de elasticitate
k12π32
Gl
× dp4
× 1.805 106´=:= Nm
Pentru a determina lungimile liniei de arbori se adoptă următoarele mărimi:
-lungimea unui cot între mijloacele a 2 paliere adiacente:
y 1.365:= m -calculat în proiectul anterior
y2 y 1.365=:= m
y3 0.75 y× 1.024=:= m
y4 y 1.365=:= m
y5 6.5y 8.873=:= m
my6 5.3y 7.234=:=
Calculăm următoarele valori:
lv y2y2
- y3+ 1.706=:= m
le y5 y6+ 16.107=:= m
k6vπ32
Glv
× dp4
× 2.106 105´=:= Nm
kvelπ32
Gle
× del4
× 5.375 105´=:= Nm
J2 J1 2.086 103´=:= Nms2 k12 1.805 106
´= Nm
J3 J1 2.086 103´=:= Nms2 k23 k12 1.805 106
´=:= Nm
J4 J1 2.086 103´=:= Nms2 k34 k12 1.805 106
´=:= Nm
J5 J1 2.086 103´=:= Nms2 k45 k12 1.805 106
´=:= Nm
J6 J1 2.086 103´=:= Nms2 k56 k12 1.805 106
´=:= Nm
87
În aceste condiții, sunt cunoscute toate caracteristicile dinamiceale sistemului oscilant echivalent și pot fideterminate pulsațiile sale proprii. În acest scop, se utilizează metodologia bazată pe așa numita metodăa resturilor și organizată în cadrul tabelei lui Holzer. În calculele de mai jos sunt cuprinse valorile estimateale pulsațiilor proprii, împreună cu valorile corespunzătoare ale momentelor reziduale, pentru fiecarepulsație în parte. Adoptăm valori ale pulsației până avem valoriminime ale momentului rezidual.
ω0I 13.8846694:=
J1 2.086 103´= a1 1:= M12 a1 J1× ω0I
2× 4.022 105
´=:=
J2 2.086 103´= a2 a1
M12k12
- 0.777=:= M23 M12 J2 ω0I2
× a2×+ 7.147 105´=:=
a3 a2M23k23
- 0.381=:=J3 2.086 103´= M34 M23 J3 ω0I
2× a3×+ 8.681 105
´=:=
J4 2.086 103´= a4 a3
M34k34
- 0.099-=:= M45 M34 J4 ω0I2
× a4×+ 8.281 105´=:=
J5 2.086 103´= a5 a4
M45k45
- 0.558-=:= M56 M45 J5 ω0I2
× a5×+ 6.037 105´=:=
J6 2.086 103´= a6 a5
M56k56
- 0.892-=:= M6v M56 J6 ω0I2
× a6×+ 2.447 105´=:=
Jv 291.839= av a6M6vk6v
- 2.054-=:= Mvel M6v Jv ω0I2
× av×+ 1.292 105´=:=
Jel 291.954= ael avMvelkvel
- 2.295-=:= MrI Mvel Jel ω0I2
× ael×+ 0.015=:=
88
ω0II 28.0172893:=
J1 2.086 103´= a1 1:= M12 a1 J1× ω0II
2× 1.637 106
´=:=
J2 2.086 103´= a2 a1
M12k12
- 0.093=:= M23 M12 J2 ω0II2
× a2×+ 1.79 106´=:=
a3 a2M23k23
- 0.898-=:=J3 2.086 103´= M34 M23 J3 ω0II
2× a3×+ 3.188 105
´=:=
J4 2.086 103´= a4 a3
M34k34
- 1.075-=:= M45 M34 J4 ω0II2
× a4×+ 1.441- 106´=:=
J5 2.086 103´= a5 a4
M45k45
- 0.277-=:= M56 M45 J5 ω0II2
× a5×+ 1.894- 106´=:=
J6 2.086 103´= a6 a5
M56k56
- 0.773=:= M67 M56 J6 ω0II2
× a6×+ 6.291- 105´=:=
Jv 291.839= av a6M6vk6v
- 0.389-=:= Mvel M6v Jv ω0II2
× av×+ 1.556 105´=:=
Jel 291.954= ael avMvelkvel
- 0.679-=:= MrII Mvel Jel ω0II2
× ael×+ 0.016-=:=
ω0III 34.7876587:=
J1 2.086 103´= a1 1:= M12 a1 J1× ω0III
2× 2.524 106
´=:=
J2 2.086 103´= a2 a1
M12k12
- 0.398-=:= M23 M12 J2 ω0III2
× a2×+ 1.519 106´=:=
a3 a2M23k23
- 1.24-=:=J3 2.086 103´= M34 M23 J3 ω0III
2× a3×+ 1.61- 106
´=:=
J4 2.086 103´= a4 a3
M34k34
- 0.348-=:= M45 M34 J4 ω0III2
× a4×+ 2.488- 106´=:=
J5 2.086 103´= a5 a4
M45k45
- 1.03=:= M56 M45 J5 ω0III2
× a5×+ 1.133 105´=:=
J6 2.086 103´= a6 a5
M56k56
- 0.968=:= M67 M56 J6 ω0III2
× a6×+ 2.556 106´=:=
Jv 291.839= av a6M6vk6v
- 0.194-=:= Mvel M6v Jv ω0II2
× av×+ 2.003 105´=:=
Jel 291.954= ael avMvelkvel
- 0.567-=:= MrIII Mvel Jel ω0III2
× ael×+ 0.015=:=
89
Conform datelor din calculele anterioare, pentru valorile minime ale momentelorreziduale (resturilor), pot fi stabilite pulsațiile proprii ale sistemului oscilant echivalent:
ω0I 13.885= s 1- ω0II 28.017= s 1- ω0III 34.788= s 1-
In încheierea acestui calcul preliminar la vibrațiilor torsionale, s-a realizat stabilirea regimurilor derezonanță, respectiv precizarea turațiilor critice corespunzătoare acestor regimuri. În acest scop, înfigura de mai jos a fost trasată diagrama pulsație-turație. Astfel, pulsațiile proprii fiind independente deturație, ele sunt reprezentate prin drepte paralele cu abcisa diagramei. Pe de altă parte, pulsațiileexcitației de diverse ordine v variază liniar cu turația.În continuare se stabilesc valorile turațiilor criticepentru diverse grade γ ale pulsației și diverse ordine v ale excitației. Se folosesc următoarele expresii:
v 1 30..:=
nmin 0.5 n× 1.15 103´=:=
nmax 1.1 n× 2.53 103´=:=
ncrI v( )30 ω0I×
v π×:= min 1-
ncrII v( )30 ω0II×
v π×:= min 1-
ncrIII v( )30 ω0III×
v π×:= min 1-
Turaţile critice pentru diverse grade (I, II, III) ale pulsaţiei şi diverse ordine v ale excitaţiei.
Calculul valorii momentului motor produs de forta de presiune a gazelor
Momentul motor mediu
Mpm 106Rπ D2
×4
×172
×
1
72
j
pTjå=
× 189.99=:= Nm
D1 1.4:= D5 0.3:= D9 0.7:=
D2 1.3:=daN
cm2D6 0.2:=
daN
cm2D10 0.7:=
daN
cm2D3 0.9:= D7 1.5:= D11 0.6:=
D4 0.6:= D8 1.1:= D12 0.4:=
Mp1 106π D2× R×
D14
× 1.203 103´=:= Mp9 106π D2
× R×D14
× 1.203 103´=:=
Mp2 106 D2× R×
D24
× 355.469=:= Mp10 106 D2× R×
D24
× 355.469=:=
Nm
Mp3 106 D2× R×
D34
× 246.094=:= Mp11 106 D2× R×
D34
× 246.094=:=
Nm
Mp4 106 D2× R×
D44
× 164.063=:= Mp12 106 D2× R×
D44
× 164.063=:=
Mp5 106π D2× R×
D14
× 1.203 103´=:=
Mp6 106 D2× R×
D24
× 355.469=:= Nm
Mp7 106 D2× R×
D34
× 246.094=:=
Mp8 106 D2× R×
D44
× 164.063=:=
Deci valoarea momentului motor este
Mp Mpm Mp1 sin 1 ω× τ× θ0+( )×+ Mp2 sin 2 ω× τ× θ0+( )×+ Mp3 sin 3 ω× τ× θ0+( )×+ Mp4 sin 4 ω× τ× θ0+( )×+ Mp5 sin×+:=
Mp 189.99= Nm
Calculul momentului produs de forta de inertie a maselor cu miscare de translatie
Valorile coeficientilor armonici se determina astfel
E114λd×
116
λd3
×+15512
λd5
×+ 0.064=:=
114
E21-
2132
λd4
×-132
λd6
×- 0.5-=:=
E33-
4λd×
9-32
λd3
×+81-
512λd
5×+ 0.192-=:=
E41-
4λd
2×
1-8
λd4
×+1-
16λd
6×+ 0.016-=:=
E5532
λd3
×75512
λd5
×+ 2.584 10 3-´=:=
E6332
λd4
×332
λd6
×+ 3.891 10 4-´=:=
Amplitudinea momentului produs de forta de inertie a maselor cu miscare de translatie
Mit1 mit R2× ω2
× E1× 221.528=:= Nm
Mit2 mit R2× ω2
× E2× 1.745- 103´=:= Nm
Mit3 mit R2× ω2
× E3× 669.932-=:= Nm
Mit4 mit R2× ω2
× E4× 56.262-=:= Nm
Mit5 mit R2× ω2
× E5× 9.015=:= Nm
Mit6 mit R2× ω2
× E6× 1.357=:= Nm
Calculul momentului total de excitatie
Ψ1 0:= Ψ5 0:= Ψ6 0:= Ψ2 π:= Ψ3 π:= Ψ4 π:=
M1 Mp1 sin θ1( )× Mit1 sin Ψ1( )×+( )2 Mp1 cos θ1( )× Mit1 cos Ψ1( )×+( )2+ 1.342 103
´=:=
NmM2 Mp2 sin θ2( )× Mit2 sin Ψ2( )×+( )2 Mp2 cos θ2( )× Mit2 cos Ψ2( )×+( )2+ 1.839 103
´=:=
M3 Mp3 sin θ3( )× Mit3 sin Ψ3( )×+( )2 Mp3 cos θ3( )× Mit3 cos Ψ3( )×+( )2+ 746.741=:=
M4 Mp4 sin θ4( )× Mit4 sin Ψ4( )×+( )2 Mp4 cos θ4( )× Mit4 cos Ψ4( )×+( )2+ 200.066=:=
NmM5 Mp5 sin θ5( )× Mit5 sin Ψ5( )×+( )2 Mp5 cos θ5( )× Mit5 cos Ψ5( )×+( )2+ 1.212 103
´=:=
M6 Mp6 sin θ6( )× Mit6 sin Ψ6( )×+( )2 Mp6 cos θ6( )× Mit6 cos Ψ6( )×+( )2+ 356.299=:=
M7 Mp7 sin θ7( )×( )2 Mp7 cos θ7( )×( )2+ 246.094=:=
NmM8 Mp8 sin θ8( )×( )2 Mp8 cos θ8( )×( )2+ 164.063=:=
115
M9 Mp9 sin θ9( )×( )2 Mp9 cos θ9( )×( )2+ 1.203 103
´=:=
M10 Mp10 sin θ10( )×( )2 Mp10 cos θ10( )×( )2+ 355.469=:=
M11 Mp11 sin θ11( )×( )2 Mp11 cos θ11( )×( )2+ 246.094=:= Nm
M12 Mp12 sin θ12( )×( )2 Mp12 cos θ12( )×( )2+ 164.063=:=
Avand in vedere ordinea de aprindere a motorului, 1-5-3-6-2-4, avem rangul aprinderii din cilindrii.
j1 1:= j5 2:=
j2 5:= j6 4:=
j3 3:=
j4 6:=
Φc 720:=
RACsω
180 2300×30
1.38 104´=:=
Determinam valorile momentului total de excitatie
Mv1 Mv sin v ω× τ×( )×:= Nm Mv2 Mv sin v ω τ× j2 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
×:= Nm
M11 M1 sin 1 ω× τ×( )× 712.652-=:=M12 M1 sin 1 ω τ× j2 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 230.215=:=
M21 M2 sin 2 ω× τ×( )× 1.548 103´=:=
M22 M2 sin 2 ω τ× j2 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 802.38=:=M31 M3 sin 3 ω× τ×( )× 727.683-=:=
M32 M3 sin 3 ω τ× j2 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.426 103´=:=M41 M4 sin 4 ω× τ×( )× 1.089- 103
´=:=
M51 M5 sin 5 ω× τ×( )× 1.735 103´=:= M42 M4 sin 4 ω τ× j2 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 1.782 103´=:=
M61 M6 sin 6 ω× τ×( )× 768.403-=:= M52 M5 sin 5 ω τ× j2 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.681 103´=:=
M71 M7 sin 7 ω× τ×( )× 74.387-=:= M62 M6 sin 6 ω τ× j2 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.144 103´=:=
M81 M8 sin 8 ω× τ×( )× 323.338-=:=M72 M7 sin 7 ω τ× j2 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 384.274=:=M91 M9 sin 9 ω× τ×( )× 900.818=:=
M82 M8 sin 8 ω τ× j2 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 313.545-=:=M101 M10 sin 10 ω× τ×( )× 261.45-=:=
116
M111 M11 sin 11 ω× τ×( )× 1.157- 103´=:= M92 M9 sin 9 ω τ× j2 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 719.368-=:=
M121 M12 sin 12 ω× τ×( )× 1.736 103´=:= M102 M10 sin 10 ω τ× j2 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 743.855-=:=
M112 M11 sin 11 ω τ× j2 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 459.145-=:=
M122 M12 sin 12 ω τ× j2 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 20.083-=:=
Mv3 Mv sin v ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
×:= Nm Mv4 Mv sin v ω τ× j4 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
×:= Nm
M13 M1 sin 1 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 740.43-=:= M14 M1 sin 1 ω τ× j4 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 688.204=:=
M23 M2 sin 2 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.492- 103´=:= M24 M2 sin 2 ω τ× j4 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 1.581 103´=:=
M33 M3 sin 3 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 416.149-=:= M34 M3 sin 3 ω τ× j4 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 970.771=:=
M43 M4 sin 4 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.437 103´=:= M44 M4 sin 4 ω τ× j4 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 764.366-=:=
M53 M5 sin 5 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.715 103´=:= M54 M5 sin 5 ω τ× j4 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 1.659- 103´=:=
M63 M6 sin 6 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 460.899=:= M64 M6 sin 6 ω τ× j4 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 963.839-=:=
M73 M7 sin 7 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 212.171-=:= M74 M7 sin 7 ω τ× j4 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 47.915-=:=
M83 M8 sin 8 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 360.313=:= M84 M8 sin 8 ω τ× j4 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 248.167-=:=
M93 M9 sin 9 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 587.48=:= M94 M9 sin 9 ω τ× j4 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.006- 103´=:=
M103 M10 sin 10 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 510.611-=:= M104 M10 sin 10 ω τ× j4 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 852.21-=:=
M113 M11 sin 11 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.659- 103´=:= M114 M11 sin 11 ω τ× j4 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 436.72=:=
M123 M12 sin 12 ω τ× j3 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.286- 103´=:= M124 M12 sin 12 ω τ× j4 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 1.585 103´=:=
Mv5 Mv sin v ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
×:= Nm Mv6 Mv sin v ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
×:= Nm
117
M15 M1 sin 1 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 892.358-=:= M16 M1 sin 1 ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 313.33-=:=
M25 M2 sin 2 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 86.264=:= M26 M2 sin 2 ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.058- 103´=:=
M35 M3 sin 3 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 2.016 103´=:= M36 M3 sin 3 ω τ× j6 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 1.78- 103´=:=
M45 M4 sin 4 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 220.669-=:= M46 M4 sin 4 ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 2.043- 103´=:=
M55 M5 sin 5 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.727- 103´=:= M56 M5 sin 5 ω τ× j6 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 1.699- 103´=:=
M65 M6 sin 6 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 183.248=:= M66 M6 sin 6 ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 956.672-=:=
M75 M7 sin 7 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 390.156=:= M76 M7 sin 7 ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 234.323-=:=
M85 M8 sin 8 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 76.696=:= M86 M8 sin 8 ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 97.009=:=
M95 M9 sin 9 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 979.47=:= M96 M9 sin 9 ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 86.797-=:=
M105 M10 sin 10 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 387.543-=:= M106 M10 sin 10 ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 629.691-=:=
M115 M11 sin 11 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.633- 103´=:= M116 M11 sin 11 ω τ× j6 1-( )
Φci
×-éêë
ùúû
×éêë
ùúû
× 1.229- 10´=:=
M125 M12 sin 12 ω τ× j5 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 551.557=:= M126 M12 sin 12 ω τ× j6 1-( )Φci
×-éêë
ùúû
×éêë
ùúû
× 1.601- 10´=:=
Decalajul dintre aprinderi
118
δvj v- j 1-( )×Φci
×:=
δ1j1 1- j1 1-( )×Φci
× 0=:=
δ1j2 1- j2 1-( )×Φci
× 480-=:=
δ1j3 1- j3 1-( )×Φci
× 240-=:=
δ1j4 1- j4 1-( )×Φci
× 600-=:=
δ1j5 1- j5 1-( )×Φci
× 120-=:=
δ1j6 1- j6 1-( )×Φci
× 360-=:=
VERIFICAREA LA REZONANTA A ARBORELUI COTIT
Pentru determinarea amplitudinii momentului de torsiune la rezonanta se pleaca de la regimurile derezonanta, respectiv turatiile critice corespunzatoare acestor regimuri. Obesrvam ca regimurile derezonanta nu se inregistreaza pentru nici o pulsatie.Totusi pentru calcul vom considera armonica 1pentru pulsatia de grad 2 si armonica 1 pentru pulsatia de grad 3.n 2300:=
Calculam amortizarea a, factorul de amplificare A si valorile diferentei de faza a elongatiei miscarii fortateβv [rad]
ωπ n×30
240.855=:=radsec
ordinul 1-II ordinul 1-III
discul nr. 1,3,5 discul nr. 1,3,5
a1ξ1ξ1cr
0.034=:= a'1ξ2ξ1cr
0.04=:=
A1'1
11 ω×ω0II
æçè
ö÷ø
2-
éêêë
ùúúû
2
4 a12
×1 ω×ω0II
æçè
ö÷ø
2×+
:= A'1'1
14 ω×ω0III
æçè
ö÷ø
2-
éêêë
ùúúû
2
4 a'12
×1 ω×ω0III
æçè
ö÷ø
2×+
:=
A1' 0.014= A'1' 1.305 10 3-´=
β1 atan
2 a1× 1 ω×
ω0II
11 ω×ω0II
æçè
ö÷ø
2-
éêêêêêë
ùúúúúúû
8.058- 10 3-´=:= β'2 atan
2 a'1× 1 ω×
ω0III
11 ω×ω0III
æçè
ö÷ø
2-
éêêêêêë
ùúúúúúû
0.012-=:=
discul nr. 2,4,6 discul nr. 2,4,6
122
a1ξ1ξ2cr
0.075=:= a'1ξ2ξ2cr
0.088=:=
A1''1
11 ω×ω0II
æçè
ö÷ø
2-
éêêë
ùúúû
2
4 a12
×1 ω×ω0II
æçè
ö÷ø
2×+
:= A'1''1
11 ω×ω0III
æçè
ö÷ø
2-
éêêë
ùúúû
2
4 a'12
×1 ω×ω0III
æçè
ö÷ø
2×+
:=
A1'' 0.014= A'1'' 0.021=
β1 atan
2 a1× 1 ω×
ω0II
11 ω×ω0II
æçè
ö÷ø
2-
éêêêêêë
ùúúúúúû
0.018-=:= β'2 atan
2 a'1× 1 ω×
ω0III
11 ω×ω0III
æçè
ö÷ø
2-
éêêêêêë
ùúúúúúû
0.026-=:=
volant volant
a1ξ1ξvcr
0.244=:= a'1ξ2ξvcr
0.288=:=
A11
11 ω×ω0II
æçè
ö÷ø
2-
éêêë
ùúúû
2
4 a12
×1 ω×ω0II
æçè
ö÷ø
2×+
:= A'11
11 ω×ω0III
æçè
ö÷ø
2-
éêêë
ùúúû
2
4 a'12
×1 ω×ω0III
æçè
ö÷ø
2×+
:=
A1 0.014= A'1 0.021=
β1 atan
2 a1× 1 ω×
ω0II
11 ω×ω0II
æçè
ö÷ø
2-
éêêêêêë
ùúúúúúû
0.058-=:= β'2 atan
2 a'1× 1 ω×
ω0III
11 ω×ω0III
æçè
ö÷ø
2-
éêêêêêë
ùúúúúúû
0.085-=:=
eliceelice
a'1ξ2ξelcr
0.288=:=a1ξ1ξelcr
0.244=:=
A'11
12 ω×ω0III
æçè
ö÷ø
2-
éêêë
ùúúû
2
4 a'12
×1 ω×ω0III
æçè
ö÷ø
2×+
:=A11
11 ω×ω0II
æçè
ö÷ø
2-
éêêë
ùúúû
2
4 a12
×1 ω×ω0II
æçè
ö÷ø
2×+
:=
A'1 5.239 10 3-´=A1 0.014=
β1 atan
2 a1× 1 ω×
ω0II
11 ω×ω0II
æçè
ö÷ø
2-
éêêêêêë
ùúúúúúû
0.058-=:= β'1 atan
2 a'1× 1 ω×
ω0III
11 ω×ω0III
æçè
ö÷ø
2-
éêêêêêë
ùúúúúúû
0.085-=:=
123
Calculam deformatia unghiulara statica, amplitudinea vibratiei, faza initiala a vibratiei de grad γ si ordinζ, diferenta de faza a vibratiei de grad γ si ordin v, functia caracteristica a vibratiei de grad γ si ordinv, deformatia unghiulara pentru discul ζ si ordinul v, faza initiala a vibratiei de grad γ si ordin v ,diferenta de faza a vibratiei de grad γ si ordin v, functia caracteristica a vibratiei de grad γ si ordin v.