梁柱端板连接节点优化设计智能系统构建及参数识别方法与流程

1.本发明涉及梁柱端板连接节点优化设计智能系统构建及参数识别方法，属 于梁柱连接节点技术领域。


背景技术：

2.对节点结构进行优化设计势必要将结构各个构成参数进行综合评估，采用 试验研究，可以得到建立适合物理模型所要求的较完整信息，从而可以作为提 出设计模型的基础。然而，试验方法既不经济又费时费力，虽然根据工程设计 的需要，可以作一些近似简化，但一般试验并不能完全考虑影响节点性能的重 要因素。随着计算机应用技术的发展，有限元模型仿真技术也得到了迅速地发 展，可以提供更全面的分析数据。目前对结构参数评估采用的基本方法为有限 元分析和灵敏度计算。
3.有限元分析软件中，ansys的pds模块具有有限元分析与灵敏度计算功能。 nsys有限元pds分析采用蒙特卡罗法，该方法基于大数定理，当在进行大量抽 样时，有限元分析过程就会随着抽样次数的增加而增加，要想获得更精确的结 构参数灵敏度计算结果，要进行大量的有限元分析。对于简单结构而言，导致 计算时间过长；复杂结构仿真模型的分析都是基于初次参数化建模，当结构模 型的参数变化范围增大时，模型建立过程中会出现错误的网格划分机制导致计 算分析终止，因此利用有限元软件进行概率设计时，参数的变化范围往往只能 局限于一个极小的范围内，由此得到的灵敏度计算结果也是基于小范围参数变 化。
4.基于ansys有限元分析和灵敏度计算的结构参数重要性识别量化方法存在 以下问题：
5.1对端板连接节点结构进行ansys有限元pds分析，当参数变化范围超过 2％时，ansys报错不能运行，参数变化范围小于2％得到的灵敏度计算结果对于 结构的优化设计没有实际意义；
6.2梁柱节点结构优化设计进行参数识别采用灵敏度分析方法，该方法得到 的是局部灵敏度，无法度量全局的参数敏感性，且掩盖了不同类型物理参数量 级差别很大时分散度对结构的全局影响程度。


技术实现要素：

7.针对现有技术的上述问题，本发明提供了一种梁柱端板连接节点优化设计 智能系统构建及参数识别方法，用以解决梁柱连接节点这一类有限元模型较复 杂的结构在结构优化时对重要参数识别量化采用anasys的局限性问题。
8.为解决上述问题，本发明采用了如下技术方案：
9.梁柱端板连接节点优化设计智能系统构建及参数识别方法，包括如下步骤：
10.步骤一：参数初始化，设定元粒族群规模g，权重因子w，竞技因子m1、 m2，元粒群优化算法的最大竞技代数t，混沌映射的迭代次数q，设置当前竞技 代数t＝1，能量函数e，能
量函数变化次数δ＝0，元粒群优化算法能量变化次数 限值δ
max
，适应度函数u，元粒移动位移函数s，可行域中随机产生两类元粒， 分别表示神经网络的权值和阀值；其中，适应度函数u用以早熟判断，元粒移 动位移函数s用以稳态判断；
11.步骤二：根据上述参数构建能量函数e，能量函数e用适应度函数u和元 粒移动位移函数s表示，评价每个元粒的初始适应度值，选出初始化元粒群的全 局最好地点矢量g
e
，即能量函数最低值，设p
e
为当前的元粒的最好地点矢量；
12.步骤三：更新元粒的地点矢量和速度，产生新的元粒，评价每个元粒的能 量值，如果第i个元粒的能量函数值比适应度函数u低或元粒移动位移函数s低， 则计算产生q个混沌列点，选取其中能量函数最低的点作为元粒的新地点矢量； 如果第i个元粒的能量函数值比适应度函数u高且比元粒移动位移函数s高，则 直接进入步骤四；
13.步骤四：评价新元粒的能量值，如果比更新前的能量函数值低，则更新p
e
的值及其对应的地点矢量；如果所有元粒中最低的能量值低于g
e
，则更新g
e
和 相应的地点矢量，且δ＝0,否则，δ＝δ 1；
14.步骤五：如果δ≥δ
max
，则使用误差反向传播算法在g
e
附近进行局部搜索； 如果搜索结果比g
e
好，则用此搜索结果代替p
e
同时更新对应的地点矢量，否则 用此搜索结果代替p
e
中性能最差的个体，δ＝0，t＝t 1；当t≥t时，寻优结束， 否则返回步骤三；
15.步骤六：保存结果，梁柱端板连接节点优化设计智能系统构建完成；
16.步骤七：对智能系统结构优化参数识别量化计算。
17.进一步地，所述步骤三中，更新元粒的地点矢量和速度，产生新的元粒的 具体步骤如下：
18.神经网络采用matlab实现，计算过程如式1所示：
[0019][0020]
输入矢量为n维变量，隐层有β个神经元，输出矢量为m维变量，输入层 到隐层的权值为w
ih
，阀值为b
h
，隐层到输出层权值为w
hl’，阀值为b
l’，f为传递 函数；d为输入神经元值，dn’为隐层神经元值；采用levenberg
‑
marquardt法 对神经网络进行训练，确定搜索方向的表达式为(2)
[0021][0022]
其中：d(x
(ι)
)为第ι次迭代的搜索方向，x
(ι)
为第ι次迭代输入值，λ
(ι)
为第 ι次迭代梯度值，为梯度算子；h为hessian矩阵，可通过雅可比矩阵进行计 算，i为单位矩阵；
[0023]
元粒竞技算法是在v维搜索空间中，由k个元粒组成族群g＝(g
(1)
,g
(2)
,
…
,g
(k)
)
t
， 第i个元粒的地点矢量为x
i
＝(x
i,1
,x
i,2
,
…
,x
i,v
)
t
，速度为v
i
＝(v
i,1
,v
i,2
,..,.v
i,v
)
t；
个体极值 为p
i
＝(p
i,1
,p
i,2
,
…
,p
i,v
)
t
，族群的全局极值为p
g
＝(p
g,1
,p
g,2
,
…
,p
g,k
)
t
；元粒群寻优过 程用以下函数(3)和(4)实现：
[0024]
[0025][0026]
其中，χ为区间(0,1)的随机数，m1、m2为竞技因子，为元粒i在第j次迭 代中第v维的速度，为元粒i在第j次迭代中第v维的地点矢量，为元粒i 在第v维的个体极值地点矢量，为群体在第v维的全局极值地点矢量。
[0027]
进一步地，所述步骤三中，更新元粒的地点矢量和速度，产生新的元粒的 具体步骤还包括如下：
[0028]
在元粒群运动过程引入权重因子w以调整元粒的速度；带有权重因子的元 粒群算法为式(5)和(6)：
[0029][0030]
w(j)＝exp[
‑
u(j)/u(j
‑
1)]
ꢀꢀ
(6)
[0031]
式中，为元粒i在第v维的最好地点矢量，为群体在第v维的全局最 好地点矢量，w(j)为第j次迭代的权重因子，u(j)为元粒第j次迭代的适应度函数 值，u(j
‑
1)为元粒第j
‑
1次迭代的适应度函数值。
[0032]
进一步地，所述步骤三中，适应度函数值计算方法如下：
[0033]
以目标函数作为适应度函数判断元粒的状态，将神经网络输出值与实际响 应值的相对误差作为适应度函数，建立函数(7)
[0034][0035]
其中，ζ
i
为神经网络输出值，κ(x
i
)为结构实际响应值；
[0036]
通过元粒当前移动距离s来判别稳定状态判别，建立函数(8)
[0037]
s＝|x
i
(t)
‑
x
i
(t
‑
1)|
ꢀꢀ
(8)
[0038]
式中，t为元粒竞技代数，当s或u(x)小于一定限值时即认为元粒群算法达 到局部最优解p
g
＝(p
g,1
,p
g,2
,
…
p
g,v
)
t
，需利用混沌搜索将变量混沌化。
[0039]
进一步地，所述步骤三中，计算产生q个混沌列点的具体方法如下：
[0040]
采用logistic混沌系统，其映射为式(9)
[0041]
x
n 1
＝1
‑
λx
n2
，n＝1,2,...
ꢀꢀ
(9)
[0042]
1)将元粒x
i
的每一维x
i,v
映射到区间[
‑
1,1]；
[0043][0044]
其中[a
v
,b
v
]为第v维x
i,v
的定义域，μ为一极小值。
[0045]
2)由式c
i(j 1)
＝1
‑
2(c
i(j)
)2迭代j次产生混沌序列
[0046]
3)根据混沌列点的相应位置，按照式(10)有：
[0047]
将混沌序列中的点映射到原搜索空间，每 一维混沌序列点为
[0048]
进一步地，所述步骤七具体为：
[0049]
结构优化参数识别量化是对结构各个构成参数进行灵敏度数值的计算，该 方法基于智能系统的预测结果，根据概率密度函数对输入变量进行蒙特卡罗随 机抽样，考察对应于随机变量取值范围内的输出响应的变化量，从而判定各随 机变量的灵敏度；综合考虑各随机变量对结构响应的影响程度，采用重要性度 量参数pearson相关系数度量随机变量与结构响应的相关关系，为消除各随机变 量的量纲和测量尺度的影响，在计算过程中引入随机变量的变异系数γ
r
，具体 计算公式为式(11)和式(12)
[0050][0051][0052]
式中，r
r
为第r变量对结构的pearson相关系数，s
r
为第r变量对结构的灵 敏度数值，ξ为随机变量总维数，α为蒙特卡罗抽样次数，x
i
为随机变量，ζ
i
为 在输入随机变量x
i
条件下的网络输出值，γ
r
为蒙特卡罗抽取随机变量样本点对应 的变异系数，s
r
的绝对值反应了随机变量对结构响应的影响程度，正负反应了随 机变量与结构响应的相关程度；该方法计算得到的灵敏度数值为全局随机灵敏 度。
[0053]
进一步地，所述步骤七之后，还包括步骤八，即精度判断过程，具体如下：
[0054]
精度判断过程是根据智能系统计算值和实际值的偏差进行判别，其中：
[0055]
(1)r2[0056][0057]
其中，f(a
i
)为智能系统的输出响应值，y(a
i
)为结构的真实响应值，为 m个真实响应值的平均值；如果r2的值越接近1，说明智能系统对结构的模拟 越精确；
[0058]
(2)r
aae
[0059][0060]
其中，s
td
代表了m个结构真实响应值的标准差，r
aae
的值越小表示智能系 统值越精确；
[0061]
(3)r
mae
[0062][0063]
r
mae
反应了智能系统与实体模型在局部某一区域的最大误差，r
mae
的值越 小表示智能系统的精度越高；r2和r
aae
是对智能系统的整体评价，r
mae
是对 智能系统的局部评价。
[0064]
进一步地，所述步骤八之后，还包括步骤九，即将matlab与ansys对接， 具体对接方法如下：
[0065]
(1)将均匀设计表保存成文件；
[0066]
(2)在matlab中循环读取均匀设计表中的参数，并保存至相应矩阵中；
[0067]
(3)建立利用ansys安装地址并建立空文档使得在matlab中打开ansys分 析软件，并将ansys分析过程写入空文档；
[0068]
(4)在ansys中对节点模型应用参数化建模，设计参数采用变量形式，设 立与设计参数对应的中间存储变量，将设计参数的中间存储变量读入matlab保 存的数据，进行有限元分析，同理建立与结构响应值对应的中间存储变量；
[0069]
(5)matlab读取有限元分析结果保存变量。
[0070]
进一步地，所述步骤九中的均匀设计表中的均匀设计是使每个因素的每个 水平仅做一次试验，试验点在试验范围内均匀分布，通过最少的试验来获得最 多的信息；这种设计方法适合于多因素多水平的试验和系统模型完全未知的情 况；
[0071]
均匀设计表采用好格子点法设计样本点，构造方法为：1)给定实验数ε，寻 找比ε小的整数o，且使ε和o的最大公约数为1，符合这些条件的正整数组成 向量o；2)均匀设计表的第j列由a
ij
＝io
j
[modε]生成，modε为同余运算，若io
j
超过试验数ε，则用它减去一个ε的适当倍数，使其差落在[1,ε]之中；a
ij
可以递 推生成：
[0072][0073]
本发明取得了如下技术效果：
[0074]
1)构建智能系统，采用混沌元粒群竞技智能算法优化神经网络参数，提高了 计算结果的精度和效率。
[0075]
2)运用智能系统解决了节点结构ansys pds模块分析参数变化范围局限 于某一较小数值的问题，该方法可以根据节点结构自身特点设定参数范围。
[0076]
3)将全局灵敏度计算方法与结构的随机性相结合，改善了各个参变量的量纲 和取值范围分散度较大时对结构灵敏度值的真实影响。
附图说明
[0077]
为了更清楚的介绍本发明的方案，下面对本方案的所需要的附图进行简单 介绍：
[0078]
图1为本发明智能系统构建流程；
[0079]
图2为matlab与ansys对接流程；
[0080]
图3为输入参数变化范围为2％时的节点智能系统与有限元模型相对误差；
[0081]
图4为增大输入参数变化范围时的节点智能系统与有限元模型相对误差。
具体实施方式
[0082]
下面结合附图1
‑
4以对本发明做进一步说明。
[0083]
如图1所示，本发明的梁柱端板连接节点优化设计智能系统构建及参数识 别方法，包括如下步骤：
[0084]
步骤一：参数初始化，设定元粒族群规模g，权重因子w，竞技因子m1、 m2，元粒群优化算法的最大竞技代数t，混沌映射的迭代次数q，设置当前竞技 代数t＝1，能量函数e，能量函数变化次数δ＝0，元粒群优化算法能量变化次数 限值δ
max
，适应度函数u，元粒移动位移函数s，可行域中随机产生两类元粒， 分别表示神经网络的权值和阀值；其中，适应度函数u用以早熟判断，元粒移 动位移函数s用以稳态判断；
[0085]
步骤二：根据上述参数构建能量函数e，能量函数e用适应度函数u和元 粒移动位移函数s表示，评价每个元粒的初始适应度值，选出初始化元粒群的全 局最好地点矢量g
e
，即能量函数最低值，设p
e
为当前的元粒的最好地点矢量；
[0086]
步骤三：更新元粒的地点矢量和速度，产生新的元粒；具体步骤如下：
[0087]
神经网络采用matlab实现，计算过程如式1所示：
[0088][0089]
输入矢量为n维变量，隐层有β个神经元，输出矢量为m维变量，输入层 到隐层的权值为w
ih
，阀值为b
h
，隐层到输出层权值为w
hl’，阀值为b
l’，f为传递 函数；d为输入神经元值，dn’为隐层神经元值；采用levenberg
‑
marquardt法 对神经网络进行训练，确定搜索方向的表达式为(2)
[0090][0091]
其中：d(x
(ι)
)为第ι次迭代的搜索方向，x
(ι)
为第ι次迭代输入值，λ
(ι)
为第 ι次迭代梯度值，为梯度算子；h为hessian矩阵，可通过雅可比矩阵进行计 算，i为单位矩阵；
[0092]
元粒竞技算法是在v维搜索空间中，由k个元粒组成族群g＝(g
(1)
,g
(2)
,
…
,g
(k)
)
t
， 第i个元粒的地点矢量为x
i
＝(x
i,1
,x
i,2
,
…
,x
i,v
)
t
，速度为v
i
＝(v
i,1
,v
i,2
,..,.v
i,v
)
t；
个体极值 为p
i
＝(p
i,1
,p
i,2
,
…
,p
i,v
)
t
，族群的全局极值为p
g
＝(p
g,1
,p
g,2
,
…
,p
g,k
)
t
；元粒群寻优过 程用以下函数(3)和(4)实现：
[0093][0094][0095]
其中，χ为区间(0,1)的随机数，m1、m2为竞技因子，为元粒i在第j次迭 代中第v维的速度，为元粒i在第j次迭代中第v维的地点矢量，为元粒i 在第v维的个体极值地点矢量，为群体在第v维的全局极值地点矢量。
[0096]
在元粒群运动过程引入权重因子w以调整元粒的速度；带有权重因子的元 粒群算法为式(5)和(6)：
[0097]
[0098]
w(j)＝exp[
‑
u(j)/u(j
‑
1)]
ꢀꢀ
(6)
[0099]
式中，为元粒i在第v维的最好地点矢量，为群体在第v维的全局最 好地点矢量，w(j)为第j次迭代的权重因子，u(j)为元粒第j次迭代的适应度函数 值，u(j
‑
1)为元粒第j
‑
1次迭代的适应度函数值。
[0100]
适应度函数值计算方法如下：
[0101]
以目标函数作为适应度函数判断元粒的状态，将神经网络输出值与实际响 应值的相对误差作为适应度函数，建立函数(7)
[0102][0103]
其中，ζ
i
为神经网络输出值，κ(x
i
)为结构实际响应值；
[0104]
通过元粒当前移动距离s来判别稳定状态判别，建立函数(8)
[0105]
s＝|x
i
(t)
‑
x
i
(t
‑
1)|
ꢀꢀ
(8)
[0106]
式中，t为元粒竞技代数，当s或u(x)小于一定限值时即认为元粒群算法达 到局部最优解p
g
＝(p
g,1
,p
g,2
,
…
p
g,v
)
t
，需利用混沌搜索将变量混沌化。
[0107]
评价每个元粒的能量值，如果第i个元粒的能量函数值比适应度函数u低或 元粒移动位移函数s低，则计算产生q个混沌列点，计算产生q个混沌列点的具 体方法如下：
[0108]
采用logistic混沌系统，其映射为式(9)
[0109]
x
n 1
＝1
‑
λx
n2
，n＝1,2,...
ꢀꢀ
(9)
[0110]
1)将元粒x
i
的每一维x
i,v
映射到区间[
‑
1,1]；
[0111][0112]
其中[a
v
,b
v
]为第v维x
i,v
的定义域，μ为一极小值。
[0113]
2)由式c
i(j 1)
＝1
‑
2(c
i(j)
)2迭代j次产生混沌序列
[0114]
3)根据混沌列点的相应位置，按照式(10)有：
[0115]
将混沌序列中的点映射到原搜索空间，每 一维混沌序列点为
[0116]
选取其中能量函数最低的点作为元粒的新地点矢量；如果第i个元粒的能量 函数值比适应度函数u高且比元粒移动位移函数s高，则直接进入步骤四；
[0117]
步骤四：评价新元粒的能量值，如果比更新前的能量函数值低，则更新p
e
的值及其对应的地点矢量；如果所有元粒中最低的能量值低于g
e
，则更新g
e
和 相应的地点矢量，且δ＝0,否则，δ＝δ 1；
[0118]
步骤五：如果δ≥δ
max
，则使用误差反向传播算法在g
e
附近进行局部搜索； 如果搜索结果比g
e
好，则用此搜索结果代替p
e
同时更新对应的地点矢量，否则 用此搜索结果代替p
e
中性能最差的个体，δ＝0，t＝t 1；当t≥t时，寻优结束， 否则返回步骤三；
[0119]
步骤六：保存结果，梁柱端板连接节点优化设计智能系统构建完成；
[0120]
步骤七：对智能系统结构优化参数识别量化计算。具体为：
[0121]
结构优化参数识别量化是对结构各个构成参数进行灵敏度数值的计算，该 方法基于智能系统的预测结果，根据概率密度函数对输入变量进行蒙特卡罗随 机抽样，考察对应于随机变量取值范围内的输出响应的变化量，从而判定各随 机变量的灵敏度；综合考虑各随机变量对结构响应的影响程度，采用重要性度 量参数pearson相关系数度量随机变量与结构响应的相关关系，为消除各随机变 量的量纲和测量尺度的影响，在计算过程中引入随机变量的变异系数γ
r
，具体 计算公式为式(11)和式(12)
[0122][0123][0124]
式中，r
r
为第r变量对结构的pearson相关系数，s
r
为第r变量对结构的灵 敏度数值，ξ为随机变量总维数，α为蒙特卡罗抽样次数，x
i
为随机变量，ζ
i
为 在输入随机变量x
i
条件下的网络输出值，γ
r
为蒙特卡罗抽取随机变量样本点对应 的变异系数，s
r
的绝对值反应了随机变量对结构响应的影响程度，正负反应了随 机变量与结构响应的相关程度；该方法计算得到的灵敏度数值为全局随机灵敏 度。
[0125]
步骤七之后，还包括步骤八，即精度判断过程，具体如下：
[0126]
精度判断过程是根据智能系统计算值和实际值的偏差进行判别，其中：
[0127]
(1)r2[0128][0129]
其中，f(a
i
)为智能系统的输出响应值，y(a
i
)为结构的真实响应值，为 m个真实响应值的平均值；如果r2的值越接近1，说明智能系统对结构的模拟 越精确；
[0130]
(2)r
aae
[0131][0132]
其中，s
td
代表了m个结构真实响应值的标准差，r
aae
的值越小表示智能系 统值越精确；
[0133]
(3)r
mae
[0134][0135]
r
mae
反应了智能系统与实体模型在局部某一区域的最大误差，r
mae
的值越 小表示智能系统的精度越高；r2和r
aae
是对智能系统的整体评价，r
mae
是对 智能系统的局部评价。
[0136]
步骤八之后，还包括步骤九，即将matlab与ansys对接，如图2所示，具 体对接方法
如下：
[0137]
(1)将均匀设计表保存成文件；
[0138]
(2)在matlab中循环读取均匀设计表中的参数，并保存至相应矩阵中；
[0139]
(3)建立利用ansys安装地址并建立空文档使得在matlab中打开ansys分 析软件，并将ansys分析过程写入空文档；
[0140]
(4)在ansys中对节点模型应用参数化建模，设计参数采用变量形式，设 立与设计参数对应的中间存储变量，将设计参数的中间存储变量读入matlab保 存的数据，进行有限元分析，同理建立与结构响应值对应的中间存储变量；
[0141]
(5)matlab读取有限元分析结果保存变量。
[0142]
其中，步骤九中的均匀设计表中的均匀设计是使每个因素的每个水平仅做 一次试验，试验点在试验范围内均匀分布，通过最少的试验来获得最多的信息； 这种设计方法适合于多因素多水平的试验和系统模型完全未知的情况；
[0143]
均匀设计表采用好格子点法设计样本点，构造方法为：1)给定实验数ε，寻 找比ε小的整数o，且使ε和o的最大公约数为1，符合这些条件的正整数组成 向量o；2)均匀设计表的第j列由a
ij
＝io
j
[modε]生成，modε为同余运算，若io
j
超过试验数ε，则用它减去一个ε的适当倍数，使其差落在[1,ε]之中；a
ij
可以递 推生成：
[0144][0145]
实施例1
[0146]
外伸端板连接节点，梁和柱分别采用ipe300和he340m型钢，螺栓采用等 级为8.8级手工加预应力拧紧的m20全螺纹线高强螺栓。
[0147]
(1)输入参数变化范围为2％
[0148]
节点参数变化区间取值为输入参数初始值的2％，如表1所示。随机抽取50 个样本点进行智能系统构建分析得出相应的结构响应值即半刚性节点的转角 值，并将其与有限元分析的结果进行对比，如图3所示，最大误差小于1％，说 明智能系统的计算结果可应用来进行节点结构重要参数识别量化。依赖有限元 分析软件进行200次仿真模拟概率分析用时约26.67小时，利用智能系统用时约 6.8小时，计算效率得到了显著提高。智能系统基本参数：最大迭代次数为100， 学习因子分别为1.3和2.7，混沌算法的最大迭代次数为30，神经网络的输入神 经元为8，输出神经元为1，隐藏层神经元为10，训练函数采用trainlm函数， 由输入层到隐藏层神经元的传递函数采用logsig函数，由隐藏层神经元到输出 层神经元的传递函数采用purelin函数。
[0149]
表1输入参数取值区间
[0150][0151]
随机分析抽样次数为100000次，计算各个输入参数对输出参数的灵敏度值， 灵敏度计算结果如表2所示。
[0152]
表2节点相关系数
[0153][0154]
(2)增大输入参数变化范围
[0155]
如果要得到在输入参数较大范围内的弯节点结构参数重要性量化指标(灵 敏度数值)，目前常用的ansys有限元分析软件不能达到计算目的。原因主要是 当参数变化范围较大时，采用ansys有限元分析软件划分网格时依旧采用元模型 的建模方法，导致网格划分出错计算终止。
[0156]
采用智能系统。增大输入参数变化范围，智能系统预测值与有限元值的相 对误差如图4所示。节点参数变化范围和识别量化指标如表3所示。
[0157]
表3节点参数表
[0158][0159]
(注：表中各尺寸参数均为本领域通用公知参数，h，b，t，t
w
，t
h
，h
mid
，h
up
，v，h
bt
分别为 梁高、梁宽、钢板厚度、腹板厚度、端板厚度、中层螺栓至端板中轴线距离，上层螺栓至端 板中轴线距离，螺栓纵轴线至端板纵轴线距离，底层螺栓至端板中轴线距离)
[0160]
综上，尽管已经对本发明的实施例进行描述，可以理解的是，上述实施例 是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本 发明的原理和宗旨的情况下可以对上述实施例进行变化、修改、替换和变型。