EI、Scopus 收录
中文核心期刊

留言板

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

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

2021年  第53卷  第6期

研究综述
基于高温气体效应的磁流体流动控制研究进展
罗凯, 汪球, 李逸翔, 李进平, 赵伟

高超声速飞行器强激波后高温气体形成具有导电性的等离子体流场,电离气体为磁场应用提供了直接工作环境.磁流体控制技术利用外加磁场影响激波后的离子或电子运动规律,可有效地改善高超声速飞行器气动特性,在飞行器气动力操控和热环境管理等方面均具有广阔的应用前景; 同时,超导材料及电磁技术的发展又重新推动了这一领域的研究热潮.虽然国内外在高超声速磁流体流动控制领域已开展了一些研究工作,但其实验研究依然极具挑战, 且由于实验条件及测量技术等限制,其压力、热流等参数的测量并没有得出较为系统的结论,因此需要对影响脱体激波距离、热流、压力变化的规律及机理进行深入研究; 同时,数值模拟方法和理论分析也亟待可靠的实验数据来对其进行验证.本综述调研和讨论了基于高温真实气体效应的磁流体流动控制技术研究,主要针对磁流体流动控制的试验技术、数值模拟、理论方法以及流动控制的主要研究方向等进行了总结,并对其发展趋势进行了讨论和展望.

2021, 53(6): 1515-1531. doi: 10.6052/0459-1879-21-067
流体力学
基于组合神经网络的雷诺平均湍流模型 多次修正方法
张珍, 叶舒然, 岳杰顺, 王一伟, 黄晨光
求解雷诺平均(Reynolds-averaged Navier-Stokes, RANS)方程依然是工程应用中有效且实用的方法, 但对雷诺应力建模的不确定性会导致该方法的预测精度具有很大差异. 随着人工智能的发展, 湍流闭合模型结合机器学习元素的数据驱动方法被认为是提高RANS模型预测性能的有效手段, 然而这种数据驱动方法的稳定性和预测精度仍有待进一步提高. 本文通过构建一个全连接神经网络对RANS方程中的涡黏系数进行预测以实现雷诺应力的隐式求解,该神经网络记作涡黏系数神经网络(eddy viscosity neural network, EVNN). 此外, 也使用张量基神经网络(tensor basis neural network, TBNN)预测未封闭量与解析量之间的高阶涡黏关系, 并利用基张量保证伽利略不变性. 最后, 采用多次修正的策略实现修正模型对流场预测的精度闭环. 上述方法使用大涡模拟(large eddy simulation, LES)方法产生的高保真数据, 以及RANS模拟获得的基线数据对由EVNN和TBNN组合的神经网络进行训练, 然后用训练好的模型预测新的RANS模拟的流场. 通过与高保真LES结果进行对比, 结果表明, 相比于原始RANS模型, 修正模型对后验速度场、下壁面平均压力系数和摩擦力系数的预测精度均有较大提升. 可以发现对雷诺应力线性部分的隐式处理可以增强数值求解的稳定性, 对雷诺应力非线性部分的修正可以提升模型对流场各向异性特征预测的性能, 并且多次修正后的模型表现出更高的预测精度. 因此, 该算法在数据驱动湍流建模和工程应用中具有很大的应用潜力.
2021, 53(6): 1532-1542. doi: 10.6052/0459-1879-21-073
基于数据驱动的流场控制方程的稀疏识别
江昊, 王伯福, 卢志明
利用有限数据建立系统的非线性动力学模型是具有挑战性的重要课题. 数据驱动的稀疏识别方法是近年来发展的从数据识别动力系统控制方程的有效方法. 本文基于数据驱动稀疏识别方法对不同流场的控制方程进行了识别. 采用非线性动力学偏微分方程函数识别(partial differential equations functional identification of nonlinear dynamics, PDE-FIND)方法和最小绝对收缩和选择算子(least absolute shrinkage and selection operator, LASSO)方法对二维圆柱绕流、顶盖驱动方腔流、Rayleigh-Bénard (RB)对流和三维槽道湍流的控制方程进行了识别. 在稀疏识别过程中, 采用直接数值模拟得到的流场数据来计算过完备候选库中的每一项, 候选库中变量最高保留到二次, 变量导数最高保留到二阶, 非线性项最高保留到四阶. 结果发现PDE-FIND方法和LASSO方法对于不含有非线性项的控制方程, 如涡量输运方程、热输运方程和连续性方程, 都能准确识别. 对于含有强非线性项的控制方程, 如Navier-Stokes方程的识别, PDE-FIND方法正确地识别出了控制方程及流场的Rayleigh数和Reynolds数, 而LASSO方法识别结果不正确, 这是因为候选库中的项之间存在分组效应, LASSO方法通常只取分组中的一项. 本文还发现选择流动结构丰富的区域的数据进行控制方程的稀疏识别可以提高识别的准确性.
2021, 53(6): 1543-1551. doi: 10.6052/0459-1879-21-052
串列布置三圆柱涡激振动频谱特性研究
涂佳黄, 胡刚, 谭潇玲, 梁经群, 张平
对串列三圆柱体双自由度涡激振动问题进行了数值计算, 并分析了雷诺数、固有频率比和约化速度对串列三圆柱体结构动力响应及频谱特性的影响. 研究发现: 雷诺数、频率比对上游圆柱的振幅和流体力系数的影响较小. 中游圆柱频率锁定区域随着雷诺数的增大而增大, 其动力响应受上游圆柱尾流的影响较大, 但频率比的影响较小. 同时, 流体力系数在约化速度较小时受雷诺数和频率比的影响较大. 另外, 下游圆柱的振幅和流体力系数受雷诺数及频率比的影响较大. 雷诺数、频率比和约化速度对圆柱流体力系数能量谱密度(PSD)曲线中主峰幅值、频谱成分及波动性的影响较大. 流体力系数PSD曲线波动性的增强, 导致圆柱运动轨迹会从"8"字形转变成不规则形状. 当频率比为2.0时, 上游圆柱尾流出现P$+$S模式, 导致其发生非对称运动, 且升、阻力系数PSD曲线主峰重合. 最后, 激励荷载平均功率值随约化速度的变化趋势与对应的结构动力响应的变化类似. 在同一约化速度区间内, 结构振动响应的强弱与位移的平均功率值成正比. 对不同约化速度区间内的升力系数功率谱密度分析时, 振动频率比($f_{s}/f_{n, y})$对结构振动响应的影响更大.
2021, 53(6): 1552-1568. doi: 10.6052/0459-1879-21-036
颗粒群碰撞搜索及CFD-DEM耦合分域 求解的推进算法研究
刘巨保, 王明, 王雪飞, 姚利明, 杨明, 岳欠杯
在采用计算流体力学-离散元耦合方法(computational fluiddynamics-discrete element method, CFD-DEM)进行固液两相耦合分析时, 颗粒计算时间步的选取直接影响到耦合计算精度和计算效率. 为此, 本文选取每个目标颗粒为研究对象, 引入插值函数计算时间步的运动位移, 构建可变空间搜索网格; 通过筛选可能碰撞颗粒建立搜索列表, 采用逆向搜索方式判断碰撞颗粒, 从而提出一种改进的DEM方法(modified discreteelement method, MDEM). 该算法在颗粒群与流体耦合计算中, 颗粒计算初始时间步选取不受颗粒碰撞时间限制, 通过自动调整和修正实现大步长, 由颗粒和流体耦合条件实时更新流体计算时间步, 使颗粒计算时间步选取过小导致计算效率低、选取过大导致颗粒碰撞漏判的问题得以解决, 为颗粒与流体耦合的数值模拟提供了行之有效的计算方法. 通过两个颗粒和多个颗粒的数值模拟, 得到的颗粒间碰撞力、碰撞位置及次数, 与理论计算结果的相对误差均低于2%, 与传统的DEM碰撞搜索算法相比, 在选取的3种计算时间步均不会影响计算精度, 且有较高的计算效率. 通过多个颗粒与流体的耦合数值模拟, 采用传统的CFD-DEM方法, 只有颗粒计算时间步选取10$^{-6}$ s或更小才能得到精确解, 而采用本文方法取10$^{-4}$ s也能够得到精确解, 避免了颗粒碰撞随时间步增大而出现的漏判问题, 且计算耗时降低了16.7%.
2021, 53(6): 1569-1585. doi: 10.6052/0459-1879-21-002
一种基于块雅可比迭代的高阶FR格式隐式方法
于要杰, 刘锋, 高超, 冯毅
最近, 基于非结构网格的高阶通量重构格式(flux reconstruction, FR)因其构造简单且通用性强而受到越来越多人的关注. 但将FR格式应用于大规模复杂流动的模拟时仍面临计算开销大、求解时间长等问题. 因此, 亟需发展与之相适应的高效隐式求解方法和并行计算技术. 本文提出了一种基于块Jacobi迭代的高阶FR格式求解定常二维欧拉方程的单GPU隐式时间推进方法. 由于直接求解FR格式空间和隐式时间离散后的全局线性方程组效率低下并且内存占用很大. 而通过块雅可比迭代的方式, 能够改变全局线性方程组左端矩阵的特征, 克服影响求解并行性的相邻单元依赖问题, 使得只需要存储和计算对角块矩阵. 最终将求解全局线性方程组转化为求解一系列局部单元线性方程组, 进而又可利用LU分解法在GPU上并行求解这些小型局部线性方程组. 通过二维无黏Bump流动和NACA0012无黏绕流两个数值实验表明, 该隐式方法计算收敛所用的迭代步数和计算时间均远小于使用多重网格加速的显式Runge-Kutta格式, 且在计算效率方面至少有一个量级的提升.
2021, 53(6): 1586-1598. doi: 10.6052/0459-1879-20-404
基于LBM的铝微滴斜柱沉积水平偏移研究
任彦霖, 刘赵淼, 逄燕, 王翔
金属微滴沉积制造技术采用逐点堆砌方式成型, 为斜柱沉积提供无支撑制造方式, 具有高度灵活性. 本文针对铝液滴斜柱连续沉积过程, 建立格子玻尔兹曼模型进行数值模拟, 研究液滴在凝固表面上的水平偏移运动. 根据表面能充放过程, 沉积运动被划分为下落、快速扩张、慢速扩张、回弹4个阶段, 其受力状态由表面能、重力势能、动能和黏性耗散趋势得到. 液滴内部流动在扩张阶段中表现为滑动状态, 而在回弹阶段中表现为滚动状态. 液滴偏移运动的加速阶段主要发生在扩张阶段, 而偏移距离则在回弹阶段中产生. 扩张阶段的受力状态表明偏移运动的主要推动力为重力和毛细力. 随着液滴轴线距离的增大, 扩张阶段中的加速段时间缩短、速度峰值提高, 使水平偏移距离呈先增大后减小的趋势, 这种阶段化特征源于加速段时长和速度极大值的竞争关系. 不同沉积高度和固液浸润度下, 偏移距离均保持相同的演化趋势. 在相同的轴线距离下, 偏移距离随固液浸润度的增大而减小, 随沉积高度的增大而减小. 通过拟合水平偏移距离演化规律、优化扫描步距, 能够实现斜柱的均匀沉积, 并使倾角与理论结果一致.
2021, 53(6): 1599-1608. doi: 10.6052/0459-1879-21-022
固体力学
磁场力及膜曲率对磁敏感薄膜-基底界面 黏附性能的影响与调控
韩明杰, 彭志龙, 姚寅, 张博, 陈少华
界面黏附和脱黏的可调控在攀爬装置、黏附开关、机械抓手等方面具有重要的应用需求. 针对磁敏感薄膜-基底界面, 开展了薄膜初始曲率及外加磁场对界面黏附性能影响机制的研究. 首先实验制备了具有初始曲率的磁敏感薄膜, 分别开展了具有初始曲率的磁敏感薄膜-基底界面撕脱实验及理论研究, 研究了薄膜初始曲率、弯曲刚度和外加磁场强度对界面黏附性能的影响规律. 实验和理论结果一致表明: 具有初始曲率的磁敏感薄膜-基底界面黏附力随薄膜初始曲率的增大而减小, 而外加磁场能够有效提高界面黏附力;相比于初始零曲率薄膜-基底界面稳态撕脱力与薄膜弯曲刚度无关, 薄膜弯曲刚度减弱了具有初始曲率薄膜-基底界面的稳态撕脱力. 进一步从能量角度分析了界面等效黏附性能, 揭示了薄膜弯曲能、磁场势能、界面黏附能的相互竞争机制. 最后, 基于本文的实验及理论结果, 提出了一种磁场和薄膜初始曲率协同调控的简易机械抓手, 可连续实现物体的拾取、搬运和释放功能. 本文结果不仅有助于理解多场调控的界面可逆黏附机制, 对界面黏附可控的功能器件设计亦提供了一种新方法.
2021, 53(6): 1609-1621. doi: 10.6052/0459-1879-21-091
颗粒增强复合材料压缩行为的位错动力学模拟
丁一凡, 魏德安, 陆宋江, 刘金铃, 康国政, 张旭
颗粒增强铜基复合材料因具有极高的强度和弹性模量, 优异的导电、导热性能和抗磨损能力, 被广泛应用于航天航空、汽车、电子工业等领域. 第二相强化是其主要的强化方式, 其通过合金中弥散的微粒阻碍位错运动, 可有效提高金属材料的力学性能, 提高其服役安全. 针对该问题本文采用三维离散位错动力学(three-dimensional discrete dislocation dynamics, 3D-DDD)方法, 对微尺度颗粒增强铜基复合材料进行了微柱压缩模拟, 分析了位错与第二相颗粒交互作用对材料力学响应的影响, 揭示第二相颗粒强化的微观机理. 本研究将第二相颗粒视为位错不可穿透的球形微粒, 采用位错绕过机制模拟颗粒与位错的交互作用过程. 通过调控滑移面相对于第二相颗粒中心的距离发现: 屈服应力和应变硬化率均随距离的增大而减小. 研究也发现Schmid因子越高的滑移系, 屈服应力越低, 后续应变硬化率越低. 多位错与颗粒交互作用的模拟发现, 同一滑移面中位错间的反应和不同滑移系中位错的交互作用可能是导致屈服应力和应变硬化率降低的关键.
2021, 53(6): 1622-1633. doi: 10.6052/0459-1879-21-028
考虑颗粒转矩的接触网络诱发各向异性分析
王怡舒, 沈超敏, 刘斯宏, 陈静涛
颗粒材料的宏观力学行为与接触网络的组构各向异性密切相关, 根据接触点的滑动与否、转动与否和强弱力情况, 可以将颗粒间的接触系统分为不同的子接触网络. 一般而言, 不同的子接触网络在颗粒体系中的传力机制不同, 对宏观力学响应的贡献也有不同. 采用离散单元法(discrete element method, DEM)模拟了不同抗转动系数$\mu_r$下颗粒材料三轴剪切试验, 分析了剪切过程中不同子接触网络的组构张量的演变规律, 并探究了颗粒抗转动效应对子接触网络各向异性指标演变规律的影响. 研究发现: 剪切过程中转动、非转动接触的组构张量变化不是独立的, 受到颗粒间滑动与否的影响; 非滑动、强接触网络是颗粒间的主要传力结构, 非滑动接触网络的接触法向和法向接触力各向异性均随$\mu_r$的增大而增大, 其对宏观应力的贡献程度随$\mu_r$的增大而减小;强接触网络的接触法向各向异性随$\mu_r$的增大而增大, 但法向接触力各向异性随$\mu_r$的增大无明显变化, 强接触网络对宏观应力的贡献程度在不同$\mu_r$情况下均相同.
2021, 53(6): 1634-1646. doi: 10.6052/0459-1879-21-090
完整岩石拉压应力状态下的张裂破坏准则
吕爱钟, 刘宜杰, 尹崇林
一点的应力状态可由3个主应力$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$来表示, 当规定主应力以压为正时, 沿最大主应力$\sigma_{1}$方向将产生收缩变形, 若中间主应力$\sigma_{2}$和最小主应力$\sigma_{3}$都远小于$\sigma_{1}$, 则沿$\sigma_{2}$和$\sigma_{3}$方向会产生横向扩张变形, 当横向扩张变形达到一定极限时, 将会在平行于$\sigma _{1}$的方向产生张裂破坏. 如何建立这种张裂破坏的强度准则目前尚缺乏研究, 最大拉应变理论(第二强度理论)有时被用来解释张裂破坏, 但最大拉应变理论难以应用于三向受力状态. 本文分别用$\varepsilon_{1}$, $\varepsilon_{2}$表示最大张应变和次大张应变, 则最大拉应变理论认为当$\varepsilon_{1}$达到单向拉伸屈服应变时, 材料将产生破坏. 而本文将根据$\varepsilon_{1}+\varepsilon_{2}$之和达到极限值$\varepsilon_u$来建立张裂破坏准则. 可以证明$\varepsilon_{1} +\varepsilon_{2}$所表示的是$\sigma_{1}$主平面的面积增长率. 当$\sigma_{3}<\sigma_{2} \ll \sigma_{1}$时, 大部分岩石都具有脆性破坏的特点, 所以可将破坏前的岩石视为满足广义胡克定律的线弹性材料, 这样用$\varepsilon_{1}$, $\varepsilon_{2}$表示的强度准则可通过$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$来表示. 在这个过程中还可考虑岩石在拉伸和压缩时具有不同弹性参数和强度的特点, 并可通过单向拉伸和单向压缩的破坏状态来确定$\varepsilon_u$. 不管$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$是压应力, 还是拉应力, 或者$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$中有拉有压的情形, 基于$\varepsilon_{1} +\varepsilon_{2} =\varepsilon_u$都可建立相应的强度准则. 所建立的准则可以反映中间应力$\sigma_{2}$对强度的影响规律, 通过建立的强度准则还可以证明: 静水拉力能引起屈服, 而静水压力不能产生屈服; 压缩破坏能使塑性体积增大, 其结果比Mohr-Coulomb准则更能反映实际情形. 并通过拉压应力状态下的试验数据验证了所建立的强度准则, 所得理论计算结果和已有的试验数据吻合得很好. 通过提出的强度准则和圆盘劈裂的试验结果, 可获得更为可靠的岩石单轴抗拉强度.
2021, 53(6): 1647-1657. doi: 10.6052/0459-1879-21-026
铁木辛柯梁中的卸载弯曲波及二次断裂
龙龙, 郑宇轩, 周风华, 任会兰, 宁建国
半无限长梁承受恒定弯矩作用后, 如果自由端的初始弯矩突然释放, 将在梁中激发出一列卸载弯曲应力波. 采用铁木辛柯梁和瑞利梁来研究突然卸载所激发出的弯曲波的传播特征. 利用拉普拉斯变换方法进行分析, 首先推导出铁木辛柯梁和瑞利梁中的卸载弯曲波的像函数解析解, 采用数值反变换方法给出了时域上波传播的响应解, 并研究了梁中各点的横向位移、弯矩和剪力随时间的变化规律. 计算结果表明: 与简化的欧拉梁不同, 旋转惯性的引入使铁木辛柯梁和瑞利梁中的弯曲波传播具有强烈的局部化效应, 特别是梁中各点经历的弯矩变化, 和其距离自由端的位置相关, 不同时刻的弯矩峰值大小不同;瑞利梁中离自由端不同距离各点的峰值弯矩先增大后降低, 最后达到一个渐近值;铁木辛柯梁中各点的峰值弯矩总体上随着时间单调增大到同一个渐近值, 该渐近值与欧拉梁中的峰值弯矩值相同, 均为1.43.切应力效应的引入进一步降低了铁木辛柯梁中卸载弯曲波的波速, 同时也使得铁木辛柯梁中弯矩峰值的最大值小于瑞利梁中的最大值. 对于脆性细长梁的纯弯曲断裂, 铁木辛柯梁可以较好地预测二次断裂的发生位置, 相应的碎片尺寸约为7倍梁横截面厚度.
2021, 53(6): 1658-1670. doi: 10.6052/0459-1879-21-106
半结晶聚合物损伤演化的实验表征与数值模拟
张毅, 薛世峰, 韩丽美, 周博, 刘建林, 贾朋
损伤本构模型对研究材料的断裂失效行为有重要意义, 但聚合物材料损伤演化的定量表征实验研究相对匮乏. 通过4种高密度聚乙烯(high density polythylene, HDPE)缺口圆棒试样的单轴拉伸实验获得了各类试样的载荷-位移曲线和真应力-应变曲线, 采用实验和有限元模拟相结合的方法确定了HDPE材料不同应力状态下的本构关系, 并建立了缺口半径与应力三轴度之间的关系;采用两阶段实验法定量描述了4种HDPE试样单轴拉伸过程中的弹性模量变化, 并建立了基于弹性模量衰减的损伤演化方程, 结合中断实验和扫描电子显微镜分析了应力状态对HDPE材料微观结构演化的影响. 结果表明缺口半径越小, 应力三轴度越大, 损伤起始越早、演化越快; 微观表现为: 高应力三轴度促进孔洞的萌生和发展, 但抑制纤维状结构的产生;基于实验和有限元模拟获得的断裂应变、应力三轴度、损伤演化方程等信息提出了一种适用于聚合物的损伤模型参数确定方法, 最后将本文获得的本构关系和损伤模型用于HDPE平板的冲压成形模拟, 模拟结果与实验结果吻合良好.
2021, 53(6): 1671-1683. doi: 10.6052/0459-1879-21-101
基于人工弹簧模型的周期结构带隙计算方法研究
冯青松, 杨舟, 郭文杰, 陆建飞, 梁玉雄
能量法具有将求解微分方程边值问题转化为泛函极值问题的优点,故而在结构动力学分析中被广泛使用, 近年来也被引入到周期结构带隙计算中. 然而,由于周期结构边界条件相对复杂,采用传统能量法(如Rayleigh-Ritz法)分析时位移函数构造难度大;且由于位移函数中包含波数项,扫描波数计算带隙的过程中质量、刚度矩阵需不断重算, 导致计算量较大. 鉴于此,本文对传统能量法进行改进,通过引入人工弹簧来模拟包含周期边界在内的各类边界条件,可将边界约束转化为人工弹簧的弹性势能,故而各能量分部中仅有周期边界弹性势能包含波数项,扫描波数时仅需重新计算与其对应的刚度矩阵,其余的质量、刚度矩阵只需要计算一次, 继而显著降低了计算量. 研究结果表明,本文方法准确、可靠, 且相较于传统能量法, 本文方法的计算效率更高,随着结构质量、刚度矩阵的维度增大, 或者扫描波数点数的增多,本文方法计算效率优势更加明显. 此外, 人工弹簧模型使用灵活、便捷,可进一步地拓展到更为复杂的周期性组合结构带隙分析中.
2021, 53(6): 1684-1697. doi: 10.6052/0459-1879-21-007
动力学与控制
张拉整体结构的动力学等效建模与实验验证
陈占魁, 罗凯, 田强
为了实现张拉整体结构高效动力学计算, 并考虑其大范围运动中柔性杆局部动态屈曲, 提出了一种受压细长杆动力学降阶模型, 采用五节点弹/扭簧集中质量离散模型等效连续杆的静力学和动力学特性. 首先, 通过静力学等效分析推导了弹簧拉压刚度和扭簧弯曲刚度表达式, 可准确预测杆件受压屈曲和近似预测其后屈曲行为. 第二, 通过动能等效分析推导了集中质量表达式, 可准确预测杆在线速度场下的运动. 第三, 通过弯曲振动固有模态等效分析确定弯曲刚度和节点质量的分布参数, 合适的分布参数取值组合可将降阶模型前两阶固有频率相对误差均降低至1%以内. 第四, 在全局坐标系下建立张拉整体结构瞬态动力学方程, 并利用静力凝聚法实现方程高效迭代求解. 最后, 分别对球形张拉整体结构准静态压缩、模态分析和碰撞动力学进行仿真和实验对比分析, 证明了提出的动力学降阶模型可有效预测张拉整体结构的静力学行为、固有振动特性及瞬态动力学响应, 并分析了结构参数变化对其力学特性的影响规律. 本文提出的动力学等效建模与计算方法, 可望用于软着陆行星探测器、大型可展开空间结构及点阵材料等复杂张拉整体系统的动力学分析与控制.
2021, 53(6): 1698-1711. doi: 10.6052/0459-1879-21-056
Hamel 框架下几何精确梁的离散动量守恒律
高山, 史东华, 郭永新
Hamel场变分积分子是一种研究场论的数值方法, 可以通过使用活动标架规避几何非线性带来的计算复杂度, 同时数值上具有良好的长时间数值表现和保能动量性质. 本文在一维场论框架下, 以几何精确梁为例, 从理论上探究Hamel场变分积分子的保动量性质. 具体内容包括: 利用活动标架法对几何精确梁建立动力学模型, 通过变分原理得到其动力学方程, 利用其动力学方程及Noether定理得到系统动量守恒律; 将几何精确梁模型离散化, 通过变分原理得到其Hamel场变分积分子, 利用Hamel场变分积分子和离散Noether定理得到离散动量守恒律, 并给出离散动量的一阶近似表达式; Hamel场变分积分子可在计算中利用系统对称性消除系统运动带来的非线性问题, 但此框架中离散对流速度、离散对流 应变及位形均不共点, 而这种错位导致离散动量中出现级数项, 本文对几何精确梁的离散动量与连续形式的关系及其应 用进行了讨论, 并通过算例验证了结论. 上述证明方法也同样适用一般经典场论场景下的Hamel场变分积分子. Hamel场变分积分子的动量守恒为进一步研究其保结构性质提供了参考依据.
2021, 53(6): 1712-1719. doi: 10.6052/0459-1879-21-092
负刚度时滞反馈控制动力吸振器的等峰优化
代晗, 赵艳影
相比于传统动力吸振器, 负刚度动力吸振器同时具有更好的减振能力和更宽的有效减振频带宽度, 为了进一步降低共振峰幅值, 在负刚度吸振器系统耦合时滞反馈控制. 对负刚度时滞反馈控制动力吸振器系统进行等峰优化设计, 优化设计的准则是:第一和第二共振峰的峰值相等; 同时兼顾两个目标, 一个目标是在优化时的最大共振峰幅值小于被动负刚度吸振器系统的反共振峰幅值, 另一目标是在优化时共振峰幅值与反共振峰幅值差小于被动吸振器系统. 接着, 通过设计和调节负刚度系数、吸振器阻尼系数和时滞反馈控制系数对控制系统进行等峰优化设计. 最后, 在降低幅值的同时, 分析结构参数对有效减振频带宽度的影响. 经过等峰优化之后, 选择本文的一组结构参数与两个典型的模型进行对比. 为了定量比较不同模型的降幅效果, 定义了减幅百分比, 研究发现在有效减振频带区间内减幅百分比超过40%以上. 结果表明, 通过等峰优化准则对结构参数进行优化设计和调节增益系数和时滞量, 共振峰幅值的减幅百分比也近似达到40%, 也可以调节增益系数和时滞量, 使得幅频响应曲线具有较宽的有效减振频带和较低的共振峰幅值与反共振峰幅值的差值.
2021, 53(6): 1720-1732. doi: 10.6052/0459-1879-21-074
电磁场下神经元模型中混合簇的分岔机制
冀文超, 段利霞, 齐会如
Pre-B?tzinger复合体是新生哺乳动物呼吸节律起源的关键部位, 是呼吸节律产生的中枢. 忆阻器的功能类似于神经元突触的可塑性, 可用其模拟磁通量.本文在Butera动力学模型的基础上引入刺激电流和磁通控制忆阻器, 分别研究这两个因素对单个pre-B?tzinger复合体神经元中混合簇放电模式的影响.通过无量纲化的方法对变量进行时间尺度分析, 结果表明, 模型包含3个不同的时间尺度.通过快慢分解和分岔分析研究了神经元混合簇放电产生和转迁的动力学机制.电流和磁通量都可以影响混合簇中胞体簇的个数, 减小电流和磁通量的值, 混合簇中胞体簇的个数也会相应减少, 并使簇的类型由"fold/homoclinic"型簇放电转迁为经由"fold/homoclinic"滞后环的"Hopf/Hopf"型簇放电.双参数分岔分析表明, 随着钙离子浓度的逐渐增加, 全系统轨线在鞍结分岔曲线和同宿轨分岔曲线之间来回跃迁, 是混合簇的产生分岔机制.全系统轨线在鞍结分岔曲线和同宿轨分岔曲线之间跃迁的次数, 与混合簇中胞体簇的个数相对应.
2021, 53(6): 1733-1746. doi: 10.6052/0459-1879-21-071
生物、工程及交叉力学
激波风洞高低压段钢膜片破裂特性研究
聂少军, 汪运鹏, 薛晓鹏, 姜宗林
激波风洞是用于高超声速飞行器气动外形设计和优化的常用地面试验装置,基于爆轰驱动技术,激波风洞能够在短时间(毫秒级)内产生高温、高压的驱动气体来模拟高超声速试验气流.主膜片位于激波风洞中的爆轰驱动段和激波管段之间,试验时膜片在爆轰脉冲压力下打开,膜片的打开状态和脱落情况对激波风洞气流品质有很大的影响. 同时,膜片也是形成激波的先决条件. 传统的风洞采用铝质膜片进行试验,在激波风洞中需要承压能力更强的膜片, 此时铝质膜片不再适用, 需要采用钢质膜片.因此, 对激波风洞中的钢膜片破裂特性进行研究很有必要.将数值计算结果与试验结果进行比较, 发现数值计算结果与试验结果吻合得比较理想,计算结果具有可靠性. 基于膜片的应力-应变模型, 建立了膜片打开的动力学模型,根据CJ爆轰理论, 采用有限元软件计算模拟了膜片破裂的过程,分析总结了膜片破裂的机制和力学特性规律.采用控制变量法对不同厚度和凹槽长度的膜片进行分析研究,得到了膜片破膜压力和有效破膜时间的变化规律. 在激波风洞试验中,根据膜片总破膜时间设计了适用于JF-12复现风洞的膜片参数.
2021, 53(6): 1747-1757. doi: 10.6052/0459-1879-20-341
考虑界面力学性能的组件及结构的协同优化
马晶, 赵明宣, 王浩淼, 刘湃, 亢战
包含多个内嵌功能组件及支撑结构的多组件结构因其轻量化、多功能等优良特性被广泛应用于航空航天等领域.已有的多组件结构拓扑优化研究大多基于理想界面假设,忽略了材料连接界面可能发生的破坏. 本文针对包含多个内嵌式功能性组件的结构,考虑连接界面的力学性能, 对组件的形状、布局及支撑材料的拓扑进行协同优化,以实现多组件结构的优良承载性能. 首先,基于超椭圆模型对内嵌组件的形状及布局进行显式的参数化描述,并构造其水平集函数表达; 进而,结合组件及支撑材料的水平集拓扑描述、内聚力模型及扩展有限元方法(extended finite element method, XFEM), 在固定网格下对随优化迭代不断演化的结构拓扑及连接界面的力学性能进行准确描述;进一步, 建立水平集法框架下考虑界面力学性能的多组件结构拓扑优化列式,基于伴随变量法推导解析的灵敏度并采用梯度优化算法求解优化问题.本文采用该优化框架分别对内嵌组件的悬臂梁和MBB梁进行协同优化, 在优化过程中,发现组件的初始布局对最终设计有很大的影响, 并且可能导致不良结构.为了避免此情况, 本文提出了两个阶段的优化策略,即首先对组件布局和形状进行优化, 再进行结构和内嵌组件的协同优化.数值结果显示, 在优化结果中功能性组件及界面通常分布于结构受压应力作用的区域,且连接界面最优形状呈现为曲率较小的光滑曲线,该设计避免界面发生拉伸及剪切破坏, 有效提高了结构的承载力,同时也表明了本文所提出考虑连接界面力学性能拓扑优化方法的有效性.
2021, 53(6): 1758-1768. doi: 10.6052/0459-1879-21-010
接地惯容式减振器对悬臂输流管稳定性 和动态响应的影响研究
郭梓龙, 王琳, 倪樵, 贾青青, 杨文正
输流管道广泛应用于机械、航空、核电和石油等重要工程领域.为防止管道结构因流致振动破坏造成的损失, 很有必要对其稳定性、动力学响应及其调控进行深入研究.本文提出一种由惯容器、弹簧和阻尼器并联组成的减振器模型, 研究了这种接地惯容减振器对悬臂输流管稳定性和非线性振动的影响. 首先, 基于哈密顿原理给出了带有接地惯容减振器非保守系统的非线性动力学模型; 然后, 利用高阶伽辽金方法对非线性方程进行离散化; 最后, 分别从线性和非线性角度分析了不同减振器参数下输流管道的被动控制效果, 着重讨论了惯容系数和减振器安装位置对悬臂管稳定性和动态响应的影响机制.线性理论模型的研究结果显示, 接地惯容减振器可显著影响悬臂管的失稳临界流速, 故通过调节减振器参数能有效提高输流管道的稳定性;惯容系数和弹簧刚度对系统稳定性的控制效果还与减振器的安装位置密切相关.非线性理论模型的分析结果显示, 惯容系数和减振器位置对输流管的非线性动态响应也有显著影响, 且这种影响还依赖于管道的流速取值; 在某些参数条件下, 减振器还可使输流管道由周期运动演化为复杂的混沌行为. 本文研究结果表明, 通过设计合理的惯容式减振器参数, 可提升悬臂输流管道的稳定性并有效抑制其颤振幅值.
2021, 53(6): 1769-1780. doi: 10.6052/0459-1879-21-105
地形和土-结相互作用效应对三维跨峡谷 桥梁地震响应的影响分析
陈少林, 伍锐, 张娇, 谷音
评估跨峡谷桥梁的地震性态需要考虑地形效应、行波效应以及土-结相互作用效应.将峡谷-桥梁系统在地震波输入下的反应分析看作波动散射问题,即桥梁及其邻近非规则区域对峡谷场地"自由场"的扰动. 基于此思想,本文发展了一套跨峡谷桥梁地震反应分析方法,通过二维模型分析得到峡谷场地的"自由场", 结合人工边界输入到峡谷-桥梁体系,采用土-结相互作用分区并行方法对其进行分析, 并编制了相应的分析程序.该方法可在自由场分析时考虑非垂直入射地震波, 计入行波效应,因此可综合考虑行波效应、地形效应和土-结相互作用效应. 通过峡谷场地分析算例,验证了自由场和人工边界实施的正确性; 并以马水河大桥为对象,通过5种计算模型结果的比较,分析了地形效应和土-结相互作用效应对跨峡谷桥梁地震反应的影响, 算例结果表明,地形效应对墩底剪力、弯矩和轴力有明显影响,对位移的影响要比对剪力、弯矩的影响小; 土-结相互作用对桥梁反应的影响较大,较大地减小了桥梁反应.
2021, 53(6): 1781-1794. doi: 10.6052/0459-1879-21-039
TA2钛合金开口柱壳外爆碎片分布研究
吴文苍, 董新龙, 庞振, 周风华
金属柱壳爆炸膨胀断裂机制及其对碎片分布、特征尺寸的影响是应用物理、力学、兵器工程等领域共同关心的重要课题, 但目前除数值模拟外, 考虑断裂机制的简单二维碎裂模型尚未出现.开展TA2钛合金开口柱壳在不同装药条件下的碎裂实验研究, 通过对软回收碎片的统计及微观分析, 探讨金属柱壳外爆断裂模式及二维碎片分布规律, 结果显示:(1) TA2钛合金柱壳在实验爆压(7 $\sim$ 25 GPa)下宏观断口均为剪切断裂模式, 但机制不同, 在较高爆压下柱壳剪切断裂由多重绝热剪切带破坏控制, 在较低压力下为剪切破坏;(2) 与一维拉伸碎裂相比, 柱壳爆炸碎裂不充分, 碎片质量更符合$\beta=1$ (或更接近1)的指数分布; 爆炸碎裂越充分, 碎片越小并趋于均匀, $\beta$趋于较小的值, 趋向Mott和Linfoot提出的泊松统计分布形式;(3) Rayleigh分布可以较好描述柱壳碎片的宽度分布规律, 不同爆压下柱壳碎片宽度归一化尺寸分布具有相似性, 呈现"量子化"特性, 即存在最小的特征尺寸; (4) TA2柱壳碎片特征尺寸远大于G-K剪切断裂公式预测的尺寸, G-K剪切式描述的是多重绝热剪切带间距.本研究为金属柱壳碎片特征、分布规律及其模型分析提供了重要参考.
2021, 53(6): 1795-1806. doi: 10.6052/0459-1879-21-017