一种健身车骑行阻力模拟方法与流程

本发明属于健身设备
技术领域：
，特别涉及一种健身车骑行阻力模拟方法。
背景技术：
：申请号为cn02144446.3的专利文件，公开了“一种健身车与时实场景播放互动的方法，其特征在于有如下步骤：(a)由安装在健身车的脚踏板的转轴处的脉冲发生和电磁控装置，设置一个基本阻力，该阻力对应于模拟场景的一般路面行使，阻力设为v；(b)对应上述的阻力v，不同路段上程序预先设定不同的阻力值，进入爬坡路段阻力值会加大，大于v而在下坡的路段阻力值会减小，小于v或变为负值；(c)当一定力量施加于健身车使运动装置达到某一速度时，由脉冲发生和电磁控装置引发模拟信号产生，并采集该脚踏板转轴的转动速度的模拟信息；(d)该模拟信号通过i/o板，将该模拟信号转换为数字信号；(e)该数字信号通过健身车的接口单元以无线或有线方式发送给数据计算处理接口单元；(f)数据计算处理单元中：将场景影象的三维模型和控制三维模型转化为可视对象的程序存入非易失性存储器；临时存储器读取非易失性存储器中的各种数据，和微处理器一起运行，将三维模型转化成的一系列连续的可视图象复制到帧存贮器；帧存贮器将可视图象输出到接收显示单元上；(g)随着模拟场景的变换，出现不同路段，数据计算处理单元分析场景变化与之对应的阻力值，分析数字信号；(h)该数字信号通过接口单元传递给i/o板；(i)i/o板将数字信号变为模拟信号传递给脉冲发生和电磁控装置；(j)脉冲发生和电磁控装置产生相应的阻力，影响脚踏板的转轴的转动；(k)返回步骤(c)，继续采集运动模拟信号，循环上述步骤”。然而，该文件并没有公开健身车如何逼真模拟各模拟场景中的骑行阻力。技术实现要素：一种健身车骑行阻力模拟方法，将骑行阻力的输出设定分为三种模式，包括固定阻力模式、路线模式和比赛模式，固定阻力模式中，根据坡度和路面参数，设定多个阻力级别，通过阻力控制器输出骑行阻力，路线模式中，根据路线情况、模拟车辆信息、骑行人数据、踏频反馈数据、调速档综合计算骑行阻力输出，比赛模式在阻力控制上与路线模式一致，但设定多车比赛同时开始，初速度为0，必须等发令枪响后才可以骑行。所述的骑行阻力包括滚动阻力、坡度阻力、空气阻力和惯性阻力，滚动阻力ff为路面对轮胎的法向作用力为fy，车轮半径为r，滚动阻力矩为tf，δ为滚动阻力系数，轮胎与地面的法向作用力fz前移距离b，坡度阻力fθ＝mgsinθ，其中θ为爬坡角度，空气阻力cd为阻力常数，a为迎风面积，v为相对速度km/h，滚动阻力系数与轮胎类型、路面类型和速度的关系如下：轮胎类型路面类型速度0.015f1h200.016t1s300.020x2n40fy＝(m1+m2+m+m)gδcosθm1，m2为模拟自行车车型的前后轮的质量，m为车架的质量，m为骑行者的质量，g为重力常数θ为路面倾角，δ为滚动阻力常数。对于空气阻力，迎风面积a骑行者身高cm骑行者体重kg0.40170600.42175650.50175900.3517555相对速度为v＝vt+vd，其中，骑行速度vt，风速vd。坡度阻力fθ＝(m1+m2+m+m)gsinθ≈(m1+m2+m+m)gθm1，m2为模拟自行车车型的前、后轮的质量，m为车架的质量，m为骑行者的质量，g为重力常数θ为路面倾角。惯性阻力其中m1，m2为模拟自行车车型的前、后轮的质量，m为车架的质量，m为骑行者的质量，dv/dt为加速度，i1，i2为前后轮的转动惯量。经过测试，本发明的健身车骑行阻力模拟方法能够较为精确地模拟输出健身车骑行阻力输出。附图说明图1是本发明实施例中滚动阻力分析原理图。图2是本发明实施例中自行车动力学模型原理图。图3是本发明实施例中骑行驱动力分析原理图。图4是本发明实施例中健身车力学模型原理图。图5是本发明实施例中坡度阻力计算原理图。图6是本发明实施例中路线模式示意图。图7是本发明实施例中路线模式计算示意图。具体实施方式在骑行过程中，人动力和骑行阻力的共同作用决定了自行车的运动的过程。骑行阻力决定了运动者体能的消耗，骑速的大小，它是健身车模拟现实路况中最主要的因素。模拟自行车的路面骑行，关键是最大限度的模拟路面骑行环境条件，比如温度，湿度，大气压力，噪音，行驶阻力等等。这其中最重要的是骑行阻力，这是健身车运动力学分析的基础。骑行阻力一般包括滚动阻力、重力坡面阻力、空气阻力以及惯性阻力。车轮是自行车的运动部件，车轮滚动时，轮胎与地面的接触会有法向，切向的相互作用力，而且路面和支撑地面会有应变，在硬的路面运动时，轮胎变形产生的弹性迟滞，形成阻碍车轮前行的阻力矩，是法向作用力fz前移距离b，如图1所示。图中，路面对轮胎的法向作用力为fy，车轮半径为r，滚动阻力矩为tf，则滚动阻力ff，上述等式中，δ为滚动阻力系数，该系数与路面种类骑行速度，轮胎的构造，材料等因素有关。低速度时，滚动阻力系数随车速变化不大。滚动阻力在受力图上无法表现出来。如果地面切向作用力fx小于轮轴的推力fp，则会产生打滑现象。fx的极限值称为附着力fψ，在硬的路面上，它与车轮的法向反作用力fy成正比，fψ＝ψfy，比例系数ψ称为附着系数，由路面和轮胎决定，通常比滚阻系数δ要大很多。坡度是以坡高于底长比值表示的，根据目前的公路工程技术标准，平原丘陵地区i级路面最大坡度为4％，山岭区域路面坡度最大9％。坡度阻力fθ＝mgsinθ，其中θ为爬坡角度(下坡时为负值)，弧度单位。一般路面坡度较小的情况下sinθ＝tanθ＝θ所以fθ＝mgsinθ＝mgtanθ＝mgθ空气阻力又包括压力阻力和摩擦阻力两部分。空气阻力可以描述为：cd为阻力常数，a为迎风面积(平方米)，v为相对速度km/h。可以直接简化为作用于车子质心处的作用力。惯性阻力：惯性阻力包括车身总质量平移的惯性，也包含车轮转动的惯性，一般将车轮转动的惯性简化为平移的惯性。本发明建立的自行车的动力学模型如图2所示，自行车车体及骑行者的重力(m+m)g，车体重心为c，距离前后车轮轴心c1，c2分别为l1，l2，从动轮、驱动轮的质量为m1，m2，转动惯量为i1，i2，骑行加速度为dv/dt，相应的车轮角加速度为dw/dt，行驶中空气阻力fw，坡度阻力为fθ。w1，w2为从动轮和驱动轮上的载荷，fy1，fy2，为路面对从动轮、驱动轮的法向作用力，fp1，fp2为从动轮驱动轮轴作用于车轮平行于路面的力，tf1，tf2为从动驱动轮滚动阻力矩，fx1，fx2为作用在从动轮驱动轮上的地面切向反作用力，tm1，tm2为从动轮驱动轮转动惯性力矩，t为人力驱动子驱动轮上的转矩。从动轮在车轮轴传来的推力fp1，地面切向反作用力fx1，滚动阻力矩tf1，载荷w1，以及重力m1g的共同作用下，由平衡条件可得：所以：其中驱动轮在车轮轴传来的推力fp2，地面切向反作用力fx2、滚动阻力矩tf2、驱动力矩t以及载荷w2,，重力m2g的共同作用下由平衡条件得到：人作用于驱动轮的驱动力f＝t/r，根据以上可得：其中对车身进行分析，在从动轮、驱动轮对其反作用力f’p1，f’p2，支撑力w’1，w’2，以及重力(m+m)g，空气阻力fw的共同作用下：由于f’p1，f’p2与fp1，fp2为作用力与反作用力，fy＝fy1+fy2＝(m1+m2+m+m)gδcosθfθ＝(m1+m2+m+m)gsinθ≈(m1+m2+m+m)gθ也就是驱动力等于滚动阻力、空气阻力、坡度阻力与惯性阻力之和。骑行的驱动力如图3所示，踏力为f，链条拉力为f1＝f2,，阻力器张力为fm，形成推力为fp。车轮的半径为r。踏柄的长度为l，齿数为n，飞轮齿数为n。齿轮齿数与半径的比例系数为k。踏力力矩为tf＝fltf＝f1nk，车轮的扭力矩tl＝f2nk，f1＝f2。得到：所以也就是这个推力要等于前面提到的阻力之和。自行车飞轮的转速采用霍尔传感器测量，飞轮每转动一圈，发生三个脉冲。设飞轮的转动频率为f0，轮毂半径为r(m)，得到的车速为：v＝2πrf0m/s，换算成公里小时v＝7.2πrfkm/h。示例：轮毂尺寸，轮胎型号决定骑行轮胎转动一周的距离26*1.5寸的是2.026米26*1.6寸的是2.051米26*1.75寸的是2.070米26*1.9寸的是2.089米26*2.0寸的是2.114米26*2.125寸的是2.133米齿数比为52/16，踏频为1r/s。则行车速度为3.6*2.026*1*52/16＝23.7km/h(参考1.5寸轮胎计算)。健身车的力学模型如图4所示，踏力为f，链条拉力为f1＝f2,，阻力器张力为fm，力矩为tm。轮子的转动惯量为i。踏柄的长度为l，齿数为n，飞轮齿数为n，齿轮齿数与半径的比例系数为k，飞轮的半径为r0，踏力力矩为tf＝fl，惯性力矩为tf＝f1nk，飞轮的扭力矩tl＝f2nk，f1＝f2。得到：fp为坡度阻力，空气阻力，惯性阻力，滚动阻力之和。m为健身车飞轮的质量。r为一般自行车的轮毂半径，r0为健身车飞轮的半径。其中：fy＝fy1+fy2＝(m1+m2+m+m)gδcosθfθ＝(m1+m2+m+m)gsinθ≈(m1+m2+m+m)gθ从上面的力学模型分析可以看出，健身车需要模拟出滚动阻力，坡度阻力，空气阻力和惯性阻力四个部分。滚动阻力与滚动系数δ有关，而阻力滚动系数与轮胎种类，路面关系最为密切。同时和骑行速度也有关。因此需要构建滚动系数与轮胎准类、路面、速度与之间的关系。车型：计算数据需要：重量，变速齿轮比例，轮毂尺寸。从展现角度可能还需要：外形图片，3d，制造商，型号等等(需要仔细考虑)。路面：路面类型，土路，柏油，水泥，级别等等。建立滚动阻力系数与这些参数之间的关系。δ轮胎类型路面类型速度0.015f1h200.016t1s300.020x2n40fy＝fy1+fy2＝(m1+m2+m+m)gδcosθ滚动阻力的计算公式如上。m1，m2为模拟自行车车型的前后轮的质量，m为车架的质量，m为骑行者的质量，g为重力常数θ为路面倾角，δ为上表中的滚动阻力常数。滚动阻力系数取值范围参考表上表在低速下测得的数据。空气阻力的计算如下面的公式：cd为阻力常数，a为迎风面积(平方米)，v为相对速度km/h。空气阻力系数cd与迎风面的光滑程度，材料，造型等等有关。我们这里设置为固定值。专用的骑行服就是为了减小风阻系数。一般风阻系数在0.5到0.7左右。迎风面积a的大小与骑行姿势，车辆，骑行者的身高体重有较大关系。这里可以模拟建立体重、身高跟面积的关系。(数值不准确，仅作为示例)迎风面积a身高cm体重kg0.40170600.42175650.50175900.3517555相对速度为骑行速度vt与风速vd的和。顺风vd为负值。v＝vt+vd坡度阻力fθ＝(m1+m2+m+m)gsinθ≈(m1+m2+m+m)gθ滚动阻力的计算公式如上。m1，m2为模拟自行车车型的前后轮的质量，m为车架的质量，m为骑行者的质量，g为重力常数θ为路面倾角。下坡时候θ为负值。计算如下所示：惯性阻力：滚动阻力的计算公式如上。m1，m2为模拟自行车车型的前后轮的质量，m为车架的质量，m为骑行者的质量，dv/dt为加速度。i1，i2为前后轮的转动惯量。因为车轮的质量分布主要集中于轮毂。所以：i1≈m1r2i2≈m2r2εm＝2m1+2m2+m+m本发明的阻力输出数据模型包括三种模式：a.固定阻力训练模式b.路线模式c.比赛模式固定阻力模式，可以通过控制台控制阻力级别，比如陡坡，平坡，平地。也可以考虑路面类型，比如柏油路，水泥路，砂石路等等。实际健身车可以模拟256级，当然控制台只需要几个级别即可。后续可以根据市场需求情况进行优化调整。也就是控制台只需要向阻力控制器发送1-256之间的一个数字即可。计算出以上所有的阻力之和，对照0-256级的阻力表，得到阻力级别发送到阻力控制器。路线模式下，通过前面的分解，可以看到，需要根据路线情况、模拟车辆信息、骑行人数据、踏频反馈数据、调速档综合计算出电磁阻力矩tm的大小。比赛模式在阻力控制上与路线模式一致。但需要考虑比赛同时开始，初速度为0，必须等发令枪响后才可以。骑行者能量消耗及上述的所有阻力的和与骑行距离的乘积的积分单位为焦耳，l为骑行距离，f为骑行阻力。使用求和：fn为当前阻力和，an为当次行进距离。单位焦耳，可以换算成大卡，卡路里等人们更熟悉的单位。因为实际骑行的惯性比较大，而模拟骑行过程中惯性较小，会降低用户的部分体验。需要通过技术手段改善。基本技术原理即：改变骑行者停止踩踏和开始踩踏时候，降低阻力，这样更加接近与实际骑行过程。当前第1页12