EI、Scopus 收录
中文核心期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

预出版

预出版栏目展示本刊经同行评议确定正式录用的文章,这些文章目前处在编校过程,尚未确定卷期及页码,但可以根据DOI进行引用。
显示方式:
Gappy POD方法重构湍流数据的研究
李天一, BuzzicottiMichele, BiferaleLuca, 万敏平, 陈十一
Gappy POD 是一种基于本征正交分解(proper orthogonal decomposition, POD)的数据重构方法. 本文研究了gappy POD在湍流数据重构中的应用, 主要关注了以下两个因素的影响: 第一, 数据本身的复杂程度, 即构成流场的POD模态数量; 第二, 破损区域的面积大小和几何形状. 考虑到上述因素, 本文重新严格地表述了gappy POD的重构过程, 并推导出gappy POD重构误差的公式. 论文选取旋转湍流数据为案例进行了gappy POD重构的研究, 并解释了构成gappy POD重构误差的两个部分. 第一部分来自流场POD展开的截断误差, 该截断误差会被POD基函数在已知点上的值组成的矩阵的最小特征值放大. 这部分误差主要取决于流场的复杂程度, 当流场复杂程度较低时, 相应误差随采用的POD模态数目增大而减小. 当流场复杂程度较高时, 很小的POD截断误差也会导致很大的重构误差, 此时需要采用流场所有的POD模态进行重构以消除截断误差. 重构误差的第二部分来自POD基函数在已知点上的值组成的矩阵的非列满秩性, 它主要取决于破损区域的面积大小和几何形状. 破损区域的面积越大, 或者破损面积相同时, 破损区域内信息所包含的相关性越大, 第二部分的重构误差越大.
, doi: 10.6052/0459-1879-21-464
基于尾流时程目标识别的流场参数选择研究
战庆亮, 葛耀君, 白春锦
浸入流场中的固体壁面会形成高度复杂且具有一定特征的尾流流场, 利用尾流所包含的信息对物体的外形特征进行识别具有重要的应用价值. 然而, 在较高雷诺数情况下尾流流场形态及其时序特征复杂, 难以通过传统的数学物理方法对流场信号进行特征的识别与提取. 本文提出了基于尾流时程数据深度学习的流场特征提取与分析方法, 实现了基于一点的物理量时程进行流场中物体外形的识别; 同时, 对流场中不同物理参数时程的识别精度与识别结果进行分析与研究, 得到适用于目标识别的最优物理量参数. 通过对圆柱和方柱的尾流数据研究结果表明, 本文提出的基于卷积神经网络的模型具有好的训练收敛性和高的预测精度, 能够识别并提取得到时程数据中包含的流场特征, 采用流场横向速度时程作为物体外形识别信号的模型准确率高. 证明了本方法用于浸入流场中物体外形识别的可行性, 是一种目标识别的高精度方法.
, doi: 10.6052/0459-1879-21-332
变后掠翼的参变气动弹性建模与分析
张立启, 岳承宇, 赵永辉
参变气动弹性建模是可变后掠翼气动弹性研究中的难点之一. 当地建模技术是一种构造线性参变 (LPV) 模型的实用方法, 但一直缺乏有效途径来解决当地气动弹性模型的不一致问题. 气动弹性模型的不一致性体现在当地结构动力学模型和非定常气动力模型随参数变化的不连续性. 本文提出了一种自下而上的方法, 对变后掠翼不一致的当地气动弹性模型进行了一致性处理. 首先, 采用匈牙利算法跟踪结构模态并按模态分支进行排序, 使得匹配后的模态能够保证结构动力学模型的一致性; 其次, 对有理函数拟合表达式中的系数矩阵进行缩放处理, 解决了空气动力系数矩阵的不一致问题. 采取上述两项措施后, 最终生成了一致的当地状态空间气动弹性模型. 这样, 对一致的状态空间模型进行插值, 就可快速生成任意后掠角下的气动弹性模型, 使得系统的稳定性分析和慢参变响应计算得以高效进行. 数值仿真结果验证了一致性处理的必要性:未经一致性处理的原始模型经插值后得到的模型会出现严重的模型误差. 本文为变后掠翼的参变气动弹性系统建模提供了一种实用、准确和高效的建模方法.
, doi: 10.6052/0459-1879-21-275
高马赫数燃烧强化的激波风洞试验研究
张旭, 张启帆, 岳连捷, 孟东东, 罗苇航, 于江鹏, 张晓源, 李进平, 陈宏, 李飞
基于中国科学院力学研究所的JF-24激波风洞, 通过开展高马赫数超燃冲压发动机的直连试验, 研究了高马赫数燃烧的强化方法以及燃料类型对燃烧的影响. 试验段是采用凹腔结构的圆截面燃烧室, 喷孔布置在隔离段, 燃料分别是氢气和乙烯, 当量比均为0.7. 燃料喷注分别采用无支板和小支板两种构型, 后者部分喷孔位于小支板顶部. 两种构型均设置了流向近距双排喷孔, 可分别进行单环和双环喷注. 试验结果论证了飞行马赫数10.0条件下氢气和乙烯在超高速气流中的稳定燃烧性能. 并且, 相比于单环喷注, 双环喷注以及补充小支板可以强化燃烧. 推测其原因是双环射流和激波/分离结构的近距离交互作用很可能改善掺混, 而补充小支板顶部喷注还能利用更多空气组织掺混. 在同样采用双环耦合小支板顶部喷注的强化措施下, 氢气与乙烯燃烧效率接近, 但氢推力性能更优. 这是因为较高热值氢的释热更多. 此外, 试验还证明了在当前来流条件下, 释热受控于掺混, 且高温离解效应限制释热上限. 这是由于释热降低流速且提高静温, 使高温离解的吸热效应更加显著.
, doi: 10.6052/0459-1879-21-348
基于绝对节点坐标法的平面梁有限变形下变形重构
吴懋琦, 谭述君, 高飞雄
现有的对有限变形条件下柔性结构变形重构的研究往往单纯基于曲率与应变间的几何关系, 同时忽略了被测体的纵向变形及其与弯曲变形的耦合效应. 为得到一种更加精确且能借助现有的力学工具进行应用方向扩展的变形重构方法, 以平面梁为对象, 借鉴变形重构逆有限元法的思想, 将平面梁的变形重构问题视作一类最优化问题. 首先, 通过引入绝对节点坐标法(absolute nodal coordinate formulation, ANCF)对柔性结构大变形下非线性的平面梁应变−位移关系进行精确描述, 构造了一种逆梯度缩减ANCF平面索梁单元. 然后, 对此逆ANCF单元进行改进, 在简化节点自由度的同时通过引入罚函数确保单元节点处的曲率连续性, 既保证了本问题的适定性, 也提升了最终解的精确性. 最后, 基于该单元利用Newton法构造了平面梁有限变形下变形重构问题的两种求解算法, 即逐单元算法和多单元整体算法, 以实现不同需求下的稳定求解. 数值仿真结果表明, 本方法在大变形条件下的变形重构误差小于1%, 而且在测点较少的情况下依然保持较高的精度, 同时验证了本方法的收敛性与计算效率.
, doi: 10.6052/0459-1879-21-338
周期性非均匀介质中气相爆轰波演变模式研究
陈达, 宁建国, 李健
气相爆轰波在周期性非均匀介质中的起爆, 稳态传播和失效机制都极为复杂, 很多物理机制尚不明确, 是当前爆轰物理领域研究的热点和难点. 本文使用反应欧拉方程和两步化学反应模型对爆轰波在非均匀介质中的传播机理进行了数值模拟研究, 非均匀性由横向周期性分布的温度扰动体现, 重点分析不同波长、不同幅度的温度扰动对波阵面波系结构的影响. 计算结果表明, ZND爆轰波在温度扰动下向胞格爆轰波的转变主要受制于两种竞争性因素: 一是爆轰波内在的不稳定性; 二是温度扰动的波长和幅度, 前者是内因, 后者是外因. 温度扰动的存在抑制横波的发展, 延迟了ZND爆轰波向胞格爆轰波的演化, 并且内在不稳定性的增加可以减慢这种延迟现象. 这说明, 温度扰动可以在一定的范围内抑制胞格不稳定性的发展, 但是不能够终止这一过程. 温度的不连续性使得爆轰波阵面更为扭曲, 并在横波附近存在较弱的三波点结构, 即温度扰动可增加爆轰波固有的不稳定性, 改变爆轰波阵面的传播机理. 幅值较大的人工温度扰动可抑制爆轰波的传播和爆轰波自身的不稳定性. 爆轰波阵面胞格结构的形成取决于温度扰动与其自身的不稳定性.
, doi: 10.6052/0459-1879-21-069
磁力耦合道路能量收集设计与动力学分析
邹鸿翔, 郭丁华, 甘崇早, 唐曙光, 袁俊, 魏克湘, 张文明
通过在交通环境布置无线传感器等小型机电系统, 实现交通状况监测、交通系统管控、交通设施健康状态监测等, 可以使交通系统更加安全、有序、高效地运行. 但是, 如何为这些广泛分布的小型机电系统供能?本文提出了一种磁力耦合道路能量收集设计, 用以收集车辆滚压能量并转换成电能. 通过磁力耦合进行无接触能量传递, 减小了装置受到的冲击并使得装置具有良好密封性, 从而提升装置的鲁棒性. 通过升频齿轮机构、棘轮机构将车辆滚压激励转换为高速单向旋转, 并且通过换向齿轮机构能够继续收集复位弹性势能, 提高了收集装置的输出功率. 基于磁力耦合道路能量收集系统的工作原理建立了机电耦合动力学模型. 数值仿真探究了减速带限位距离和复位弹簧刚度等关键设计参数对能量采集系统动力学和电学性能的影响. 能量采集系统在车速为50 km/h时最大输出电压为76.28 V, 最大功率为59.94 W. 磁力耦合道路能量收集装置可以成为未来智慧交通系统的重要组成部分, 俘获交通环境能量为交通环境中小型机电系统提供可持续的绿色无碳电力.
, doi: 10.6052/0459-1879-21-374
亚毫米球体撞击液滴过程实验研究
左子文, 蒋鹏, 王军锋, 王林, 霍元平
球体撞击液面是自然界和工业过程中的普遍现象. 目前相关研究主要关注毫米级及更大尺寸球体撞击水平液面. 对于亚毫米球体撞击过程的动力学特性及液滴弯曲液面对撞击行为的影响仍需深入研究. 本研究基于高速显微数码摄像技术, 开展了不同速度亚毫米球体撞击液滴不同位置的实验研究. 撞击行为呈现振荡和浸入两种模式. 由于液滴弯曲液面的存在, 撞击现象与水平液面不同, 润湿过程中三相接触线 (three phase contact line, TPCL)固定点方位角与撞击角度线性正相关, 非轴对称的空穴在TPCL固定点较高一侧率先形成且曲率半径较大. 揭示了撞击过程中球体主导力变化和能量转化机制. 分析了撞击速度和撞击角度对撞击行为的影响, 并给出了撞击模式图. 结果表明: 冲击阶段形状阻力主导撞击行为, 球体动能损耗量与撞击速度正相关. 空穴发展阶段由表面张力主导, 球体动能转化为维持空穴形状的表面能. 振荡模式空穴长度与韦伯数We正相关, 空穴发展速度差异较小. 根据量纲分析及实验结果拟合得到临界浸入韦伯数Wecr与撞击角度α关系式$We_{cr}^{1/2} $ = α/40.
, doi: 10.6052/0459-1879-21-351
复杂流动中的智能颗粒游动策略研究进展
邱敬然, 赵立豪
智能颗粒定义为可以根据环境的变化而主动调整自身在流场中运动的颗粒, 一般用于描述微小的游动物体, 如微生物、浮游生物和微纳机器人. 由于颗粒运动学特性和流场环境的复杂性, 探索智能颗粒的运动策略是一项具有现实意义与挑战性的研究. 近期强化学习算法被应用于智能颗粒的运动策略研究, 并取得一定进展. 本文将讨论强化学习在智能颗粒研究中的应用, 并介绍浮游生物运动策略的相关研究进展, 包括适用于海洋浮游生物的游动颗粒模型, 以及强化学习的基本原理和此类问题的研究框架. 海洋浮游生物的昼夜垂直迁移对其生存和繁衍至关重要, 生物学研究表明某些浮游生物可感知周围的流场信息, 但能否利用这些信息提高垂直方向游动效率仍是一个未知的问题. 基于这一背景, 相关研究考虑了重力沉降和颗粒形状对浮游生物垂直游动策略的影响. 研究发现细长的颗粒能够更加高效地垂直向上运动, 而重力沉降则导致智能游动策略发生较大变化. 在此基础上, 后续研究进一步考虑了局部流体信号的作用, 并讨论了基于局部信息在全局坐标系中定向运动的可能性. 当颗粒只能感受到局部流体信号时, 必须破坏动力学系统的对称性, 否则颗粒无法学到有效的垂向游动策略. 此外, 研究还发现强化学习能够利用流体信号背后的机制, 得到颗粒在二维定常流动和三维非定常湍流中高效垂直运动的策略. 这些智能游动策略依赖于本质的物理规律, 因此这些策略在更加复杂或真实的流动中也可能有优秀的表现.
, doi: 10.6052/0459-1879-21-402
增材制造微结构演化及疲劳分散性计算
易敏, 常珂, 梁晨光, 周留成, 杨阳祎玮, 易新, 胥柏香
为了预测增材制造中工艺参数−微结构−力学性能之间的关联规律, 提出了集成离散元、相场模拟、晶体塑性有限元和极值概率理论的计算方法, 揭示了激光扫描速度对微结构演化、屈服应力和疲劳分散性的影响. 首先, 采用离散元实现了重力作用下粉床在已凝固层表面上的逐层铺设; 其次, 通过热−熔体−微结构耦合的非等温相场模拟, 获得了熔体、气孔、晶界、晶粒取向等的时空演化以及最终形成的多晶微结构; 然后, 应用晶体塑性有限元计算了增材制造多晶微结构的宏观力学响应, 并得到表征疲劳裂纹萌生驱动力的疲劳指示参数(FIP); 最后, 采用极值概率理论分析了增材制造多晶微结构的FIP极值分布规律及疲劳分散性. 以316L不锈钢选区激光熔化增材制造为例的计算结果表明: 增材制造微结构的宏观屈服强度随激光扫描速度的增加而降低, 且呈各向异性; FIP极值符合Gumbel极值分布规律, 激光扫描速度增加可降低增材制造微结构疲劳分散性, 但会导致FIP极值升高, 使得疲劳裂纹萌生驱动力增加, 疲劳寿命降低.
, doi: 10.6052/0459-1879-21-389
双自由面溶质−热毛细液层的不稳定性
赵诚卓, 胡开鑫
溶质−热毛细对流是流体界面的浓度和温度分布不均导致的表面张力梯度驱动的流动, 它主要存在于空间微重力环境、小尺度流动等表面张力占主导的情况中, 例如晶体生长、微流控、合金浇筑凝固、有机薄液膜生长等. 对其流动进行稳定性分析具有重要意义. 本文采用线性稳定性理论研究了双自由面溶质−热毛细液层对流的不稳定性, 得到了两种负毛细力比(η)下的临界Marangoni数与Prandtl数(Pr)的函数关系, 并分析了临界模态的流场和能量机制. 研究发现: 溶质−热毛细对流和纯热毛细对流的临界模态有较大的差别, 前者是同向流向波、逆向流向波、展向稳态模态和逆向斜波, 后者是逆向斜波和逆向流向波. 在Pr较大时, Pr增加会降低流动稳定性; 在其他参数下, Pr增加会增强流动稳定性. 在中低Pr, 溶质毛细力使流动更加不稳定; 在大Pr时, 溶质毛细力的出现可能使流动更加稳定; 在其他参数下, 溶质毛细力会减弱流动稳定性. 流动稳定性不随η单调变化. 在多数情况下, 扰动浓度场与扰动温度场都是相似的. 能量分析表明: 扰动动能的主要能量来源是表面张力做功, 但其中溶质毛细力和热毛细力做功的正负性与参数有关.
, doi: 10.6052/0459-1879-21-148
微重力下成一定夹角平板间的表面张力驱动流动的研究
陈上通, 吴笛, 王佳, 段俐, 康琦
空间微重力环境中, 由于重力基本消失, 表面张力等次级力发挥主要作用, 流体行为与地面迥异, 因此有必要深入探究微重力环境中的流体行为规律和特征. 板式贮箱利用板式组件在微重力环境中对流体进行管理, 从而为推力器提供不夹气的推进剂, 这对航天器精确进行姿态控制、轨道调整具有重要意义. 板式组件中常包含成一定夹角的平板结构, 比如蓄液叶片之间. 本文研究了微重力环境中成一定角度平板间的表面张力驱动流动问题, 考虑了液体与壁面的动态接触角、对流引起的压力损失、黏滞阻力、液池内弯曲的液面等因素的影响, 推导出了表面张力驱动流动中液体爬升高度的二阶微分方程. 该方程可用四阶Runge−Kutta方法求解. 通过同时考虑两个主导力, 可将流动过程分为三个阶段, 并得到了不同阶段内的爬升高度的近似方程. 本研究建立了6个不同尺寸的计算模型、选用3种不同型号的硅油, 利用有限体积法开展仿真工作, 仿真结果与理论结果吻合良好, 验证了理论解的正确性. 本文的研究结果可为板式贮箱的研制和空间流体管理提供理论依据和数据支撑.
, doi: 10.6052/0459-1879-21-261
激光定向能量沉积的粉末尺度多物理场数值模拟
黄辰阳, 陈嘉伟, 朱言言, 廉艳平
激光定向能量沉积技术作为一种同轴送粉式金属增材制造技术, 以其制造效率高、成形尺寸大等优势在航空、航天、交通等领域具有广阔的应用前景. 然而, 该技术在金属零件的尺寸精度和形状精度控制方面存在诸如尺寸偏差大、表面不平整等控形问题, 亟需发展高效高精度预测熔覆层成形尺寸形貌的数值模拟方法. 针对该问题, 本文建立了考虑激光-粉末-熔池交互过程的高保真多物理场数值模型. 其中, 采用Gauss面热源等效激光光束, 采用拉格朗日质点法求解粉末输送及其与激光交互的过程, 进一步结合有限体积法和流体体积法求解粉末-熔池的交互及其流动凝固过程, 并通过TC17合金单道熔覆层实验结果进行了验证. 基于该模型, 本文首先预测了不同工艺参数下单道熔覆层形貌尺寸, 并对熔覆层形貌的变化趋势及其内在的物理机理进行了深入分析. 结果表明, 依赖于工艺参数的粉末温度分布和粉末基板能量分配比例对熔池流场和熔覆层尺寸有显著的影响. 本文所建立的数值模型可辅助激光定向能量沉积增材制造技术控形工艺参数优化, 所得结论可为成形件尺寸和形状精度控制提供理论指导.
, doi: 10.6052/0459-1879-21-420
基于扭摆的微冲量测量方法及实验研究
杨超, 贺建武, 章楚, 康琦, 段俐
超高微重力水平的卫星平台在空间引力波探测、地球重力场测量中发挥着重要的作用, 脉冲微推力器可以帮助微重力卫星实现姿态控制. 微冲量是评价脉冲微推力器性能的重要指标之一, 常用的基于扭摆的微冲量测量方法有两种, 方法一是根据单个冲量元瞬间作用于无阻尼扭摆后, 扭摆转动最大角位移计算冲量, 方法二是根据高固定频率的连续脉冲作用于有阻尼扭摆后, 扭摆转动的平均角位移计算冲量. 为了在地面实现对脉冲微推力器的微冲量测量, 利用已有的基于扭摆的亚微牛级推力测量系统, 进行了实验研究. 利用静电梳齿产生的标准静电力标定已有的扭摆推力测量系统, 通过电容式位移传感器测量扭摆角位移, 进而得到推力与角位移的关系, 以及其他扭摆系统参数; 然后, 根据两种冲量测量方法, 再以电磁螺线圈与永磁体分别产生瞬时磁力与固定频率的磁力作用于扭摆, 研究推力测量系统微冲量测量性能. 实验结果表明: 使用方法一时, 推力测量系统冲量测量范围为0.05 ~ 220 μN·s, 分辨力可达到0.02 μN·s; 而利用方法二测量微冲量相比于方法一而言, 能够扩大冲量测量范围, 提高冲量分辨能力.
, doi: 10.6052/0459-1879-21-191
面向增材制造的熔池凝固组织演变的相场研究
肖文甲, 许宇翔, 宋立军
激光增材制造(laser additive manufacturing, LAM)技术极适合复杂整体构件的近净成形和高附值损伤件的快速修复. 然而, 激光增材制造熔池内部复杂的动态凝固过程显著影响成形件的终态组织, 进而制约其服役性能. 本文针对激光直接能量沉积(direct energy deposition by laser, DED-L) Inconel 718过程, 构建宏观传热传质与多相场耦合的多尺度数学模型, 解决了熔池宏−微观温度场的直接耦合, 并基于MPI并行程序设计实现了熔池二维的全域定量模拟, 研究了凝固过程中的晶粒演变过程. 结果表明, 模拟的熔池尺寸、凝固界面与实验结果吻合较好. 熔池凝固界面形态和晶体择优取向是影响晶粒演变的重要因素. 在熔池横截面上, 凝固过程主要受温度梯度方向的驱使, 取向与温度梯度方向夹角越小的晶粒占优生长. 在纵截面上, 晶粒的生长表现出弯曲生长以及“上三角”的晶粒特征, 温度梯度方向的渐变导致了晶粒弯曲, 相邻晶粒的竞争行为决定了晶粒形貌. 本文阐明了金属激光增材制造晶粒演变的机理, 有助于厘清增材制造热物理、化学、冶金过程, 为凝固组织的预测和调控提供理论指导. 此外, 该多尺度数学模型也适用于其他金属材料的激光增材制造过程.
, doi: 10.6052/0459-1879-21-364
基于半解析VOF-DEM的激光直接沉积多尺度过程模拟
王泽坤, 刘谋斌
与传统铸造技术相比, 基于金属粉末的增材制造技术因其生产周期短、可操作性强而在航空航天、生物医学等领域具有很好的优越性. 尤其是激光直接沉积技术, 因其自由度高, 在复杂构件制造、部件修复中有着广泛的运用. 但是该激光直接沉积过程涉及多物理场、跨尺度、极端高温高压环境和相变问题, 仅靠实验不能很好地研究其中的机理. 已有数值模拟技术一般通过预设或者射入拉格朗日点作为颗粒输入, 不能做到同时考虑环境气体、颗粒碰撞和相变过程. 本文在近期发展的基于核函数近似背景流场的半解析CFD-DEM耦合方法中引入了流体体积分数法(VOF), 发展了可以同时模拟含热、刚体颗粒、相变和自由液面及相变界面的半解析VOF-DEM (或半解析CFD-DEM-VOF)方法, 从而首次实现了真实物理环境下激光直接沉积技术的数值模拟. 其中, VOF中的气相为环境气体, 液相为熔融和凝固的金属相, 界面通过iso-Advector重构, DEM为未熔化的金属粉末, 且流体网格可解析离散元颗粒形状. 这一模拟框架可以有效复现颗粒之间的碰撞、粘结、熔化、融合, 以及熔池熔道的形成, 为激光直接沉积技术的数值模拟提供了开拓性的范式, 并可以运用到其他带相变的颗粒系统中.
, doi: 10.6052/0459-1879-21-361
翼型颤振压电俘能器的输出特性研究
田海港, 单小彪, 张居彬, 隋广东, 谢涛
压电俘能器能够为自然界中低功率的微机电系统持续供能. 为了模拟机翼的沉浮−俯仰二自由度运动和有效俘获气动弹性振动能量, 本文提出一种新颖的翼型颤振压电俘能器. 基于非定常气动力模型, 推导翼型颤振压电俘能器流−固−电耦合场的数学模型. 建立有限元模型, 模拟机翼的沉浮−俯仰二自由度运动, 获得机翼附近的涡旋脱落和流场特性. 搭建风洞实验系统, 制作压电俘能器样机. 利用实验验证理论和仿真模型的正确性, 仿真分析压电俘能器结构参数对其气动弹性振动响应和俘获性能的影响. 结果表明: 理论分析、仿真模拟和实验研究获得的输出电压具有较好的一致性, 验证建立数学和仿真模型的正确性. 仿真分析获得机翼附近的压力场变化云图, 表明交替的压力差驱动机翼发生二自由度沉浮−俯仰运动. 当风速超过颤振起始速度时, 压电俘能器发生颤振, 并表现为极限环振荡. 当偏心距为0.3和风速为16 m/s时, 可获得最大输出电压为17.88 V和输出功率为1.278 mW. 功率密度为7.99 mW/cm3, 相比较于其他压电俘能器, 能实现优越的俘获性能. 研究结果对设计更高效的翼型颤振压电俘能器提供重要的指导意义.
, doi: 10.6052/0459-1879-21-377
机翼尺度效应对等离子体分离流动控制特性的影响
阳鹏宇, 张鑫, 赖庆仁, 车兵辉, 陈磊
等离子体流动控制技术是一种以等离子体气动激励为控制手段的主动流动控制技术. 为了进一步提高等离子体激励器可控机翼尺度, 以超临界机翼SC(2)-0714大迎角分离流为研究对象, 以对称布局介质阻挡放电等离子体为控制方式, 以测力、粒子图像测速仪(PIV: Particle Image Velocimetry)为研究手段, 从等离子体激励器特性研究出发, 深入开展了机翼尺度效应对等离子体控制的影响研究, 提出了适用于分离流控制的能效比系数, 探索了分离流等离子体控制机理, 掌握了机翼尺度对分离流控制的影响规律. 结果表明, (1)随着机翼尺度的增大, 布置到机翼上的激励器电极长度会相应增加; 在本文的参数研究范围内, 激励器的平均消耗功率不会随电极长度的增加而线性增大; 当电极长度达到一定阈值时, 激励器的平均消耗功率趋于定值; (2)在固定雷诺数的情况下, 随着机翼尺度的增大, 等离子体的控制效果并未降低, 激励器能效比系数提高; (3)等离子体在主流区诱导的大尺度展向涡与在壁面附近产生的一系列拟序结构成为分离流控制的关键. 研究结果为实现真实飞机的等离子体分离流控制, 推动等离子体流动控制技术工程化应用提供了技术支撑.
, doi: 10.6052/0459-1879-21-379
多个椭圆柱波浪力的一种解析解
赵密, 龙彭振, 王丕光, 张超, 杜修力
波浪在大尺寸结构表面会产生不可忽略的散射波, 该散射波会在多柱体体系中继续传播, 并在同体系中的其他柱体上产生高次散射波. 本文基于椭圆坐标系和绕射波理论首先推导了波浪作用下椭圆单柱体产生的散射波压力公式, 随后考虑将该散射波在多柱体系中的传播, 将其视为第二次入射波, 推导出柱体上第二次散射波压力公式, 同理可以推导出高次散射波压力公式, 最后得到椭圆多柱体波浪力解析解, 并用数值解验证了本文解析方法的正确性; 本文以双柱体和四柱体体系为例, 分析了不同参数(波数, 净距, 波浪入射角度等)下, 高次散射波对柱体上波浪作用的影响. 结果表明: 波数较大的情况下, 高次散射波引起柱体上的波浪力不能忽略; 结构间距较大的情况下, 虽然高次波的作用有减小的趋势但仍然明显; 高次散射波来自多个柱体对入射波的散射, 柱体数目的增加后, 高次波的影响会增加, 结构所受的高次波作用因参数变化而起的波动会变剧烈; 高次波对上游柱体波浪力的贡献较对下游柱体的贡献大.
, doi: 10.6052/0459-1879-21-318
带附加质量块的压电圆板能量采集器振动分析
孟莹, 丁虎, 陈立群
基于圆板的压电能量采集技术在取代化学电池为低功耗电子器件提供能源方面具有巨大的潜能. 本文通过理论建模和数值仿真研究了考虑附加质量接触面积的压电圆板能量采集器的采集性能. 首先, 基于基尔霍夫薄板理论, 用广义哈密顿原理推导了带附加质量块的压电圆板能量采集器的机电耦合方程, 并用伽辽金法对方程近似离散, 通过离散方程得到电压、功率输出和最优负载阻抗的闭合解. 用有限元仿真对所提出的理论模型进行了验证, 结果表明该理论模型可以成功地预测压电圆板能量采集器输出电压和功率. 最后, 基于闭合解探讨了负载阻抗、附加质量块、压电圆板的内外半径等相关参数对压电圆板能量采集器固有频率、输出电压和功率的影响. 结果表明, 当质量块与复合板的接触半径足够小(本文中接触半径小于板半径的1/14)时, 质量块与复合圆板的接触面积可以忽略; 相较于无孔的压电片, 内径位于2.5 ~ 4 mm范围内的压电片可以提高能量采集器的采集性能; 附加质量、压电片外径和负载阻抗的合理选择既可以降低压电圆板的固有频率, 还可以提高其采集性能.
, doi: 10.6052/0459-1879-21-441
压电与摩擦电复合型旋转能量采集动力学协同调控机制研究
赵林川, 邹鸿翔, 刘丰瑞, 魏克湘, 张文明
低转速激励下能量采集性能差是目前制约旋转能量采集技术应用的瓶颈问题. 本文提出了动力学协同调控机制, 并用于调控系统的动力学行为, 可以使器件在低转速激励下有效工作, 提高了旋转能量采集系统的电学性能. 旋转刚度软化、非线性磁力、几何边界的协同调控既可以增加系统在低速下的振动位移以及压电材料的形变, 也可调控系统的最大位移, 使其振动可控并限制位移过大提高可靠性. 此外, 几何边界可以方便地集成摩擦纳米发电机, 实现压电与摩擦两种机电转换机制在振动和碰撞过程中协同发电, 有效利用空间和提高输出电能. 基于哈密顿原理建立了系统的机电耦合动力学模型并进行了实验验证. 实验结果表明系统能够在0~250 r/min的低转速范围内有效工作, 在转速为250 r/min时, 压电单元和摩擦纳米发电机的最大峰峰值电压分别为132 V和1128 V, 总平均功率为1426 μW. 本文提出的动力学协同调控机制为能量采集系统动力学和电学性能改进提供新的途径, 有益于促进自供能物联网技术的发展与应用.
, doi: 10.6052/0459-1879-21-410
气相爆轰驱动二级轻气炮内弹道数值模拟
尚甲豪, 邢好运, 汪球, 李进平, 赵伟, 魏炳忱
二级轻气炮是超高速弹丸驱动技术中使用最广泛的技术之一, 它在超高速气动物理现象及材料高速碰撞下力学性能的实验研究和验证方面起着不可或缺的作用. 中国科学院力学研究所基于爆轰驱动方法研制了一座大型二级轻气炮, 可弥补高压气体驱动能力有限和火药使用受限的不足. 本文基于经过实验验证的准一维数值模拟方法, 详细研究了该设备的内弹道动力学参数及发射性能, 并探讨了不同发射方法及装填参数对设备性能的影响规律和机理. 研究结果表明, 氢氧爆轰驱动相比于高压气体驱动具有明显优势; 不同爆轰驱动方式对弹丸发射性能影响较小, 但其影响到整个设备的强度设计; 对装填运行参数的研究表明增大爆轰段充气压力可以有效加强轻气炮发射性能, 而活塞质量变化对发射速度的影响较为复杂, 轻气炮实际运行中受设备设计指标及模型材料性能的限制, 优化过程中需要同时调整3种参数以达到轻气炮最佳性能.
, doi: 10.6052/0459-1879-21-437
正常和早期膝骨关节炎的软骨生物力学研究
林伟健, 李俊言, 陈瑱贤, 靳忠民
膝骨关节炎是导致膝关节疼痛和慢性残疾的一种常见的关节疾病. 膝关节的软骨生物力学是评价膝骨关节炎程度的重要指标. 然而, 早期膝骨关节炎的软骨生物力学依然有待研究, 正常、内侧和内外侧早期膝骨关节炎的软骨生物力学差异尚不清楚. 本文基于固−液双相纤维增强的软骨有限元建模方法, 分别建立了正常膝关节模型、内侧早期膝骨关节炎模型和内外侧早期膝骨关节炎模型, 在步态周期中最大载荷时刻和最大屈曲角度时刻下分别对比分析了正常、内侧和内外侧早期膝骨关节炎3种情况下的软骨生物力学差异. 结果表明, 与正常膝关节相比, 内侧早期膝骨关节炎模型的内侧软骨的流体压力减少, 固相等效应力减少以及应变增大; 外侧软骨的结果基本没有差异. 然而, 内外侧早期膝骨关节炎模型的内外侧软骨的流体压力都减少, 固相等效应力都减少以及应变都增大. 早期膝骨关节炎中退变软骨的属性变化会导致软骨的承载能力下降以及变形增大, 从而增加软骨进一步退变的风险. 本文提出的基于双相纤维增强软骨模型的膝关节有限元模型有效预测了正常和关节炎情况下的软骨生物力学差异, 该模型也可以推广应用于髋、踝和脊柱等其他关节生物力学的研究.
, doi: 10.6052/0459-1879-21-390
基于人工神经网络的亚格子应力建模
吴磊, 肖左利
亚格子(SGS)应力建模在湍流大涡模拟(LES)中有着极为重要的作用. 传统亚格子应力模型存在相对误差较大、耗散过强等问题. 近年来, 计算机技术的发展使得人工神经网络(ANN)等机器学习方法逐渐成为亚格子应力建模型的新研究范式. 本文着重考虑滤波宽度及雷诺数影响, 在不可压缩槽道湍流中建立了亚格子应力的ANN模型. 该模型以滤波后的直接数值模拟(fDNS)流场物理量及滤波尺度为输入信息, 相应滤波尺度下的亚格子应力为输出量. 通过对不同滤波尺度及不同雷诺数数据的训练, ANN模型能够给出与直接数值模拟(DNS)高度吻合的亚格子应力. 此外, 模型在亚格子耗散等非ANN建模量上也有着优异的预测性能, 与基于DNS获得的对应物理量的相关系数大都在0.9以上, 较梯度模型及Smagorinsky模型有明显提升. 在后验测试中, ANN模型对流向平均速度剖面的预测同样优于梯度模型、Smagorinsky模型及隐式大涡模拟(ILES)等传统LES模型. 在脉动速度均方根预测方面, 除了某些法向位置外ANN模型的性能整体上相对其他3个模型有所提升. 然而, 随着网格尺度的增大ANN模型预测的结果与fDNS结果的偏差逐渐增大. 总之, ANN方法在发展高精度亚格子应力模型上具有很大的潜力.
, doi: 10.6052/0459-1879-21-356
基于SBFEM和深度学习的裂纹状缺陷反演模型
江守燕, 万晨, 孙立国, 杜成斌
结构内部缺陷的识别是结构健康监测的重要研究内容, 而当前以无损检测为主的结构安全检测多以定性分析为主, 定量识别缺陷的尺度较困难. 本文将比例边界有限元法(scaled boundary finite element methods, SBFEM)和深度学习相结合, 提出了基于Lamb波在结构中传播时的反馈信号定量识别结构内部裂纹状缺陷的反演模型. 通过随机生成缺陷信息(位置、大小), 采用SBFEM模拟Lamb波在含不同缺陷信息的结构中的信号传播过程, SBFEM仅需对结构边界离散可最小化网格重划分过程, 大大提高了计算效率. Lamb波在含裂纹状缺陷结构中传播时观测点的反馈信号包含大量的裂纹信息, 基于这一特性可为深度学习模型提供足够多的反映问题特性的训练数据. 建议的缺陷反演模型规避了传统反分析问题的目标函数极小化迭代过程, 在保证计算精度的前提下大大减少了计算成本. 对含单裂纹和多裂纹板的数值算例进行分析, 结果表明: 建立的缺陷识别模型能够准确地量化结构内部的缺陷, 对浅表裂纹亦有很好的识别效果, 且对于含噪信号模型仍具有较好的鲁棒性.
, doi: 10.6052/0459-1879-21-360
物理增强的流场深度学习建模与模拟方法
金晓威, 赖马树金, 李惠
流体运动理论上可用Navier−Stokes方程描述, 但由于对流项带来的非线性, 仅在少数情况可求得方程解析解. 对于复杂工程流动问题, 数值模拟难以高效精准计算高雷诺数流场, 实验或现场测量难以获得流场丰富细节. 近年来, 人工智能技术快速发展, 深度学习等数据驱动技术可利用灵活网络结构, 借助高效优化算法, 获得对高维、非线性问题的强大逼近能力, 为研究流体力学计算方法带来新机遇. 有别于传统图像识别、自然语言处理等典型人工智能任务, 深度学习模型预测的流场需满足流体物理规律, 如Navier−Stokes方程、典型能谱等. 近期, 物理增强的流场深度学习建模与模拟方法快速发展, 正逐渐成为流体力学全新研究范式: 根据流体物理规律选取网络输入特征或设计网络架构的方法称为物理启发的深度学习方法, 直接将流体物理规律显式融入网络损失函数或网络架构的方法称为物理融合的深度学习方法. 研究内容涵盖流体力学降阶模型、流动控制方程求解领域.
, doi: 10.6052/0459-1879-21-373
水中开孔腔流激振荡控制实验研究
章文文, 徐荣武, 何琳, 潘龙德, 赵佳锡
水中开孔腔流激振荡是水下航行器的一类突出噪声源. 为探究有效抑制水中开孔腔流激振荡的控制方法和作用特性, 首先以水下航行器的表面开孔结构为对象设计了开孔腔模型, 并提出一种基于来流边界层分流原理的流动控制装置——前缘分流体, 借助循环水洞装置对水中开孔腔的流激振荡特性及其控制进行了实验研究. 通过沿流向和展向安装于腔底的动态压力传感器测量腔内脉动压力, 分别从腔内脉动压力的频谱特性和空间分布特性两方面, 探讨了水中开孔腔在不同流速下的流激振荡特性和前缘分流体对水中开孔腔流激振荡的控制效果, 并对前缘分流体的主要作用机理进行了分析. 研究结果表明: 水中开孔腔流激振荡形式以剪切层自持振荡为主, 在流速较低时, 如2.4 m/s, 就会产生稳定的自持振荡, 且具有随流速升高而急剧增大的趋势; 前缘分流体对水中开孔腔绕流自持振荡具有良好的抑制效果, 且抑制效果随流速增加而显著提升, 对腔内脉动压力频谱峰值和总级的最大抑制量分别达到25.3 dB和15.6 dB; 此外, 前缘分流体对开孔腔流激振荡具有低频频移作用, 有益于避免发生流激空腔共振; 脉动压力空间分布特性表明, 前缘分流体对水中开孔腔流激振荡抑制机理主要在于破坏了腔内流场受到的周期调制作用.
, doi: 10.6052/0459-1879-21-143
机器学习在力学模拟与控制中的应用专题序
王建春, 晋国栋
为了改进机器学习方法在力学问题中的应用效果, 并提高该方法的预测精度、计算效率和泛化能力, 需要回答一些基本问题, 包括: (1)如何从数据中系统地学习抽象的力学规律和本构关系? (2)如何将物理原理和机器学习有效地结合在一起? (3)如何提升模型的鲁棒性和泛化能力? (4)如何大幅度地提升复杂力学系统的模拟效率? (5)针对多尺度或高维度的力学系统, 如何构造可靠、高效的降阶模型? (6)如何系统地提升对大型设计空间的快速搜索能力? 围绕上述这些重要的问题, 《力学学报》组织了《机器学习在力学模拟与控制中的应用》这一专题. 由于篇幅限制, 该专题包含了3篇综述论文和7篇研究论文, 从一个侧面反映了国内科研人员在该方向上的一部分最新研究进展, 供读者参考.
, doi: 10.6052/0459-1879-21-501
基于离散单元法和人工神经网络的近壁颗粒动力学特征研究
段总样, 赵云华, 徐璋
颗粒与壁面的相互作用往往对颗粒流动具有显著影响. 为研究颗粒与壁面作用机理, 对滚筒内颗粒流动过程进行离散单元法(DEM)数值模拟. 基于模拟结果统计分析靠近壁面处颗粒的运动特征, 结果表明, 小摩擦系数时颗粒平动和旋转速度均近似满足正态分布, 但由于壁面影响, 摩擦系数增大时颗粒沿滚筒轴向的旋转速度偏离正态分布, 颗粒动力学理论推导壁面边界条件时应考虑速度正态分布的修正及速度脉动的各向异性. 采用人工神经网络(ANN)构建了颗粒无因次旋转温度、滑移速度和平动温度之间的函数模型, 进而可以在常规双流模型壁面边界条件中考虑颗粒旋转的影响. 基于DEM模拟及结果分析可以为壁面边界条件的理论构造和半经验修正提供基础数据和封闭模型.
, doi: 10.6052/0459-1879-21-313
均匀电场中气泡上升特性的实验研究
陈烁, 王太, 苏硕, 谢英柏, 刘春涛
电场中的气泡对于强化传热有显著的作用, 对于电场中气泡动力学特性的研究对增强换热器效率, 提高能源利用率有重要意义. 为了获得外加电场作用下气泡的动力学特性, 设计与搭建了可视化实验平台. 采用50 kV高压直流电源构建均匀电场, 高清摄像机拍摄实验图像. 引入电场强度、气泡体积与溶液介电常数作为变量, 探究其对于气泡动力学特性的影响. 观测了竖直与水平均匀电场中气泡的上升过程, 分析了不同变量下气泡变形状况与上升速度的变化. 引入气泡长宽比L/D用于表示气泡拉伸变形程度, 截取单个气泡上升过程分时段图像展示形态变化过程. 研究结果表明, 气泡沿电场方向伸长, 且电场强度越大, 变形越明显; 竖直电场中气泡伸长导致上升速度增大, 而水平电场中气泡上升速度减小. 气泡尺寸增大, 浮升力作用增强, 气泡上升速度增大. 溶液介电常数增加, 电场力作用明显增加, 气泡变形更加明显.
, doi: 10.6052/0459-1879-21-352
基于能力评估的空间翻滚目标抓捕策略优化
许若男, 罗建军, 王明明
由于目标的翻滚运动, 空间双臂机器人对动态目标的抓捕相比于静态目标更具有挑战性. 对抓捕策略进行优化可以提高空间双臂机器人对翻滚目标的操作能力以保证任务的成功. 本文提出了一种基于能力评估的抓捕策略优选方法. 空间双臂机器人捕获目标时, 双臂末端执行器与目标同时接触形成闭链系统, 闭链约束的引入使操作能力的评估更加复杂. 在对双臂空间机器人协调操作翻滚目标的运动学与动力学分析基础上, 建立了考虑闭链约束的协调工作空间, 并分析了基于任务兼容度的消旋能力评估指标. 建立的协调工作空间同时包含位置和姿态信息, 可以用于灵巧度的计算. 接着, 基于协调工作空间的全局灵巧度指标确定机械臂末端执行器对目标的最优抓捕点, 以及考虑相机视角约束和末端执行器对目标速度跟踪约束下的力任务兼容度指标确定空间双臂机器人捕获翻滚目标时的最优抓捕构型. 利用能力评估确定抓捕策略可以充分利用双臂的协调性以增加对动态目标的操作能力, 通过仿真验证了所提抓捕策略的可行性和有效性.
, doi: 10.6052/0459-1879-21-436
考虑混凝土塑性耗散的CDM-XFEM裂缝计算方法
金浩, 余朔
混凝土结构在服役期间受外界荷载的影响容易产生裂缝, 导致结构刚度降低、构件承载性能衰退, 而采用准确的计算方法预测混凝土裂缝的发展是治理裂缝的基本前提, 也是保障结构安全的重要手段. CDM方法能够描述微裂缝的扩展过程, 但不能表示离散的开裂面, 且存在网格诱导偏差及虚假应力传递的弊端, XFEM方法能够描述宏观裂纹的扩展过程, 但不能反映微裂缝的动态扩展, 两者计算出的裂纹分布与实际差异均较大. 现有的CDM-XFEM方法已经能够模拟混凝土微裂缝及宏观裂缝发展的整个过程, 但忽略了宏观裂缝出现时混凝土产生的塑性应变, CDM与XFEM的能量转化过程欠缺平衡性. 因此, 本文重点考虑能量转化时的塑性耗散, 选取指数型函数为粘结裂缝的牵引-分离模式, 基于能量及应力等效的条件重新构建了CDM与XFEM之间的能量转化方程. 采用广义逆最小二乘法求解能量转化系数, 确定能量转化时的临界位移, 并给出了裂缝面水平集的更新算法及整体计算方法的程序流程. 以双切口混凝土受剪拉开裂试验为例, 采用多种裂缝计算方法与试验进行了对比. 结果表明, 采用考虑混凝土塑性耗散的CDM-XFEM方法算出的裂缝分布及拉力-张开位移曲线与试验结果差异最小, 说明采用考虑混凝土塑性耗散的CDM-XFEM计算方法能够更好地计算混凝土裂缝.
, doi: 10.6052/0459-1879-21-272
页岩储层纳微米孔隙CO2/CH4吸附及驱替特性研究进展
邓佳, 吕子健, 张奇, 宋付权, 李久江, 赵广杰
利用CO2开采页岩气不仅能够提高页岩气采收率, 还能够节省水资源并且对CO2进行地质封存, 有助于实现页岩气开采过程的碳中和. 富有机质页岩储层纳微米孔隙中气体运移机制不同于常规储层, CO2在储层中具有超临界特性, 致使开采机理复杂, 无法得到CO2开采页岩气微观机理的准确认识, 所以研究CH4、CO2及其二元混合物在页岩储层纳微米孔隙中的吸附及驱替特性对准确评估和高效开采页岩气至关重要. 本文从实验、理论以及模拟方面对页岩储层纳微米孔隙中CH4的吸附特性、CO2/CH4二元混合物竞争吸附特性以及驱替特性进行了综合分析, 对气体在纳微米孔隙中吸附及驱替特性的基础研究及关键问题进行讨论分析并提出了展望. 研究表明CH4在页岩储层中表现为物理吸附, 有机质特征(丰度、成熟度、类型)、孔隙结构、无机矿物组成、温度和压力、含水率对页岩的CH4吸附能力均有一定程度的影响. 在相同条件下, CO2比CH4更易被页岩储层吸附, 在页岩储层中注入CO2可以促进CH4的解吸, 并有利于CO2的地质埋存. 开采方案的部署可采用井网形式的注采方式, 可以通过调整注入井的位置、数量以及CO2注入速率对开采方案进行优化.
, doi: 10.6052/0459-1879-21-292
斜爆轰发动机的推力性能理论分析
杨鹏飞, 张子健, 杨瑞鑫, 滕宏辉, 姜宗林
爆轰燃烧具有释热快、循环热效率高的特点. 斜爆轰发动机利用斜爆轰波进行燃烧组织, 在高超声速吸气式推进系统中具有重要地位. 以往研究主要关注斜爆轰波的起爆、驻定以及波系结构等, 缺少从整体层面出发对斜爆轰发动机开展推力性能分析. 本文将斜爆轰发动机内的流动和燃烧过程分解成进气压缩、燃料掺混、燃烧释热和排气膨胀4个基本模块并分别进行理论求解, 建立了斜爆轰发动机推力性能的理论分析模型. 在斜爆轰波系研究成果的基础上, 选取了过驱动斜爆轰、Chapman−Jouguet斜爆轰、过驱动正爆轰和斜激波诱导等容燃烧等4种燃烧模式来描述燃烧室内的燃烧释热过程, 并对比分析了不同燃烧模式对发动机比冲性能的影响. 此外, 还获得了不同来流参数、燃烧室参数和进排气参数等对发动机推力的影响规律, 发现来流马赫数和尾喷管的膨胀面积比是发动机理论燃料比冲的主要影响因素. 最后, 结合以往关于受限空间内斜爆轰波驻定特性等方面的研究成果, 提出了斜爆轰发动机燃烧室的设计方向.
, doi: 10.6052/0459-1879-21-206
倍频激励双转子振动同步机理与实验研究
邹敏, 方潘, 侯勇俊, 彭欢, 王德金
振动筛作为第一级钻井液净化设备, 其筛分性能直接决定后续固控系统的生产效率. 针对在筛分工程中由相同驱动频率激励的振动系统所表现出的单一振动特性难以与筛面上的物料粒度相匹配, 从而引发筛孔堵塞的问题, 本文提出了倍频激励双转子自同步振动系统. 首先根据数学模型运用广义Lagrange公式建立多自由度振动系统的运动微分方程, 并引入复变函数求解箱体的稳态幅值响应; 其次探明了振动系统实现倍频同步的前提条件和稳定性评估标准; 然后通过数值计算定量地讨论转子间的倍频动力学特征与系统结构参数之间的关系, 并结合Runge−Kutta算法建立振动系统的多自由度机电耦合动力学仿真模型, 详细讨论转子与箱体间的倍频同步机理; 最后设计实验样机, 针对系统在不同工况下的动态特性和同步运动状态展开实验测试, 进一步验证理论研究和计算机模拟的正确性. 研究表明, 系统的倍频同步能力指数随安装距离的不断增加而接近于零值附近, 此时高频转子与低频转子获得稳定同步振动的可能性越来越高; 电机同步状态几乎不受弹簧刚度的影响, 但稳态相位差值在单周期内会随激振器安装倾角和距离的变化而逐渐减小并趋于在一个恒定值. 研究成果不仅对石油工业新型振动筛系统的研制具有重要的参考价值, 也将促进其他同步振动机械的发展.
, doi: 10.6052/0459-1879-21-359
一种改进的颗粒移动床μ(I)拟流体模型及应用
吴坤, 刘向军, 戴椰凌
颗粒移动床在工业领域应用广泛, 发展实用可靠的颗粒移动床模型具有理论和应用价值. 本文基于颗粒流μ(I)模型, 补充局部颗粒体积分数与颗粒局部压力和局部颗粒流密度的关系式, 将移动床内密集颗粒处理成可压缩拟流体, 建立了颗粒流单相可压缩流μ(I)模型, 并建立了颗粒流−壁面摩擦条件, 在计算中对颗粒流拟黏度和拟压力项进行正则化处理. 采用上述模型与方法对3种典型散料在移动床缩口料仓内的流动进行模拟, 与实验对比, 得到了玻璃珠、刚玉球和粗沙的μ(I)模型参数, 分析了3种不同散料在料仓内的颗粒速度、体积分数等分布特性, 模拟结果较好地揭示了料仓内不同物料的整体流和漏斗流特性; 进而以玻璃珠为例, 对移动床颗粒单管绕流流动进行了模拟, 所得结果合理揭示了管流附近的流动特性. 计算结果表明, 对于本文的计算工况, 颗粒体积分数变化最大范围为0.510 ~ 0.461, 绝大部分区域流动惯性数小于0.1, 改进的单相μ(I)模型能合理预测出密集颗粒流移动床内的流动特性, 方法可行且较多相流算法能明显减小计算量.
, doi: 10.6052/0459-1879-21-320
基于基因表达式编程的数据驱动湍流建模
赵耀民, 徐晓伟
计算流体动力学是湍流研究的重要手段, 其中雷诺平均模拟在航空航天等实际工程中得到了广泛应用. 雷诺平均模拟的结果很大程度上依赖于湍流模型的预测精度, 而实际工程应用中常用的模型往往精度有限. 近年来, 数据驱动的湍流建模方法得到越来越多的关注. 本文介绍了基于基因表达式编程 (gene-expression programming, GEP) 方法的湍流建模相关进展. 本文首先讨论基因表达式编程应用于湍流建模的具体方法, 包括基本算法、显式代数应力模型和湍流传热两种建模框架、模型测试方法以及损失函数设置等. 在此基础上, 基因表达式编程方法被应用于涡轮叶栅尾流混合、竖直平板间自然对流、三维横向流中的射流等问题. 结果表明, GEP可以有效提升常用模型对于尾流混合损失、壁面热通量等关键参数的预测精度. 基因表达式编程方法可以显式给出模型方程, 因此模型具有可解释性强等特点. 基于双向耦合方法得到的模型还被证明具有较好的后验测试精度和鲁棒性. 基因表达式编程方法还被初步应用于大涡模拟亚格子应力和边界层转捩等问题的建模, 在不同湍流建模领域表现出很大的潜力.
, doi: 10.6052/0459-1879-21-391
基于深度强化学习算法的颗粒材料应力−应变关系数据驱动模拟研究
狄少丞, 冯云田, 瞿同明, 于海龙
颗粒材料的宏观力学行为受颗粒组分等材料参数, 孔隙率、配位数等状态参数的影响, 同时又具备复杂的加载路径和加载历史相关性, 建立包含多个内变量以及各变量间相互关联的颗粒材料本构模型是一个重要的科学难题. 不同于传统的基于屈服面、流动法则和硬化函数框架下的唯象本构模型, 本文基于颗粒物质力学的研究基础, 以颗粒材料平均孔隙率、细观组构参数和弹性刚度参数作为内变量, 结合深度学习方法建立以有向图表征的数据本构模型. 有向图中以不同的链接网络表示不同的内变量信息流动方向, 各个内变量间的映射关系采用循环神经网络来建立, 将各个神经网络相互组合, 形成包含不同内变量且具有不同预测能力的本构模型. 该本构模型的建立过程等价于在众多可能的内变量链接关系空间中寻找最能描述实际材料宏观应力应变行为的优化问题. 因此, 可将有向图本构模型的建立过程看作“马尔可夫决策过程”, 采用深度强化学习算法构建有向图的内变量链接组合优化过程, 具体采用AlphaGo Zero算法自动寻找最优的颗粒材料数据驱动本构模型建模路径. 研究结果表明, 采用有向图和深度强化学习算法可建立起完全依靠“数据驱动”的颗粒材料应力−应变关系. 此外, 本方法提供了一种将不同理论模型从数据角度统一起来, 且基于人工智能算法发展更优模型的研究思路, 可为相似问题的研究提供借鉴.
, doi: 10.6052/0459-1879-21-312
基于人工神经网络的非结构网格尺度控制方法
王年华, 鲁鹏, 常兴华, 张来平, 邓小刚
网格自动化生成和自适应是制约计算流体力学发展的瓶颈问题之一, 网格生成质量、效率、灵活性、自动化程度和鲁棒性是非结构网格生成的关键问题. 在非结构网格生成中, 网格空间尺度分布控制至关重要, 直接影响网格生成质量、效率和求解精度. 采用传统的背景网格法进行空间尺度分布控制需要在背景网格上求解微分方程得到背景网格上的尺度分布, 再将网格尺度从背景网格插值到真实空间点, 过程十分繁琐且耗时. 本文从效率和自动化角度提出两种网格尺度控制方法, 首先发展了基于径向基函数(RBF)插值的网格尺度控制方法, 通过贪婪算法实现边界参考点序列的精简, 提高了RBF插值的效率. 同时, 还采用人工神经网络进行网格尺度控制, 初步引入相对壁面距离和相对网格尺度作为神经网络输入输出参数, 建立人工神经网络训练模型, 采用商业软件生成二维圆柱和二维翼型非结构三角形网格作为训练样本, 通过训练和学习建立起相对壁面距离和相对网格尺度的神经网络关系. 进一步实现了二维圆柱、不同的二维翼型的尺度预测, RBF方法和神经网络方法的效率与传统背景网格法相比提高了5~10倍, 有助于提高网格生成的效率. 最后, 将方法推广应用于各向异性混合网格尺度预测, 得到的网格质量满足要求.
, doi: 10.6052/0459-1879-21-334
时间尺度上Lagrange系统的Hojman守恒量
张毅, 田雪, 翟相华, 宋传静
利用对称性和守恒律, 可以简化动力学问题甚至求解力学系统的精确解, 更好地理解其动力学行为. 时间尺度分析将连续和离散动力学模型统一并拓展到时间尺度框架, 既避免了重复研究又可揭示两者之区别和联系. 因此, 通过对称性来探寻在时间尺度的框架下新的守恒定律很有必要. 文章首先建立了时间尺度上Lagrange方程, 利用时间尺度微积分性质导出了时间尺度上Lagrange系统的两个重要关系式; 其次, 依据微分方程在单参数Lie变换群下的不变性, 建立了时间尺度上Lie对称性的定义和确定方程; 最后, 建立了时间尺度上Lie对称性定理并利用上述关系式给出了证明, 得到了时间尺度上Lagrange系统的新守恒量. 当时间尺度取为实数集时, 该守恒量退化为著名的Hojman守恒量. 文末考察了一个两自由度时间尺度Lagrange系统, 在三种不同时间尺度情形下得到了该系统的Hojman守恒量, 数值计算结果验证了定理的正确性.
, doi: 10.6052/0459-1879-21-413
K0超固结结构性黏土的本构模型
万征, 曹伟, 刘媛媛, 张红芬
K0固结黏土在自然界广泛分布, 其通常同时具有超固结性与天然结构性, 而K0超固结性又与K0正常固结性质存在很大差异. 为了有效的描述K0超固结性质, 在结构性模型基础上, 做了如下三点改进, 使得原模型拓展为同时考虑K0超固结特性与天然结构性影响的本构模型. (1)引入相对应力比来描述屈服面, 并引入初始各向异性转轴参量ξ来表达初始各向异性对屈服面在p-q空间的位置影响. (2)基于给定的屈服面方程, 推导得到变相应力比参量, 并将变相应力比引入到统一硬化参数中, 利用统一硬化参数可以有效描述初始各向异性固结黏土在剪切加载下的剪缩与剪胀, 应变硬化及软化现象. (3)引入反映结构性胶结强度性质的胶结参量pe, 并给出pe随塑性偏应变的衰减演化方程, 利用胶结参量可描述结构性黏土的剪胀特性. 预测与试验结果对比表明, 所提的K0超固结结构性模型可有效描述K0超固结黏土的刚度提高效应, 黏土的包辛格效应, 结构性黏土胶结强度的损失现象以及结构性黏土的应变软化现象. 证明了所提模型的适用性以及合理性.
, doi: 10.6052/0459-1879-21-265