一种软体机器人仿真方法、装置、电子设备及存储介质与流程

1.本发明属于机器人领域，特别涉及一种软体机器人仿真方法、装置、电子设备及存储介质。


背景技术：

2.在设计软体机器人时，设计人员需要首先给出一个设计方案，该方案可以从仿生学或简化模型分析获得。之后将该设计方案的几何与材料参数输入到仿真场景中，通过基于连续介质假设的离散化方法(常采用有限元法)求解偏微分方程的数值解，得到仿真结果。在对柔性体
‑
柔性体交互情形进行仿真时，在发生接触时需要针对离散化网格的位置关系进行调整，以防止网格穿透与畸形网格现象。在对柔性体
‑
流体交互情形进行仿真时，在同一个更新步内，需要对柔性体及流体分别更新，流体的更新采用光滑粒子流体动力学方法、有限元法、有限体积法等，均存在将流体从仿真对象转变为边界条件的步骤。因此，现有的软体机器人仿真尚无法提供简单、统一的处理手段来解决复杂场景中的柔性体
‑
柔性体及柔性体
‑
流体交互仿真问题。
3.现有技术中，粒子
‑
有限元法(cremonesi et al.,“a state of the art review of the particle finite element method(pfem)”,2020)是一种以有限元法为基础的改进方法。这种方法将网格单元的顶点作为可以自由移动的粒子，在每次根据有限元求解结果进行单元移动后，依据新的顶点分布重新划分单元并用于下一次有限元求解。
4.上述方法在仿真软体机器人时，一定程度上解决了网格穿透的问题，但柔性体与流体依然需要异步求解，且每一次求解之前的网格重划分需要额外的算法设计。


技术实现要素：

5.本说明书实施例的目的是提供一种软体机器人仿真方法、装置、电子设备及存储介质，该方法追踪软体机器人可以通过设置少量标记的物质点完成，无需通过额外的算法解决不同运算步之间网格边界变化的问题，灵活性高。
6.为解决上述技术问题，本说明书实施例通过以下方式实现的：
7.第一方面，提供一种软体机器人仿真方法，该方法包括：
8.确定第一时刻软体机器人物质点的第一运动学参数；
9.通过插值函数确定物质点的第一运动学参数对应的第一网格动力学参数；
10.基于第一网格动力学参数，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数，第二时刻为第一时刻的下一时刻；
11.通过插值函数确定第二网格动力学参数对应的物质点的第二运动学参数。
12.在其中一个实施例中，基于第一网格动力学参数，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数，包括：
13.基于第一网格动力学参数且满足施加于网格的边界条件，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数。
14.在其中一个实施例中，施加于网格的边界条件包括驱动软体机器人的驱动源的施加压力。
15.在其中一个实施例中，确定驱动源的施加压力，包括：
16.标记软体机器人处于初始受压表面的物质点；
17.确定包含物质点在内的若干个邻域物质点组成的邻域；
18.在邻域内任意两个物质点之间构造邻接线；
19.基于邻接线确定物质点组成的邻域的面积、法向量；
20.根据邻域的面积、法向量，确定物质点受压力的大小及方向。
21.在其中一个实施例中，第一网格动力学参数包括网格的速度、应力、质量；
22.基于第一网格动力学参数，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数，包括：
23.根据网格的速度及应力，构造网格位置上力的平衡方程；
24.根据平衡方法，构造非线性方程组；非线性方程组包括所有有质量的网格的位置上力的平衡方程；
25.采用拟牛顿法求解非线性方程组，得到第二网格动力学参数。
26.在其中一个实施例中，第一运动学参数包括物质点的速度、质量、应力、角动量中至少一者。
27.在其中一个实施例中，第一时刻与第二时刻之间间隔预设时间步长。
28.第二方面，提供一种软体机器人仿真装置，该装置包括：
29.第一确定模块，用于确定第一时刻软体机器人物质点的第一运动学参数；
30.第二确定模块，用于通过插值函数确定物质点的第一运动学参数对应的第一网格动力学参数；
31.第三确定模块，用于基于第一网格动力学参数，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数，第二时刻为第一时刻的下一时刻；
32.第四确定模块，用于通过插值函数确定第二网格动力学参数对应的物质点的第二运动学参数。
33.第三方面，提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行程序时实现如第一方面的软体机器人仿真方法。
34.第四方面，提供一种可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如第一方面的软体机器人仿真方法。
35.由以上本说明书实施例提供的技术方案可见，
36.(1)本发明提出一种软体机器人仿真方法、装置、电子设备及存储介质，该方案追踪软体机器人边界可以通过设置少量标记的物质点完成，无需pfem中通过额外的算法解决不同运算步之间网格边界变化的问题，具有更高的灵活性。
37.(2)本发明可以在统一的更新框架内同时完成对不同物体的仿真，同时结合后向欧拉法进行求解，保证了较高的精度。
附图说明
38.为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或
现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
39.图1为本技术实施例提供的软体机器人仿真方法的流程示意图；
40.图2为本技术实施例提供的软体机器人仿真方法的更新步骤示意图；
41.图3为本技术实施例提供的对输入的几何模型进行预处理的示意图；
42.图4为本技术实施例提供的软体机器人仿真装置的结构示意图；
43.图5为本技术实施例提供的电子设备的结构示意图。
具体实施方式
44.下面结合附图和实施例对本技术作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释相关发明，而非对该发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与发明相关的部分。
45.为了使本技术领域的人员更好地理解本技术方案，下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本技术保护的范围。
46.本技术的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”、“第四”等(如果存在)是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便描述的本技术的实施例能够以除了在这里图示或描述的那些以外的顺序实施。
47.此外，术语“包括”和“具有”以及它们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或模块的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或模块，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或模块。
48.需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本技术。
49.现有技术中，软体机器人仿真时，采用例如粒子
‑
有限元法，一种以有限元法为基础的改进方法。这种方法将网格单元的顶点作为可以自由移动的粒子，在每次根据有限元求解结果进行单元移动后，依据新的顶点分布重新划分单元并用于下一次有限元求解。该方法在一定程度上解决了网格穿透的问题，但柔性体与流体依然需要异步求解，且每一次求解之前的网格重划分需要额外的算法设计。
50.基于上述缺陷，本技术实施例提出了一种软体机器人仿真方法，无需通过额外的算法解决不同运算步之间网格边界变化的问题，具有更高的灵活性。
51.参照图1、2，其示出了适用于本技术实施例一种软体机器人仿真方法的流程示意图。本技术中软体机器人仿真方法可以适用于流体驱动软体机器人、线驱动软体机器人等，可以理解的，生物能源驱动软体机器人、化学反应驱动软体机器人等可以通过将其转化为力，采用本技术实施例提出的仿真方法进行仿真。
52.说明：本技术实施例中，x表示位置；v表示速度；m表示质量；f表示力；v表示体积；f
表示形变梯度；上标int表示基于应力表示的内力；上标ext表示由边界条件产生的外力(如重力、电磁力等)；下标i表示网格变量；下标p表示物质点变量；上标n表示第n时刻，上标n 1表示第n 1时刻，第n时刻和第n 1时刻之间间隔为预设时间步长δt。
53.如图1、2所示，一种软体机器人仿真方法，可以包括：
54.s110、确定第一时刻软体机器人物质点的第一运动学参数；
55.s120、通过插值函数确定物质点的第一运动学参数对应的第一网格动力学参数；
56.s130、基于第一网格动力学参数，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数，第二时刻为第一时刻的下一时刻；
57.s140、通过插值函数确定第二网格动力学参数对应的物质点的第二运动学参数。
58.具体的，本实施例中将软体机器人的材料区域离散为一组相对于背景网格运动的物质点，每个物质点均代表一块材料区域并携带有所以物质的物理信息(即运动学参数)，因此，所有物质点的集合代表了软体机器人的整个材料区域。可选的，运动学参数可以包括质量m、速度v、应力f、应变、角动量、形变梯度f、体积v中至少一者。可以理解的，更新时，每个物质点的质量、体积等是定量，而速度、应力、应变等是变量。其中，应力用于流体驱动软体机器人仿真。可以由材料模型(第一pk应力张量形式ψ为应变能密度函数，由材料模型给出)计算应力：
[0059][0060]
其中，
[0061]
其中，第一时刻可以是任意时刻，当第一时刻是初始时刻时，即软体机器人处于静止态或平衡态，此时，第一运动学参数可以是预先设置的初始值，例如，速度和应力都可以是0，可以理解的，初始速度和应力也可以不为0，例如设置的一些边界条件等。当第一时刻不是初始时刻时，第一运动学参数是根据前一时刻的运动学参数更新得到的。
[0062]
其中，第二时刻的前一时刻即为第一时刻。可选的，第一时刻(例如第n时刻)与第二时刻(例如第n 1时刻)之间间隔可以为预设时间步长δt，即每δt更新一次物质点的运动学参数。可以理解的，预设时间步长δt可以根据实际需求进行设定。
[0063]
软体机器人中的网格是指在空间中固定的，起到积分点的作用，用来计算空间导数和求解动量方程。
[0064]
第一网格动力学参数是指物质点的第一运动学参数通过插值函数转移到网格上后对应的动力学参数。
[0065]
示例性的，插值函数可以采用：
[0066][0067]
在进行质量与动量的转移时，物质点p对网格i的贡献度由权重α
ip
求得，其中，α
ip
为：
[0068][0069]
其中，h表示相邻网格之间的距离，物质点p的位置为x
p
＝(x
p
,y
p
,z
p
)，网格i的位置为x
i
＝(x
i
,y
i
,z
i
)。
[0070]
示例性的，将物质点的质量与动量转移到网格上为：
[0071][0072][0073]
其中，为对局部角动量耗散的修正项。
[0074]
第二网格动力学参数是指第一网格动力学参数在预设时间步长δt后，根据隐式时间积分方法更新的网格动力学参数。
[0075]
然后通过插值函数将第二网格动力学参数转换为物质点的第二运动学参数。
[0076]
示例性的，将位置、速度、形变梯度、仿射速度场c的网格动力学参数转换到物质点为：
[0077][0078][0079][0080][0081]
其中，
[0082]
在一个实施例中，基于第一网格动力学参数，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数，可以包括：
[0083]
基于第一网格动力学参数且满足施加于网格的边界条件，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数。
[0084]
具体的，边界条件可以是第一类边界条件、第二类边界条件等。示例性的，边界条件可以包括驱动软体机器人的驱动源的施加压力，还可以包括对软体机器人固定端的约束、在某一个时刻施加的额外速度、外力及设置初始速度等。例如，设置速度v
i
和(或)外力其中，上标ext表示由边界条件产生的外力，下标i表示网格变量。
[0085]
可以理解的，软体机器人的一大常见类型以流体作为驱动源。在运动时，由软体机器人内部的流体提供压强使软体机器人受迫运动。由于软体机器人的形变程度较大，受压表面的形状与位置都随时间发生改变，因而难以给出简单的几何表示方法。
[0086]
在一个实施例中，确定驱动软体机器人的驱动源的施加压力，可以包括：
[0087]
标记软体机器人处于初始受压表面的物质点；
[0088]
确定包含物质点在内的若干个邻域物质点组成的邻域；
[0089]
在邻域内任意两个物质点之间构造邻接线；
[0090]
基于邻接线确定物质点组成的邻域的面积、法向量；
[0091]
根据邻域的面积、法向量，确定物质点受压力的大小及方向。
[0092]
具体的，假设流体驱动软体机器人的初始受压表面上的点始终处于任意时刻的受压表面。
[0093]
在进行算法的初始化时，如图3所示，对输入的几何模型进行预处理：首先，标记出软体机器人处于初始受压表面的物质点；其次，找出包含物质点自身在内的b个邻域物质点组成的邻域并构造邻接线(该邻接线为在邻域内任意两个物质点之间构造的)。在每次更新时，通过邻接线计算某个物质点的b 1个点组成的邻域的面积、法向量，以确定该物质点受驱动力的大小与方向。通过计算标记出的每一个物质点的受驱动力大小和方向，可以近似压强情况。
[0094]
可以理解的，b个邻域物质点可以根据实际需求进行设置，例如3个邻域物质点、5个邻域物质点等。
[0095]
在每次更新时，通过邻接线计算某个物质点的b 1个点组成的邻域的面积、法向量，以确定该物质点受驱动力的大小与方向，可以包括：
[0096]
设共有q个在初始受压表面时被标记出的物质点，暂时记为(x
p1
…
x
pm
…
x
pn
)，m∈[1,q]，进行如下操作：
[0097]
第一步：对于第m个被标记出的物质点x
pm
，找出其k个邻域物质点；
[0098]
第二步：将第一步中b 1个物质点的坐标作为列向量，排列为3
×
(b 1)的矩阵g’，求该矩阵g’的每一行的平均值组成列向量h，将g’的每一行减去列向量h，完成零均值化，得到矩阵g；
[0099]
第三步，计算gg
t
，其中t表示矩阵的转置，得到3
×
3的方阵l，对该方阵做svd分解，得到3个特征值与3个对应的特征向量。取最小的特征值向量v’所对应的特征向量作为该物质点上所受驱动力所在的直线，并对其作归一化得到单位向量v＝v’/||v’||，||v’||表示向量v’的模；
[0100]
第四步，记w(x
pm
,v)表示如下操作：将x
pm
的坐标按照v的方向移动1个单位的距离，取该移动后坐标最近的网格的质量m
i
；完成两次这样的操作m
i1
＝w(x
pm
,v)，m
i2
＝w(x
pm
,
‑
v)，若m
i1
>m
i2
，则该物质点上所受驱动力的方向为v，若m
i1
<m
i2
，则该物质点上所受驱动力的方向为
‑
v；
[0101]
第五步，将其余两个特征向量也分别做归一化得到两个正交的单位向量，组成单位正交基，每个向量作为行向量，组成一个2
×
3的矩阵s，并右乘k 1个物质点的坐标，将三维坐标投影至二维的平面，得到这k 1个物质点在该平面的二维坐标；
[0102]
第六步，求取该b 1个物质点在该平面的外包凸多边形，并计算该凸多边形的表面积，使用该面积乘以此刻驱动源所提供的压强，即得到一个力的大小，该值即该物质点x
pm
上所受驱动力的大小f
pext
。
[0103]
可以理解的，将该力从物质点转移到网格上参与动力学更新，可以通过f
iext
＝σ
p
α
ip
f
pext
完成。
[0104]
前向欧拉法的显式时间积分会造成运算不稳定，精度低等问题。后向欧拉法具有
运算稳定、精度高的优点。但需要将后向欧拉法与更新框架相结合。由于物质点的数量多，作为自变量时的维数高，而网格的数量较少，且通过物质点到网格的物理量转移，网格上已经具有动力学信息，如力和速度。因而，将网格作为隐式时间积分的变量是合理的选择。
[0105]
在一个实施例中，基于第一网格动力学参数，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数，可以包括：
[0106]
根据网格的速度及应力，构造网格位置上力的平衡方程；
[0107]
根据平衡方法，构造非线性方程组；非线性方程组包括所有有质量的网格的位置上力的平衡方程；
[0108]
采用拟牛顿法求解非线性方程组，得到第二网格动力学参数。
[0109]
具体的，网格的加速度a
in 1
可转化为网格位置的函数a
i
(x
in 1
)，而网格上的力f
in 1
也是网格位置的函数f
i
(x
in 1
)，因而可以构造在网格更新后位置上的力的平衡方程m
i
a
i
(x
in 1
)
‑
f
i
(x
in 1
)＝r。通过下述方式求使得r＝0。
[0110]
其中，网格的速度为：
[0111][0112]
网格的加速度用网格的位移和当前时刻的速度的形式可表示为：
[0113][0114]
网格的加速度用网格上合力的形式可表示为：
[0115][0116]
其中，和均与δx
i
有关，记
[0117]
则
[0118]
当为真实值的时候，两个表达式的加速度应相等，即
[0119][0120]
由于并非网格的真实位移，而是进行时间积分的中间变量，在每个时间步之间重置为0，因而上式在未进行迭代求解的真实值之前不为0，即
[0121][0122]
此处，的上标0表示第0次迭代求解的结果，即的初始值为0，类似的表示第k次迭代求解的结果，最终目标是求得一个使得
[0123]
是一个关于的非线性方程组，需要通过迭代的方式求得。
[0124]
由于网格数量多，运算复杂，的雅可比矩阵难以计算，因此考虑采用broyden形式的拟牛顿法进行迭代求解，该方法只要求未知量初始值时的雅可比矩阵令b0＝j0，作为对雅可比矩阵的初始估计，之后每一次迭代求解的雅可比矩阵估计值b
k
可以通过迭代得到而不用耗费大量计算资源计算j
k
。
[0125]
使用broyden形式的拟牛顿法的具体步骤为
[0126]
1)构造线性方程组(第一次迭代时为其中p
k
为该线性方程组的未知量，也是在第k步的增量，解该线性方程组得到p
k
；
[0127]
2)通过线搜索算法确定增量p
k
的系数α，并更新未知量
[0128]
3)计算更新b
k 1
需要的中间变量s
k
和y
k
；
[0129]
其中，s
k
＝αp
k
[0130]
而y
k
为r的增量，即
[0131]
4)更新下一次迭代用的对雅可比矩阵的估计b
k 1
：
[0132][0133]
其中为s
k
的转置。
[0134]
至此，完成一次broyden形式的拟牛顿法的迭代求解。
[0135]
当满足条件，即小于设定的阈值时，停止迭代，令
[0136]
则网格的位移为速度为随后通过网格动力学参数转换到物质点的公式将其转移到物质点上。
[0137]
参照图4，其示出了根据本技术一个实施例描述的软体机器人仿真装置的结构示意图。
[0138]
如图4所示，软体机器人仿真装置，可以包括：
[0139]
第一确定模块410，用于确定第一时刻软体机器人物质点的第一运动学参数；
[0140]
第二确定模块420，用于通过插值函数确定物质点的第一运动学参数对应的第一网格动力学参数；
[0141]
第三确定模块430，用于基于第一网格动力学参数，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数，第二时刻为第一时刻的下一时刻；
[0142]
第四确定模块440，用于通过插值函数确定第二网格动力学参数对应的物质点的第二运动学参数。
[0143]
可选的，第三确定模块430还用于：
[0144]
基于第一网格动力学参数且满足施加于网格的边界条件，通过隐式时间积分方法，确定第二时刻的第二网格动力学参数。
[0145]
可选的，施加于网格的边界条件包括驱动软体机器人的驱动源的施加压力。
[0146]
可选的，第三确定模块430还用于：
[0147]
标记软体机器人处于初始受压表面的物质点；
[0148]
确定包含物质点在内的若干个邻域物质点组成的邻域；
[0149]
在邻域内任意两个物质点之间构造邻接线；
[0150]
基于邻接线确定物质点组成的邻域的面积、法向量；
[0151]
根据邻域的面积、法向量，确定物质点受压力的大小及方向。
[0152]
可选的，第一网格动力学参数包括网格的速度、应力、质量；第三确定模块430还用于：
[0153]
根据网格的速度及应力，构造网格位置上力的平衡方程；
[0154]
根据平衡方法，构造非线性方程组；非线性方程组包括所有有质量的网格的位置上力的平衡方程；
[0155]
采用拟牛顿法求解非线性方程组，得到第二网格动力学参数。
[0156]
可选的，第一运动学参数包括物质点的速度、质量、应力、角动量中至少一者。
[0157]
可选的，第一时刻与第二时刻之间间隔预设时间步长。
[0158]
本实施例提供的一种软体机器人仿真装置，可以执行上述方法的实施例，其实现原理和技术效果类似，在此不再赘述。
[0159]
图5为本发明实施例提供的一种电子设备的结构示意图。如图5所示，示出了适于用来实现本技术实施例的电子设备500的结构示意图。
[0160]
如图5所示，电子设备500包括中央处理单元(cpu)501，其可以根据存储在只读存储器(rom)502中的程序或者从存储部分508加载到随机访问存储器(ram)503中的程序而执行各种适当的动作和处理。在ram 503中，还存储有设备500操作所需的各种程序和数据。cpu 501、rom 502以及ram 503通过总线504彼此相连。输入/输出(i/o)接口505也连接至总线504。
[0161]
以下部件连接至i/o接口505：包括键盘、鼠标等的输入部分506；包括诸如阴极射线管(crt)、液晶显示器(lcd)等以及扬声器等的输出部分507；包括硬盘等的存储部分508；以及包括诸如lan卡、调制解调器等的网络接口卡的通信部分509。通信部分509经由诸如因特网的网络执行通信处理。驱动器510也根据需要连接至i/o接口506。可拆卸介质511，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器510上，以便于从其上读出的计算机程序根据需要被安装入存储部分508。
[0162]
特别地，根据本公开的实施例，上文参考图1描述的过程可以被实现为计算机软件程序。例如，本公开的实施例包括一种计算机程序产品，其包括有形地包含在机器可读介质上的计算机程序，计算机程序包含用于执行上述软体机器人仿真方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分509从网络上被下载和安装，和/或从可拆卸介质511被安装。
[0163]
附图中的流程图和框图，图示了按照本发明各种实施例的系统、方法和计算机程
序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分，前述模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个接连地表示的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
[0164]
描述于本技术实施例中所涉及到的单元或模块可以通过软件的方式实现，也可以通过硬件的方式来实现。所描述的单元或模块也可以设置在处理器中。这些单元或模块的名称在某种情况下并不构成对该单元或模块本身的限定。
[0165]
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、笔记本电脑、行动电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
[0166]
作为另一方面，本技术还提供了一种存储介质，该存储介质可以是上述实施例中前述装置中所包含的存储介质；也可以是单独存在，未装配入设备中的存储介质。存储介质存储有一个或者一个以上程序，前述程序被一个或者一个以上的处理器用来执行描述于本技术的软体机器人仿真方法。
[0167]
存储介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd
‑
rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
[0168]
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
[0169]
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。