Диссертация (1173093), страница 20
Текст из файла (страница 20)
– NewYork : Springer, 2011. – 237 p.127. Blundell, M. The Multibody Systems Approach to Vehicle Dynamics[Text] / M. Blundell, D. Harty. – 2-d ed. – Boston, MA : Elsevier, 2014. – 768 p.128. Datta, J. Synthesis and Investigation of Glycolysates and Obtained Polyurethane Elastomers [Text] / J. Datta // J. of Elastomers and Plastics. – 2010 – V. 42. – P.117-127.129. Dixon, J.C. The Shock Absorber Handbook [Text] / J.C. Dixon.
– 2-d ed. –Chichester : John Wiley, 2007. – 432 p.130. Maity, M. Polyblend Systems of Polyurethane Rubber and Silicone Rubberin the Presence of Silane Grafting Agent [Text] / M. Maity // J. of Elastomers and Plastics. – 2001. – V. 33. – P. 211-224.131. Manfanovskii, S.B. Internally Sprung Wheels [Text] / S.B. Manfanovskii,A.A. Enaev // Russian Engineering Research – 2016. – Vol. 36, No. 1 – pp. 6-9.132.
Pat. JP200621637 Japan, IPC B60B9/10, B60B9/26. Caster / IshikawaTomohito, Maki Kenji, Imamura Tadashi; 2004-2011319; заявл. 08.07.04 ; опубл.26.01.06.133. Pat. US2011/0030861 United States, IPC B60B9/06. Internal wheel suspension and shock absorbing system / David L. Wichern; 12/924.958; заявл. 08.10.10 ;опубл. 10.02.11.149134. Pat. US2014/0034198 United States, IPC B60B9/02. Internal wheel suspension system with shock absorption / David L. Wichern; 13/136.930; заявл. 15.08.11; опубл. 06.02.14.135.
Pat. US6.698.480 United States, IPC B60B9/08, B60B9/18. Nonpneumatic tire and wheel system / Maurice H. Cornellier; 10/290.840; заявл. 08.11.02 ;опубл. 02.03.04.136. Polynomial Hyperelastic Model [Electronic resource] // Solved withCOMSOLMultiphysics5.0.–P.1-14.–URL:https://www.comsol.com/model/download/197707/models.nsm.polynomial_hyperelastic.pdf.137.
Rek, V. Thermal Degradation of Polyurethane Elastomers: Determinationof Kinetic Parameters [Text] / V. Rek // J. of Elastomers and Plastics. – 2003. – V. 35. –P. 311-323.138. Thermoplastic Polyurethane Elastomers (TPU). Elastollan – Material Properties [Electronic resource] // Technical Information. BASF The Chemical Company. –2011.November.–44p.–URL:http://www.polyurethanes.basf.de/pu/solutions/us/function/conversions:/publish/content/group/Arbeitsgebiete_und_Produkte/Thermoplastische_Spezialelastomere/Infomaterial/elastollan_material_uk.pdf.139.
Vorobyov, Y. Numerical simulation of laminated plastics pulse rivetingprocess [Electronic resource] / Y. Vorobyov, I. Pechenizkiy, V. Garin, Y. Tsegelnyk //Авиационно-космическая техника и технология. – 2007. №3 (39). – P. 47-51. –URL:http://www.khai.edu/csp/nauchportal/Arhiv/AKTT/2007/AKTT307/Vorobyev.pdf.150ПРИЛОЖЕНИЕПриложение АПрограмма расчета колебаний масс транспортного средства, оснащенного:КВП взамен подвески; стандартной подвеской и традиционными колесамиПриведение системы уравнений к форме Коши.Ввод дополнительных переменных: dξ и dψ.dξ:= dξψ:= dψdtddtПреобразование исходных уравнений.ddt( dξ ) := −1⋅ 2 ⋅ η B ⋅ ( dξ − dψ ) + 2 ⋅ Cb ⋅ ( ξ − ψ )m 1d( dψ ) := −⋅ η ⋅ ( dψ − dq ( t ) ) + Cs ⋅ ( ψ − q ( t ) ) − η B ⋅ ( dξ − dψ ) − Cb ⋅ ( ξ − ψ )m0 shdtq ( t ) – входное воздействие;dq ( t ) – производная входного воздействия.Ввод «машинных переменных».x0 := ξx1 := ψx2 := dξx3 := dψЗапись системы уравнений в «машинных переменных».ddtx0 := x2151ddtddtddtx1 := x3x2 := −x3 := −⋅ 2 ⋅ η B ⋅ ( x2 − x3) + 2 ⋅ Cb ⋅ ( x0 − x1)m 1⋅ η ⋅ ( x3 − dq ( t ) ) + Cs ⋅ ( x1 − q ( t) ) − η B ⋅ ( x2 − x3) − Cb ⋅ ( x0 − x1)m0 sh1Ввод параметров колебательной системы.m := 550Подрессоренная масса, кгm0 := 10Масса обода колеса, кгŋsh := 1260Коэффициент демпфирования шины, Н·с/мŋВ := 1035Коэффициент демпфирования упругого элемента, Н·с/мCb := 330000Коэффициент нормальной жесткости упругого элемента, Н/мCS := 94000Коэффициент нормальной жесткости шины, Н/мQ := 0,03Величина входного ступенчатого воздействия, мts := 0,01Время наступления перепада ступенчатого воздействия, cτ := 0,1Длительность импульса, аппроксимирующего производнуюступенчатого воздействия, сq0 := 0,03Половина высоты неровности, мω := 35·2·πЧастота возмущающего воздействия, рад / cvar := 0Вариант возмущающего воздействия:0 - ступенчатый сигнал; 1 - смещенный гармоническийn1 := 10000Число точек расчета для временного интервала 1сt0 := 0Начальное время, сtkon := 5Конечное время, сПолное число точек расчета.n := n1 ⋅ tkonФункция, задающая гармоническое возмущающее воздействие.()q_g t , q0 , ω := q0 ⋅ ( 1 − cos ( ω ⋅ t ) )Функция, задающая производную гармонического возмущающего воздействия.()dq_g t , q0 , ω := q0 ⋅ ω ⋅ sin( ω ⋅ t )Функция, задающая ступенчатое возмущающее воздействие.q_s ( t , Q , ts) :=q ← 0 if t < tsq ← Q otherwiseq152Функция, задающая производную ступенчатого возмущающего воздействия, как импульс конечной длины и амплитуды.dq_s ( t , ts , τ , dQ) :=dq ← 0 if t < tsotherwisedq ← 0 if t ≥ ts + τdq ← dQ otherwisedqАмплитуда импульса.QdQ :=τПроизводная возмущающего ступенчатого воздействия амплитудой Q моделируется какимпульс конечной длительностью τ.
Таким образом, выполняется условие dQ· τ= Q, или дляслучая единичной ступенчатой функции dQ· τ= 1.В зависимости от заданного вида возмущающего воздействия выбирается первая иливторая строка функций q(t) и dq(t).q ( t ) :=q_s ( t , Q , ts) if var < 0.5(q_g t , q0 , ωdq ( t ) :=)otherwisedq_s ( t , ts , τ , dQ ) if var < 0.5(dq_g t , q0 , ω)otherwisex2x31D( t , x) := − ⋅ 2 ⋅ η B ⋅ x2 − 2 ⋅ η B ⋅ x3 + 2 ⋅ Cb ⋅ ( x0 − x1)m − 1 ⋅ η ⋅ x − dq ( t) + C ⋅ x − q ( t) − η ⋅ x − x − C ⋅ x − x ) s (1) B ( 2 3) b ( 0 1) m0 sh ( 3Начальные условия.00y0 := 00Решение системы дифференциальных уравнений.
Извлечение переменных из решения.y := AdamsBDF( y0 , t0 , tkon , n , D)〈0〉t := yξ〈1〉:= y153ψ〈2〉:= y〈3〉dξ := y〈4〉dψ := yВычисление производных. dξ1 i ξ i dψ1 iψi := Dt , ddξ i i dξ i ddψ i dψ i 154Приложение БПрограмма расчета колебаний масс транспортного средства,оснащенного традиционными колесами без подвескиПриведение уравнений к форме Коши.Ввод дополнительной переменной: dξ .dξdt:= dξПреобразование исходного уравнения:M⋅ddtddt( dξ ) + 2 ⋅ η sh ⋅ dξ −( dξ ) := −M 1ddtq ( t ) + 2 ⋅ Csh ⋅ ( ξ − q ) := 0⋅ 2 ⋅ η sh ⋅ dξ −ddtq ( t ) + 2 ⋅ Csh ⋅ ( ξ − q )q ( t ) - входное воздействие;dq ( t ) - производная входного воздействия.Ввод «машинных переменных».x0 := ξx1 := dξЗапись системы уравнений в «машинных переменных»:ddtddtx0 := x1( x1):= −M 1⋅ 2 ⋅ η sh ⋅ x1 −ddtq ( t ) + 2 ⋅ Csh ⋅ ( x0 − q )Ввод параметров колебательной системы.М := 650Подрессоренная масса, кгŋsh := 1260Коэффициент демпфирования шины, Н·с/мCS := 130000Коэффициент нормальной жесткости шины, Н/м155Вариант возмущающего воздействия:var := 00 - ступенчатый сигнал; 1 - смещенный гармонический;2 - косоугольный сигнал; 3 - «наезд на неровность».n1 := 10000Число точек расчета для временного интервала 1сt0 := 0Начальное время, сtkon := 5Конечное время, сПараметры, задаваемые для ступенчатого воздействия:Q := 0,025Величина входного ступенчатого воздействия, мts := 0,01Время наступления перепада ступенчатого воздействия, cτ := 0,1Длительность импульса, аппроксимирующего производнуюступенчатого воздействия, сПараметры, задаваемые для смещенного гармонического воздействия:q0 := 0,03Половина высоты неровности, мω := 35·2·πЧастота возмущающего воздействия, рад / cПараметры, задаваемые для косоугольного воздействия:Q1 = 0,5Величина входного ступенчатого воздействия, мvQ = 10Скорость изменения входного воздействия при косоугольном входномсигнале, м/cВремя наступления перепада ступенчатого воздействия, cts1 = 0,1Параметры, задаваемые для воздействия "наезд на неровность":l = 0,01Длина неровности, мh = 0,05Высота неровности, мV= 5∙Wt1 = 0,11000Z3600Скорость движения, м /сВремя наезда на неровность, сПолное число точек расчета.n := n1 ⋅ tkonФункция, задающая гармоническое возмущающее воздействие()q_g t , q0 , ω := q0 ⋅ ( 1 − cos ( ω ⋅ t ) )Функция, задающая производную гармонического возмущающего воздействия()dq_g t , q0 , ω := q0 ⋅ ω ⋅ sin( ω ⋅ t )156Функция, задающая ступенчатое возмущающее воздействиеq_s ( t , Q , ts) :=q ← 0 if t < tsq ← Q otherwiseqФункция, задающая производную ступенчатого возмущающего воздействия, как импульс конечной длины и амплитуды.dq_s ( t , ts , τ , dQ) :=dq ← 0 if t < tsotherwisedq ← 0 if t ≥ ts + τdq ← dQ otherwisedqАмплитуда импульсаdQ :=QτПроизводная возмущающего ступенчатого воздействия амплитудой Q моделируется какимпульс конечной длительностью τ Амплитуда этого импульса dQ= Q/τ.
Таким образом, выполняется условие dQ·τ= Q, или для случая единичной ступенчатой функции dQ·τ= 1.Функция, задающая косоугольное возмущающее воздействие()q_v t , Q1 , vQ , ts :=q ← 0 if t < tsq ← vQ ⋅ ( t − ts) otherwiseq ← Q1 if q ≥ Q1qФункция, задающая производную косоугольного возмущающего воздействия()dq_v t , Q1 , vQ , ts :=if t < tsq ←0dq ← 0otherwiseq ← vQ ⋅ ( t − ts)dq ← vQif q ≥ Q1q ← Q1dq ← 0dq157Функция, задающая возмущающее воздействие "наезд на препятствие"q_n ( t , t1 , l , h , V) :=τ1←lVq ←0 π ⋅ V ⋅ ( t − t1) if ( t > t1) ∧ [ t < ( t1 + τ1 ) ] l q ← h ⋅ sinqФункция, задающая производную возмущающего воздействия "наезд на препятствие"dq_n ( t , t1 , l , h , V) :=τ1←lVq ←0 π ⋅ V ⋅ h ⋅ cos π ⋅ V ⋅ ( t − t1) if ( t > t1) ∧ [ t < ( t1 + τ1 ) ] l l q ←qВ зависимости от заданного вида возмущающего воздействия выбирается первая, вторая,третья или четвертая строка функций q(t) и dq(t)q ( t ) :=q_s ( t , Q , ts) if var < 0.5() if (var ≥ 0.5) ∧ (var < 1.5)q_v ( t , Q , vQ , ts) if ( var ≥ 1.5) ∧ ( var < 2.5)q_g t , q0 , ωq_n ( t , t1 , l , h , V) otherwisedq ( t ) :=dq_s ( t , ts , τ , dQ) if var < 0.5() if (var ≥ 0.5) ∧ (var < 1.5)dq_v ( t , Q , vQ , ts) if ( var ≥ 1.5) ∧ ( var < 2.5)dq_g t , q0 , ωdq_n ( t , t1 , l , h , V) otherwiseddtddtx0 := x1( x1):= −M 1⋅ 2 ⋅ η sh ⋅ x1 −ddtq ( t ) + 2 ⋅ Csh ⋅ ( x0 − q )x1D( t , x) := − 1 ⋅ 2 ⋅ η sh ⋅ ( x1 − dq ( t) ) + 2 ⋅ Cs ⋅ ( x0 − q ( t) ) M158Начальные условия00y0 := y := AdamsBDF( y0 , t0 , tkon , n , D)y := rkfixed( y0 , t0 , tkon , n , D)〈〉t := y 0ξ〈〉:= y 1〈〉dξ := y 2i := 0 ..