Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

MATHEMATICAL MODEL OF ELASTIC-PLASTIC DEFORMATIONS OF VEHICLES FOR CARRYING OUT AUTOTECHNICAL EXPERTISE OF ROAD TRANSPORT ACCIDENTS

Bessarabov E.N. 1 Bessarabov A.N. 1 Gasanov B.G. 1 Efimov A.D. 1 Nefedov V.V. 1
1 Platov South-Russian State Polytechnic University (NPI)
The absolute majority of cars sold on the territory of Russia are the anti-lock braking system (ABS). If you do not have a guarantee order, then you can contact us by email. Therefore, it is necessary to improve or develop new methods of autotechnical examination (ATE), using the laws of mechanics, the theory of operational properties of vehicles (ATS), continuum mechanics and the theory of plasticity. The article presents the principles and possible directions of application of the relevant procedures when conducting an auto-technical examination of an accident. The resulting mathematical expressions can be used in determining the speeds and directions of movement of the means at the moment of collision. Proposals for the principles of creating a mathematical model of an accident, in accordance with the conduct of an autotechnical examination of an accident, in particular, when arriving at a fixed obstacle (standing vehicle) and collisions. It has been substantiated that in order to carry out ATE in collisions and on arrivals, it is necessary to apply a mathematical model based on the laws of conservation of energy and momentum; deformation theory of plasticity and dynamics of deletions to determine the stopping distance, taking into account road conditions; describing analytical dependencies for estimating the stress-strain state based on finite element methods.
autotechnical examination
road traffic accident
vehicle
finite element method
law of conservation of momentum and energy

При проведении АТЭ ДТП следственные органы перед экспертами чаще всего ставят следующие вопросы [1]:

1. С какими скоростями двигались ТС в моменты возникновения опасной ситуации на дороге и столкновения (наезда)?

2. Имели ли водители техническую возможность предотвратить ДТП путем экстренного торможения или маневра?

3. Соответствовали (или не соответствовали) действия участников ДТП требованиям правил дорожного движения, и какие несоответствия послужили причиной ДТП?

Как известно [2], для ответа на указанные вопросы при АТЭ эксперты в основном определяют остановочный (So) и тормозной путь (SТ) по следующим формулам:

bess01.wmf;

bess02.wmf, (1)

где V0 – начальная скорость автомобиля, км/ч;

tc, tн – время запаздывания тормозной системы и время нарастания замедления соответственно;

tуст – установившееся замедление;

A – коэффициент, характеризующий время срабатывания тормозной системы [3].

В формулах (1) начальная скорость неизвестна. В случае экстренного торможения из выражений (1) получают [4]:

bess03.wmf; (2)

bess04.wmf, (3)

где j – замедление ТС за время нарастания тормозных сил до максимального значения (tЗ);

Sю – длина пути юза;

Vпс – скорость ТС в момент первого столкновения (наезда).

Анализ рынка сбыта в России за последние три года показывает, что абсолютное большинство легковых автомобилей, продаваемых на территории России, оборудованы АБС [5]. Поскольку ТС с АБС при экстренном торможении практически не оставляют следы юза, то формулы (2) и (3) не могут быть использованы для ответа на поставленные вопросы. Поэтому нужно совершенствовать или разработать новые методы проведения АТЭ, используя законы механики, теорию эксплуатационных свойств АТС, механики сплошной среды и теории пластичности.

Целью настоящей работы является обобщение математических моделей для экспертизы ДТП, позволяющих ответить на главные вопросы, возникающие при АТЭ.

Материалы и методы исследования

Изменение кинетической энергии (работа) ТС при ДТП по времени равно мощности результирующей силы (F) [6]. Например, при наезде на недеформируемое препятствие

bess05.wmf

где hдеф – глубина деформации ТС;

tз – интенсивность замедления ТС за период удара.

Все силы, действующие на ТС при ДТП, можно разделить на внешние и внутренние. Пренебрегая работой сил инерции контактирующих деталей ТС, допускаем, что работа внешних сил равна работе сил взаимодействия. К внутренним силам относятся потенциальные и диссипативные. Следовательно, при столкновениях ТС

bess06.wmf

где

bess07a.wmf

bess07aa.wmf;

bess08.wmf – скорость ТС после столкновения (наезда).

Приращение работы потенциальных сил равно сумме приращения работы упругой и пластической деформации:

bess09.wmf

К работе диссипативных сил можно отнести работу трения скольжения, тогда

bess10.wmf (4)

где Fтр, μК – сила и коэффициент контактного трения;

dSК – путь скольжения каждого элемента в зоне взаимодействия ТС;

PК – нормальное давление.

Все ДТП по механике деформации условно можно разделить на три типа: абсолютно упругие, абсолютно неупругие (идеально пластическое тело) и упруго-пластические. Это означает, что в первом случае пренебрегаем пластической (остаточной) деформацией ТС при ДТП. Во втором случае вся кинетическая энергия автомобиля при ДТП затрачивается на пластическую деформацию. Наиболее достоверен третий вариант, когда при ДТП работа деформации включает и пластическую, и упругую составляющие. Для аналитического описания лучше использовать законы сохранения импульса и энергии.

Неупругое столкновение (упругопластическое деформирование объектов взаимодействия при ДТП), в результате которого внутренняя энергия bess11.wmf после столкновения одного или обоих ТС изменяется, и изменяется направление их движения (хотя бы одного).

Запишем импульс каждого ТС после столкновения [2]

bess12.wmf. (5)

Для определения возможных направлений движения ТС после столкновения построим векторную диаграмму импульсов (рисунок) в системе отсчета Эйлера.

bessaar1.tif

Векторная диаграмма импульсов в случае неупругого столкновения

Импульс ТС с массой m1 P1 = AB делят точкой О на две части, пропорциональные массам ТС bess13.wmf. Из точки O проводим окружность радиусом P′1, определяемым по формуле (5). Эта окружность является геометрическим местом точек возможных положений вершины треугольника импульсов АВ, стороны АС и СВ которого равны импульсам соответствующих ТС после столкновения.

Покажем вариант оценки предельной скорости или кинетической энергии, при котором кинетическая энергия пойдет только на деформацию, т.е. увеличение внутренней энергии, и после столкновения объекты взаимодействия остаются на месте ДТП или двигаются как одно целое. В системе отсчета Лагранжа это условие можно записать в следующем виде:

bess14.wmf

Рассмотрим случай наезда на стоящее ТС, когда кинетическая энергия одного из ТС близка к нулю. Если после столкновения оба ТС не двигались, то их суммарный импульс равен импульсу P1 перед ДТП первого ТС, а кинетическая энергия системы

bess15.wmf; bess16.wmf

Если после ДТП оба ТС переместились на определенное расстояние вместе, то минимальная кинетическая энергия одного из них, достаточная для их движения, определяется по следующей формуле:

bess17.wmf;

или

bess18.wmf (6)

Поскольку

bess19.wmf (7)

то из формул (6) и (7) находим предельную скорость одного из ТС до столкновения, выше которой они могут двигаться как одно целое, или, наоборот, ниже которой они остаются на месте ДТП. Для этого определим bess20.wmf по следующей формуле:

bess21.wmf (8)

В формуле (8) bess22.wmf – это приращение внутренней энергии системы, т.е. энергии, затраченной на работу деформации ТС.

Рассмотрим методы определения работы деформации элементов кузова ТС при ДТП.

Приращение работы деформации конструкции ТС при ДТП можно определить, если известно их напряженно-деформированное состояние [7]:

bess23.wmf

где σij – компоненты тензора напряжения;

bess24.wmf, bess25.wmf – приращение пластической и упругой деформаций каждого отдельного элемента ТС в процессе ДТП.

Связь между напряжениями и деформацией для упругого тела установлена на основе законов Гука и Пуассона:

bess26.wmf;

bess27.wmf (9)

где σср, εср – компоненты шаровых тензоров
(bess28.wmf; bess29.wmf);

δij – символ Кронекера.

В частности, из выражения (9) следует, что

bess30.wmf; bess31.wmf,

где Е – модуль Юнга;

v – коэффициент Пуассона;

G – модуль сдвига (G = E/2(1 + v).

Для малых упругих деформаций [6]

bess32.wmf

где dUi и dUj – приращение вектора перемещения выделенного элемента (например, прямоугольника на поверхности деформированной детали кузова).

Рассмотрим методы определения приращения работы пластической деформации деталей кузова (кабины) и других конструкций ТС.

Согласно деформационной теории пластичности [7] приращение работы пластической деформации пропорционально приращению (изменению) объема материала и изменения формы изделия:

bess33.wmf

где T – девиатор напряжения (интенсивность касательных напряжений);

Г – интенсивность деформации сдвига.

Если при ДТП не образовались трещины и разрывы, то изменением объема при пластической деформации (в результате увеличения плотности дефектов) dθ можно пренебречь. Тогда

bess34.wmf (10)

Если известны компоненты тензора деформации εij, то Г определяют по следующей формуле:

bess35.wmf (11)

Интегрированием формулы (10) работу пластической деформации не можем определить, так как Т зависит от Г. Следовательно, нужно использовать дополнительные условия или различные теории пластичности. В частности, если девиатор напряжений Т = τs, bess36.wmf равен пределу текучести материала (стали, алюминиевых сплавов), то вместо выражения (10) можно использовать формулу Генке [7]:

bess37.wmf (12)

где Sij – компоненты тензора девиатора напряжений, характеризующие изменение формы (конфигурации) рассматриваемого элемента детали,

bess38.wmf

Нужно отметить, что выражение (12) не взаимно однозначно, т.е. если известны значения εij, то можно найти значения Sij, так как по формуле (11) определяется и девиатор деформации. Если заданы Sij, то невозможно найти εij, поскольку неизвестен девиатор деформации Г.

В необратимых процессах деформации ТС при ДТП очень сложно оценить деформированное состояние, еще сложнее определить напряженное состояние, соответственно, и работу деформации. В теории обработки давлением установлено, что девиатор напряжений пропорционален девиатору скоростей деформации.

Компоненты тензора скоростей деформации ηij определяются дифференцированием dεij:

bess39.wmf

Тогда имеем

bess40.wmf

где dVi и dVj – приращение скоростей сдвига за кульминационный период ДТП отдельного геометрического элемента детали ТС.

В этом случае уравнение теории течения для развитой пластической деформации (уравнение Сен-Венана – Мизеса) имеет следующий вид [7]:

bess41.wmf (13)

где H – характеризует интенсивность скорости деформации сдвига

bess42.wmf; bess43.wmf.

Результаты исследования и их обсуждение

Для определения Sij по формулам (12) и (13) эксперту следует последовательно описывать процесс геометрической деформации элементов автомобиля при ДТП до получения конечной картины в этой зоне. Для этого, как правило, применяют различные математические методы, наиболее известным из которых является метод конечных элементов (МКЭ) [8, 9]. Так, при исследовании ДТП поверхность транспортного средства можно представить как сочетание прямоугольников или треугольников. В подвижной системе координат изменение положения каждого узла можно задать как dUxy и dUyx или как dVxy и dVyx. Соответственно:

bess44.wmf; bess45.wmf; bess46.wmf

или

bess47.wmf; bess48.wmf; bess49.wmf

При наличии скользящего удара нужно иметь в виду, что в формуле (4) коэффициент контактного трения скольжения не будет соответствовать коэффициенту трения скольжения в кинематических парах. Допускается изменение его значения в пределах от 0,2–0,5, в зависимости от ряда внешних факторов [7]. Контактное давление (PК) при этом определяется напряженным состоянием в зоне деформации и площадью контакта поверхностей.

Выводы

1. Обосновано, что для проведения АТЭ при столкновениях и наездах необходимо применять математическую модель, основанную на законах сохранения энергии и импульса; деформационной теории пластичности и динамики удара для определения тормозного пути с учетом дорожных условий; описывающая аналитические зависимости для оценки напряженно-деформированного состояния на основе методов конечных элементов.

2. Эксперт, которому поручено провести АТЭ, обязан принять во внимание все введенные идеализации и допущения, а также оценить степень значимости всех факторов, влияющих на результаты экспертизы.