Расчет 1 (1259510), страница 2
Текст из файла (страница 2)
NN(AMaa := М прДВ_ном Aϕ iiAϕ := s ( i)i)cMaa := cspline Aϕ , AMaa()М прДВ_ном ( ϕ) := interp cMaa , Aϕ , AMaa , ϕs1 := 4πi := 0 .. NNAMbb := М прДВ_хх A1ϕ iiA1ϕ := s1 ( i)i(cMbb := cspline A1ϕ , AMbb()М прДВ_хх ( ϕ) := interp cMbb , A1ϕ , AMbb , ϕ)Решение задачи динамики при номинально нагрузкеϕ := 0 , 1deg .. 720degРабота внешних сил:AДВ_ном( ϕ) := AC_ном( ϕ) := M прС_ном ϕN := 200s := 4ϕМ прДВ_ном ( ϕ) dϕ0πi := 0 ..
NNi(AMzbb := AДВ_ном Az1 ϕ iiAz1 ϕ := s ( i ))cMzbb := cspline Az1ϕ , AMzbb()AДВ_ном( ϕ) := interp cMzbb , Az1 ϕ , AMzbb , ϕ3AC_ном( f ) = -1.47471 10AΣ_ном ( ϕ) := AДВ_ном ( ϕ) + AC_ном( ϕ)AΣ_ном ( 4π) = -0.025883AΣ_ном ( f ) = 1.84601 1042104AΣ_ном( ϕ)110- AC_ном( ϕ)AДВ_ном( ϕ)0AC_ном( ϕ)4- 1104- 2100200400ϕdeg600800Регулирование движения по методу Мерцалова:ωср_хх :=n 1хх 2π60= 73.30383TII_хх( ϕ) := Jпр2( ϕ) ωср_хх2TII_хх( f ) = 408.650972500400TII_хх( ϕ) 3002001000200400600ϕdegTнач_хх := Jпр2( f ) ωср_хх22= 408.65097ΔTI_хх( ϕ) := AΣ_хх ( ϕ) - TII_хх( ϕ)ΔTI_хх( f ) = 898.6766931.5 1031 10ΔTI_хх( ϕ)500AΣ_хх( ϕ)0- 5000200400600800ϕdegРассчитаем наибольшее изменение кинетическ ой энергии первой группы звеньев:t1 := 180deg = 3.14159Given(t1min)t1min := Minimize ΔTI_хх , t1(deg= 3.61383)ΔTI_хх t1min = -135.86778t2 := 330 deg = 5.75959Given(t1max)t1max := Maximize ΔTI_хх , t2(deg= 87.33946)ΔTI_хх t1max = 899.65476()()ΔTInb_хх := ΔTI_хх t1max - ΔTI_хх t1minОпределим момент инерции первой группы звеньев:δ :=130Jпр1_хх :=ΔTInb_хх2ωср_хх δJпрΣ( ϕ) := Jпр2( ϕ) + Jпр1_хх3ΔTInb_хх = 1.03552 10Jпр1_хх = 5.78133Колебания обобщенной координаты:(ΔTI_хх(t1max) + ΔTI_хх(t1min))ΔTI_хх( ϕ) 2Δωхх( ϕ) :=ωср_хх Jпр1_ххωхх( ϕ) := ωср_хх + Δωхх( ϕ)Δωхх( f ) = 1.21942ωхх( 0 ) = 72.106557574ωхх( ϕ)ωср_хх73720200400600800ϕdegM Σпр_хх( ϕ) - dJ ( ϕ) ε хх( ϕ) :=ωхх( ϕ)22ε хх( f ) = -6.85887JпрΣ( ϕ)300200100εхх( ϕ)0- 100- 200- 3000200400ϕdeg600800Габаритные размеры и масса маховика:J10 := 0.0036 9.81 = 0.03532J0пр = 7.2594JM := Jпр1_хх - J10 = 5.74601Используем маховик в виде диск а:5ДиаметрD := 0.366 JM = 0.51922Ширинаb := 0.2 D = 0.10384Массаm := 1230 D = 172.174263m :=8 JM2= 170.5DРазмеры и масса маховика в виде обода со спицами и втулкой:5DMXc := 0.437 JM = 0.61995b MX := 0.2D = 0.10384DMX1 := 0.8 DMXc = 0.4959622mMXc := 6123 DMXc - DMX1 b MX = 88G1y := -mMXc g = -863.03731Силовой расчет:Определение ускорений:2aAx( ϕ) := aqAx( ϕ) ω1 ( ϕ) + ε 1 ( ϕ) VqAx( ϕ)2aAy( ϕ) := aqAy( ϕ) ω1 ( ϕ) + ε 1 ( ϕ) VqAy( ϕ)aA( ϕ) :=2aAx( ϕ) + aAy( ϕ)222aA1y( ϕ) := aqA1y( ϕ) ω1 ( ϕ) + ε 1 ( ϕ) VqA1y( ϕ)2aA1x( ϕ) + aA1y( ϕ)aAy( f ) = -0.88322aA( f ) = 2.10681 10aA1x( ϕ) := aqA1x( ϕ) ω1( ϕ) + ε1 ( ϕ) VqA1x( ϕ)aA1( ϕ) :=3aAx( f ) = -2.10681 10233aA1x( f ) = 2.10681 10aA1y( f ) = 0.883223aA1( f ) = 2.10681 102aS2_14x( ϕ) := aqS2_14x( ϕ) ω1 ( ϕ) + VqS2_14x( ϕ) ε 1 ( ϕ)3aS2_14x( f ) = -1.5169 102aS2_14y( ϕ) := aqS2_14y( ϕ) ω1 ( ϕ) + VqS2_14y( ϕ) ε 1 ( ϕ)aS2_14( ϕ) :=2aS2_14x( ϕ) + aS2_14y( ϕ)aS2_14y( f ) = 162.805732aS2_14( f ) = 1.52562 102aS2_23x( ϕ) := aqS2_23x( ϕ) ω1 ( ϕ) + VqS2_23x( ϕ) ε 1 ( ϕ)aS2_23x( f ) = 1.5169 10332aS2_23y( ϕ) := aqS2_23y( ϕ) ω1 ( ϕ) + VqS2_23y( ϕ) ε 1 ( ϕ)aS2_23( ϕ) :=2aS2_23x( ϕ) + aS2_23y( ϕ)22aBy( ϕ) := aqBy( ϕ) ω1 ( ϕ) + VqBy( ϕ) ε1 ( ϕ)aB( ϕ) :=aBy( ϕ)aS2_23( f ) = 1.52581 10aBy( f ) = 583.720152aB( f ) = 583.720152aB1y( ϕ) := aqB1y( ϕ) ω1 ( ϕ) + VqB1y( ϕ) ε1 ( ϕ)aB1( ϕ) :=aS2_23y( f ) = 164.57216aB1y( ϕ)2aB1y( f ) = 585.48658aB1( f ) = 585.486582ε 2_14( ϕ) := ε q2_14( ϕ) ω1 ( ϕ) + ωq2_14( ϕ) ε 1 ( ϕ)2ε 2_23( ϕ) := ε q2_23( ϕ) ω1 ( ϕ) + ωq2_23( ϕ) ε 1 ( ϕ)3ε 2_14( f ) = -7.79471 103ε 2_23( f ) = 7.79471 103Звено 1:Jпр1_хх = 5.78133Φ 1x := 0Φ 1y := 0G1y = -863.03731()M Φ1 ( ϕ) := - -Jпр1_хх - J0пр ε 1 ( ϕ)M Φ1 ( f ) = 154Звено 2:Φ 2_14x( φ) := -m2 aS2_14x( φ)Φ 2_14x( f ) = 5400Φ 2_14y( φ) := -m2 aS2_14y( φ)Φ 2_14y( f ) = -580Φ 2_23x( φ) := -m2 aS2_23x( φ)Φ 2_23x( f ) = -5400Φ 2_23y( φ) := -m2 aS2_23y( φ)Φ 2_23y( f ) = -586G2_14y := G2y = -34.9236G2_23y := G2y = -34.9236M Φ2_14( φ) := -J2S ε 2_14( φ)M Φ2_14( f ) = 141.5M Φ2_23( φ) := -J2S ε 2_23( φ)M Φ2_23( f ) = -141.5Звено 3Φ 3_14x( φ) := 0Φ 3_14y( φ) := -m3 aBy( φ)Φ 3_14y( f ) = -1868G3_14y := -m3 gG3_14y = -31.4M Φ3_14 := 0F1y_ном( f ) = -1.81905 10F1y_ном( ϕ)43F4y_ном( f ) = 1.64065 10F4y_ном( ϕ)M Φ3_23 := 0Φ 3_23x( φ) := 0Φ 3_23y( φ) := -m3 aB1y( φ)G3_23y := -m3 gF2y_ном( ϕ)Φ 3_23y( f ) = -1874F3y_ном( ϕ)G3_23y = -31.43F2y_ном( f ) = -319.99263F3y_ном( f ) = -1.64065 100Φ(ϕ)+G+F(ϕ)( 3_14y3_14y1y_ном )Φ2_14x( ϕ)T1 ( ϕ) := Φ(ϕ)+G( 2_14y2_14y)0-Φ2_14x( ϕ) yS2_14( ϕ) + ( Φ 2_14y( ϕ) + G2_14y) xS2_14( ϕ) + MΦ2_14( ϕ) 04 2.00898 10 3T1 ( f ) = -5.40018 10 614.51198 0 300.40944 1 0 0K1( ϕ) := 0 -yB( ϕ) 0 1 00K1( f ) = 0 -0.27029 0010000100-10100-1 01 -0.27029 000 0.270290000001010000100-1010 0-101 -yB( ϕ) xB00 yB( ϕ) -xB -yA( ϕ) xA( ϕ) 00 0 1 00.075 DR1( ϕ) := lsolve( K1( ϕ) , T1 ( ϕ) )T(3DR1( f ) = 4.63363 10R130x( ϕ) := DR1( ϕ)0R132x( ϕ) := DR1( ϕ)2R121x( ϕ) := DR1( ϕ)4R112x( ϕ) := -R121x( ϕ)30 -4.63363 1042.00898 104-1.00338 10M130( ϕ) := DR1( ϕ)1R132y( ϕ) := DR1( ϕ)3R121y( ϕ) := DR1( ϕ)5R112y( ϕ) := -R121y( ϕ)'2.07043 104)0Φ(ϕ)+G+F(ϕ)( 3_23y3_23y2y_ном )Φ2_23x( ϕ)T2 ( ϕ) := (Φ2_23y( ϕ) + G2_23y)0-Φ2_23x( ϕ) yS2_23( ϕ) + ( Φ 2_23y( ϕ) + G2_23y) xS2_23( ϕ) + MΦ2_23( ϕ) 03 2.22494 10 3T2 ( f ) = 5.40018 10 620.80049 0 -300.74902 1 00K2( f ) = 0 -0.27029 0 1 0 0K2( ϕ) := 0 -yB1( ϕ) 0010000100-10100-1 01 -0.27029 000 0.270290000001010000100-1010 0-101 -yB1( ϕ) xB100 yB1( ϕ) -xB1 -yA1( ϕ) xA1( ϕ) 00 0 1 0-0.075 DR2( ϕ) := lsolve( K2( ϕ) , T2 ( ϕ) )T(3DR2( f ) = 323.05749 0 -323.05749 2.22494 10R230x( ϕ) := DR2( ϕ)0R232x( ϕ) := DR2( ϕ)2R221x( ϕ) := DR2( ϕ)4R212x( ϕ) := -R221x( ϕ)35.07712 10R232y( ϕ) := DR2( ϕ)3R221y( ϕ) := DR2( ϕ)5R212y( ϕ) := -R221y( ϕ)32.84574 10)M230( ϕ) := DR2( ϕ)10-( Φ 3_23y( ϕ) + G3_23y + F3y_ном( ϕ) )-Φ2_23x( ϕ)T3 ( ϕ) :=-( Φ2_23y( ϕ) + G2_23y)0--Φ2_23x( ϕ) yS2_23( ϕ) + ( Φ 2_23y( ϕ) + G2_23y) xS2_23( ϕ) + MΦ2_23( ϕ) 1 0 0K3( ϕ) := 0 -yB1( ϕ) 0 1 00K3( f ) = 0 -0.27029 000001010000100-1010 0-101 -yB1( ϕ) xB100 yB1( ϕ) -xB 1 -yA1( ϕ) xA1( ϕ) 010000100-10100-1 01 -0.27029 000 0.270290000 1 0-0.075 0DR3( ϕ) := lsolve( K3( ϕ) , T3 ( ϕ) )T(3DR3( f ) = -43.40198 0 43.40198 3.5456 10R330x( ϕ) := DR3( ϕ)0R332x( ϕ) := DR3( ϕ)2R321x( ϕ) := DR3( ϕ)4R312x( ϕ) := -R321x( ϕ)5.44358 103M330( ϕ) := DR3( ϕ)1R332y( ϕ) := DR3( ϕ)3R321y( ϕ) := DR3( ϕ)5R312y( ϕ) := -R321y( ϕ)34.1664 10)0-( Φ 3_14y( ϕ) + G3_14y + F4y_ном( ϕ) )-Φ2_14x( ϕ)T4 ( ϕ) :=-( Φ2_14y( ϕ) + G2_14y)0--Φ2_14x( ϕ) yS2_14( ϕ) + ( Φ 2_14y( ϕ) + G2_14y) xS2_14( ϕ) + MΦ2_14( ϕ) 1 0 0K4( ϕ) := 0 -yB( ϕ) 00258.64574 -5.40018 103 T4 ( f ) = 614.51198 0 300.40944 1 00K4( f ) = 0 -0.27029 0010000100-10100-1 01 -0.27029 000 0.270290000001010000100-1010 0-101 -yB( ϕ) xB00 yB( ϕ) -xB -yA( ϕ) xA( ϕ) 00 0 1 00.075 DR4( ϕ) := lsolve( K4( ϕ) , T4 ( ϕ) )T(3DR4( f ) = -869.15873 0 869.15873 258.64574 -4.53102 10R430x( ϕ) := DR4( ϕ)0M430( ϕ) := DR4( ϕ)R432x( ϕ) := DR4( ϕ)2R432y( ϕ) := DR4( ϕ)3R421x( ϕ) := DR4( ϕ)4R412x( ϕ) := -R421x( ϕ)1R421y( ϕ) := DR4( ϕ)5R412y( ϕ) := -R421y( ϕ)873.15772)R112x( ϕ) + R212x( ϕ) + R312x( ϕ) + R412x( ϕ)R112y( ϕ) + R212y( ϕ) + R312y( ϕ) + R412y( ϕ) + G1yT5 ( ϕ) := -( -yA( ϕ) R112x( ϕ) + xA( ϕ) R112y( ϕ) ) + -yA( ϕ) R412x( ϕ) + xA( ϕ) R412y( ϕ) ... + -yA1( ϕ) R212x( ϕ) + xA1( ϕ) R212y( ϕ) ... + -yA1( ϕ) R312x( ϕ) + xA1( ϕ) R312y( ϕ) + MΦ1 ( ϕ) -4.04413 103 T5 ( f ) = 4 2.94526 10 938.82792 1 0 0 K5( ϕ) := 0 1 0 0 0 1 DR5( ϕ) := lsolve( K5( ϕ) , T5 ( ϕ) )T(3DR5( f ) = -4.04413 10R 10x( ϕ) := DR5( ϕ)42.94526 10938.82792R 10y( ϕ) := DR5( ϕ)04R 10y( f ) = 2.94526 10M 1 ( f ) = 938.82792M 1 ( ϕ) := DR5( ϕ)13R 10x( f ) = -4.04413 102M 1 ( 10deg) = 938.82792M прС_ном + M 1( f )M прС_ном= 7.6637 10)M 1 ( 210deg) = 938.82792M прС_ном = -938.82792- 12938.8279207938.82792065M 1938.8279206( ϕ)938.82792055938.82792050510ϕ15R130( ϕ) :=R130x( ϕ)R230( ϕ) :=2R230x( ϕ)R330( ϕ) :=32R330x( ϕ)2R130( f ) = 4.634 10R130x( f ) = 4.63363 10R230( f ) = 323.057R230x( f ) = 323.057493R330( f ) = 43.402R330x( f ) = -43.40198R430( ϕ) :=fi := 0 , 4fideg=π24R430x( ϕ)2R430( f ) = 869.159..
4 πR130x( fi)200=R230x( fi) =2.62031-0.0131-197.6460538.04969-291.328576031.64193323.057499023.16817753.1143412015.63891-1.29648·1031506.13656-2.620311800.01317.60991·103210-3.383736.32836·103240-4.06554.63363·103270-0.216973.12779·1033006.45981.22732·1033307.333152.62031360-0.0131-676.75282-9.52644-813.10693-10.36843-43.38721-4.345871.29198·1030.146261.46665·1035100.80287-2.620315400.0131-1.9053·103570-0.98831-2.07371·103600-1.45673-869.184396301.6154229.259516603.76601160.57408690-6.482132.62031720-0.0131030390420450480R430x( f ) = -869.15873R330x( fi) =2.62031R430x( fi) =-2.62031-676.73956-1.90523·103-813.09364-2.07363·103-43.40198-869.158731.29192·10329.238711.4666·103160.56215-2.620312.62031-1.90525·103-197.65213-2.07366·103-291.33517-869.16591323.0646729.2471753.1493160.56894-1.29645·1032.62031-2.62031-197.659317.60988·103-291.341866.32833·103323.072264.63362·103753.171843.12779·103-1.29642·1031.22732·103-2.620312.620317.60986·103-676.755446.32831·103-813.111424.63361·103-43.376323.1278·1031.29201·1031.22732·1031.46665·1032.62031-2.620312R132( ϕ) :=R132x( ϕ) + R132y( ϕ)(2φ132( ϕ) := angle R132x( ϕ) , R132y( ϕ)3R132( f ) = 20.617 10)fi := 0 , 4π24φ132( f )..