EI、Scopus 收录
中文核心期刊

留言板

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

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

2012年  第44卷  第3期

显示方式:
研究论文
钝头体高超声速绕流底部失稳特征数值模拟
朱德华, 沈清, 王强, 袁湘江
利用数值模拟方法对高超声速钝锥及Apollo返回舱底部尾迹流场进行了研究, 分析尾迹流动的失稳过程. 对钝锥模型, 在M=6, Re=1.71× 106(Re以球头半径为参考长度)条件下观察到了底部流动的不稳定性. 不添加任何扰动, 数值模拟首先得到的流动是稳定解, 在底部发展出一个主分离区和一个二次分离区, 流动是轴对称状态. 继续进行计算, 发现二次分离线率先变形, 底部流场发展出非定常周期流动. 对Apollo返回舱模型, 在相同条件下 (Re以前面圆弧半径为参考长度), 数值模拟首先得到的流动同样是稳定解, 出现以二次分离线率先变形为起始的结构失稳, 演化出周期性过程, 但持续时间较短, 很快出现了非周期非对称状态. 研究表明, 高超声速钝锥及Apollo返回舱底部流场均存在不稳定性问题, Apollo返回舱的底部流场更加不稳定.
2012, (3): 465-472. doi: 10.6052/0459-1879-2012-3-20120301
柱形汇聚激波冲击球形重气体界面的实验研究
王显圣, 司廷, 罗喜胜, 杨基明
采用高速纹影法实验研究了柱形汇聚激波与球形重气体界面相互作用的 Richtmyer-Meshkov不稳定性问题. 激波管实验段基于激波动力学理论设计, 将马赫数为1.2 的平面激波转化为柱形汇聚激波, 气体界面由肥皂膜分隔六氟化硫(内)和空气(外)得到. 采用高速摄影机在单次实验中拍摄激波运动的全过程, 对柱形激波的形成进行了实验验证, 并进一步观测了汇聚激波与球形气体界面相互作用过程中的波系发展和气体界面变形以及反射激波同已变形界面二次作用的流场演化. 结果表明: 当柱形汇聚激波穿过气泡界面以后, 气泡左侧界面极点沿激波传播方向保持匀速运动, 气泡右侧界面发展成为射流结构, 气泡主体发展成为涡环结构; 在反射激波的二次作用下, 流场中无序运动显著增强并很快进入湍流混合阶段.
2012, (3): 473-480. doi: 10.6052/0459-1879-2012-3-20120302
磁场对双扩散液层热毛细对流的影响
李炜, 姜燕妮, 颜君毅, 陈启生
通过数值模拟的方法对磁场作用下的双扩散液层热毛细对流进行了研究, 模型中同时考虑了热毛细效应和溶质毛细效应的存在. 研究结果显示, 外部磁场能够有效削弱液层内热毛细对流的强度, 改变热毛细对流的对流结构; 随着磁场强度的增大, 液层内热毛细对流的对流强度逐渐减小, 热质传递过程中扩散效应逐渐得到增强; 最终, 溶质浓度沿水平方向呈梯度分布. 因此, 当磁场强度足够大时能够实现晶体生长中所需的纯扩散条件.
2012, (3): 481-486. doi: 10.6052/0459-1879-2012-3-20120303
LES/RANS方法混合函数特性研究
陈逖, 刘卫东, 孙明波, 范晓樯, 梁剑寒
将两方程k-ω SST湍流模型和Sagaut的混合尺度亚格子模型通过一个混合函数相结合, 构造出一种混合大涡/雷诺平均N-S方程模拟方法(hybird large eddy simulation/reynolds-averaged navier-stokes, Hybrid LES/RANS), 采用这种混合模拟方法结合5阶WENO格式对Ma=2.8平板湍流边界层进行了数值模拟, 并在计算区域上游入口处采用“回收/调节”方法生成湍流脉动边界条件, 通过考查RANS区域向LES区域的过渡参数及网格分辨率对这种混合模拟方法进行了评价. 计算结果表明: 该文采用的混合模拟方法可以捕捉到湍流边界层中的大尺度结构且入口边界层平均参数不会发生漂移, 混合函数应当将RANS区域和LES区域的过渡点设置在对数律层和尾迹律层的交界处, 而过渡应当迅速以获得正确的雷诺剪切应力分布, 在该文采用的模型及数值方法的条件下, 流向及展向的网格小至与Escudier混合长相当时, 能够获得可以接受的脉动速度的单点-二阶统计值.
2012, (3): 487-493. doi: 10.6052/0459-1879-2012-3-20120304
流固耦合分析的一种改进CBS有限元算法
周岱, 何涛, 涂佳黄
针对流固耦合问题, 发展了一种基于任意拉格朗日-欧拉(ALE)描述有限元法的弱耦合分区算法. 运用半隐式特征线分裂算法求解Navier-Stokes方程, 在压力Poisson 方程中引入质量源项以满足几何守恒律; 运用子块移动技术更新动态网格, 并配以光滑处理防止网格质量下降; 采用Newmark-β 法求解结构运动方程. 为保持流体-结构界面处速度和动量守恒, 利用修正结合界面边界条件方法求解界面处速度通量和动量通量. 运用本方法分别模拟了不同雷诺数下单圆柱横向和两向流致振动、串列双圆柱两向流致振动. 计算表明, 本文方法计算效率高, 计算结果与已有实验和数值计算数据吻合.
2012, (3): 494-504. doi: 10.6052/0459-1879-2012-3-20120305
用物理黏性构建高阶不振荡对流扩散差分格式
高智
利用数值摄动算法, 通过扩散格式数值摄动重构把对流扩散方程的2阶中心差分格式(2-CDS)重构为高精度高分辨率格式, 解析分析和模型方程计算证实了新格式的高精度不振荡性质. 新格式是把物理黏性使流动光滑化的扩散运动规律引入2-CDS 中的结果. 该法显然与构建高级离散格式的常见方法不同. 证实: 数值摄动重构中引入扩散运动规律的结果格式与引入对流运动规律(下游不影响上游的规律)的结果格式一致, 说明对离散方程的数值摄动运算, 在维持原格式结构形式不动的条件下, 不仅能提高格式精度和稳健性, 且可揭示对流离散运动规律与扩散离散运动规律之间的内在关联;同时证实, 文中提出和使用的上、下游分裂方法是构建高精度不振荡离散格式的一个有效方法.
2012, (3): 505-512. doi: 10.6052/0459-1879-2012-3-20120306
可压缩流场中气泡脉动数值模拟
王诗平, 孙士丽, 张阿漫, 陈玉
在应用边界元方法对气泡动力学的研究中, 绝大多数模型是建立在不压缩势流理论基础之上, 针对可压缩流场中气泡运动特性的研究很少. 从波动方程出发, 分别在气泡运动前期和后期对波动方程进行简化, 得到气泡运动局部和全局简化方程, 采用双渐进方法对简化方程进行匹配, 提出了考虑流场可压缩性的非球状气泡运动模型. 该模型的计算结果与Prospertti 等的解析结果吻合很好, 气泡脉动最大半径和内部最大压力随气泡脉动逐渐减小. 基于该模型对比了自由场中药包爆炸考虑可压缩性与不考虑可压缩性的计算结果, 发现考虑可压缩性气泡射流速度较小, 随后基于该模型计算了刚性边界下气泡的运动特性.
2012, (3): 513-519. doi: 10.6052/0459-1879-2012-3-20120307
水力机械非定常流动的三维涡方法计算
王龙步, 祝宝山, 王宏, 曹树良
利用拉格朗日型涡方法计算了混流式水轮机转轮内部的非定常流动. 采用高精度的边界元法计算域内势流速度以满足固体边界不可穿透条件, 通过在固体边界上导入新生涡满足边界的无滑移条件. 利用拉格朗日法计算涡元的运动, 而离散涡元之间的诱导速度则通过毕奥-萨伐尔公式结合快速多级子算法和自适应树结构算法获得. 通过对某混流式转轮水轮机在设计工况及非设计工况下的数值模拟, 以及模拟结果与实验测量的对比, 验证了该拉格朗日涡方法在计算转轮内部高雷诺数非定常流动中的有效性.
2012, (3): 520-527. doi: 10.6052/0459-1879-2012-3-20120308
基于拓扑优化技术的集中力扩散结构设计
牛飞, 王博, 程耿东
对于运载飞行器一类由多个部段装配而成的复杂结构, 外载荷往往被转化为多点集中载荷传递到结构某一部段, 为此, 结构上需要采取使集中力扩散的措施. 基于连续体拓扑优化, 以结构最小柔顺性为目标函数, 同时考虑设计域内材料用量约束和考查区域(承载结构)内力均匀性约束, 提出了集中力扩散结构优化设计的理论模型, 并给出了平面和三维两种不同的优化算例, 均得到了合理的优化设计结果. 最后, 针对运载火箭燃料贮箱短壳, 考虑工程实际的受力及约束情况, 给出了可行的概念性设计方案.
2012, 44(3): 528-536. doi: 10.6052/0459-1879-2012-2-20120309
热力耦合结构的弹性支撑分析与拓扑优化设计
杨军刚, 张卫红, 王丹, 蔡守宇
弹性支撑是保证热力耦合载荷作用下结构有效承载的一种支撑设计方案, 可以在满足结构刚度设计要求的同时有效降低结构内的应力集中, 保证结构的热变形协调. 从理论分析与数值计算两个方面, 研究了热力耦合载荷作用下结构的弹性支撑优化设计. 首先以弹性支撑梁模型为例, 通过给出热力耦合载荷下应力计算的解析表达式, 从理论上说明弹性支撑对结构应力峰值的影响, 阐述了弹性支撑的承载与热变形协调的双重效果. 在此基础上, 提出了解决该类问题的弹性支撑通用优化数学模型和拓扑构型优化方法, 通过数值算例与优化结果展示了方法的有效性.
2012, (3): 537-545. doi: 10.6052/0459-1879-2012-3-20120310
模糊分布参数条件下结构系统的近似效应分析
王攀, 吕震宙, 唐樟春
针对结构系统中随机变量的分布参数具有不确定性的问题, 采用模糊数来描述分布参数的主观不确定性, 并且利用隶属度来衡量主观不确定性的大小. 在不同隶属度水平下, 分别建立了输出响应和失效概率的主效应和总效应模型, 研究了两种模型和隶属度水平之间的关系, 并且给出分布参数不确定性的平均效应指标. 针对所提指标计算量大的问题, 通过对输出响应的各阶矩进行一阶泰勒展开, 建立了两种指标的计算方法. 算例结果验证了所提指标的合理性和所提方法的高效性.
2012, (3): 546-556. doi: 10.6052/0459-1879-2012-3-20120311
多层地基条带基础动力刚度矩阵的精细积分算法
林皋, 韩泽军, 李伟东, 李建波
提出应用精细积分算法计算多层地基的动力刚度问题. 精细积分是计算层状介质中波传播的高效而精确的数值方法. 利用傅里叶积分变换将层状地基的波动方程转换为频率-波数域内的两点边值问题的常微分方程组, 运用精细积分方法求解格林函数, 最后再将得到的频率-波数域内地基表面的动力刚度矩阵转换到频率-空间域内, 进而得到刚性条带基础频率域的动力柔度或刚度矩阵. 所建议的精细积分算法, 可以避免一般传递矩阵计算中的指数溢出问题, 对各种情况有广泛的适应性, 计算稳定, 在高频段可以保障收敛性, 并能达到较高的计算精度.
2012, (3): 557-567. doi: 10.6052/0459-1879-2012-3-20120312
弹性力学的杂交自然单元法
董轶, 马永其, 冯伟
为了解决自然单元法无法直接求解节点应力以及应力解精度不高的问题, 将应力杂交的思想引入自然单元法中, 与弹性问题的Hellinger-Reissner变分原理结合, 提出了弹性问题的杂交自然单元法. 算例表明: 杂交自然单元法的计算结果与解析解吻合, 证明该方法是可行的; 在求解应力方面, 杂交自然单元法比自然单元法有着更高的计算精度, 而且可以直接求解出节点的应力.
2012, (3): 568-575. doi: 10.6052/0459-1879-2012-3-20120313
快慢耦合振子的张驰簇发及其非光滑分岔机制
张晓芳, 陈小可, 毕勤胜
通过引入适当的参数值, 得到了两时间尺度下的快慢耦合振子, 分析了耦合系统及子系统的平衡点及其性质, 进而利用微分包含理论, 探讨了非光滑分界面上的奇异性, 指出在适当的参数条件下, 系统轨迹在穿越分界面时会产生由Hopf分岔和Fold分岔组合的非常规分岔. 给出了不同参数条件下的周期簇发行为, 分析了簇发过程的振荡特性, 指出激发态的频率取决于快子系统在非光滑分界面上的Hopf分岔频率, 而慢子系统的固有频率影响了簇发行为的振荡周期, 并进一步揭示了由非光滑分岔引起的不同周期簇发的分岔机制.
2012, (3): 576-583. doi: 10.6052/0459-1879-2012-3-20120314
SD振子的设计及非线性特性实验研究初探
陈恩利, 曹庆杰, 冯明, 田瑞兰
具有光滑与不连续转迁特征的SD振子发现和提出以来, 引起了广泛关注. 基于双稳系统大位移特征的测量法困难, SD振子的实验研究还未见报道. 该文提出并设计了具有SD振子系统光滑特征的非线性实验装置, 用实验的方法揭示由几何关系产生的强非线性系统的非线性动力学行为. 设计的非线性实验装置基本振动参数均有良好的可调性和可测量性, 对SD振子在不同频率及幅值的简谐激励作用下的非线性动力学响应进行了实验研究. 为克服大位移测量难题, 研究采用高速摄像机采集振子振动视频信号并进行分析. 结果表明, SD振子系统在一定的参数条件下会产生周期振动、周期5振动及混沌运动等复杂非线性动力学现象, 在相同实验参数条件下进行了数值仿真, 仿真结果与实验结果一致.
2012, (3): 584-590. doi: 10.6052/0459-1879-2012-3-20120315
脑Willis环的一维血流动力学及氧输运特性的数值研究
姬长金, 贺缨
Willis环是大脑侧枝循环的重要组成部分, 研究其血流动力学特性以及氧输运规律对脑缺血疾病的认知和预防有着非常重要的作用. 该文旨在利用一维血流动力学模型模拟整个Willis环的流量变化和压力分布, 并建立动脉内氧输运的一维模型以模拟Willis 环内氧分压的变化规律, 为研究脑组织内血液流动和氧输运打下基础. 首先, 基于弹性圆管内的一维非线性流动方程和状态方程建立血流动力学模型, 在一维对流扩散方程的基础上, 考虑由管腔向壁面的扩散和壁面细胞的新陈代谢消耗推导出氧输运特性方程. 通过 Lax-Wendroff两步法对血流动力学方程进行离散, 而在进行对流扩散方程的离散时, 则运用迎风格式. 通过数值计算得到了正常情况下Willis环各个血管任意位置的流量、压力和氧分压的变化曲线, 正常情况下各个位置的氧分压处于稳定的平衡状态. 最后, 还通过此模型进一步模拟了右侧颈内动脉狭窄对各个血管内流动的影响. 当狭窄程度达到80%时, 中脑动脉的流量和压力会明显下降, 造成其供应区域的血流减少. 同时, Willis环右侧血管内的氧分压会大大降低, 而左侧血管的氧分压会出现上升趋势, 但幅度要小于右侧血管降低的幅度.
2012, (3): 591-599. doi: 10.6052/0459-1879-2012-3-20120316
生物芯片中周期性电渗驱动液体薄膜的流动特性
唐文跃, 胡国辉
研究了二维周期性电渗驱动液体薄膜的流动特性. 以Debye-Hückel 假设近似下线性化的Poisson-Boltzmann方程描述双电层电动势分布和电荷密度的分布关系, 与黏性不可压缩流体Navier-Stokes方程相耦合, 得到流体在自由面与固壁之间的周期电渗流流场的精确解. 结果显示, 薄膜内速度振幅与流体黏性密切相关, 雷诺数越大, 速度振幅就越小. 该文还细致分析了雷诺数和自由面ζ电势对自由面的流速振幅和薄膜内速度相位差的影响.
2012, (3): 600-606. doi: 10.6052/0459-1879-2012-3-20120317
鼻腔加温功能特征及其与气流场关系的研究
苏英锋, 孙秀珍, 刘迎曦, 于申, 王吉喆, 高飞
为研究鼻腔加温功能特征及其与气流场之间的关系, 选用1例健康国人的鼻腔进行CT扫描. 据CT数据对鼻腔气道进行表面三维重建, 运用计算流体动力学方法分析通气量为12L/min时吸气相0.15s, 0.45s, 0.75s的鼻腔气流场与温度场. 结果显示吸气相0.15s, 0.45s, 0.75s鼻腔气流场主要表现为双侧气流量分布不对称, 其中气流主要流经左侧; 双侧均为总鼻道中、下部气流量较多, 嗅裂、中鼻道和下鼻道气流量少. 吸气相0.15s, 0.45s, 0.75s温度场均表现为温度自鼻腔前端至鼻咽部逐渐增高, 其中温度主要上升区域为鼻内孔-下鼻甲前端-中鼻甲前端对应气道, 且在吸气速度和流量增大后, 这一主要加温区域无明显向后延长征象; 双侧鼻腔及单侧鼻腔不同部位气道气流分布差别较大, 但双侧温度场基本对称, 单侧鼻腔不同气道部位温度差值亦较小, 幅度均小于1℃.
2012, (3): 607-613. doi: 10.6052/0459-1879-2012-3-20120318
血红细胞力学性能的纳米压痕实验和有限元模拟
李卫卫, 杨庆生, 刘志远
采用纳米压痕技术和有限元方法研究了血红细胞的生物力学性能. 进行了血红细胞的纳米压痕实验, 得到了血红细胞的材料参数和变形形貌; 在实验基础上, 建立了血红细胞的三维有限元模型, 模拟了血红细胞的压痕载荷-位移曲线, 并考虑了参数效应. 数值模拟结果和实验数据符合很好. 通过改变压头与材料之间的摩擦系数和压头曲率半径等参数, 比较了载荷-位移曲线的变化情况. 研究表明摩擦系数对压痕载荷-位移曲线和应力分布影响很小, 而压头曲率对载荷-位移曲线的影响明显.
2012, (3): 614-621. doi: 10.6052/0459-1879-2012-3-20120319
蜗窗激励评价的有限元计算模型研究
王学林, 胡于进
蜗窗驱动是实现中耳听力装置与耳蜗耦合的新途径. 利用外耳道、中耳和耳蜗集成有限元模型, 分别模拟外耳道声激励经听骨链、前庭窗传入内耳的正向传递过程和蜗窗机械激励在耳中逆向传递过程, 计算获得中耳和耳蜗的传递函数. 比较蜗窗激励和外耳道激励下耳蜗基底膜的振动, 提出了以基底膜最佳反应部位位移相等为准则的蜗窗等效激励力计算方法. 理论计算获得的蜗窗等效激励力与相关文献根据实验数据预测结果一致. 计算结果还表明耳蜗窗正、逆向激励时, 基底膜上最佳反应部位无变化, 但在蜗窗逆向激励耳蜗时, 低频下驱动基底膜运动的效率比高频时低. 所获得的理论结果可为蜗窗驱动的听力装置设计和现有装置的应用提供参考.
2012, (3): 622-630. doi: 10.6052/0459-1879-2012-3-20120320
ACL单束!/!双束重建对多屈曲角位姿的胫股关节力学特性影响
徐强, 黄荣瑛, 许勇刚, 郭云飞, 郑红光
进行膝关节的轴移仿真, 在证实了单例个体膝关节的重建模型能反映不同重建术后膝关节力学特性有效性的基础上, 应用SONATA MAESTRO 1.5T扫描得到的屈曲角位姿为0°/25°/60°/80° 的膝关节MRI图像数据, 建立了对应屈曲角的正常/ACL(anterior cruciate ligament)单束/双束重建胫股关节3D模型, 通过对该系列关节模型施加轴向压力与压扭载荷, 来分析多屈曲角度位姿的ACL单束/双束重建法对胫股关节软骨、半月板和韧带应力分布影响及韧带张力特性影响. 结果表明: (1) ACL单束/双束重建关节的软骨与半月板上的应力分布改变跟屈曲角度相关, 某些角度位姿下软骨应力分布改变显著; (2) 单束重建关节的软骨和半月板上的最大等效应力有较明显的增大, 最大增幅达40%左右; 双束重建关节在各屈曲角位姿下软骨和半月板上的最大等效应力值更接近于正常关节; (3) ACL单束重建虽降低了PCL(posterior cruciate ligament)的最大等效应力值, 却使高屈曲角位姿的内侧/外侧副韧带等效应力明显增大; 双束重建后MCL(medial collateral ligament)上的最大等效应力随屈曲角的变化明显, 但LCL(lateral collateral ligament)和PCL的最大等效应力随屈曲角度的变化趋势与正常关节相一致性; (4) 单束/双束重建后MCL上平均张力要高于另外两条韧带, 而双束重建的韧带张力特性比单束重建更接近于正常关节. 总之, 综合ACL重建后软骨、半月板和3条韧带随屈曲角度的等效应力分布、张力变化等多种特性表明: ACL双束重建的胫股关节力学特性比单束重建更接近于正常关节, 且无论ACL单束还是双束重建, 都引起术后关节软骨与韧带上应力分布变化与等效应力峰值增大, 这将是诱发术后关节慢性退变与膝关节并发症的根源.
2012, (3): 631-637. doi: 10.6052/0459-1879-2012-3-20120321
研究简报
基于压缩蠕变试验的流变参数理论解析反演
杨文东, 张强勇, 陈芳, 李术才, 李文纲, 王建洪, 贺如平, 曾纪全
大坝坝基的长期稳定性是大坝运营安全的重要保障, 而坝基岩体中的软弱夹层是影响其变形和稳定的重要因素. 为研究大岗山水电站坝基中辉绿岩脉软弱夹层在长期载荷作用下的变形机制, 在坝基边坡的试验平硐内垂直于软弱夹层进行了现场大型圆形刚性承压板压缩蠕变试验. 辨识得到可较准确表示其蠕变特性的5参量广义Kelvin模型, 克服了3 参量广义Kelvin模型收敛过快的缺陷. 基于弹性力学中的布辛涅斯克问题, 通过黏弹性理论中的拉普拉斯变换及逆变换, 推导出了刚性承压板下部岩体的5参量广义Kelvin模型表示的黏弹性变形公式, 并以此为基础反演得到流变参数.
2012, (3): 638-642. doi: 10.6052/0459-1879-2012-3-20120322
轴向移动局部浸液单向板的1:3内共振分析
李红影, 郭星辉, 王延庆, 李健, 谢里阳, 陈波
考虑单向板的轴向速度、轴向张力、流固耦合作用以及阻尼等因素, 基于由 von Kármán薄板大挠度方程得到的轴向移动局部浸液单向板的非线性振动方程, 研究了外激励作用下单向板在1:3内共振情况时的非线性振动特性. 首先利用Galerkin法对非线性振动方程离散化, 然后分别应用数值法和近似解析法对离散后模态方程组进行求解, 获得了系统内共振情况下复杂的幅频特性曲线, 并讨论了周期解的稳定性. 最后研究了1:3内共振系统平均方程组的运动分岔现象.
2012, (3): 643-647. doi: 10.6052/0459-1879-2012-3-20120323
基于重要样本法的结构动力学系统的首次穿越
任丽梅, 徐伟, 肖玉柱, 王文杰
基于Gisranov定理, 提出一种估计稳态高斯白噪声激励的结构动力学系统首穿失效概率的重要样本法. 文章重点是构造控制函数, 控制函数促使随机响应尽量集中在样本空间中最易导致首次穿越发生的部分. 利用设计点构造控制函数, 在线性系统场合, 结合时不变系统的结构可靠性理论, 通过解有约束的优化问题得到设计点; 在非线性系统场合, 利用Heonsang Koo提出的设计点激励, 通过镜像法得到设计点. 最后给出例子, 将所提方法与原始蒙特卡罗法相比较, 模拟结果显示方法的正确性与有效性.
2012, (3): 648-652. doi: 10.6052/0459-1879-2012-3-20120324