考虑空气阻力下的自行车挡泥板覆盖范围计算.docxVIP

PAGE \* MERGEFORMAT13 考虑空气阻力下的 自行车挡泥板覆盖范围计算 作者：尹达恒 所属院系：物联网工程学院物联网工程系 学号：1030616134  内容摘要 自行车前后挡泥板在自行车上起到阻挡车轮飞溅泥水打湿骑车人的作用。挡泥板的设计制造是自行车设计制造的重要一环。本文在考虑空气阻力的条件下，对不同车速以及飞溅水滴直径时自行车挡泥板最小覆盖范围进行了研究,主要内容如下： 使用常微分方程组建立了空气相对地面静止（无风），车轮相对地面不打滑滚动时的车轮-水滴飞溅模型。 利用上述模型计算出了在相同的水滴直径下两种不同种类自行车在正常车速情况下挡泥板的覆盖范围。 利用上述模型进一步讨论了相同的水滴直径下一般两种不同种类自行车挡泥板覆盖范围与车速之间的关系。 关键词 自行车挡泥板 数学模型 常微分方程组 最小值  目录 内容摘要 II 目录 III 一、绪论 1 二、建模 2 （一）基本假设 2 （二）模型分析 2 1、水滴受力分析 2 2、水滴初始状态分析 3 （三）建模 3 1、水滴运动方程 3 2、水滴初始状态方程 3 （四）模型验证 4 三、编写计算函数 5 四、计算 6 （一）竞速型山地自行车挡泥板位置计算(以GIANT? ATX850-M型为例) 6 1、数据估计 6 2、前轮计算 6 3、后轮计算 7 （二）竞速型山地自行车挡泥板位置与车速关系计算(以GIANT? ATX850-M型为例) 7 1、前轮 7 2、后轮 8 （三）竞速型山地自行车挡泥板位置计算(以GIANT? ATX850-M型为例) 9 1、数据估计 9 2、前轮计算 9 3、后轮计算 10 （四）竞速型山地自行车挡泥板位置与车速关系计算(以GIANT? ATX850-M型为例) 10 1、前轮 10 2、后轮 11 五、结论 12 六、参考文献 13考虑空气阻力下的自行车挡泥板覆盖范围计算 绪论 自行车在我们的日常生活中起着举足轻重的作用，尤其在大学校园，自行车可谓是师生出行的首选交通工具。但是由于骑车时自行车车轮与骑车人在同一平面上，所以在下雨天，车轮转动带起的泥水很容易飞溅到骑车人身上。为此人们为自行车设计了一许多不同形状不同大小的挡泥板，安装于前后车轮上方，起到阻挡泥水的作用。这就带来一个问题:一般情况下自行车挡泥板至少需要覆盖车轮多大的范围，才能正好挡住可能飞溅到骑车人身上的泥水?文献[1]使用maple大致计算了不考虑空气阻力的情况下自行车挡泥板的覆盖范围，但计算出的前轮挡泥板形状与现实生活中的前轮挡泥板有很大出入。 对于空中飞行的细小水滴而言，空气阻力会对水滴的飞行轨迹产生较大影响，忽略空气阻力的近似计算会使计算结果产生较大误差。针对这个问题，本文建立了考虑空气阻力的车轮-水滴飞溅模型并使用matlab计算了考虑空气阻力的情况下不同自行车挡泥板的覆盖范围，获得了与实际情况比较接近的结果。  建模 基本假设 为了讨论简单起见，现做如下假设： 假定自行车匀速直线运动。 假定人骑车时身体与自行车在同一平面内。 水滴脱离车轮时与车轮有相同的线速度。 空中飞行的水滴视为球体，直径3mm 地面水平。 空气相对地面静止（无风）。 车轮不打滑，车轮相对空气的平动速度等于其转速。 骑车人等效为一个矩形区域，落入此区域内的雨滴皆视为溅到骑车人。 覆盖范围以车轮挡泥板在车轮上方的覆盖角度θmin和θmax表示, 以车轮圆心为极点，水平向右为极轴，精度0.001rad，范围-π~π。 模型分析 水滴受力分析 以在空中飞行的水滴为原点，水平面为x轴，铅直线向上为y轴建立平面直角坐标系，如右图所示，对水滴进行受力分析 按照文献[2]给出的有空气阻力的情况下抛体的运动方程，可进一步得到有风时水滴在空中运动的微分方程组： dxdt=vx dydt=vy v=(vx-vwind)2+vy2 mdvxdt=-kv2cosα=-kv(vx-vwind) mdvydt=-kv2sinα-mg=-kvvy-mg k=12CρS (其中k为阻力系数，v表示水滴相对空气运动的速度，vwind表示风速，vx和vy分别表示水滴相对自行车的运动速度在水平和竖直方向上的分量) 水滴飞行初始状态分析 以车轮圆心为原点，平面为x轴，铅直线向上为y轴建立平面直角坐标系，如右图所示： 可知水滴飞行初始状态为 x0=Rcosθ y0=Rsinθ vx0=vcosπ2+θ=-vsinθ vy0=vsinπ2+θ=vcosθ 建模 水滴运动方程 (二)-1中给出的常微分方程组在matlab中可以转化为以下函数： function dot=drdt(~,r) m=(pi*r(6)^3)/6000000; %m=10^3*(4/3*