大变形扑翼动态模型设计及气动网格生成技术.pdfVIP

第 1期 朱志东，等 ：大变形扑翼动态模型设计及气动 网格生成技术 3 内翼极 限位 置 时 ，外翼 相 对 内翼 折弯 (￡)一 q- 。 t (4) 内翼不在极 限位置 时 ，外 翼跟 随 内翼 一 起扑动 ，即 内翼 下扑 时 ： (￡)一 (￡) (5) 内翼 上提 时 ： (￡)一J8(￡)-+-∞… t… (6) B 对 于外翼 的扭 转角 ，观 察 录像 可 以发 现扭转 角沿 着 展 fa)二维 (b)三维 向是变化 的 ，翼 梢 部分 扭 转 角最 大 ，不妨 假 设 扭 转 角 沿 展 向是 线性 的。假设 扭转 角 只发 生在外端机 翼部 分 ，由于 内 图3 位置系数和节点关系 段翼 扭转 角很 小 ，这个 假 设 是 合 理 的 。 同时注意 到 ，机 翼 z p 一 ∑ zY3DY (11) 扭转 与机 翼 扑动 存在 相 位 差 ，机翼 下 扑 时外 翼迎 角减 小 ， (5)计 算下一个 时 间步 ，网格 点 P的新 位置 。 机翼 上提 时外 翼迎 角 增 大 。那 么某 个 截 面扭 转 角 的规 律 如 下 所 示 ： 3．2 扑 翼模型 的程 序 实现 a(￡)一 一 X 垃… × sin( h) (7) 程序 实现机 翼运 动 模 型主要 思想 是 区分 内外 翼 和 整 其 中 为该截 面 与 内外 翼 连接处 距离 的无 量纲 表示 形式 。 个运 动周 期 的 4个 阶段 ，通过 网格 点到 根弦 线 的距离来 区 这样 通 过变 化各 个参 数 (如 ~Odo ， ， … ， … ，a～ ， 分 内外翼 。对 于 内外 翼 和 4个 阶段 分别在 各 自的运 动 规 )就 可 以较好 地模 拟 鸟类 翅膀 的各种 扑动 规律 。 律 下计算 。运 动模 块 的流程 图见 图 4。 3 大 变 形 扑 翼 动 网格 生成 l读取 3．1 基于 Delaunay图 映射 的网格 变形 原理 由于 扑 翼 MAV 的运 动 幅 度 很 大 ，在 大 变 形 的情 况 N 啦匦卜呼 ： l 下 ，怎么生成 扑翼 MAV 气 动计算 所需 的合 格 网格成 为研 究 中的重点和难 点 。常用 的方法有弹簧模 型 的方法和基 口塑]————— — 于 Deiaunay图映射 的方 法 。 如 果 直接使 用 弹簧 模 型 的方 法 计 算 大 变形 扑 翼 动 网 格 ，由于 网格 数 量多 ，变形 较 大 ，在 计算 时会 发生 网格 重叠 图 4 运 动模 块 流程 图 现 象 ，网格质量 下 降，甚 至 网格不能计算 。如果采用基于 如 图 4所示 ，在 扑 翼 飞 行 的一 个 周期 内 的某 个 时 刻 Delaunay图映射 的方 法 ，对 稀 疏 的背 景 网格 采 用 弹簧 模 时 ，程 序读取 当前 网格 中某 点 的坐 标 值 z， ， ，计算 网格 型 的计 算 方法 ，则 背景 网格 质 量较 好 ，无 需迭 代 计 算