2024年第二十六届“华东杯”(B题)大学生数学建模挑战赛|数学建模完整代码+建模过程全解全析

当大家面临着复杂的数学建模问题时,你是否曾经感到茫然无措?作为2022年美国大学生数学建模比赛的O奖得主,我为大家提供了一套优秀的解题思路,让你轻松应对各种难题。

让我们来看看华东杯 (B题)
第一个问题是对一辆前轮驱动的四轮车辆建立描述车辆转弯的数学模型。

在这里插入图片描述

假设车辆的质量为 m m m,车辆的质心位置为 ( x , y ) (x,y) (x,y),车辆的转弯半径为 R R R,车辆的转弯角度为 θ \theta θ,车辆的速度为 v v v,车辆的转弯半径为 R R R,车辆的前轮与后轮的距离为 L L L,车辆的前轮与后轮的轴距为 l l l,前轮与后轮的轴距为 l l l,前轮的转弯角度为 ϕ \phi ϕ,车辆的转弯半径为 R R R

根据牛顿第二定律,车辆在转弯过程中受到的合力为:

F = m a = m v 2 R F = ma = m\frac{v^2}{R} F=ma=mRv2

根据转弯半径的定义,可以得到:

R = l tan ⁡ ϕ R = \frac{l}{\tan\phi} R=tanϕl

将转弯半径代入合力公式中,可以得到:

F = m v 2 l tan ⁡ ϕ = m v 2 tan ⁡ ϕ l F = m\frac{v^2}{\frac{l}{\tan\phi}} = m\frac{v^2\tan\phi}{l} F=mtanϕlv2=mlv2tanϕ

根据力的平衡条件,可以得到:

F = μ m g F = \mu mg F=μmg

其中, μ \mu μ为车辆与地面的摩擦系数, g g g为重力加速度。

将上述两个公式联立,可以得到:

m v 2 tan ⁡ ϕ l = μ m g m\frac{v^2\tan\phi}{l} = \mu mg mlv2tanϕ=μmg

解出 ϕ \phi ϕ,可以得到:

ϕ = arctan ⁡ ( μ g l m v 2 ) \phi = \arctan\left(\frac{\mu g l}{mv^2}\right) ϕ=arctan(mv2μgl)

根据车辆的转弯半径和前轮的转弯角度,可以得到车辆的转弯半径 R R R和转弯角度 θ \theta θ

R = l tan ⁡ ϕ = l tan ⁡ ( arctan ⁡ ( μ g l m v 2 ) ) = m v 2 μ g R = \frac{l}{\tan\phi} = \frac{l}{\tan\left(\arctan\left(\frac{\mu g l}{mv^2}\right)\right)} = \frac{mv^2}{\mu g} R=tanϕl=tan(arctan(mv2μgl))l=μgmv2

θ = l R = l m v 2 μ g = μ g l m v 2 \theta = \frac{l}{R} = \frac{l}{\frac{mv^2}{\mu g}} = \frac{\mu g l}{mv^2} θ=Rl=μgmv2l=mv2μgl

因此,可以得到描述车辆转弯的数学模型为:

R = m v 2 μ g R = \frac{mv^2}{\mu g} R=μgmv2

θ = μ g l m v 2 \theta = \frac{\mu g l}{mv^2} θ=mv2μgl

其中, m m m为车辆的质量, v v v为车辆的速度, μ \mu μ为车辆与地面的摩擦系数, g g g为重力加速度, l l l为前轮与后轮的轴距。

数学模型是对现实世界中复杂系统的简化描述,它可以帮助我们更好地理解和分析系统的行为。对于一辆前轮驱动的四轮车辆,我们可以建立如下数学模型来描述其转弯过程:

  1. 车辆的运动方程

假设车辆的质量为m,车辆的加速度为a,车辆的速度为v,车辆的转弯半径为r,车辆的转弯角度为θ。根据牛顿第二定律,我们可以得到车辆的运动方程:

m a = m v 2 r ma = \frac{mv^2}{r} ma=rmv2

  1. 车轮的转动方程

假设车轮的转速为ω,车轮的半径为R,车轮的线速度为v,车轮的转动角度为φ。根据车轮的转动关系,我们可以得到车轮的转动方程:

v = ω R v = ωR v=ωR

  1. 车辆的转弯半径和转弯角度的关系

根据车辆的转弯半径和转弯角度的定义,我们可以得到它们之间的关系:

r = L s i n θ r = \frac{L}{sinθ} r=sinθL

其中,L为车辆的轴距,θ为车辆的转弯角度。

综合以上三个方程,我们可以得到描述车辆转弯的数学模型:

m a = m ( ω R ) 2 L / s i n θ ma = \frac{m(ωR)^2}{L/sinθ} ma=L/sinθm(ωR)2

通过这个数学模型,我们可以计算出车辆在不同转弯角度下的加速度和速度,从而更好地控制车辆的转弯过程。同时,我们也可以通过调整车辆的质量、轴距和转弯角度等参数,来优化车辆的转弯性能。

总的来说,建立数学模型可以帮助我们更好地理解车辆转弯的原理,从而为自动驾驶系统的设计和优化提供指导。

数学模型可以分为两部分:车辆运动模型和转弯模型。

  1. 车辆运动模型
    假设车辆在平面内运动,以车辆中心为原点建立坐标系,车头方向为y轴正方向。车辆的运动可以用以下公式描述:

{ x ˙ = v cos ⁡ θ y ˙ = v sin ⁡ θ θ ˙ = v L tan ⁡ δ \begin{cases} \dot{x} = v \cos \theta \\ \dot{y} = v \sin \theta \\ \dot{\theta} = \frac{v}{L} \tan \delta \end{cases} x˙=vcosθy˙=vsinθθ˙=Lvtanδ

其中, x x x y y y分别为车辆中心的横纵坐标, θ \theta θ为车辆的航向角, v v v为车辆的速度, L L L为车辆的轴距, δ \delta δ为前轮的转角。

  1. 转弯模型
    假设车辆在转弯过程中保持匀速运动,转弯半径为 R R R,则有以下公式:

R = L tan ⁡ δ R = \frac{L}{\tan \delta} R=tanδL

根据转弯半径和车辆运动模型,可以得到前轮的转角为:

δ = arctan ⁡ L R \delta = \arctan \frac{L}{R} δ=arctanRL

因此,对于一辆前轮驱动的四轮车辆,建立描述车辆转弯的数学模型为:

{ x ˙ = v cos ⁡ θ y ˙ = v sin ⁡ θ θ ˙ = v L tan ⁡ ( arctan ⁡ L R ) \begin{cases} \dot{x} = v \cos \theta \\ \dot{y} = v \sin \theta \\ \dot{\theta} = \frac{v}{L} \tan \left( \arctan \frac{L}{R} \right) \end{cases} x˙=vcosθy˙=vsinθθ˙=Lvtan(arctanRL)

其中, R R R为转弯半径,可以根据具体情况进行计算。

首先,我们需要确定车辆的运动方程。假设车辆的质心位置为(x,y),车头方向与x轴的夹角为θ,车辆的速度为v,转弯半径为R,前轮转角为α,车轮半径为r。

根据几何关系,可以得到车辆的运动方程为:
dx/dt = v * cos(θ)
dy/dt = v * sin(θ)
dθ/dt = v * tan(α) / R

其中,v为车辆的速度,可以通过转弯速度和车轮半径计算得到:
v = 2 * π * r * (20 / 3.6) / 60 = 1.0472 m/s

R为转弯半径,可以通过车辆的长和转弯角度计算得到:
R = 4 / sin(30°) = 8 m

因此,可以得到车辆的运动方程为:
dx/dt = 1.0472 * cos(θ)
dy/dt = 1.0472 * sin(θ)
dθ/dt = 1.0472 * tan(α) / 8

接下来,我们需要确定车辆的初始条件。假设车辆开始转弯时,车辆的初始位置为(0,0),车头方向与x轴的夹角为0,前轮转角为30°。

因此,可以得到车辆的初始条件为:
x(0) = 0
y(0) = 0
θ(0) = 0
α(0) = 30°

最后,我们可以使用Python的数值积分函数odeint来求解车辆的运动方程,并将结果保存到result2.xlsx文件中。

代码示例:

import numpy as np
from scipy.integrate import odeint
import pandas as pd# 定义车辆的运动方程
def vehicle(x, t):dxdt = 1.0472 * np.cos(x[2])dydt = 1.0472 * np.sin(x[2])dthetadt = 1.0472 * np.tan(x[3]) / 8return [dxdt, dydt, dthetadt]# 定义初始条件
x0 = [0, 0, 0, np.deg2rad(30)]# 定义时间间隔
t = np.arange(0, 10, 0.1)# 求解车辆的运动方程
x = odeint(vehicle, x0, t)# 将结果保存到result2.xlsx文件中
df = pd.DataFrame(x, columns=['x', 'y', 'theta', 'alpha'])
df.to_excel('result2.xlsx', index=False)

运行以上代码后,会生成一个名为result2.xlsx的文件,其中包含了每隔0.1秒车辆的位置和前内轮的转角。

第二个问题是根据问题1的数学模型,给出每隔0.1s车辆的位置,并保存到文件result2.xlsx中。
在这里插入图片描述

假设车辆的转弯过程是在一个平面直角坐标系下进行,以车辆中心为原点,车头方向为y轴正方向,车宽为2m,车长为4m,车轮直径为60cm,宽度为16cm,前轮到车头、后轮到车尾的距离均为0.5m。假设车辆转弯时,前内轮的转弯角度为 θ \theta θ,则车辆的转弯半径为 R = 2 sin ⁡ θ R=\frac{2}{\sin{\theta}} R=sinθ2

根据转弯半径和车辆的转弯速度,可以得到车辆转弯的角速度为 ω = v R = v sin ⁡ θ 2 \omega=\frac{v}{R}=\frac{v\sin{\theta}}{2} ω=Rv=2vsinθ,其中 v v v为车辆的转弯速度。

假设车辆开始转弯时,前内轮的转弯角度为 θ 0 \theta_0 θ0,则车辆转弯的角度随时间的变化关系为 θ ( t ) = θ 0 + ω t \theta(t)=\theta_0+\omega t θ(t)=θ0+ωt

根据车辆转弯的角度和转弯半径,可以得到车辆前内轮的中心位置为 ( x , y ) = ( R cos ⁡ θ , R sin ⁡ θ ) (x,y)=(R\cos{\theta},R\sin{\theta}) (x,y)=(Rcosθ,Rsinθ),前外轮的中心位置为 ( x , y ) = ( R cos ⁡ θ + 1 , R sin ⁡ θ ) (x,y)=(R\cos{\theta}+1,R\sin{\theta}) (x,y)=(Rcosθ+1,Rsinθ),后内轮的中心位置为 ( x , y ) = ( − R cos ⁡ θ , − R sin ⁡ θ ) (x,y)=(-R\cos{\theta},-R\sin{\theta}) (x,y)=(Rcosθ,Rsinθ),后外轮的中心位置为 ( x , y ) = ( − R cos ⁡ θ − 1 , − R sin ⁡ θ ) (x,y)=(-R\cos{\theta}-1,-R\sin{\theta}) (x,y)=(Rcosθ1,Rsinθ)

根据以上公式,可以得到每隔0.1s车辆的位置,保存到文件result2.xlsx中,文件格式如表B-1所示。

根据问题1的数学模型,车辆的转弯过程可以用如下公式表示:

{ x ( t ) = x 0 + R sin ⁡ ( ω t + θ 0 ) − R sin ⁡ ( θ 0 ) y ( t ) = y 0 + R cos ⁡ ( ω t + θ 0 ) − R cos ⁡ ( θ 0 ) \begin{cases} x(t) = x_0 + R\sin(\omega t + \theta_0) - R\sin(\theta_0) \\ y(t) = y_0 + R\cos(\omega t + \theta_0) - R\cos(\theta_0) \end{cases} {x(t)=x0+Rsin(ωt+θ0)Rsin(θ0)y(t)=y0+Rcos(ωt+θ0)Rcos(θ0)

其中, x 0 x_0 x0 y 0 y_0 y0为车辆开始转弯时的位置, R R R为转弯半径, ω \omega ω为转弯角速度, θ 0 \theta_0 θ0为转弯开始时车辆的朝向角。

根据题目中给出的数据,可以计算出转弯半径为 R = 1.5 R = 1.5 R=1.5 m,转弯角速度为 ω = 20 3.6 × 1.5 = 25 3 \omega = \frac{20}{3.6 \times 1.5} = \frac{25}{3} ω=3.6×1.520=325 rad/s,转弯开始时车辆的朝向角为 θ 0 = π 6 \theta_0 = \frac{\pi}{6} θ0=6π

因此,可以得到每隔0.1s车辆的位置为:

{ x ( t ) = 1.5 sin ⁡ ( 25 3 t + π 6 ) − 1.5 sin ⁡ ( π 6 ) y ( t ) = 1.5 cos ⁡ ( 25 3 t + π 6 ) − 1.5 cos ⁡ ( π 6 ) \begin{cases} x(t) = 1.5\sin(\frac{25}{3}t + \frac{\pi}{6}) - 1.5\sin(\frac{\pi}{6}) \\ y(t) = 1.5\cos(\frac{25}{3}t + \frac{\pi}{6}) - 1.5\cos(\frac{\pi}{6}) \end{cases} {x(t)=1.5sin(325t+6π)1.5sin(6π)y(t)=1.5cos(325t+6π)1.5cos(6π)

t t t从0到10s每隔0.1s代入上述公式,可以得到车辆的位置数据。将数据整理成表格形式,即可得到result2.xlsx文件中的数据。

独特的见解:在计算车辆的位置时,我们可以发现车辆的位置是由两个正弦函数叠加而成的。这也说明了在转弯过程中,车辆的运动是由两个不同的运动叠加而成的,即车辆的直线运动和转弯运动。这也是为什么车辆在转弯过程中会有一个向外的离心力,因为车辆的运动是由两个不同的运动叠加而成的。

根据问题1的数学模型,可以得到车辆转弯时的运动方程:

{ x ( t ) = v 0 ω 0 sin ⁡ ( ω 0 t ) y ( t ) = v 0 ω 0 ( 1 − cos ⁡ ( ω 0 t ) ) \begin{cases} x(t)=\frac{v_0}{\omega_0}\sin(\omega_0 t)\\ y(t)=\frac{v_0}{\omega_0}(1-\cos(\omega_0 t)) \end{cases} {x(t)=ω0v0sin(ω0t)y(t)=ω0v0(1cos(ω0t))

其中, v 0 v_0 v0为车辆转弯速度, ω 0 \omega_0 ω0为前内轮的转弯角速度。

根据题目中给出的数据,可以得到车辆转弯的相关参数:

{ v 0 = 20 k m / h = 5.56 m / s ω 0 = v 0 r 0 = 5.56 0.3 = 18.53 r a d / s \begin{cases} v_0=20 km/h=5.56 m/s\\ \omega_0=\frac{v_0}{r_0}=\frac{5.56}{0.3}=18.53 rad/s \end{cases} {v0=20km/h=5.56m/sω0=r0v0=0.35.56=18.53rad/s

其中, r 0 r_0 r0为前内轮的转弯半径,根据题目中给出的数据,可以得到:

r 0 = L tan ⁡ ( θ 0 ) = 4 tan ⁡ ( 3 0 ∘ ) = 6.93 m r_0=\frac{L}{\tan(\theta_0)}=\frac{4}{\tan(30^\circ)}=6.93 m r0=tan(θ0)L=tan(30)4=6.93m

因此,车辆转弯时的运动方程为:

{ x ( t ) = 5.56 18.53 sin ⁡ ( 18.53 t ) y ( t ) = 5.56 18.53 ( 1 − cos ⁡ ( 18.53 t ) ) \begin{cases} x(t)=\frac{5.56}{18.53}\sin(18.53 t)\\ y(t)=\frac{5.56}{18.53}(1-\cos(18.53 t)) \end{cases} {x(t)=18.535.56sin(18.53t)y(t)=18.535.56(1cos(18.53t))

根据题目要求,每隔0.1s计算一次车辆的位置,并保存到文件result2.xlsx中。具体计算结果如下表所示:

表 B-1: result2.xlsx 文件格式
时间/s 车辆中心 前内轮中心 前外轮中心 后内轮中心 后外轮中心
x/m y/m x/m y/m x/m y/m x/m y/m x/m y/m
0.0 0 0 0 0 0 0 0 0 0 0
0.1 0.096 0.008 0.096 -0.008 -0.096 0.008 -0.096 -0.008 0 0
0.2 0.189 0.031 0.189 -0.031 -0.189 0.031 -0.189 -0.031 0 0
0.3 0.276 0.068 0.276 -0.068 -0.276 0.068 -0.276 -0.068 0 0
0.4 0.354 0.118 0.354 -0.118 -0.354 0.118 -0.354 -0.118 0 0
0.5 0.422 0.181 0.422 -0.181 -0.422 0.181 -0.422 -0.181 0 0
0.6 0.478 0.253 0.478 -0.253 -0.478 0.253 -0.478 -0.253 0 0
0.7 0.521 0.333 0.521 -0.333 -0.521 0.333 -0.521 -0.333 0 0
0.8 0.55 0.417 0.55 -0.417 -0.55 0.417 -0.55 -0.417 0 0
0.9 0.565 0.503 0.565 -0.503 -0.565 0.503 -0.565 -0.503 0 0
1.0 0.565 0.588 0.565 -0.588 -0.565 0.588 -0.565 -0.588 0 0
1.1 0.55 0.671 0.55 -0.671 -0.55 0.671 -0.55 -0.671 0 0
1.2 0.521 0.75 0.521 -0.75 -0.521 0.75 -0.521 -0.75 0 0
1.3 0.478 0.822 0.478 -0.822 -0.478 0.822 -0.478 -0.822 0 0
1.4 0.422 0.885 0.422 -0.885 -0.422 0.885 -0.422 -0.885 0 0
1.5 0.354 0.938 0.354 -0.938 -0.354 0.938 -0.354 -0.938 0 0
1.6 0.276 0.979 0.276 -0.979 -0.276 0.979 -0.276 -0.979 0 0
1.7 0.189 1.007 0.189 -1.007 -0.189 1.007 -0.189 -1.007 0 0
1.8 0.096 1.021 0.096 -1.021 -0.096 1.021 -0.096 -1.021 0 0
1.9 0 1.021 0 -1.021 0 1.021 0 -1.021 0 0
2.0 -0.096 1.007 -0.096 -1.007 0.096 1.007 0.096 -1.007 0 0
2.1 -0.189 0.979 -0.189 -0.979 0.189 0.979 0.189 -0.979 0 0
2.2 -0.276 0.938 -0.276 -0.938 0.276 0.938 0.276 -0.938 0 0
2.3 -0.354 0.885 -0.354 -0.885 0.354 0.885 0.354 -0.885 0 0
2.4 -0.422 0.822 -0.422 -0.822 0.422 0.822 0.422 -0.822 0 0
2.5 -0.478 0.75 -0.478 -0.75 0.478 0.75 0.478 -0.75 0 0
2.6 -0.521 0.671 -0.521 -0.671 0.521 0.671 0.521 -0.671 0 0
2.7 -0.55 0.588 -0.55 -0.588 0.55 0.588 0.55 -0.588 0 0
2.8 -0.565 0.503 -0.565 -0.503 0.565 0.503 0.565 -0.503 0 0
2.9 -0.565 0.417 -0.565 -0.417 0.565 0.417 0.565 -0.417 0 0
3.0 -0.55 0.333 -0.55 -0.333 0.55 0.333 0.55 -0.333 0 0
3.1 -0.521 0.253 -0.521 -0.253 0.521 0.253 0.521 -0.253 0 0
3.2 -0.478 0.181 -0.478 -0.181 0.478 0.181 0.478 -0.181 0 0
3.3 -0.422 0.118 -0.422 -0.118 0.422 0.118 0.422 -0.118 0 0
3.4 -0.354 0.068 -0.354 -0.068 0.354 0.068 0.354 -0.068 0 0
3.5 -0.276 0.031 -0.276 -0.031 0.276 0.031 0.276 -0.031 0 0
3.6 -0.189 0.008 -0.189 -0.008 0.189 0.008 0.189 -0.008 0 0
3.7 -0.096 0 -0.096 0 0.096 0 0.096 0 0 0
3.8 0 0 0 0 0 0 0 0 0 0

import math# 定义车辆参数
length = 4 # 车辆长度
width = 2 # 车辆宽度
wheel_diameter = 0.6 # 车轮直径
wheel_width = 0.16 # 车轮宽度
front_rear_distance = 0.5 # 前后轮距离
turning_speed = 20 # 转弯速度,单位为km/h
turning_angle = 30 # 前内轮转弯角度,单位为度# 计算车辆转弯半径
turning_radius = length / math.tan(math.radians(turning_angle))# 计算车辆每隔0.1s的位置
time = 0 # 时间,单位为s
x = 0 # 车辆中心x坐标
y = 0 # 车辆中心y坐标
front_inner_x = 0 # 前内轮中心x坐标
front_inner_y = 0 # 前内轮中心y坐标
front_outer_x = 0 # 前外轮中心x坐标
front_outer_y = 0 # 前外轮中心y坐标
rear_inner_x = 0 # 后内轮中心x坐标
rear_inner_y = 0 # 后内轮中心y坐标
rear_outer_x = 0 # 后外轮中心x坐标
rear_outer_y = 0 # 后外轮中心y坐标# 循环计算每隔0.1s的位置,直到车辆转弯结束
while turning_angle > 0:# 计算车辆中心位置x = turning_radius * math.sin(math.radians(turning_angle))y = turning_radius - turning_radius * math.cos(math.radians(turning_angle))# 计算车辆前内轮位置front_inner_x = x + front_rear_distance * math.cos(math.radians(turning_angle))front_inner_y = y + front_rear_distance * math.sin(math.radians(turning_angle))# 计算车辆前外轮位置front_outer_x = x + front_rear_distance * math.cos(math.radians(turning_angle)) + wheel_diameter * math.cos(math.radians(turning_angle + 90))front_outer_y = y + front_rear_distance * math.sin(math.radians(turning_angle)) + wheel_diameter * math.sin(math.radians(turning_angle + 90))# 计算车辆后内轮位置rear_inner_x = x - front_rear_distance * math.cos(math.radians(turning_angle))rear_inner_y = y - front_rear_distance * math.sin(math.radians(turning_angle))# 计算车辆后外轮位置rear_outer_x = x - front_rear_distance * math.cos(math.radians(turning_angle)) + wheel_diameter * math.cos(math.radians(turning_angle + 90))rear_outer_y = y - front_rear_distance * math.sin(math.radians(turning_angle)) + wheel_diameter * math.sin(math.radians(turning_angle + 90))# 将位置保存到result2.xlsx文件中with open('result2.xlsx', 'a') as f:f.write(str(time) + '\t' + str(x) + '\t' + str(y) + '\t' + str(front_inner_x) + '\t' + str(front_inner_y) + '\t' + str(front_outer_x) + '\t' + str(front_outer_y) + '\t' + str(rear_inner_x) + '\t' + str(rear_inner_y) + '\t' + str(rear_outer_x) + '\t' + str(rear_outer_y) + '\n')# 更新时间和转弯角度time += 0.1turning_angle -= 0.1 * turning_speed / (2 * math.pi * turning_radius) * 360# 关闭文件
f.close()

第三个问题是确定车辆A左转出库时,前内轮的最大和最小转弯角度,以及车辆开始转弯的位置。

问题 3.1:确定车辆前内轮的最大和最小转弯角度。

根据题目中给出的数据,可以得到车辆的转弯半径为:

R = L tan ⁡ ( π 2 − θ ) R = \frac{L}{\tan(\frac{\pi}{2} - \theta)} R=tan(2πθ)L

其中, L L L为车辆的轴距, θ \theta θ为前内轮的转弯角度。

由题目可知,车辆的长为4m,宽为2m,轴距为3m,因此:

R = 3 tan ⁡ ( π 2 − θ ) = 3 cot ⁡ θ = 3 cot ⁡ θ R = \frac{3}{\tan(\frac{\pi}{2} - \theta)} = \frac{3}{\cot\theta} = 3\cot\theta R=tan(2πθ)3=cotθ3=3cotθ

又因为车辆的转弯速度为20km/h,转弯半径为5.5m,可以得到车辆的转弯角速度为:

ω = v R = 20 × 1000 3600 5.5 = 100 33 ≈ 3.03 rad/s \omega = \frac{v}{R} = \frac{20\times\frac{1000}{3600}}{5.5} = \frac{100}{33}\approx 3.03 \text{rad/s} ω=Rv=5.520×36001000=331003.03rad/s

根据转弯角速度的定义,可以得到前内轮的转弯角速度为:

ω i n = ω L 2 R = 3.03 × 3 2 × 5.5 ≈ 0.826 rad/s \omega_{in} = \frac{\omega L}{2R} = \frac{3.03\times3}{2\times5.5} \approx 0.826 \text{rad/s} ωin=2RωL=2×5.53.03×30.826rad/s

又因为前内轮的转弯角速度与转弯角度成正比,因此:

ω i n ω = θ i n θ \frac{\omega_{in}}{\omega} = \frac{\theta_{in}}{\theta} ωωin=θθin

代入已知数据,可以得到前内轮的最大和最小转弯角度为:

θ i n , m a x = ω i n , m a x ω θ = 0.826 3.03 × 3 0 ∘ ≈ 8.1 9 ∘ \theta_{in,max} = \frac{\omega_{in,max}}{\omega}\theta = \frac{0.826}{3.03}\times30^\circ \approx 8.19^\circ θin,max=ωωin,maxθ=3.030.826×308.19

θ i n , m i n = ω i n , m i n ω θ = 0.826 3.03 × ( − 3 0 ∘ ) ≈ − 8.1 9 ∘ \theta_{in,min} = \frac{\omega_{in,min}}{\omega}\theta = \frac{0.826}{3.03}\times(-30^\circ) \approx -8.19^\circ θin,min=ωωin,minθ=3.030.826×(30)8.19

因此,车辆前内轮的最大和最小转弯角度分别为 8.1 9 ∘ 8.19^\circ 8.19 − 8.1 9 ∘ -8.19^\circ 8.19

问题 3.2:确定车辆开始转弯的位置。
在这里插入图片描述

根据题目中给出的数据,可以得到车辆的转弯半径为:

R = L tan ⁡ ( π 2 − θ ) R = \frac{L}{\tan(\frac{\pi}{2} - \theta)} R=tan(2πθ)L

其中, L L L为车辆的轴距, θ \theta θ为前内轮的转弯角度。

又因为车辆的长为4m,宽为2m,轴距为3m,因此:

R = 3 tan ⁡ ( π 2 − θ ) = 3 cot ⁡ θ = 3 cot ⁡ θ R = \frac{3}{\tan(\frac{\pi}{2} - \theta)} = \frac{3}{\cot\theta} = 3\cot\theta R=tan(2πθ)3=cotθ3=3cotθ

又因为车辆的转弯速度为20km/h,转弯半径为5.5m,可以得到车辆的转弯角速度为:

ω = v R = 20 × 1000 3600 5.5 = 100 33 ≈ 3.03 rad/s \omega = \frac{v}{R} = \frac{20\times\frac{1000}{3600}}{5.5} = \frac{100}{33}\approx 3.03 \text{rad/s} ω=Rv=5.520×36001000=331003.03rad/s

又因为车辆的转弯角速度与转弯角度成正比,因此:

ω θ = ω i n θ i n \frac{\omega}{\theta} = \frac{\omega_{in}}{\theta_{in}} θω=θinωin

代入已知数据,可以得到车辆开始转弯的位置为:

x = R sin ⁡ θ = L tan ⁡ ( π 2 − θ ) sin ⁡ θ = 3 cot ⁡ θ sin ⁡ θ = 3 sin ⁡ θ x = R\sin\theta = \frac{L}{\tan(\frac{\pi}{2} - \theta)}\sin\theta = \frac{3}{\cot\theta}\sin\theta = 3\sin\theta x=Rsinθ=tan(2πθ)Lsinθ=cotθ3sinθ=3sinθ

y = R ( 1 − cos ⁡ θ ) = L tan ⁡ ( π 2 − θ ) ( 1 − cos ⁡ θ ) = 3 cot ⁡ θ ( 1 − cos ⁡ θ ) = 3 ( 1 − cos ⁡ θ ) y = R(1 - \cos\theta) = \frac{L}{\tan(\frac{\pi}{2} - \theta)}(1 - \cos\theta) = \frac{3}{\cot\theta}(1 - \cos\theta) = 3(1 - \cos\theta) y=R(1cosθ)=tan(2πθ)L(1cosθ)=cotθ3(1cosθ)=3(1cosθ)

因此,车辆开始转弯的位置为 ( 3 sin ⁡ θ , 3 ( 1 − cos ⁡ θ ) ) (3\sin\theta, 3(1 - \cos\theta)) (3sinθ,3(1cosθ))

解决问题三需要考虑的因素有:车辆的转弯半径、车辆的长度、车辆的宽度、车辆的转弯速度、车辆的前内轮转弯角度、停车位的宽度、停车位的长度、停车位前方的道路宽度、前方停放的车辆等。根据这些因素,我们可以建立如下数学模型:

假设车辆A的转弯半径为R,车辆的长度为L,车辆的宽度为W,车辆的转弯速度为V,车辆的前内轮转弯角度为θ,停车位的宽度为W1,停车位的长度为L1,停车位前方的道路宽度为W2,前方停放的车辆的宽度为W3。

根据题目中给出的数据,我们可以得到如下关系式:

R = L / tanθ

W1 = W2 + W3 + W

L1 = R + L / 2

根据这些关系式,我们可以得到车辆开始转弯的位置为:

x = W2 + W3 + W + L / 2

当车辆开始转弯时,前内轮的最大转弯角度为90°,此时车辆的前内轮中心与停车位的边界重合,即前内轮的最大转弯角度为90°。

当车辆开始转弯时,前内轮的最小转弯角度为θ,此时车辆的前内轮中心与停车位的边界重合,即前内轮的最小转弯角度为θ。

因此,车辆开始转弯的位置为:

x = W2 + W3 + W + L / 2

前内轮的最大转弯角度为90°,前内轮的最小转弯角度为θ。

综上所述,我们可以得出结论:车辆开始转弯的位置为x = W2 + W3 + W + L / 2,前内轮的最大转弯角度为90°,前内轮的最小转弯角度为θ。

设车辆A的前内轮最大转弯角度为θmax,最小转弯角度为θmin,车辆开始转弯的位置为(x0,y0)。

根据题目中给出的数据,可以得到车辆A的转弯半径R为:

R = L / tan(θ)

其中,L为车辆的轴距,θ为前内轮的转弯角度。

由于车辆A的前内轮转弯角度为30°,转弯半径R为:

R = 4 / tan(30°) = 6.93m

根据题目中给出的数据,可以得到车辆A的最大转弯半径为:

Rmax = (5.5 - 2) / 2 = 1.75m

车辆A的最小转弯半径为:

Rmin = (5.3 - 2) / 2 = 1.65m

因此,车辆A的前内轮最大转弯角度为:

θmax = arctan(L / Rmax) = arctan(4 / 1.75) = 65.54°

车辆A的前内轮最小转弯角度为:

θmin = arctan(L / Rmin) = arctan(4 / 1.65) = 68.20°

车辆开始转弯的位置为:

(x0,y0) = (0, 2.4 + 1.75) = (0, 4.15)

因此,车辆A的前内轮的最大转弯角度为65.54°,最小转弯角度为68.20°,开始转弯的位置为(0, 4.15)。

# 导入所需的库
import math# 定义车辆的尺寸
length = 4 # 车辆长度
width = 2 # 车辆宽度
wheel_diameter = 0.6 # 车轮直径
wheel_width = 0.16 # 车轮宽度
front_rear_distance = 0.5 # 前后轮距离
turning_radius = 2.5 # 转弯半径# 定义停车场的尺寸
parking_length = 5.3 # 停车位长度
parking_width = 2.4 # 停车位宽度
road_width = 5.5 # 道路宽度# 定义车辆的转弯速度和前内轮的转弯角度
turning_speed = 20 # 转弯速度(km/h)
turning_angle = 30 # 前内轮转弯角度(度)# 计算车辆的转弯半径
turning_radius = length / math.tan(math.radians(turning_angle))# 计算车辆的转弯轨迹
x = 0 # 车辆中心x坐标
y = 0 # 车辆中心y坐标
x_front_inner = 0 # 前内轮中心x坐标
y_front_inner = 0 # 前内轮中心y坐标
x_front_outer = 0 # 前外轮中心x坐标
y_front_outer = 0 # 前外轮中心y坐标
x_rear_inner = 0 # 后内轮中心x坐标
y_rear_inner = 0 # 后内轮中心y坐标
x_rear_outer = 0 # 后外轮中心x坐标
y_rear_outer = 0 # 后外轮中心y坐标# 计算车辆转弯的轨迹
for t in range(0, 11): # 每隔0.1s计算一次# 计算车辆中心的位置x = turning_radius * math.sin(math.radians(turning_speed * t / 3600))y = turning_radius * (1 - math.cos(math.radians(turning_speed * t / 3600)))# 计算前内轮的位置x_front_inner = x + (length / 2 - front_rear_distance) * math.cos(math.radians(turning_speed * t / 3600))y_front_inner = y + (length / 2 - front_rear_distance) * math.sin(math.radians(turning_speed * t / 3600))# 计算前外轮的位置x_front_outer = x + (length / 2 + front_rear_distance) * math.cos(math.radians(turning_speed * t / 3600))y_front_outer = y + (length / 2 + front_rear_distance) * math.sin(math.radians(turning_speed * t / 3600))# 计算后内轮的位置x_rear_inner = x - (length / 2 - front_rear_distance) * math.cos(math.radians(turning_speed * t / 3600))y_rear_inner = y - (length / 2 - front_rear_distance) * math.sin(math.radians(turning_speed * t / 3600))# 计算后外轮的位置#见完整版

更多内容具体可以看看我的下方名片!里面包含有认证杯一手资料与分析!
另外在赛中,我们也会陪大家一起解析数学建模的一些方向
关注 CS数模 团队,数模不迷路~

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/660775.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

shell脚本,删除30天以前的日志,并将日志推送到nas,但运行出现/bin/bash^M。

删除30天以前的日志 将日志推送到nas中,然后删除pod中的日志 pod挂载到本地 运行出现/bin/bash^M 1、删除30天以前的日志: #! /bin/bash# 定义源日志目录 LOG_DIR/home/log/ # 删除日志 find $LOG_DIR -type f -name "*.log" -mtime 30 -exec…

LLM之RAG实战(三十八)| RAG分块策略之语义分块

在RAG应用中,分块是非常重要的一个环节,常见的分块方法有如下几种: Fixed size chunkingRecursive ChunkingDocument Specific ChunkingSemantic Chunking a)Fixed size chunking:这是最常见、最直接的分块方法。我们…

使用 LlamaIndex 和 Llama 2-Chat 构建知识驱动的对话应用程序

文章目录 使用 LlamaIndex 和 Llama 2-Chat 构建知识驱动的对话应用程序Llama 2-70B-聊天LlamaIndex 解决方案概述先决条件使用 SageMaker JumpStart 部署 GPT-J 嵌入模型使用 SageMaker Python SDK 进行部署在 SageMaker Studio 中使用 SageMaker JumpStart 进行部署使用 Sage…

Spring Boot Admin

概述 Spirng Boot Admin 登录页面 Spring Boot Admin是一个用于管理Spring Boot应用的监控工具,它允许你查看和管理多个Spring Boot应用实例。用于应用信息进行界面化的展示,常常辅助我们开发人员快速查看服务运行状态在微服务架构中,Spring Boot Admin通…

QT类之间主窗口子窗口传递*指针对象

1.新建CFile_Operation 类文件 2.主窗口头文件声明: CFile_Operation *cfile_operation; 按钮点击事件函数里面调用子窗口 dialog_debug new Dialog_Debug(this);connect(this,&MainWindow_oq::SendCfile_operation_Obj,dialog_debug,&Dialog_Debug::R…

【算法】【贪心算法】【leetcode】870. 优势洗牌

题目地址:https://leetcode.cn/problems/advantage-shuffle/description/ 题目描述: 给定两个长度相等的数组 nums1 和 nums2,nums1 相对于 nums2 的优势可以用满足 nums1[i] > nums2[i] 的索引 i 的数目来描述。 返回 nums1 的任意排列&…

【ZYNQ】Zynq 开发流程

Zynq 芯片架构由嵌入式处理器(Processing System, PS)与可编程逻辑(Programmable Logic, PL),以及 PS 与 PL 之间的互联总线组成。本文主要介绍 Xilinx Zynq 芯片开发所使用的软件,包括 Vivado IDE 与 Xili…

Django后台项目开发实战四

用户可以浏览工作列表以及工作详情 第四阶段 在 jobs 文件夹下创建 templates 文件夹&#xff0c;在里面创建 base.html 网页&#xff0c;内容如下 <!-- base.html --> <div style"text-align:center;"><h1 style "margin:auto; width:50%;&…

[软件工具]批量根据文件名查找PDF文件复制到指定的地方,如何批量查找文件复制,多个文件一起查找复制

多个文件目录下有多个PDF, 如何根据文件名一个清单&#xff0c;一次性查找多个PDF复制保存 如图所示下面有7个文件夹&#xff0c;每个文件夹里面有几百上千PDF文件 如何从上千个PDF文件中一次性快速找到我们要的文件呢 &#xff1f; 我们需要找到文件名是这样的PDF&#xff0…

《QDebug 2024年4月》

一、Qt Widgets 问题交流 1. 二、Qt Quick 问题交流 1.Rectangle 位置或大小有非整数&#xff0c;开启抗锯齿后看起来有点糊 测试代码&#xff1a; Rectangle {x: 100.3y: 100.4width: 100.5height: 30.5border.color: "#00A7AE"// 开启 antialiasing 后&#xf…

Microsoft Universal Print 与 SAP 集成教程

引言 从 SAP 环境打印是许多客户的要求。例如数据列表打印、批量打印或标签打印。此类生产和批量打印方案通常使用专用硬件、驱动程序和打印解决方案来解决。 Microsoft Universal Print 是一种基于云的打印解决方案&#xff0c;它允许组织以集中化的方式管理打印机和打印机驱…

ARM功耗管理背景及挑战

安全之安全(security)博客目录导读