一种基于欧拉拉格朗日耦合方法的流体仿真方法.pdf

《一种基于欧拉拉格朗日耦合方法的流体仿真方法.pdf》由会员分享，可在线阅读，更多相关《一种基于欧拉拉格朗日耦合方法的流体仿真方法.pdf（12页完整版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 104268943 A (43)申请公布日 2015.01.07 CN 104268943 A (21)申请号 201410508371.0 (22)申请日 2014.09.28 G06T 17/30(2006.01) G06T 15/00(2011.01) (71)申请人 北京航空航天大学 地址 100191 北京市海淀区学院路 37 号 申请人 北京天人同达软件科技有限公司 (72)发明人 郝爱民 李帅 高阳 (74)专利代理机构 北京科迪生专利代理有限责 任公司 11251 代理人 杨学明 顾炜 (54) 发明名称 一种基于欧拉 - 拉格朗日耦合方法的流体仿 真。

2、方法 (57) 摘要 本发明公开了一种基于欧拉 - 拉格朗日耦合 方法的流体仿真方法。其以实现高效的流体仿真 为目的， 结合欧拉法和拉格朗日法求解 N-S 方程 进行流体仿真的核心思想， 在SPH和LBM的流体仿 真基础上， 围绕多方法耦合的计算方法展开研究。 首先利用基于欧拉网格的 LBM 方法进行大规模流 场建模， 构建流体仿真的主体部分， 然后利用基于 拉格朗日粒子思想的 SPH 方法仿真诸如浪花、 水 珠等流体细节， 并通过设计耦合算法整合 LBM 流 体与 SPH 流体， 实现逼真、 快速、 兼顾细节的大规 模流体仿真。 (51)Int.Cl. 权利要求书 2 页 说明书 6 页 附。

3、图 3 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书6页 附图3页 (10)申请公布号 CN 104268943 A CN 104268943 A 1/2 页 2 1. 一种基于欧拉 - 拉格朗日耦合方法的流体仿真方法， 其特征在于包括以下四个步 骤 ： 步骤 (1)、 大规模流体的物理仿真 ： 利用求解纳维 - 斯托克斯方程 (Navier-Stokes， N-S) 的 Lattice Bhatnagar-gross-krook(LBGK) 控制方程， 采用 n 维离散空间的 m 速度 (DnQm) 多维离散网格模型， 进行大规模流场建模， 这是流体。

4、仿真的主体部分 ； 步骤 (2)、 流体自由表面追踪 ： 应用 Particle Level Set Method(PLSM) 方法构建流 体自由表面 ； 输出步骤 (1) 中得到的流体网格， 用 PLSM 方法标记流体运动界面在 Lattice Boltzmann Method(LBM) 网格内部的位置， 追踪流体界面， 以此作为流体的近似表面 ； 步骤 (3)、 流体细节仿真 ： 通过 Smoothed Particle Hydrodynamics(SPH) 方法来对流 体细节进行物理仿真， 首先建立流体仿真物理模型， 包括流体求解区域和边界条件， 其中边 界条件为PLSM追踪到的流体界面。

5、 ； 然后设定SPH初始条件， 确立粒子生成位置， 初始化粒子 属性信息 ； 再然后， 将控制流体运动的 N-S 方程离散化， 用离散了的控制方程近似表达粒子 属性变量及其导数形式 ； 接下来， 根据时间步长积分形式 t 开始累加， 更新求解区域的粒 子位置、 速度、 密度信息， 计算粒子属性在一个时间步长内的变化量 ； 紧接着 t+1， 计算下 一时刻的速度变化量和新的位置信息， 更新 SPH 粒子系统 ； 最后， 用新的速度和位置信息计 算瞬时变量 ； 重复以上各步骤， 直到指定结束时刻和计算达到收敛条件， 仿真结束 ； 以此模 拟流体运动过程中产生的浪花、 水珠细节 ； 步骤 (4)、 。

6、集成仿真场景、 进行算法耦合 ； 把步骤 (1)、 (2)、 (3) 中 LBM 流体和 SPH 仿真 的流体集成到同一个场景中， 利用耦合算法进行 LBM 和 SPH 算法的参数传递， 包括密度、 速 度、 位置 ； 根据耦合算法控制流体仿真过程中 SPH 粒子的生成时机 ； 步骤 (5)、 运用快速行进算法 (Fast Marching Mehod， FMM) 对 PLSM 追踪出来的网格进 行遍历， 并用三角面片填充， SPH 渲染的流体细节直接用 OpenGL 绘制成粒子球。 2. 根据权利要求 1 所述的基于欧拉 - 拉格朗日耦合方法的流体仿真方法， 其特 征在于 ： 步骤 (1) 。

7、中所述的大规模流体的物理仿真方法， 该方法使用空间场内的均匀 网格表示流场， 利用 DnQm 网格模型表示流体在网格格点处的属性， 将纳维斯托克斯 (Navier-Stokes， N-S) 方程转化为只与时间相关的离散形式的 LBGK 方程形式求解， 具体 过程为假设流体运动具有平衡态， 分子之间的碰撞过程会促使流体分布函数趋近于平衡 态， 在流体趋于平衡态的过程中， 从空间、 时间和速度的角度将连续矢量离散成沿各个特 定的方向的标量值， 将连续的流体运动看作为两个离散步骤碰撞 (Collision) 和迁移 (Stream)， 通过离散时间 t 驱动流体在网格格点处的碰撞或迁移过程以实现流体。

8、模拟过 程。 3.根据权利要求1所述的基于欧拉-拉格朗日耦合方法的流体仿真方法， 其特征在于 ： 步骤 (3) 中所述的流体细节仿真， 该方法使用一组离散粒子表示流体， 利用平滑核函数近 似表示粒子的物理参数， 将纳维斯托克斯 (Navier-Stokes， N-S) 方程转化为只与时间相 关的离散形式的微分方程， 再利用积分方法求解 ； 描述流体运动的控制方程可以由一系列 拉格朗日型偏微分方程表示， 这些方程遵循能量守恒定律 ； N-S 方程可用时间全导数形式 表示如下 ： 权 利 要 求 书 CN 104268943 A 2 2/2 页 3 动量方程为 ： 能量方程 ： 其中， ， 表示坐。

9、标方向， 表示流体密度， t 表示时间， u 为流体速度， x 为空间位 置， 为总应力张量， 对于给定的粒子 i， 宏观密度可根据 SPH 方法用支持域内粒子密度累 加得到 ： 其中， 光滑核函数 Wij为 ： Wij W(|xi-xj|,h) W(Rij,h) Rij表示粒子 i 和 j 之间的相对距离 ； rij/h， rij是两个粒子的绝对距离 ； h 是光滑长度， 直接影响 SPH 方法的求解精度。 权 利 要 求 书 CN 104268943 A 3 1/6 页 4 一种基于欧拉 - 拉格朗日耦合方法的流体仿真方法 技术领域 0001 本发明涉及一种基于欧拉 - 拉格朗日耦合方法的。

10、流体仿真技术。 背景技术 0002 目前， 流体仿真朝着多学科， 多角度， 多平台的方向飞速发展， 但是碍于现有硬件 水平的限制和仿真算法的局限性， 对流体模型进行精确、 迅速、 逼真的还原还存在诸多困 难。近年来， 满足工程和科研要求的真实流体仿真备受国内外学者和科研机构的关注。基 于物理的真实感流体仿真研究对于从事流体相关的科研工作者和计算机图形学的研究者 来说无疑面临着巨大的挑战， 也具有巨大的吸引力。 0003 基于物理的流体仿真方法主要分为两种 ： 基于网格思想的欧拉方法和基于粒子思 想的拉格朗日方法。拉格朗日方法中 SPH 方法多用于流体的细节彰显和模拟小规模流体， 由于流体细节的。

11、精细刻画、 流体仿真的速度和流体仿真规模这三者具有不可调和的矛盾， 在进行精细流体细节渲染的同时很难保证实现大规模流体仿真。 0004 LBM 方法是基于欧拉思想的典型网格方法， 具有天然的并行性， 且边界条件容易处 理， 适合于大尺度、 低精度、 规则边界的流体仿真， 如山洪地质灾害、 泥石流、 大规模水面等。 通过精细网格划分或非规则网格划分可以拓展 LBM 流体仿真范围， 但是由此带来的稳定性 损失和效率下降将使 LBM 失去其固有优势。因此， 与其他方法相结合， 发挥 LBM 方法天然优 势的混合方法正逐渐引起越来越多学者的关注。 0005 为了解决上述问题， 本发明提出一种基于欧拉 。

12、- 拉格朗日耦合方法的流体仿真技 术， 以实现兼顾速度、 规模和细节的流体仿真为目标， 围绕基于物理的流体仿真方法展开研 究， 首先用LBM方法进行大规模的流场建模， 并利用PLSM方法提取和追踪流体表面， 然后在 已有 LBM 仿真流体的基础上， 应用 SPH 方法刻画流体诸如水花、 泡沫的细节， 通过设计耦合 算法将 SPH 流体和 LBM 流体集成到同一场景中， 实现了兼顾细节的大规模交互场景的实时 逼真绘制。 发明内容 0006 本发明解决的技术问题是 ： 克服了现有的流体仿真方法的效率问题， 提供了一种 基于欧拉 - 拉格朗日耦合方法的流体仿真技术， 解决了流体仿真规模与细节的矛盾，。

13、 满足 了对流体运动过程进行实时模拟的需要。 0007 本发明采用的技术方案为 ： 一种基于欧拉 - 拉格朗日耦合方法的流体仿真方法， 包括以下四个步骤 ： 0008 步骤(1)、 大规模流体的物理仿真 ： 利用求解N-S方程的LBGK控制方程， 采用DnQm 多维离散网格模型， 进行大规模流场建模， 这是流体仿真的主体部分。 0009 步骤(2)、 流体自由表面追踪 ： 应用PLSM方法构建流体自由表面。 输出步骤(1)中 得到的流体网格， 用PLSM方法标记流体运动界面在LBM网格内部的位置， 追踪流体界面， 以 此作为流体的近似表面。 说 明 书 CN 104268943 A 4 2/6。

14、 页 5 0010 步骤 (3)、 流体细节仿真 ： 通过 SPH 方法来对流体细节进行物理仿真， 模拟流体运 动过程中产生的诸如浪花、 水珠等细节。 0011 步骤 (4)、 集成仿真场景、 进行算法耦合。把步骤 (1)、 (2)、 (3) 中 LBM 流体和 SPH 仿真的流体集成到同一个场景中， 利用耦合算法进行 LBM 和 SPH 算法的参数传递， 例如密 度、 速度、 位置等。根据耦合算法控制流体仿真过程中 SPH 粒子的生成时机。 0012 本发明的原理在于 ： 0013 (1) 通过 LBM 方法模拟大规模流体运动， 利用空间网格对流体的运动状态进行计 算， 以此发挥 LBM 方。

15、法适合模拟大规模、 低精度流体的特点， 并能在保证流体仿真规模的同 时保证速度优势。 0014 (2) 为了快速精确的追踪流体的表面， 本发明利用 PLSM 方法通过距离函数来近似 代表流体液面的空间位置。该方法通过使用 PLSM 追踪的流体表面在网格节点中的位置， 然 后分析 PLSM 方法进行液面追踪时造成的耗散误差， 通过补偿粒子的方法弥补流体液面处 被错误平滑掉的流体信息。接着用快速行进算法得到流体液面， 采用 OpenGL 对流体液面进 行绘制， 保持了边界的锐度并且可以直接以此作为大规模流体表面的近似代表。这样得到 的流体表面， 保持了良好的真实性， 并具有非常高的效率。 0015。

16、 (3)为了对流体进行逼真的细节刻画， 本发明提出了一种基于LBM和SPH的耦合流 体仿真方法。 该方法在LBM流体的基础上， 通过SPH方法模拟诸如浪花、 泡沫等流体细节， 并 设计 LBM-SPH 的耦合算法控制细节生成位置和生成时机， 达到良好的流体真实视觉表现。 0016 本发明与现有技术相比的优点在于 ： 0017 1、 本发明提出的高效数据组织方式， 流体仿真过程中的流体参数可在各个方法间 无障碍的传递， 保证了算法的鲁棒性。 0018 2、 高效率的网格方法， 本发明利用 PLSM 对流体液面进行追踪， 由于 LBM 和 PLSM 同 属网格方法， LBM 的空间网格可以被 PL。

17、SM 方法重用， 因此有效地减少了计算时间。 0019 3、 结合不同流体仿真方法的优势， 既发挥了 LBM 方法速度快、 可仿真大规模流场 的优势， 又发挥了 SPH 方法擅长刻画流体细节的优势。 0020 4、 快速的液面追踪和绘制， 利用 PLSM 方法可对流体液面进行追踪和补偿， 提高了 算法效率。 附图说明 0021 图 1 基于欧拉 - 拉格朗日耦合方法的流体仿真流程图 ； 0022 图 2 二维 LBM 示意图 ； 0023 图 3 算法划分示意图 ； 0024 图 4 修正网格示意图 ； 0025 图 5 PLSM 与 LBM 耦合后仿真的水滴下落过程 ； 0026 图 6 仿。

18、真结果与真实照片对比 ； 0027 图 7 渲染后耦合方法的仿真结果。 具体实施方式 0028 图 1 给出了基于欧拉 - 拉格朗日耦合方法的流体仿真的总体处理流程， 下面结合 说 明 书 CN 104268943 A 5 3/6 页 6 其他附图及具体实施方式进一步说明本发明。 0029 本发明提供一种基于欧拉 - 拉格朗日耦合方法的流体仿真技术， 具体实施为 LBM-SPH 耦合的流体仿真， 主要步骤介绍如下 ： 0030 1. 大规模流体的物理仿真 0031 本发明采用 LBM 方法来模拟大规模流体的运动， LBM 由上世纪 70 年代提出并发展 的格子气自动机(Lattice Gas 。

19、Automata， LGA)发展而来， 因此可把LBM看做是描述流体运 动的一个物理模型。 LBM方法改善了传统的LGA方法中诸如噪声、 指数复杂性等问题。 在LBM 中， 用粒子分布函数替代 LGA 中粒子本身进行演化， 直接采用格子 Boltzmann 方程 (LBE) 作 为演化方程， 根据其分布函数计算流体的密度和速度。 0032 连续的 Boltzmann 方程表示为 ： 0033 0034 是粒子速度， (f) 是表示碰撞过程的碰撞因子， 是哈密顿算子。f 是连续方 程的平衡态分布函数， 平衡态分布函数保证了流体局部的质量与动量守恒。f 线性依赖于 流体密度 和速度 u。根据上式，。

20、 可通过求数值解的方法模拟流体的宏观运动。然而碰撞 项 (f) 并非分布函数的线性方程， 只与分子作用力有关， 因此求解起来具有难度， 这使得 Boltzmann方程在实际应用中受到了很大限制。 改进了的Bhatnager-Gross-Krook(BGK)模 型假设流体运动具有平衡态， 分子之间的碰撞过程会促使分布函数 f 趋近于平衡态。由于 碰撞导致的变化量和 f 与平衡态的差值成正比， 即 ： 0035 0036 c为松弛时间， 1/c称之为平均碰撞频率。 0037 从宏观运动的角度， 粒子速度 和分布函数 f 也相应的被离散到 n 维空间， 这意 味着在 t 时刻， 某个格点 x 处将有。

21、沿各个离散方向的分布函数。格子 Boltzmann 方程 (LBE) 是 BGK 方程的特殊离散格式， 从空间、 时间和速度的角度将连续矢量离散成沿各个特定的 方向的标量值。以二维流场为例， LBM 模型示意图如图 2 所示 0038 此模型中， 流体运动被看作为两个步骤碰撞 (Collision) 和迁移 (Stream)。 不含外力项的 LBGK 方程定义为 ： 0039 0040 fq是粒子速度分布函数， 其中 q 代表离散方向 (0， 1， 2， 3， q)。是粒子的平 衡态分布函数。 此方程可简化为两部分， 将碰撞过程和迁移过程分开表示， 这样的好处在于 方便计算， 并使碰撞项完全线。

22、性化， 用单松弛时间逼近近似表现碰撞过程， 使得系统趋于平 衡。 0041 碰撞 ： 0042 迁移 ： 说 明 书 CN 104268943 A 6 4/6 页 7 0043 在上式中，表示碰撞后的分布函数。平衡态分布函数为 ： 0044 0045 q表示粒子速度方向上的权重， eq为离散速度方向， 是宏观密度。c 是格子声 速， 其值为 x/t，与声速有关， x 和 t 分别代表时间步长和网格步长。在 均匀网格计算中通常取 t x。在宏观上， LBM 对粒子分部函数进行统计， 计算出宏观 压力和宏观流速等变量 ： 0046 0047 0048 b 表示离散速度数量。 0049 此外， 在 。

23、LBM 模型中时空并不是完全独立的， 而是通过粒子的离散速度 e 联系起 来 ： 0050 x eqt 0051 粒子在一个时间步长内由一个节点运动到另一个节点， 进行 “迁移” 过程， 同时其 他节点运动到此节点的物理量通过 “碰撞” 过程进行更新。 0052 应 用 于 流 体 仿 真 的 LBM 方 法 有 许 多 离 散 速 度 模 型， Lattice Bhatnager-Gross-Krook(LBGK) 模型是迄今为止应用最为广泛的模型， 本发明使用应用广 泛的 DnQm 模型。其中 n 代表空间维数， m 代表离散速度个数。不同模型的平衡态分布函数 中权系数 和格子声速 c 也。

24、不同， 本发明应用 D3Q15 模型进行大规模流体建模， 即在三维 空间内， 每个格子节点有 15 个离散的子速度。 0053 2. 基于 PLSM 的流体液面追踪 0054 a) 流体分区域仿真算法分配如下 ： 0055 1) 用 LBM 方法模拟流体主体， 发挥 LBM 方法善于计算大尺度流体的优点 ； 0056 2)在需要精细刻画和细节渲染的液面附近结合SPH和LBM方法， 渲染流体细节， 增 强显示效果 ； 0057 3) 用水平集方法对液面进行跟踪绘制流体表面， 保证流体的真实感绘制。 0058 如图 3 所示， 在进行大规模流体仿真的区域， 相对平稳的流体占据绝大部分空间， 因此采。

25、用无需刻画细节且具有速度优势的 LBM 方法建模 ； Coupling Band 表示 PLSM 作用 区域， 在此区域内应用 LBM 和 PLSM 偶和算法进行流体液面跟踪 ； 同时在液面和空气中利用 SPH 生成粒子表示液面诸如浪花、 飞沫、 气泡等细节增强真实感。 0059 b) 网格优化 0060 流体仿真根据区域的不同， 对精度的要求也不同， LBM 求解区域对精度要求较低， 速度较快， 可用粗糙网格求解， PLSM 和 SPH 求解区域对精度要求较高， 可应用精细网格求 解。基于这一思想， 本文对网格进行重新划分， 修正网格尺寸， 在保证算法效率的同时尽可 能提高算法精确度。重新划。

26、分后的网格示意图如图 4 图示。 说 明 书 CN 104268943 A 7 5/6 页 8 0061 粗细格点之间数据传递遵循以下方程 ： 0062 0063 0064 上式分别表示表示精细网格的分布函数和粗糙网格的分布函数。 其中，和 表示粗细网格的平衡态分布函数， c、 f分别表示粗糙网格和精细网格的松弛时间， 具有 如下关系 ： 0065 0066 在构建 PLS 区域的同时设置粒子和网格的耦合区域， 通过耦合算法实现粒子仿真 流体与网格仿真流体区域之间的衔接。在液面附近设置过渡区域， 区域以上的流体 ( 主要 是浪花 ) 完全采用 SPH 方法 ； 区域以下的流体主要采用 LBM 。

27、方法进行流体仿真， 由于冲撞、 挤压等物理过程产生的运动至流体内部的气泡采用 SPH 方法渲染 ； 在耦合区域的泡沫粒子 采用 SPH 方法。浪花、 气泡和泡沫这三种粒子通过不同的图形学绘制手段区分。根据粒子 类型选择参数不同的光滑核函数， 以体现不同粒子属性 ； 渲染流体粒子 ； 对 LBM 网格格点属 性进行插值， 交换格点和粒子属性， 通过耦合方程实现两种方法的数据传递。 0067 3.LBM 方法与 PLSM 的耦合 0068 PLSM 方法通过追踪流体运动构建液面网格， 而流体运动由 LBM 方法计算得到， 因 此需要构建耦合算法， 使 LBM 方法求得的流体参数准确传递给 PLSM。

28、 方法， 进而对液面进行 追踪。在本文中， 流体仿真算法 LBM 和液面跟踪技术 PLSM 同属欧拉网格方法， 因此在参数 传递、 模型合并、 方法耦合等方面具有优势。 0069 为实现模型结合后算法的正确耦合， 需要对 LBM 的算法形式加以变换， PLSM 的网 格尺寸与 LBM 网格尺寸相同， 这使得 LBM 速度场能够直接传递给 PLSM 的距离函数。在本发 明中， 应用半拉格朗日对流项传递速度参数。半拉格朗日方法的基本思想是 ： 已知给定的 Level Set方程和速度场u， 格点x处的速度每t时间进行一次更新， 液面位置通过(x) 进行跟踪 ： 0070 (x) (x-u(x)t)。

29、 0071 在液面跟踪过程中， 需要对网格进行分类 ： 节点处没有速度值的网格标记为空网 格 ( 空气网格 ) ； 节点处全部具有速度信息的网格标记为液体网格 ； 部分节点具有速度信 息， 部分节点为空的网格标记为液面网格。其中， 流体网格用 LBM 算法进行计算， 空网格所 有节点处属性置零， 因为只需对液面网格的算法进行改进即可。根据距离函数标记的网格 类型， 气体网格标记为 -1， 液体网格标记为 +1， 因此同时具有气体网格和液体网格邻居的 网格必为液面网格。液面网格的演进过程遵照 LBM 演进规则， 碰撞过程按照如下方程进行 修正 ： 0072 说 明 书 CN 104268943 。

30、A 8 6/6 页 9 0073 x 表示液面网格， i 是液面网格相邻的气体网格方向， 与 i 方向相反， A是密度， u 是 x 处的速度。 0074 算法表示的步骤为 ： 0075 1) 流体部分依据 LBM 算法模型初始化流场区域， 根据 LBGK 方程计算分布函数 ； 0076 2) 时间步长 t 开始计数， 各节点依次进行碰撞和演进两个过程， 平衡态分布函 数及分布函数同时更新 ； 0077 3) 根据流场内网格节点属性判断节点类型， 液体节点标记为 1， 空节点标记为 0， 介于液体节点和空节点之间的节点标记为液面节点 ； 0078 4) 得到液面位置信息， 修正碰撞后的分布函数。

31、 ； 0079 5) 根据修正后的分布函数计算液面节点处的速度信息， 并传递到 PLS 算法中， 重 新初始化距离函数 ； 0080 6) 根据修正的距离函数跟踪液面新位置信息 ； 0081 7) 根据得到的液面位置信息判断当前格点 x 的类型， 在下一个时间步长内反馈给 LBM 算法 ； 0082 8)LBM 算法根据当前格点类型重新初始化整个流场网格， 重新初始化分布函数 ； 0083 9) 时间步长 t+1， 进行新一轮迭代。 0084 耦合后的模型示意图如图 5 所示。 0085 4.LBM 方法与 SPH 方法的耦合 0086 网格划分过程中， 远离流体表面的区域划分为稀疏网格以节省。

32、资源， 提高运算效 率， 在液面耦合区域 ( 即 PLS 区域 ) 为了提高仿真的精确性， 实现细节的逼真刻画， 需要划 分精细的网格。区域内 LBM 和 SPH 仿真的流体之间的速度和密度传导通过以下方程实现 ： 0087 0088 0089 Xc 表示耦合位置， 通过以上两个方程， 确保了 SPH 粒子可以准确得到初始化属性 信息 ( 速度、 密度 )， 并通过局部求解将运算结果反馈给 LBM 算法， 实现了全局数值传递， 从 而达到了两种算法的耦合。 0090 5. 液面绘制 0091 最后， 运用快速行进算法(Fast Marching Mehod， FMM)对PLSM追踪出来的网格进。

33、 行遍历， 并用三角面片填充， SPH 渲染的流体细节直接用 OpenGL 绘制成粒子球。 0092 本发明未详细阐述的技术内容属于本领域技术人员的公知技术。 0093 尽管上面对本发明说明性的具体实施方式进行了描述， 以便于本技术领的技术人 员理解本发明， 但应该清楚， 本发明不限于具体实施方式的范围， 对本技术领域的普通技术 人员来讲， 只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内， 这些变 化是显而易见的， 一切利用本发明构思的发明创造均在保护之列。 说 明 书 CN 104268943 A 9 1/3 页 10 图 1 图 2 说 明 书 附 图 CN 104268943 A 10 2/3 页 11 图 3 图 4 图 5 说 明 书 附 图 CN 104268943 A 11 3/3 页 12 图 6 图 7 说 明 书 附 图 CN 104268943 A 12 。