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

news2024/9/24 13:19:20

当大家面临着复杂的数学建模问题时,你是否曾经感到茫然无措?作为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]) / 8
    return [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.1
    turning_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.coloradmin.cn/o/1644962.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

链表经典面试题02--链表的带环问题

目录 引言 环形链表 题目描述: 思路分析: 代码展示: 面试中遇到的问题: 环形链表Ⅱ 题目描述: 思路分析: 代码展示: 面试中遇到的问题: 方法二: 随机链表的复制 题目描述: 思路分析: 代码展示: 小结 引言 这个专题专门讲解链表的带环问题,并且对面试有关链表带环…

华中科技大学雷达站部署

一:项目地址 GitHub - HUSTLYRM/HUST_Radar_2023: 华中科技大学狼牙战队 RoboMaster 2023赛季 雷达站 二:安装依赖 2.1创建虚拟环境 首先是程序是基于python3.8完成,所以创建虚拟环境的时候,选择3.8的虚拟环境 conda create -…

【Linux】详解用户态和内核态内核中信号被处理的时机sigaction信号自定义处理方法

一、用户态和内核态的理解 在操作系统中,用户态和内核态是两种主要的执行模式,它们代表了不同的访问级别和权限,用于确保系统的安全和稳定性。 1.1、用户态 用户态是操作系统中用户进程的运行状态。在这种状态下,进程只能访问受限…

天锐绿盾 | 如何保护企业数据安全?

天锐绿盾为企业数据安全提供了一系列综合性的保护措施,确保数据在存储、处理、传输等各个环节的安全性。 PC地址: https://isite.baidu.com/site/wjz012xr/2eae091d-1b97-4276-90bc-6757c5dfedee 具体措施包括: 文件透明加密:采…

病毒防护的基础知识

目录 1.震网病毒Stuxnet是蠕虫病毒 2.宏病毒Macro virus 3.X卧底Phishing 4.钓鱼网站 5.病毒类型 1.震网病毒Stuxnet是蠕虫病毒 Stuxnet蠕虫病毒(超级工厂病毒)是世界上首个专门针对工业控制系统编写的破坏性病毒,能够利用对windows系统…

菜鸟跟着一起学如何批量修改文件日期时间属性 (创建、修改、访问时间)

需求:批量修改指定文件夹下的文件的创建时间、修改时间和访问时间 工具:【汇帮元数据编辑器】和【汇帮批量重命名】 文件一般有三个时间属性,分别是创建时间、修改时间和访问时间,如果你从别的地方拷贝的文件,创建时…

【手撸RPC框架】netty入门

🐼作者简介:一名大三在校生🎋 空有想法,没有实践,难成大事 专栏前言:探索RPC框架的奥秘 简介:在现代软件开发中,随着微服务架构的普及,远程过程调用(RPC&…

全网最全!场外个股期权和场内个股期权的区别的详细解析

场外个股期权和场内个股期权的区别 场外个股期权是指在沪深交易所之外交易的个股期权,其本质是一种金融衍生品,允许投资者在股票交易场所外以特定价格买进或卖出证券。场内个股期权是以单只股票作为标的资产的期权合约,其内在价值是基于标的…

编程入门(六)【Linux系统基础操作一】

读者大大们好呀!!!☀️☀️☀️ 🔥 欢迎来到我的博客 👀期待大大的关注哦❗️❗️❗️ 🚀欢迎收看我的主页文章➡️寻至善的主页 文章目录 🔥前言🚀Linux操作系统介绍与环境准备Linux操作系统介…

Java中的枚举类型介绍

一、背景及定义 情景: 枚举是在JDK1.5以后引入的。 主要用途是: 将一组常量组织起来,在这之前表示一组常量通常使用定义常量的方式: 这种定义方式实际上并不好。 例如:如果碰巧有另一个变量也是1,那么…

举个栗子!Alteryx 技巧(11):运用目录工具

你了解目录工具吗,它有什么作用?目录工具可以返回指定目录中所有文件的列表,即返回文件名和有关文件的其他信息,例如文件大小、创建日期、上次修改日期等。那么,怎样运用目录工具返回相关信息呢? 本期《举…

面对 800G 以太网设计的挑战

以太网是一种广泛使用的网络技术,用于连接局域网和广域网中的设备。它以 10 Mbps 的适度速度开始,多年来经历了大规模的创新。如今,以太网技术正朝着 800 Gbps 的速度发展,为数据中心和云计算基础设施等高需求环境提供超快的数据传…

Python-VBA函数之旅-pow函数

目录 一、pow函数的常见应用场景 二、pow函数使用注意事项 三、如何用好pow函数? 1、pow函数: 1-1、Python: 1-2、VBA: 2、推荐阅读: 个人主页:神奇夜光杯-CSDN博客 一、pow函数的常见应用场景 Py…

PyQt5中重要的概念:信号与槽

PyQt中信号与槽概念定义如下(网络上引用的): 信号(signal)和槽(slot)是Qt的核心机制,也是在PyQt编程中对象之间进行通信的机制。在创建事件循环之后,通过建立信号和槽的…

asp.net朱勇项目个人博客(3)

引文:按照书上的项目,我们最后实现管理端的三个增删改查的功能即可,相对与三个增删改查,文章,分类和留言,这里我们所需要用的的关联的一个表就是文章表,因为文章表每一个文章的增加显示和修改都需要对应的一个分类&…

OpenHarmony实战开发-如何实现动画帧

请求动画帧 请求动画帧时通过requestAnimationFrame函数逐帧回调,在调用该函数时传入一个回调函数。 runframe在调用requestAnimationFrame时传入带有timestamp参数的回调函数step,将step中的timestamp赋予起始的startTime。当timestamp与startTime的差…

分布式与一致性协议之ZAB协议(二)

ZAB协议 ZAB协议是如何实现操作地顺序性的? 如果用一句话解释ZAB协议到底是什么,我觉得它是能保证操作顺序性的、基于主备模式的原子广播协议。 接下来,还是以指令X、Y为例具体演示一下,帮助你更好地理解为什么ZAB协议能实现操作…

8086 汇编学习 Part 8

移位指令 当 C N T > 1 CNT > 1 CNT>1 时,CNT 必须是 CL 寄存器 逻辑左移 SHL OPR , CNT 将寄存器或内存单元中的数据向左移 CNT 位,最后移除的一位写入 CF,最低位用 0 补充 循环左移 ROL OPR , CNT 将寄存器中的值的最高位存…

山海鲸医疗科技:引领智慧医疗新潮流

随着科技的飞速发展,智慧医疗已经成为医疗行业创新的重要方向。在这个背景下,山海鲸智慧医疗解决方案应运而生,以其先进的技术和全面的服务,为医疗行业带来了前所未有的变革。 山海鲸智慧医疗解决方案是一套集成医疗信息化、大数…

三丰云:免费虚拟主机和云服务器的全面评测

引言: 在当前数字化时代,云计算已经成为企业和个人实现在线业务和数据存储的重要方式。然而,选择适合自己需求的云服务提供商并非易事。今天,我们将对三丰云的免费虚拟主机和云服务器进行全面评测,帮助您了解其优势和…