EI、Scopus 收录
中文核心期刊

留言板

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

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

2015年  第47卷  第5期

显示方式:
研究论文
透水床面明渠湍流时空平均特征的大涡模拟与分析
韩旭, 何国建, 方红卫, 符松
基于大涡模拟数据,研究了理想粗糙透水床面明渠湍流的时空平均特性. 考虑到空间异构性,对比分析了不同位置的时空平均流速、雷诺剪应力、构造剪应力、脉动幅度的垂线分布. 结果表明:第一,顶层床面之上,空间异构性的影响较小,不同位置的双平均流速符合类似的对数分布,但由于透水床面影响,卡门常数较不透水床面小;在床面附近,空间异构性影响较大,不同位置的双平均流速分别符合线性分布与多项式分布;在透水河床内部,靠近底层球孔的双平均流速为上部球孔双平均流速的1.55 倍. 第二,床面之上,雷诺剪应力占总剪应力的95% 以上,占有主体地位;床面附近,紊动较大,构造剪应力不能忽略,其值大约占总剪应力的15%.由于流场的各向异性,纵向与垂向的脉动幅度有所差异.
2015, 47(5): 713-721. doi: 10.6052/0459-1879-14-382
长方体腔内关于密度极值温度对称加热-冷却时冷水瑞利-贝纳德对流稳定性
胡宇鹏, 李友荣
为了解具有密度极值流体瑞利-贝纳德对流特有现象和规律,利用有限容积法对长方体腔内关于密度极值温度对称加热-冷却时冷水瑞利-贝纳德对流的分岔特性进行了三维数值模拟,得到了不同条件下的对流结构型态及其分岔序列,分析了密度极值特性、瑞利数、热边界条件以及宽深比对瑞利-贝纳德对流的影响. 结果表明:具有密度极值冷水瑞利-贝纳德对流系统较常规流体更加稳定,且流动型态及其分岔序列更加复杂;相同瑞利数下多种流型可以稳定共存,各流型在相互转变中存在滞后现象;随着宽深比的增加,流动更易失稳,对流传热能力增强;系统在导热侧壁时比绝热侧壁更加稳定,对流传热能力有所减弱;基于计算结果,采用线性回归方法,得到了热壁传热关联式.
2015, 47(5): 722-730. doi: 10.6052/0459-1879-15-171
并列双圆柱流致振动的不对称振动和对称性迟滞研究
陈威霖, 及春宁, 徐万海
对雷诺数Re = 100 间距比s/D = 2.5 和5.0 的并列双圆柱流致振动进行了数值模拟研究, 其中圆柱质量比m = 2.0, 折合流速Ur 在2.0~10.0 之间, 两圆柱仅能做横流向振动. 研究发现, 当间距比s/D = 2.5 时, 在折合流速4.4 < Ur< 4.8区间内, 两圆柱流致振动响应出现不对称振动现象, 在折合流速4.4 < Ur< 4.8 区间内, 两圆柱流致振动响应出现对称性迟滞现象; 而当间距比s/D = 2.5时, 圆柱流致振动响应与单圆柱涡激振动响应相似, 没有出现不对称振动和对称性迟滞现象. 在不对称振动区间内, 两圆柱的升、阻力参数也出现了不相等的情况. 此外, 当两圆柱不对称振动时, 圆柱间隙流稳定地偏斜向其中的一个圆柱; 相应地, 尾涡也出现了宽窄不等的模式. 窄尾流圆柱的振幅和升、阻力均较宽尾流圆柱的大. 通过对比不对称振动现象发生前后的尾涡模式, 对新现象的产生机制进行了阐述.
2015, 47(5): 731-739. doi: 10.6052/0459-1879-15-007
二维方腔热对流系统中纳米颗粒混合及凝并特性的数值模拟
徐飞彬, 周全, 卢志明
采用泰勒展开矩方法对二维瑞利-贝纳德热对流系统(1×106 ≤Ra ≤1 ×108) 中纳米颗粒群的混合和凝并特性进行了数值模拟. 结果显示颗粒群随时间演化经历了扩散阶段、混合阶段、充分混合阶段3 个阶段, 随着颗粒群混合和凝并的进行, 颗粒数目浓度减少, 颗粒群的平均体积增大; 得到了颗粒分布函数各特征量与温度相关系数以及各特征量的空间分布标准偏差在3 个阶段的不同特征; 得到了颗粒分布函数各阶矩以及平均体积长时间演化的渐近行为, 结果与零维渐近解析解一致. 最后, 本文进一步研究了无量纲数(包括瑞利数Ra, 斯密特数ScM, 达姆科勒数Da) 对颗粒群达到自保持分布时间的影响, 发现该时间随着RaScM的增大呈对数率减小, 随着Da的增大呈线性增大
2015, 47(5): 740-750. doi: 10.6052/0459-1879-15-062
基于应变梯度中厚板单元的石墨烯振动研究
徐巍, 王立峰, 蒋经农
基于应变梯度理论建立了单层石墨烯等效明德林(Mindlin) 板动力学方程,推导了四边简支明德林中厚板自由振动固有频率的解析解. 提出了一种考虑应变梯度的4 节点36 自由度明德林板单元,利用虚功原理建立了单层石墨烯的等效非局部板有限元模型. 通过对石墨烯振动问题的研究,验证了应变梯度有限元计算结果的收敛性. 运用该有限元法研究了尺寸、振动模态阶数以及非局部参数对石墨烯振动特性的影响. 研究表明,这种单元能够较好地适用于研究考虑复杂边界条件石墨烯的尺度效应问题. 基于应变梯度理论的明德林板所获得石墨烯的固有频率小于基于经典明德林板理论得到的结果. 尺寸较小、模态阶数较高的石墨烯振动尺度效应更加明显. 无论采用应变梯度理论还是经典弹性本构关系,考虑一阶剪切变形的明德林板模型预测的固有频率低于基尔霍夫(Kirchho) 板所预测的固有频率.
2015, 47(5): 751-761. doi: 10.6052/0459-1879-15-074
含固支边矩形叠层厚板的状态空间新解法
胡文锋, 刘一华
通过把固支边上的边界位移函数作为状态变量引入状态方程,得到了含固支边矩形单层与叠层厚板的精确解. 在求解过程中,将非齐次状态方程的求解变为齐次状态方程的求解,省略了求解待定常数的中间过程,使求解过程变得简单. 所得到的解能够严格满足固支边界条件,在同一材料层内不需作分层处理,因而更加精确. 此外,对固支边的应力提出了新的计算方法,能够得到更精确的边界应力. 算例表明,本文解比现有精确解收敛快,与有限元解吻合的更好,尤其是在固支边处体现得更加明显.
2015, 47(5): 762-771. doi: 10.6052/0459-1879-15-033
涂层式裂纹监测系统中基体裂纹穿越行为研究
孙洋, 王彪, 王巧云, 刘马宝
建立了智能涂层的两相模型与三相模型,基于能量准则分别用这两种模型研究了基体裂纹达到涂层界面后的穿越/偏转行为. 用有限元法分析了相对裂纹扩展长度、弹性错配参数及界面层厚度对偏转裂纹与穿越裂纹能量释放率之比的影响,结果表明当基体裂纹到达驱动层与基体界面时,能量释放率之比不仅与基体和驱动层之间的弹性错配相关,而且当驱动层较薄时对驱动层与传感层之间的弹性错配亦有较强的依赖性. 此外,随着驱动层厚度的增加,能量释放率之比对驱动层与传感层之间的弹性错配的依赖性逐渐降低. 通过与实验结果相比,建立的模型能够较好的解释基体裂纹在界面的扩展行为,可用于智能涂层裂纹传感器的优化设计.
2015, 47(5): 772-778. doi: 10.6052/0459-1879-15-142
基于圆筒模型的热障涂层安定分析
徐颖强, 孙戬, 李万钟, 张一弛, 陈亚亚
热载荷作用下,由于热障涂层(thermal barrier coatings, TBCs) 各层材料的热不匹配以及材料参数的温度相关等因素,会使热障涂层界面区域存在复杂的应力应变场,影响系统安定性,并导致涂层开裂和剥落. 将热障涂层外凸和内凹微观界面结构简化为多层圆筒模型,借助经典机动安定定理,利用特雷斯卡(Tresca) 屈服准则和增量破坏准则处理对时间的积分问题,避免了常规安定性分析的数学规划问题,建立了热障涂层安定极限分析方法,将材料屈服强度随温度变化关系简化为双线性关系,利用补偿变换的方法简化求解过程,对典型热障涂层安定性进行了研究. 结果表明,利用基于圆筒的安定极限分析方法,能够方便求解安定极限,便于工程应用;热障涂层安定极限值明显高于弹性设计值,且界面外凸区域安定极限高于内凹区域极限值,结构首先在内凹处失效;圆筒模型基体曲率和涂层厚度越大,结构安定极限越高,分析结果与试验结果一致;所建立的热障涂层安定分析方法,对进一步研究考虑蠕变因素影响的热障涂层安定性具有重要意义.
2015, 47(5): 779-788. doi: 10.6052/0459-1879-15-073
大型液体火箭姿控与跷振大回路耦合动力学建模与分析
王庆伟, 谭述君, 吴志刚, 杨云飞, 陈宇
大型液体火箭结构模态的空间化分布特征导致结构振动、姿态运动和推进系统液路脉动存在相互耦合,进而影响传统姿控回路的稳定性. 针对大型液体火箭, 充分考虑姿态控制系统对箭体姿态动力学和弹性振动的影响, 以及箭体结构弹性振动与推进系统的耦合作用(跷振(POGO)), 建立了姿控与跷振大回路耦合模型. 该模型包含了推进系统、结构系统与姿控系统之间的耦合因素, 可进行姿控-结构-推进大回路耦合机理研究. 该模型具有非奇异的优点, 可以直接用于频域分析和时域仿真. 基于该模型研究了我国某型号液体捆绑火箭推进系统参数——泵增益和蓄压器能量值对姿态运动与结构振动稳定性的影响. 研究得出, 泵增益和蓄压器能量值的变化不仅导致了结构振动的不稳定, 而且也导致了姿态运动的发散. 因此, 对于大型液体捆绑火箭, 推进系统与姿控系统之间存在不可忽略的耦合作用, 在设计姿控系统时, 有必要考虑推进系统对姿控系统稳定性的影响.
2015, 47(5): 789-798. doi: 10.6052/0459-1879-15-137
航天器集群编队最优单脉冲机动
王伟, 袁建平, 罗建军
对航天器集群编队最优单脉冲机动问题进行了研究. 针对不同的任务约束,基于非线性相对运动的周期性条件,以解析的思路分别研究了机动时刻给定和机动时刻未定情况下集群编队的最优单脉冲机动问题. 对于机动时刻给定的情况,从高斯变分方程和基于能量匹配条件的拉格朗日乘子法两个角度分别进行了探讨,将问题转化为对一元二次方程求极值或对一个单零点非线性方程求根;对于机动时刻未定的情况,将问题转化为对一个多零点非线性方程求根,通过傅里叶-贝塞尔级数展开可以得到任意高阶近似解. 对于每种情况,推导得到二范数意义下能量最省对应的最优参考长半轴,以及所施加的最优速度脉冲. 数值仿真验证了本文方法的正确性,并对仿真结果进行了解释和分析.
2015, 47(5): 799-806. doi: 10.6052/0459-1879-14-386
1:1内共振对随机振动系统可靠性的影响
王浩宇, 吴勇军
研究了二自由度耦合非线性随机振动系统在高斯白噪声激励下基于首次穿越模型的可靠性问题. 在1:1内共振情形,原始系统的运动方程经平均后化为一组关于慢变量的伊藤随机微分方程. 建立了后向柯尔莫哥洛夫方程以及庞德辽金方程,在一定的边界条件和(或) 初始条件下求解这两个偏微分方程,分别得到系统的条件可靠性函数以及平均首次穿越时间. 进而建立了无内共振情形系统的后向柯尔莫哥洛夫方程与庞德辽金方程.将无内共振情形的结果与1:1 内共振情形的结果做比较,发现1:1 内共振能显著降低系统可靠性. 用蒙特卡罗数值模拟验证了理论结果的有效性.
2015, 47(5): 807-813. doi: 10.6052/0459-1879-15-058
含摩擦与碰撞平面多刚体系统动力学线性互补算法
王晓军, 王琪
基于接触力学理论和线性互补问题的算法, 给出了一种含接触、碰撞以及库伦干摩擦, 同时具有理想定常约束(铰链约束) 和非定常约束(驱动约束) 的平面多刚体系统动力学的建模与数值计算方法. 将系统中的每个物体视为刚体, 但考虑物体接触点的局部变形, 将物体间的法向接触力表示成嵌入量与嵌入速度的非线性函数,其切向摩擦力采用库伦干摩擦模型. 利用摩擦余量和接触点的切向加速度等概念, 给出了摩擦定律的互补关系式; 并利用事件驱动法, 将接触点的黏滞-滑移状态切换的判断及黏滞状态下摩擦力的计算问题转化成线性互补问题的求解. 利用第一类拉格朗日方程和鲍姆加藤约束稳定化方法建立了系统的动力学方程, 由此可降低约束的漂移, 并可求解该系统的运动、法向接触力和切向摩擦力, 还可以求解理想铰链约束力和驱动约束力. 最后以一个类似夯机的平面多刚体系统为例, 分析了其动力学特性, 并说明了相关算法的有效性.
2015, 47(5): 814-821. doi: 10.6052/0459-1879-15-168
气体动理学统一算法的隐式方法研究
毛枚良, 江定武, 李锦, 邓小刚
目前的气体动理学统一算法(unified gas kinetic scheme, 简称UGKS) 在求解高速流动问题时的计算效率,难以满足求解复杂工程问题的需求. 为了提高该算法的计算效率, 本文对模型方程的对流项和碰撞项进行了隐式处理, 并针对UGKS 界面通量与演化时间相关的特点, 引入了演化时间平均界面通量, 通过对控制方程矩阵进行近似LU 分解(lower-upper decomposition), 实现了隐式UGKS. 不同来流马赫数的圆柱绕流算例测试表明, 只要演化时间选取得当, 隐式方法可以得到与显式方法完全相同的结果, 且计算效率可以提高1~2 个量级.
2015, 47(5): 822-829. doi: 10.6052/0459-1879-14-408
使用时间插值的特征线法在光滑粒子法拉氏边界条件中的应用
宫翔飞, 张树道, 杨基明
在计算流体动力学的实际应用中无反射边界条件是一个重要的研究课题. 文中应用时间插值的特征线方法构造了一种新式的拉氏边界条件,并应用在光滑粒子法中. 该方法与使用特征线方法的欧拉边界条件方法相比,不需要区分超声速和亚声速流的不同,并且在入流和出流中具有相同的形式,因而更加简便易行. 数值结果表明,采用时间插值特征线法的拉氏边界条件方法在稀疏波、激波以及爆轰波的模拟中都能够得到较好的无反射效果.
2015, 47(5): 830-838. doi: 10.6052/0459-1879-14-412
光滑节点插值法:计算固有频率下界值的新方法
杜超凡, 章定国
将光滑节点插值法用于悬臂梁的静力学,并首次用于旋转柔性梁的频率分析. 采用梯度光滑技术,用线性插值形函数描述梁的位移场,求解4 阶微分方程. 在静力学分析中,将该方法所得梁中各点位移与假设模态法、有限元法及解析解的结果对比,可知该方法虽用简单的线性插值形函数描述梁的位移场,但精度却很高. 进一步研究表明,采用模态高于9 阶的假设模态法会使刚度阵条件数变差,导致结果发散. 在频率分析中,与有限元法、假设模态法和解析解对比,表明该方法一个重要特性:能提供固有频率的下界值,而有限元法和假设模态法只能提供固有频率的上界值,说明该方法结合有限元法在处理无解析解的问题时可以从上下界最大程度的逼近真实解,提高精度. 光滑节点插值法具有形函数结构简单、独立变量少且能提供固有频率下界值的特性,因此,具有较高的推广及应用价值.
2015, 47(5): 839-847. doi: 10.6052/0459-1879-15-146
数字图像相关法测量局域变形场中形函数和模板尺寸的影响
徐小海, 苏勇, 蔡玉龙, 程腾, 张青川
数字图像相关法测量均匀变形场已被普遍接受, 其测量结果可与应变片测量结果比较. 然而, 在工程测量中, 针对局域变形场(应变高度集中, 如波特文-勒夏特利埃带、试件缺口附近和裂纹尖端等), 应变片受限于其尺寸, 其测量结果是接触面内的平均应变值. 此时, 采用数字图像相关法能够测量这些局域变形场. 但形函数和模板尺寸等计算参数对计算结果影响很大, 这也导致使用者很难判断计算结果的可靠性. 论文通过对合金拉伸实验获得的不同应变梯度的波特文-勒夏特利埃带和模拟生成的带的计算分析, 发掘了形函数和模板尺寸作用于计算结果的深层机制, 证明了二阶形函数比一阶形函数更适用于高度非均匀的局域变形场. 提出了在局域应变场测量中, 当一阶和二阶形函数计算结果的相对误差小于10% 时, 二阶形函数的结果是可靠的判据.
2015, 47(5): 848-862. doi: 10.6052/0459-1879-15-119
研究简报
临界产量Dupuit公式的讨论及一种新的方法
韩国锋, 陈方方, 刘曰武, 朱永峰, 马小平
临界产量的确定对底水油藏的高效开发起着重要作用. 然而,临界产量Dupuit 公式反映临界产量随着打开程度的降低而单调增加,存在不合理的地方. 本文与一些学者的研究表明,Dupuit 公式推导过程暗含的“当临界状态时水锥正好处于井底” 的假设不合理,水锥不能稳定在井底位置. 在不作水锥到达井底的假设,假设油井上部为径向流,井底周围为半球形流动的情况下,推导出了一种新的底水临界产量方程. 算例与工程实例表明,该方法比较符合实际,能够求得临界水锥高度和对打开程度进行优化.
2015, 47(5): 863-867. doi: 10.6052/0459-1879-15-043
功能梯度材料动态断裂力学的径向积分边界元法
高效伟, 郑保敬, 刘健
采用径向积分边界元法分析功能梯度材料动态断裂力学问题. 该方法使用与弹性模量无关的弹性静力学开尔文基本解作为问题的基本解,在导出的边界-域积分方程中含有由材料的非均质性和惯性项引起的域积分,通过径向积分法将域积分转化为等效的边界积分,得到只含边界积分的纯边界积分方程;从而建立只需边界离散的无内部网格边界元算法. 采用候博特方法求解关于时间二阶导数的系统离散的常微分方程组. 最后通过数值算例验证本文方法的精度和有效性.
2015, 47(5): 868-873. doi: 10.6052/0459-1879-15-150
弹性支撑圆柱绕流稳定性分析
李新涛, 张伟伟, 蒋跃文, 叶正寅
基于CFD 技术,采用系统辨识方法,建立了亚临界雷诺数(Re < 47) 下绕圆柱流动的非定常气动力模型(reduced order model, ROM). 耦合结构运动方程和降阶气动力模型,建立了弹性支撑圆柱绕流的稳定性分析模型. 算例分析了亚临界雷诺数下,结构固有频率、质量比等参数以及支撑方式对弹性系统稳定性的影响. 对于单自由度横向支撑圆柱,当结构固有频率趋近流动最不稳定模态频率时,弹性系统会在一定频率范围内失稳,这种现象最低可在Re~20 时出现. 旋转自由度的释放能够进一步降低系统的稳定性,可将临界雷诺数进一步降低至18 左右. ROM 方法不仅具有很高的效率,而且清晰地指出了弹性系统失稳的根本原因:流动模态和结构模态耦合作用导致结构模态失稳所致. 因此,失稳状态下系统振荡频率锁定于结构固有频率. 基于ROM 技术预测的失稳边界与直接CFD/CSD 仿真结果吻合,证明了该方法的正确性和精度.
2015, 47(5): 874-880. doi: 10.6052/0459-1879-14-307
科学基金
2015年度力学科学处面上项目、青年科学和地区科学基金资助情况介绍
詹世革, 张攀峰, 许向红, 孙中奎
对2015 年度国家自然科学基金委员会数理科学部力学科学处面上项目、青年科学和地区科学基金资助情况进行了简要介绍,给出了资助项目清单.
2015, 47(5): 881-898.