韩连强 陈学超 ,2 余张国 ,2 高志发 黄 岩 ,2 黄 强 ,2

双足机器人在复杂环境中的运动适应能力是其走向应用的重要指标之一.近年来双足机器人的动态运动得到大量研究,其在复杂环境中的适应能力有了很大提升[1-2].但传统全驱动位置控制机器人由于动态性限制,很难在崎岖地形中灵活运动.而欠驱动双足机器人由于本身具有的动态特性,表现出很强的地形适应性[3].这类机器人的脚踝具有1 个或0 个驱动关节,所以脚板通常会被点足替代,因此欠驱动双足机器人像人类高跷运动一样需要交替迈动双脚完成行走功能.由于灵活的运动能力通常需要连续地形环境,相比能够实现准确落脚位置[4]的全驱动位置控制双足机器人[5-7],采用关节力矩控制方法[8]的欠驱动双足机器人缺少精确的落脚点控制以保证自身平稳通过随机离散地形环境.因此使欠驱动双足机器人具备精确落脚控制能力可以扩大其应用场景,活动范围覆盖全地形环境,进而完成任务式运动需求.

欠驱动双足机器人的动态行走研究近几年得到巨大突破,但行走环境被限定为连续地形.Kim 等[9]提出的基于落脚点调节的全身运动优化控制,使机器人成功在室内进行三维行走;Luo 等[10]提出一种三维欠驱动双足机器人整体操作空间控制框架,观测质心状态更新控制周期与步长位置,实现了鲁棒平衡;Daneshmand 等[11]提出了可变尺度模型预测控制 (Model predictive control,MPC) 框架,通过实时修正落脚点实现室内三维运动;为了获得更鲁棒的控制效果,早期Westervelt 等[12]提出了混合零动力 (Hybrid zero dynamic,HZD) 的控制方法,离线构建多个步态库,通过在线调节使二维机器人在平地上进行行走,随后该方法被推广至三维运动控制中,已经实现连续起伏地形下的三维行走[13-14].而Matthew 等[15]在HZD 基础上使用基于快速指数稳定控制Lyapunov 方程的模型预测控制方法,获得了机器人二维动态行走能力.Guo 等[16]提出基于质心模型预测的步态合成规划方法,能够在线实时控制三维机器人稳定行走.Gong 等[3]在此基础上增加了实时摆动腿角度调节,使得机器人能够在草地、雪地等复杂环境行走.虽然上述方法的控制形式不同,但都并未将落脚位置作为最终控制目标,本质都是通过落脚点位置控制机器人的运动速度,因此无法实现期望步长的控制.并且,上述研究方法的测试场景都是连续不平整地面或起伏的草地,当环境出现如石头等离散立足点时,这些控制方法将出现很大局限性.

欠驱动双足机器人在随机离散地形中的运动控制在上述问题中被深入研究,主要思想是步长的实时修正.Negri 等[17]应用神经网络训练出不同步长可调节的被动步态参考生成器,使用非线性MPC跟踪参考轨迹,但只应用在被动步行中.姚渊等[18]提出了基于自适应前馈控制算法的变步长稳定运动控制策略,通过改变步长与质心参考跟踪速度,成功实现了机器人的非连续地面行走.Yao 等[19]针对人体变步长步态特征,提出基于质心状态的前馈控制策略,将变高度等效于坡度变化控制,实现了连续台阶行走仿真.由于欠驱动双足机器人与地面接触点上自由度的不可控原因,机器人在摆动中的状态不能准确预测与控制[20],并且这种随时间变化的步行轨迹很容易受外界扰动出现提前触地的情况,造成自身的不稳定.

为了解决这一问题,Grizzle 等[21]提出通过虚拟约束 (Virtual constraint,VC) 建立机器人状态与控制目标变量间关系的方法.这种方法虽然在动力学特性上与机械约束不同,但运动学特性能够满足机器人步态的要求.在此基础上,Yang 等[22]提出在全身动力学模型上设计多个动态周期步态并预先获得控制器,实现了步态库每一个元素的精准切换,尽管这种方法可实现不同步长步行,但控制器数量呈指数增长.Nguyen 等[23]建立了长度和高度地形离散信息,在HZD 基础上通过离线轨迹优化获得不同的步态库,在已知下一步离散落脚点后插值控制参数得到所需的步态.但当机器人自由度增多时,非线性优化求解困难,并且对未构建的地形适应能力减弱[24].

为避免离线优化问题并实现在线实时调节,以增强欠驱动双足机器人离散地形的鲁棒性,在前期动态运动控制的工作基础上[25-26],本文提出基于虚拟约束的变步长调节方法.其创新点如下:

1) 提出VC 参数化的尺度缩放因子实时调节方法,不仅将机器人实时状态与参考轨迹建立关系,并构建行走中的步态轨迹,还能在不同步长与高度需求间在线任意衔接;

2) 为实现步态轨迹的精确跟踪,使用反馈线性化的MPC 控制实现期望的规划步态跟随控制;

3) 采用平面欠驱动双足机器人模型进行了多种随机离散地形的稳定行走控制和算法的仿真验证.

本文结构内容安排如下: 首先建立了欠驱动双足机器人动力学模型,然后详细描述了基于非时变尺度缩放因子的步态规划设计方法,之后阐述了基于反馈线性化的MPC 轨迹跟踪控制算法,随后展示了在不同场景的离散地形下实时稳定运动的仿真结果,最后是总结与展望部分.

1 欠驱动双足机器人动力学模型

基于欠驱动双足机器人行走运动的周期性和对称性,可以将运动分为单脚支撑期和双脚碰撞冲击期两个阶段.忽略双脚支撑状态,摆动脚瞬间触地碰撞后成为下一阶段的支撑腿,而原本的支撑腿变为摆动腿,进而形成步态循环.考虑实际机器人结构,使用多连杆刚体动力学简化模型构建该过程.如图1 所示,机器人模型拥有5 个连杆,包括1 个躯干 tor、2 个大腿 fem 和2 个小腿 tib,以及2 个髋关节 hip 和2 个膝关节 knee.各连杆动力学参数有长度Lυ、质量mυ、相对关节的质心位置lυ和转动惯量Iυ,其中υ∈[tor,fem,tib],默认左右腿参数相同.

1.1 单脚支撑模型

在机器人单脚支撑期时,如图1 左图所示,脚与地面的接触认为是铰接触,机器人无滑动,即模型为固定基座的运动链.因此以接触点为原点建立x-z平面坐标系,θ1表示支撑腿小腿与坐标轴z方向的夹角,而θ2,θ3,θ4,θ5依次是下一连杆相对上一连杆的夹角,如图1 中图所示,使用虚线表示上一连杆,规定顺时针为正方向.

图1 机器人动力学模型与运动阶段Fig.1 Dynamic model and motion stage of robot

通过牛顿-欧拉方程可得到单脚支撑模型的动力学方程:

1.2 碰撞模型

欠驱动双足机器人碰撞模型被描述为摆动腿与地面的瞬间碰撞,进而导致速度产生突变而位置不变的映射模型.该模型已被大量使用在机器人运动控制中,详细的理论推导过程请参考文献[12].设是碰撞后的系统状态,S是碰撞发生时的状态判定域.状态判定域可以包括摆动腿在支撑腿前、摆动腿末端高度为零等.则当x处于S时,满足离散更新的条件:

Δ:R10×1→R10×1表示碰撞映射关系,即将碰撞前后的机器人状态进行关联.因此,平面欠驱动双足机器人可以被一个由连续时间动态方程式(3)和一个离散更新方程式(4)组成的混合非线性冲击系统完整描述.

2 基于虚拟约束的尺度缩放因子设计

VC 能够建立多个变量间的相互关系,因此可以使双足机器人步态不依赖于时间变化,而只依赖自身的运动学状态.通过尺度因子调节变量间比例关系即可快速生成不同环境需求下的步态轨迹并在运行中随时修正.基于VC 的尺度缩放可以简化运动学的规划,但会间接设计出一组零动力学.本文并不通过设计的零动力学产生周期轨道运动,但为了使拥有混杂动力学的系统轨迹相容,通过切换后的平滑处理来减小跳变的发生.

2.1 行走步态的虚拟约束设计

在欠驱动双足机器人行走运动中,上身姿态、摆动腿轨迹和支撑腿长度的设计规划非常重要[9-10,12,14-16],它们共同组成了运动步态.如图2 所示,由于平面欠驱动双足机器人只有4 个驱动自由度,因此行走步态中规划4 个控制目标即可.为了避免计算腿部逆运动学以及由于欠驱动自由度导致的不可控量,本文分别选择上身与竖直方向的夹角αtor、支撑腿髋关节的竖直高度zhip、摆动腿末端 foot 的竖直高度zfoot和水平位置xfoot进行规划设计.而为构建VC关系,需要选择机器人运动过程中单调递增的状态作为相变量δ.如图2 选择机器人髋关节与支撑点的水平距离xhip,即δ=xhip,这是因为机器人在向前或向后运动中xhip能够单调递增或递减.以上各位置均可由正运动学计算得到.接下来,首先规划平地行走时虚拟约束设计,然后在此基础上表述变步长与高度的尺度缩放因子的设计.

图2 机器人步态中的虚拟约束设计Fig.2 Virtual constraint design in robot gait

为稳定欠驱动双足机器人行走运动,通常控制上身姿态αtor为恒定的目标角度,而通过模型可知αtor=θ1+θ2+θ3与相变量δ无关,因此只需要设定控制目标:

机器人在平地行走时,稳定运动的基本条件就是周期性的步态,在实际中的体现为碰撞后的运动状态是下一运动阶段的开始.为了步态的周期对称性,摆动腿和支撑腿的相关轨迹采用一元二次方程的约束规划方式,具体操作如下: 设定固定的单步步长为Lstep(在碰撞阶段两脚的水平距离),为了使机器人在设定的步长上运动,需要满足摆动腿末端在xfoot=Lstep时zfoot=0 以保证双足碰撞阶段的发生,进而切换到下一单脚支撑阶段.同时设髋关节时刻保持在两腿末端中间以稳定质心状态,因此摆动腿末端关于髋关节与支撑点应满足对称条件:

由于上述约束,在平地行走时相变量将被限制在一个取值范围-0.5Lstep≤δ ≤0.5Lstep内.这也符合目标落脚点的要求,即机器人不应该超出目标步长,从而进一步保证机器人不会踩空.

2.2 变步长与变高度的尺度缩放设计

当下一步落脚点水平距离不等于固定步长Lstep(在碰撞阶段两脚的水平距离)时,需要改变式(8a)和式(9a),并考虑上一步和下一步步长信息,以获得平滑的虚拟约束轨迹.在上述平地行走的VC 参数设计之上,通过设计尺度缩放因子就可以适应不同步长与高度的地形行走.假设已知当前一步的参考落脚点,如图3 所示,设定上一步步长为和相对地面高度为,而当前目标步长为及相对地面高度为,这些量都是相对支撑点.需要说明的是,由于碰撞冲击阶段的转换及以支撑点为坐标原点的原因,在碰撞发生后,当前相对地面高度变成上一步高度时满足取反变换,即

图3 随机离散地面的参数示意图Fig.3 Parameter diagrams of random discrete ground

欠驱动双足机器人能够准确落脚在规划的目标步长和高度的本质原理就是在机器人摆动腿末端到达期望目标时与环境发生碰撞,从而切换进入下一单脚支撑期.由上文行走步态的VC 设计可知,摆动腿末端有2 个控制量,分别是竖直高度zfoot和水平位置xfoot.其中xfoot被实时控制在2 倍于相变量δ的目标上,因此可以通过缩放当前步长,控制竖直高度zfoot在目标步长时等于0,进而达到控制步长的目的.

因此设定改变步长的尺度缩放因子γ=Lstep/,可将式(8a)修改为:

由于机器人的机械响应限制,作为参考轨迹需要连续且没有较大跳变,但由于每一步的目标步长是离散的,因此需要对步长尺度缩放因子在发生碰撞后进行平滑过渡处理:

即当相变量大于等于0 后,调节目标步长固定不变.需要强调的是,式(11)中的平滑处理仍然只取决于相变量而不依赖时间的变化.可以想象当机器人向后倒回时,摆动腿仍然能够回退到起脚时的着地点,而不会踩到其他地方,从而进一步增强自身的稳定性.

同样,在平地变步长的基础上增加落脚高度调节时,核心思想是在达到期望步长时到达期望高度,因此调节着地时的落脚高度如下:

高度尺度缩放因子λ通过平滑函数χ(δ)和κ(δ)在上一步与下一步高度间过渡.当相变量为0 时摆动腿保持原来的抬脚高度是为了保证机器人的平稳运动,而摆动腿不会踢碰到障碍物.

由图3 所示,由于落脚高度的改变且机器人坐标原点在支撑点上的坐标系变化,因此两个支撑期间切换时虚拟约束中的支撑腿髋关节高度会产生跳变.通过调整避免上述跳变发生,设发生碰撞时髋关节相对碰撞点的位置为则下一支撑期起始时刻支撑腿髋关节参考高度应为由式(9a)可计算下一阶段VC 的髋关节高度参数:

3 基于反馈线性化的MPC 轨迹跟踪控制

在上述轨迹规划基础上,需要通过控制驱动关节的力矩使机器人达到轨迹跟踪的目的,从而控制机器人完成行走运动.在单腿支撑期时,机器人如同倒立摆一样处于不稳定状态,从而相对支撑点发生转动.这时摆动腿到达下一步落脚点,进而切换支撑腿防止机器人摔倒,重复上述过程即实现机器人的动态稳定.因此定义欠驱动双足机器人的稳定性(平衡性)为: 顺利通过离散地形而不倾倒.本文方法的算法框图如图4 所示,主要包括两部分: 当已知未来一步地面信息后,通过VC 参数化单脚支撑期轨迹并设计尺度缩放因子获得参考轨迹;然后通过状态反馈线性化的MPC 滚动优化控制实现轨迹跟踪,进而控制机器人在离散地形中的行走.接下来描述轨迹实现使用的控制方法.

图4 控制算法框图Fig.4 Control algorithm block diagrams

3.1 状态反馈线性化

为控制机器人实现第2 节设定的VC 参考轨迹,设定单脚支撑期系统输出:

对于非线性系统的输出y=y(x),由于只取决于状态中的位置变量,则其相对度为2,输出的二阶导数为:

但由于模型与实际机器人的误差,控制中需要增加一个微小的误差调节量η(y,):

3.2 MPC 滚动优化调节控制

如反馈线性化控制式(18)中所示,η(y,)表示系统输出及其导数的函数,通常被选为比例-微分(Proportional-derivative,PD) 控制器的输出.但为了能够优化能耗,使用MPC 的滚动优化特点获得这个误差调节量.为了获得能够用于MPC 控制的线性系统,将式(18)代入式(16)中,得到经过反馈线性化后的二阶系统微分方程:

设定控制周期为 ΔT,式(20)经过离散化可得离散空间状态形式:

下角标表示第k步与第k+1 步变量.设定MPC 的预测步长为NP,而控制量步长相同,通过迭代关系获得完整NP步预测模型:

由于原系统输出式(15)全部为0 时即系统达到VC 要求,所以对于预测模型式(22)的滚动优化目标即状态量趋近于0.控制量力矩为了减少能耗所以应尽量小,但状态量相对放松限制,设定优化目标函数为:

其中,W与R∈R(8×NP)×(8×NP)为权重系数的对角矩阵.将式(22)代入式(23)然后去除无关量后,得到标准的二次规划 (Quadratic programming,QP)形式:

因此建立了具有约束的QP 优化为:

其中,τmin,τmax表示误差调节量η(y,)的边界值.设实际机器人驱动关节的输出力矩边界为umin,umax,通过式(18)可知在当前状态下,物理约束的驱动力矩限制与误差调节量的边界有如下关系:

最后取优化结果U的第一组值,即是误差控制量η.

4 仿真结果

为验证本文提出的基于尺度缩放离散地形双足机器人动态行走控制算法的有效性,基于双足机器人BHR-6S 动力学模型,通过使用数学工具常微分方程 (Ordinary differential equation,ODE)积分算法构建仿真环境进行数值仿真.使用的BHR-6S机器人模型参数如图5 所示:

图5 仿真机器人模型参数Fig.5 Model parameters of simulation robot

需要说明的是模型各连杆的质心分布在几何中心处.其他控制参数如表1.虚拟约束参数中上身姿态角度限制在0°,最大抬腿高度不超过0.1 m,跨步步长不超过0.5 m,髋关节期望高度在0.6 m～0.7 m 之间.其中状态和控制权重系数给出了矩阵对角线的一个值,其他值是相同的.驱动力矩边界值同样设定为与4 个驱动关节边界值相同.

表1 仿真参数设置Table 1 Parameters setting in simulation

下面通过平地定步长(第4.1 节)、随机变步长(第4.2 节)、随机离散变步长与变高度(第4.3 节)和随机扰动(第4.4 节)的实验数据分析验证算法的有效性和鲁棒性.这里需要强调的是机器人仅已知当前一步的环境落脚点位置,对下一步环境信息是未知的,因此随机地形的变化考验着机器人实时控制策略的有效性.

4.1 机器人等步长的平地运动

首先进行了固定步长0.4 m 的平地仿真.如图6所示为仿真的动画截图.

图6 机器人平地固定步长行走仿真截图Fig.6 Simulation screenshot of robot walking with fixed step length on flat ground

如图7 所示为行走步态的极限环,横坐标是相变量,纵坐标为髋关节的前进速度.由于碰撞的跳变速度会偏离稳定点,但基于反馈线性化的MPC控制很快将其稳定在极限环轨道上,说明了本文VC设计与控制的有效性.

图7 机器人运动状态的极限环Fig.7 Limit cycle of robot motion state

行走中的单步长周期内关节驱动力矩如图8 所示,整个运动过程力矩保持在机械能力限制内,并经过快速稳定的优化控制过程后力矩实现平稳变化,进一步保证了机器人周期状态的平稳调节.

图8 机器人单脚支撑期的关节驱动力矩Fig.8 Joint driving torque of robot during one leg support period

图9 显示了单步运行过程中的VC 参考轨迹以及机器人实际状态的跟踪情况,控制目标跟踪性很好,从而保证了机器人能够准确下踩在期望的落脚点上,即实现参考步长.

图9 机器人运动过程中VC 轨迹的跟踪效果Fig.9 Tracking effect of VC trajectory in robot motion

4.2 机器人随机变步长的平地运动

图10 显示了机器人在离散变步长地形中的仿真情况,为表现机器人成功通过离散地形给出了全场景截图,但对于机器人来说仅知晓当前一步落脚点位置,控制器根据目标步长与高度实时调整轨迹步态.最终共进行了20 步的行走,随机出现的步长分别为: 0.40 m、0.20 m、0.35 m、0.45 m、0.30 m、0.20 m、0.30 m、0.45 m、0.50 m、0.20 m、0.30 m、0.20 m、0.35 m、0.20 m、0.20 m、0.30 m、0.40 m、0.30 m、0.15 m、0.30 m.

图10 机器人变步长离散地形行走仿真截图Fig.10 Simulation screenshot of robot walking on discrete terrain with variable step-size

步长尺度缩放因子与VC 参考轨迹的变化如图11 所示,实线表示摆动脚x方向的参考值.在单脚支撑期时,整个模型的原点在支撑点上,而摆动脚通常需要从后摆动到前面,所以曲线从负值连续变到正值.在当前一步完成后,经碰撞模型映射后,再次进入单脚支撑期,摆动腿参考值从正值突变到负值,所以突变表示支撑腿切换的过程.由于仅建立了单脚支撑期模型而不区分左右腿,所以之后的图像都是综合了两条腿的结果.随着步长的需求变化,γ在支撑期的前半程平滑过渡.参考步长变小时,γ变大,摆动腿水平参考位置变小;当γ变小时表示参考步长变大,γ减小到1 附近,当小于1 时表示参考步长超过Lstep,如上述参考步长中的0.45 m.这证明步长尺度缩放因子能够有效地实时控制机器人的步长变化.

图11 离散平地变步长运动中的尺度因子与VC 轨迹变化Fig.11 Scale factor and VC trajectory variation in discrete flat variable step-size motion

而图12 则展示了机器人实际步长与参考步长的误差,柱状图顶部的数字为参考步长,步长控制效果非常好,也为真实实验提供了有力理论支撑.

图12 离散平地变步长行走的步长误差Fig.12 Step error of discrete variable step-size walking on flat ground

4.3 机器人随机变步长与变高度的离散地形运动

图13 展示了机器人在随机离散地面上的行走情况,该地形是在变步长仿真的基础上更改了地面高度,高度随机分布在[-0.1 m,0.1 m]区间中.图14表示了机器人实际落脚点与参考落脚点的空间分布,实际落脚点与参考落脚点几乎重合,而精确的落脚控制保证了机器人顺利通过复杂离散地形.

图13 随机离散变步长变高度行走的仿真截图Fig.13 Simulation screenshot of random discrete variable step-size and variable height walking

图14 空间中参考落脚点与实际落脚点的示意图Fig.14 Schematic diagram of reference foothold and actual foothold in space

图15 中展示了机器人通过随机离散地形时的高度尺度缩放因子和VC 参考轨迹的变化情况.可以看出当参考高度改变时,λ先从上一阶段参考高度相对值回到0,保证机器人的抬腿高度能够避开障碍,然后被平滑至参考高度,从而控制机器人落脚在参考点.碰撞发生后为避免支撑腿参考高度与实际高度差距大而发生跳变,参数被重新修改,因此参考的支撑腿髋关节高度和摆动腿末端竖直高度发生跳变,以适应新阶段状态控制,进而保证了下一步的平稳运行.随机离散地形的实验结果验证了高度实时调节中尺度缩放因子的有效性.

图15 随机离散地形运行中的尺度缩放因子和VC 参考轨迹Fig.15 Scaling factor and VC reference trajectory in random discrete terrain

4.4 机器人一步内扰动的稳定表现

最后进行了机器人在随机离散地形中行走时受扰动的仿真验证(如图16),测试VC 轨迹在增加了尺度缩放因子后的非时变稳定性.

图16 机器人行走中受扰动的仿真截图Fig.16 Simulation screenshot of robot disturbed during walking

图17 显示了机器人通过随机离散地形时的高度尺度缩放因子和VC 参考轨迹的变化情况.当参考高度改变时,λ先被相变量改变到0,保证机器人直立时的原有抬腿高度避开障碍,按照变尺度缩放因子设计之后将再次平滑至参考高度,从而控制机器人落脚在参考点.但在行进中受到作用在髋关节位置上的500 N 水平反方向扰动,作用时间0.15 s后,相变量又再次越过0 回到负值,即机器人参考轨迹回退到阶段初始状态.这验证了依赖自身状态的非时变尺度缩放因子既可以调整机器人到目标参考点,又可以保障机器人受到大扰动后稳定落脚点变化使机器人不会 “踩空”导致倾倒,从而回退到上一步支点.本文将这种不随时间变化的平衡调节能力定义为非时变稳定性.因此本文提出的实时调节方法将机器人约束在两个已知可行支撑点上,面对更复杂的崎岖地形都能实现稳定行走.

图17 机器人行走受扰动后的尺度因子变化情况Fig.17 The change of scale factor after the robot is disturbed in walking

图18 给出了在有无外力扰动情况下,当前一步内上身躯干姿态角的变化情况.实线表示躯干角度在无扰动下的变化,暗示着在即将落脚时上身前倾的表现.虚线表示躯干角度在有扰动下的变化,当受到扰动而后退时角度偏离正常变化,但控制器仍然将其控制在期望角度附近而稳定步行.

图18 机器人行走中有无扰动下躯干姿态变化Fig.18 The posture change of the trunk with or without disturbance during robot walking

5 结论

为使欠驱动双足机器人能准确迈出特定步长及高度且稳定运动,从而通过随机离散地形,本文提出了基于实时尺度伸缩的VC 轨迹规划方法,并使用了反馈线性化的MPC 控制机器人实现设定的约束步态,最后在仿真实验中验证了算法的有效性.非时变的VC 轨迹只依赖机器人自身状态变化,这使得抗扰动能力更强,即使在运动中受到冲击仍然能够踩在预定的落脚点上或回退到上一步落脚点,从而保证运动的稳定性.尽管在本文VC 的尺度缩放调节中目标步长无法为0,但能够设置一个很小的步长在实际控制中达到相似的效果,而且通过离散地形时对步长为0 的需求很小.从仿真结果中可以看到控制效果较好,这也为实际机器人控制提供了有力保障.

在后续的研究工作中,可以将算法在真实BHR-6S 机器人平台上进行实验验证.通过为机器人左右方向增加自由度和变尺度虚拟约束,进一步扩展算法在三维模型上的应用,最终使欠驱动双足机器人适应复杂的随机离散地形环境从而走向应用.另外这种尺度因子设计还可以引入在线优化设计,例如使摆动腿能够在跨越过程中,避免碰触路面障碍,从而让机器人能够适应连续地形中的离散落脚点需求.