
Как бы то ни было, метод Эйлера, простейший численный метод решения дифференциальных уравнений, как таковой до сих пор не нашел своего применения в Минюсте, их ученым советом не апробирован и не утвержден. Хотя те, кто так говорит, сами ничего не знают.
Метод Эйлера, как и ему подобные методы, применяются в экспертизе ДТП. Поэтому эта статья полезна для автоэкспертов и адвокатов тем, что в ней просто и ясно показано, как работает модель движения автомобиля в специальных компьютерных программах реконструкции ДТП, типа PC-Crash. А это позволяет, при необходимости, сформулировать ряд вопросов для экспертного исследования или задать несколько неприятных вопросов эксперту.
Решаем простейшее дифференциальное уравнение
Адвокаты (про госавтоэкспертов точно не скажу), возможно, помнят из школьной программы, что такое производная функции. Производная – это тангенс угла наклона касательной к графику функции в некоторой точке. Это ее геометрический смысл.
Если наша функция есть зависимость пути от времени, то производная – это зависимость скорости от времени. Вторая производная от пути по времени – это зависимость ускорения от времени.
Большей информации для задачи движения автомобиля нам не надо.

Для иллюстрации метода Эйлера рассмотрим некоторую функцию y(x)=x2. Ее график показан выше слева, это – парабола. Справа на рисунке – график производной этой функции y’(x)=2x.
Видно, что, например, при x=0 производная y’=0, или касательная параллельна оси абсцисс. А при x=2 производная y’=4, или в этой точке тангенс угла наклона функции y(x) равен 4, откуда угол наклона касательной составляет около 76 градусов. Если сама функция y(x) нам не известна, а известна только ее производная, то получить эту функцию можно интегрированием производной.
Если мы проинтегрируем производную аналитически, то получим формулу искомой функции y(x), а если численно – только ее график в виде набора точек. Давайте проинтегрируем численно методом Эйлера.
Полагаем, а это – граничные условия, что искомая функция начинается с точки x=0 и y=0. Наша задача – получить правую ветвь этой параболы.
Принимаем шаг интегрирования dx=0.5. Тогда в этой точке, при x=0, значение производной тоже равно нулю 2x=2x0=0, и мы попадаем в следующую точку x=0+0.5=0.5 и y=0+2x0x0.5=0. В этой второй точке x=0.5 значение производной уже равно 2х0.5=1, поэтому далее мы попадаем в третью точку с координатами x=0.5+0.5=1 и y=0+1x1x0.5=0.5. И так далее.
Табличка со значениями аргумента x и функции y, а так же график реальной функции y=x2 (красный цвет) и ломанной (синий цвет), полученной методом Эйлера показан на рисунке ниже.

Это называется сходимостью метода – при уменьшении шага возрастает точность. Сравнивая два решения с разным шагом, можно сделать вывод о погрешности полученного решения. А сейчас, с современными компьютерами, проблем с выбором весьма малого шага интегрирования нет, так как человеку все равно, длился расчет 0.1 секунды или 0.01 секунды, лишь бы не часами.
Для примера на следующем рисунке показан результат интегрирования с шагом dx=0.01, и видно, что красная и синяя линии практически слились. Там же показана программа расчета: задаем начальную точку и шаг интегрирования, а далее в цикле для параметра i, пробегающего значения от 1 до 300 идет расчет точек ломаной по методу Эйлера.

Интегрируем дифференциальные уравнения движения автомобиля
Рассмотрим решение задачи движения полностью заторможенного автомобиля Ауди массой m=1185 кг, моментом инерции I=1566.4 кг ·м2, по горизонтальной поверхности с коэффициентом сцепления колес с поверхностью дороги f=0.7 (сухой асфальт). Скорость автомобиля в начальный момент v0=30 км/ч в направлении вправо на рисунке, скорость вращения w0=3c-1 (радиан в секунду). Такое движение может быть как в результате удара и тогда это движение после разделения автомобилей, либо занос в колее и т.п.
На рисунке ниже показано решение этой задачи в программе PC-Crash (рекомендованной Минюстом).

Начальное значение времени положим t=0, шаг интегрирования по времени dt=0.001 сек. Мы не знаем значения ускорений центра тяжести автомобиля вдоль координатных осей и угловое ускорение в начальный момент времени, поэтому их значения полаем также равными нулю: acx=0, acy=0 и e=0.
Далее будем считать, что левой переднее колесо автомобиля – это колесо 1, правое переднее – 2, правое заднее – 3 и левое заднее – 4. С помощью таблички с параметрами автомобиля определим начальные координаты его колес: x1=1.25 м, y1=0.755 м, x2=1.25 м, y2=-0.755 м, x3=-1.26 м, y3=-0.755 м, x4=1.26 м, y4=0.755 м. Для расчетов нам понадобятся длины и углы с осью x радиус-векторов всех колес.
Для передних колес это R1=1.46 м, для задних – R2=1.469 м, и углы а1=0.543 рад (31.13 град), а2=-0.543 рад (-31.13 град), а3=3.68 рад (210.9 град), а4=2.6 рад (149 град). Все исходные данные приведены в листинге ниже.





При реальном расчете эксперт может принять и иные соотношения. Так же вычисляются проекции сил сцепления на каждом колесе на оси глобальной системы координат как произведение величины указанной силы на компоненты ортов направления движения колеса со знаком «минус».
После этого (в пятой строке) проекции всех сил на оси координат суммируются. Далее вычисляется момент этих сил относительно центра тяжести как векторное произведение сил сцепления на радиус-вектор каждого колеса.


Интегрирование закончено.
Проверяем результаты
Результат интегрирования в данном случае это таблица, содержащая 1351 строку, в столбцах которой содержатся все те величины, которые мы рассчитали, то есть имеющие в листингах нижний индекс i. Таблица сама по себе суду ничего не дает, хотя некоторые адвокаты, услышав про нее, требуют представить, что я с удовольствием делаю – представляю на электронном носителе. Но лучше визуализировать результат в виде расчетного положения автомобиля в разные моменты времени и траектории его колес. В программе выше мы не вычисляли положения углов автомобиля, поэтому представим в виде графика расчетные траектории всех 4-х колес автомобиля с координатной сеткой 1м.
На рисунке ниже траектория левого переднего колеса красная, правого переднего синяя, правого заднего фиолетовая, левого заднего черная.


Обсуждение результатов
Итак, приведенный пример расчета движения заторможенного автомобиля в заносе показал, что не боги горшки обжигают, и развеял мрак тайны специальных компьютерных программ. Это немаловажно как для экспертов, так и для адвокатов, так как, надеюсь, с учетом цикла статей по алгоритму CRASH3 и уравнениям движения им стало понятно, как это все работает. Понятно тем, кто это читает. А остальным?
Для чего может понадобиться такой «ручной» расчет? Он нужен, когда положения колес автомобиля после удара существенно отличаются от их положений до удара, когда одно колесо оторвано, когда колеса после удара смотрят в разные стороны, когда одно колесо или больше спущено, когда некоторые из колес заклинены деформированными частями автомобиля, а другие свободно или со скрежетом вращаются. И это только небольшое число случаев из тех, которые могут возникнуть в реальном ДТП, и не могут быть учтены той или иной компьютерной программой.
Выше рассмотрено движение полностью заторможенного автомобиля, поэтому модель шин была простейшей – силы сцепления колес были направлены против направления движения каждого колеса. А если в реальности по другому? Какая и почему использована модель в конкретном расчете?
Резюмируя, можно сказать снова и снова – лишними знания не бывают, и теорию знать полезно.