EI、Scopus 收录
中文核心期刊

留言板

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

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

2018年  第50卷  第4期

流体力学
基于气液相界面捕捉的统一气体动理学格式
王昭, 严红
气液相界面运动的研究无论是在科学还是工程领域都是非常重要的. 其中, 非平衡流动的计算尤其受到关注. 基于此, 我们构造了捕捉气液相界面的统一气体动理学格式. 由于统一气体动理学格式将自由输运和粒子碰撞耦合起来更新宏观物理量和微观分布函数, 故而可以求解非平衡流动. 具体思路是, 通过将范德瓦尔斯状态方程所表达的非理想气体效应引入统一气体动理学格式之中来捕捉气液相界面, 两相的分离与共存通过范德瓦尔斯状态方程描述. 由于流体在椭圆区域是不稳定的, 因此气液相界面可以通过蒸发和凝结过程自动捕捉. 如此, 一个锋锐的相界面便可以通过数值耗散和相变而得到. 利用该方法得到麦克斯韦等面积律(Maxwell construction)对应的数值解, 并与其相应的理论解相比较, 二者符合良好. 而后, 通过对范德瓦尔斯状态方程所描述的液滴表面张力进行数值计算, 验证了Laplace定理. 此外, 通过模拟两个液滴的碰撞融合过程, 进一步证明了该格式的有效性. 但是, 由于范德瓦尔斯状态方程的特性, 其所构造的格式仅适用于液/气两相密度比小于5的情况.
2018, 50(4): 711-721. doi: 10.6052/0459-1879-17-364
高温气体热化学反应的DSMC微观模型分析
杨超, 孙泉华
热化学耦合的非平衡现象一直是高温气体热化学问题研究的难点, 制约了诸如爆轰波胞格结构、低温点火速率等现象的分析. 本文以高温氮气离解和氢氧燃烧中的链式置换反应为例, 从微观反应概率、振动态指定的反应速率、热力学非平衡态的宏观反应速率、碰撞后的能量再分配等角度, 分析了直接蒙特卡罗模拟中的典型化学反应模型(TCE, VFD, QK模型)的微观动力学性质. 研究发现, 无论是高活化能的高温离解反应还是低活化能的链式置换反应, 实际参与反应的分子的振动能概率分布都偏离了平衡态的Boltzmann分布, 包含较强振动能额外影响的VFD模型可以很好地模拟高温离解反应, 而TCE (VFD的一个特例)和QK模型对活化能较低的链式置换反应的预测效果相对更好. 此外, 化学反应碰撞后的能量再分配应遵循微观细致平衡原理, 细微的偏差都可能造成平动能和振动能难以达到最终的平衡状态. 直接蒙特卡罗模拟的应用评估结果表明, 化学反应的振动倾向对热化学耦合过程产生了明显的影响, 特别是由于高振动能分子更多地参与了化学反应, 气体平均振动能的下降将影响后续化学反应的进行.
2018, 50(4): 722-733. doi: 10.6052/0459-1879-18-056
单窄条控制件对横向振荡柱体尾流 2P模式 旋涡脱落的改变
曹梦圆, 金华斌, 邵传平
横向强迫振荡柱体尾流控制是柱体涡激振动控制的基础,在海洋、土木等工程中具有重要意义. 横向强迫振荡柱体尾流中存在一种锁频旋涡脱落模式,即在一个振荡周期内柱体上、下侧各脱落旋转方向相反的一对涡,称为2P模式. 本文将相对宽度b/D=0.32的窄条控制件置于横向强迫振荡柱体下游,对振幅比A/D=1.25, 无量纲振频f_e D/V_∞=0.22,雷诺数Re=1 200的2P模式旋涡脱落进行干扰,并通过改变控制件位置,研究旋涡的变化规律. 采用二维大涡模拟和实验验证方法进行研究,在控制件位置范围0.8≤X/D≤3.2, 0.4≤Y/D≤3.2内,得到了2P, 2S, P+S和另外6种新发现的旋涡脱落模式,并对各模式旋涡的形成过程作了详细描述. 在控制件位置平面上给出了各旋涡模式的存在区域,画出了旋涡脱落强度的等值线图,并发现在一个相当大的区域内,旋涡脱落强 度可减小一半以上,尾流变窄. 发现柱体大幅振荡引起的横向剪切流在旋涡生成中起关键作用. 探讨了控制件对横向剪切流的影响,分析了控制件在每种旋涡模式形成中的作用机制.
2018, 50(4): 734-750. doi: 10.6052/0459-1879-18-005
风力机叶片翼型动态试验技术研究
李国强, 张卫国, 陈立, 聂博文, 张鹏, 岳廷瑞
风力机叶片动态振荡过程往往伴随着俯仰和横摆同时进行, 以前对许多动态问题不清楚的阶段, 工程上不惜以增加叶片重量为代价而采用偏安全的设计, 通常忽略横摆振荡的影响; 大型风力机设计对获取翼型更加全面、准确的动态载荷提出了更高要求, 研究横摆振荡对翼型动态气动特性的影响规律具有重要意义. 本文首次开展翼型横摆振荡动态风洞试验研究, 采用“电子凸轮”技术代替机械凸轮实现了振荡频率和振荡角度的无级变化, 基于设计的电子外触发装置实现了对动态流场的实时测量, 实现了风洞来流、模型角位移和动态压力数据的同步采集, 分别开展了翼型静态测压、俯仰/横摆动态测压、粒子图像测速和荧光丝线等试验研究, 试验结果准度较高、规律合理; 分析了动态试验洞壁干扰影响机制. 研究表明, 横摆振荡翼型的气动曲线也存在明显迟滞效应; 随着振荡频率升高, 翼型俯仰和横摆振荡下的气动迟滞性均增强; 翼型俯仰振荡正行程的动态失速涡破裂有所延迟; 洞壁与模型端部交界处的强三维效应对翼型压力分布影响较大; 建立的横摆振荡试验技术可为风力机动态掠效应的研究提供技术支撑.
2018, 50(4): 751-765. doi: 10.6052/0459-1879-18-108
低雷诺数下串列三圆柱涡激振动中的弛振现象及其影响因素
陈威霖, 及春宁, 许栋
对间距比为1.2和雷诺数为100的串列三圆柱涡激振动进行数值模拟, 发现在某个折合流速之后, 三圆柱的响应均呈现为随着折合流速增大而增大的弛振现象, 平衡位置偏移、低频振动以及旋涡脱落与圆柱运动之间的时机三个因素共同决定了弛振现象的出现. 进一步的研究发现, 串列三圆柱的弛振现象仅出现在质量比不大于2.0和雷诺数不大于100的工况下. 当质量比较大时, 串列三圆柱的平衡位置固定不变, 且圆柱的振动不规律, 使得旋涡脱落与圆柱运动的时机处于变化之中. 当雷诺数较高时, 最上游圆柱的平衡位置在折合流速较大时回到初始位置, 不再参与对圆柱振动的调节, 使得圆柱的振动响应不再规律, 旋涡脱落与圆柱运动的时机也一直处于变化之中.
2018, 50(4): 766-775. doi: 10.6052/0459-1879-18-057
等离子体激励器诱导射流的湍流特性研究
张鑫, 黄勇, 阳鹏宇, 唐坤, 李华星
为了进一步掌握等离子体流动控制机理, 完善等离子体激励器数学模型, 提升等离子体激励器扰动能力, 采用粒子图像测速技术, 在静止空气下开展了介质阻挡放电等离子体激励器诱导射流特性研究. 实验时, 将非对称布局激励器布置在平板模型上, 随后将带有激励器的模型放置在有机玻璃箱内, 从而避免环境气流对测试结果的影响. 基于激励器诱导流场, 分析了激励电压对诱导射流特性的影响, 揭示了较高电压下诱导射流近壁区的拟序结构, 获得了卷起涡、二次涡等拟序结构的演化发展过程, 计算了卷起涡脱落频率, 阐述了卷起涡与启动涡的区别, 初步探索了卷起涡的耗散机制. 结果表明: (1)层流射流不能完全概括等离子体诱导射流特性, 激励电压是影响射流特性的重要参数. 当电压较低时, 诱导射流为层流射流; 当电压较高时, 诱导射流的雷诺数提高, 射流剪切层不稳定, 层流射流逐渐发展为湍流射流. (2)等离子体诱导湍流射流包含着卷起涡、二次涡等拟序结构; 在固定电压下, 这些涡结构存在恒定的卷起频率. (3)当激励电压较高时, 流动不稳定使得卷起涡发生了拉伸、变形, 引起了流场湍动能增大, 从而加速了卷起涡的耗散. 研究结果为全面认识激励器射流特性, 进一步挖掘激励器卷吸掺混能力, 提升激励器控制能力积累基础.
2018, 50(4): 776-786. doi: 10.6052/0459-1879-17-392
固体力学
功能梯度材料动力学问题的POD模型降阶分析
郑保敬, 梁钰, 高效伟, 朱强华, 吴泽艳
为了快速分析非均质材料结构在复杂载荷作用下的动态响应, 提出一种模型降阶方法, 只需计算结构在简单均质材料情况下的动力学问题, 进而用其计算结果对非均质材料结构进行分析. 首先, 采用结构内部任意一点处的材料参数值作为整个结构的材料参数, 利用有限元分析软件计算该均质材料结构在动态载荷作用下的位移场建立数据库, 该数据库包含计算模型各个节点(自由度为L)在某时间段内L个时刻的位移; 其次, 对数据库中的信息按照时间离散的特定方式组集成瞬像矩阵, 并利用特征正交分解方法对其进行分解, 得到该模型的H个特征正交基底, 选取其中能反应模型主要特征的H<L?N个(其中~)作为一组最优基底, 通过这组基底建立模型的低阶离散控制方程; 最后, 求解低阶离散微分方程组, 得到功能梯度材料结构在复杂载荷作用下的位移场. 文中分别给出二维和三维算例, 比较了降阶模型和全阶模型计算结果, 验证了该方法的有效性, 并且计算效率能提高1 1)2个数量级.
2018, 50(4): 787-797. doi: 10.6052/0459-1879-18-069
数字状张拉整体结构的构型设计与力学性能模拟
朱世新, 张立元, 李松雪, 张勃洋, 张清东
针对大型张拉整体结构的设计问题,选取四棱柱状张拉整体结构和截角正八面体状张拉整体结构作为基本胞元,采用节点连接节点的方式建立球柱组合式数字状张拉整体结构,并使用基于结构刚度矩阵的大变形非线性数值求解方法对其进行力学性能分析.在两类胞元满足各自的自平衡条件和稳定性条件的前提下,组合得到的数字状张拉整体结构亦处于自平衡稳定状态,搭建了实物模型进行验证.以数字8状张拉整体结构为例,模拟研究了结构承受自重等分布载荷和单轴拉压等端部载荷时的静力学响应,以及结构无阻尼振动时的固有频率和模态等动力学性能.结果表明,结构在自重作用下的变形行为受初始预应力、压杆密度和拉索刚度的影响较大,对其进行合理配置方可确保结构具有足够刚度抵抗自重;结构在单轴拉压作用下呈现非线性的载荷-位移曲线,拉伸刚度随变形量的增大而增大,压缩刚度随变形量的增大而减小;结构的固有频率随初始预应力的增大而增大,而模态振型基本不变.研究结果丰富了大型张拉整体结构的外形种类,有望推动此类结构在土木建筑、结构材料等领域的应用.
2018, 50(4): 798-809. doi: 10.6052/0459-1879-18-043
非常规态型近场动力学热黏塑性模型及其应用
王涵, 黄丹, 徐业鹏, 刘一鸣
在非常规态型近场动力学(non-ordinary state-based peridynamics, NOSB-PD) 理论框架下构建了考虑应变率效应、塑性硬化、热软化效应和材料断裂特征的非局部三维热黏塑性固体本构模型以及相应的非局部空间积分型数值算法, 并应用于金属类材料和构件在冲击载荷作用等工况下的高应变率热黏塑性变形与破坏分析. 通过对经典含初始裂纹Kalthoff-Winkler板冲击试验进行三维近场动力学模拟, 可得到裂纹的起裂角度、扩展路径、扩展速度以及裂纹扩展过程中靶板等效应力和温度分布, 所得结果与已有试验结果和其他数值方法结果吻合较好. 在此基础上, 应用该模型分析了不同冲击速度作用下金属靶板的变形与裂纹扩展过程, 结果表明: 该模型能较好地模拟不同冲击速度(应变率)情况下靶板的变形与破坏全过程. 随着冲击速度变化, 初始裂纹的起裂时间、扩展方向和扩展速度呈一定规律变化. 冲击速度越低, 起裂时间越晚(直至冲击速度低于某值时初始裂纹不扩展), 裂纹扩展速度峰值越低, 冲击过程中靶板温度峰值越低, 完全扩展所需时间越长.
2018, 50(4): 810-819. doi: 10.6052/0459-1879-18-113
PMMA膨胀环动态拉伸碎裂实验研究
李天密, 张佳, 方继松, 刘丽芝, 郑宇轩, 周风华
在强动载作用下, 脆性材料的碎裂问题是一个重要的研究课题, 而脆性材料在冲击拉伸载荷下的力学行为的实验研究相对较匮乏. 提出了一种动态拉伸断(碎)裂的液压膨胀环实验技术, 可用于准脆性/脆性材料的动态拉伸. 利用该技术对有机玻璃(PMMA)圆环试件进行了不同膨胀速度下的动态碎裂实验研究. 从回收碎片的断口形貌和碎片内部残余裂纹观察可知试件的破碎由环向拉伸应力造成, 碎片断口处发出的稀疏波会将周围的拉伸应力卸载, 从而抑制其他裂纹的进一步发展. 利用超高速相机记录了试件的膨胀碎裂过程, 利用DISAR激光速度干涉仪获得了试件外表面粒子的径向膨胀速度历史, 通过试件上的应变片获得了试件的应变历史和断裂应变. 实验结果表明: 在拉伸应变率150~500s-1范围, 材料的动态断裂应变低于准静态加载下的断裂应变, 体现出“动脆”现象; 随着加载应变率的提高, PMMA 材料的碎片尺寸减小; 无量纲化的PMMA圆环的平均碎片尺寸介于韧性碎裂模型和脆性碎裂模型的预测数值之间, 反映出材料的准脆性特性.
2018, 50(4): 820-827. doi: 10.6052/0459-1879-18-016
应力诱发界面迁移下晶内孔洞的演化
余文韬, 黄佩珍
随着微电子技术的迅猛发展, 集成电路中内连导线的失效问题引起广泛关注. 内连导线内部孔洞萌生、长大、漂移和失稳变形成狭长裂纹, 从而导致电路的开路失效. 这是内连导线失效的常见形式. 而界面迁移是导致微结构形态演化的主要机制之一. 本文基于界面迁移下微结构演化的经典理论和弱解描述, 建立了应力诱发界面迁移下微结构演化的有限单元法, 并验证了算法的可靠性. 对铜内连导线中晶内孔洞的演化进行了数值模拟, 详细分析了应力、线宽及形态比对晶内孔洞演化的影响. 研究结果表明, 椭圆形晶内孔洞存在生长和收缩两种演化分叉趋势. 通过大量数值分析得到了晶内孔洞演化的临界应力σ?c 、临界线宽h?c 和临界形态比βc . 当h??h?c, σ?h? 时, 晶内孔洞会沿长轴长大; 反之, 晶内孔洞会收缩甚至愈合. 此外, 应力β 越大、线宽σ? 越小或形态比h? 越大, 晶内孔洞越易发生长大, 且孔洞面积增大速度越快; β 越小、σ? 越大或h? 越小, 晶内孔洞越易发生收缩, 且孔洞面积减小速度越快.
2018, 50(4): 828-836. doi: 10.6052/0459-1879-18-015
基于迭代法的非线性弹性均质化研究
侯淑娟, 梁慧妍, 汪全中, 韩旭
微观结构对复合材料的宏观力学性能具有至关重要的影响, 通过合理设计复合材料微观结构可以得到期望的宏观性能. 均质化方法作为一种有效的设计方法, 它从微观结构的角度出发, 利用均匀化的概念, 实现了对复合材料宏观力学性能的预测和设计. 而当考虑非线性因素, 均质化的实现就非常困难. 本文利用双渐近展开方法, 将位移按照宏观位移和微观位移展开, 推导了非线性弹性均质化方程. 通过直接迭代法, 对非线性弹性均质化方程进行了求解, 并给出了具体的迭代方法和实现步骤. 本文基于迭代步骤和非线性弹性均质化方程编写MATLAB 程序, 对3种典型本构关系的周期性多孔材料平面问题进行了计算, 对比细致模型的应变能、最大位移和等效泊松比, 对程序及迭代方法的准确性进行了验证. 之后对一种三元橡胶基复合材料进行多尺度均质化, 将其分为芯丝尺度和层间尺度. 用线弹性的均质化方法得到了芯丝尺度的等效弹性参数, 并将其作为层间尺度的材料参数. 在层间尺度应用非线性弹性均质化方法对结构进行计算, 得到材料的宏观等效性能, 并以实验结果为基准进行评价.
2018, 50(4): 837-846. doi: 10.6052/0459-1879-18-039
多层复合壳体三维振动分析的谱--微分求积混合法
叶天贵, 靳国永, 刘志刚
对于较厚的多层复合壳体,其振动位移沿厚度方向呈锯齿形变化且层间剪切和拉、压应力呈三维耦合状态,采用传统的等效单层理论分析已不能满足精度要求. 建立不受结构厚度、铺层材料性质和铺层方式限制的三维分析方法具有重要的研究价值. 本文以独立铺层为建模对象,结合广义谱方法与微分求积技术建立了一种适用一般边界条件和铺层方式的多层复合壳体三维分析新方法——谱--微分求积混合法. 该方法应用三维弹性理论对独立铺层进行精确建模,有效克服了二维简化理论对横向变形以及层间应力估计不确切的缺点;引入微分求积技术对铺层进行数值离散,将三维偏微分问题转化为二维偏微分问题,降低了求解维度和难度;应用广义谱方法近似地表述离散计算面上的场变量,将获取的二维偏微分方程转化为以场变量谱展开系数为未知量的线性代数方程组,避免了对超越方程的求解. 数值验证结果表明该方法收敛性好,计算精度高.
2018, 50(4): 847-852. doi: 10.6052/0459-1879-18-083
基于S-R和分解定理的三维几何非线性无网格法
宋彦琦, 周涛
S-R(strain-rotation)和分解定理克服了经典有限变形理论的一些缺点, 使其可以为几何非线性数值分析提供可靠的理论基础. 对于大变形问题, 由于无网格法(element-free method)避免了对单元网格的依赖, 从而从根本上避免了有限单元法(finite element method, FEM)的单元畸变问题, 保证了求解精度. 因此, 将无网格法和S-R和分解定理结合起来势必能建立一套更加合理可靠的几何非线性数值计算方法. 目前基于S-R 定理的无网格数值方法研究较少并且只能用于二维平面问题的求解, 但实际上绝大多数问题都必须以三维模型来进行处理, 因此建立适用于三维情况的S-R无网格法是非常有必要的. 本文给出了适用于三维情况的S-R 无网格法: 采用由更新拖带坐标法和势能率原理推导出来的增量变分方程, 利用基于全局弱式的无网格Galerkin 法(EFG)得到了用于求解三维空间问题的离散格式. 利用MATLAB编制三维S-R 无网格法程序, 对受均布载荷的三维悬臂梁和四边简支矩形板结构的非线性弯曲问题进行了计算. 最后将所得的数值结果与已有文献进行了比较, 验证了本文的三维S-R无网格数值算法的合理性、有效性和准确性. 本文的三维S-R无网格数值算法可以作为一种可靠的三维几何非线性数值分析方法.
2018, 50(4): 853-862. doi: 10.6052/0459-1879-18-050
动力学与控制
多柔体系统数值分析的模型降噪方法
齐朝晖, 曹艳, 王刚
多柔体系统的动力学方程通常是一组刚性微分方程, 目前普遍采用的刚性微分方程数值解法主要通过数值阻尼滤除系统响应中的高频分量, 其求解效率难以令人满意. 为了降低多柔体系统动力学方程的刚性, 从而可采用ODE45等常规微分方程求解器进行求解, 研究了在建模过程中滤除高频振荡分量的方法. 在以当前时刻为起点的短时间内对柔性体的应力进行均匀化, 用均匀化后的应力计算柔性体的变形虚功率, 由此得到的系统动力学方程的解中不含过高频率的弹性振动, 并且可以通过调节均匀化时间区间的长度参数控制滤波的范围. 数值算例表明: 这种模型降噪方法的计算效率和精度均不低于刚性微分方程求解器, 并且在刚性微分方程求解器失效的情况下模型降噪方法仍有良好的精度和效率. 本文所提的模型降噪方法可成为求解多柔体系统动力学方程的新途径.
2018, 50(4): 863-870. doi: 10.6052/0459-1879-18-111
基于步态切换的欠驱动双足机器人控制方法
葛一敏, 袁海辉, 甘春标
由于高维、非线性、欠驱动等特点, 3-D双足机器人的稳定性控制依然是一个研究难点. 一些传统的控制方法, 如基于事件的反馈控制方法和PD控制方法, 抗扰动能力较弱, 鲁棒性较差. 通过观察, 人类受到外部扰动影响时, 会通过调整步态重新获得稳定性,相较之下仅依靠一个步态获得的稳定性是有限的. 受此启发, 本文针对上述问题提出一种基于步态切换的欠驱动3-D双足机器人控制方法. 首先, 以能耗最少为优化目标, 通过非线性优化方法预先设计多组不同步长、步速的步态作为参考步态, 以构建一个步态库; 然后, 通过综合考虑步态切换过程中的稳定性与能效, 建立了多目标步态切换函数; 最后, 将该步态切换函数作为优化目标, 并求解该最小化问题获得下一步的参考步态, 从而实现步态切换, 达到使用步态库?多轨迹方法来提高鲁棒性的目的. 在仿真实验中运用该步态切换控制方法, 欠驱动3-D双足机器人可实现相对高度在[-20, 20] mm内随机变化的不平整地面上行走, 而仅采用单步态控制策略则无法克服这样的外部扰动, 从而说明了基于步态切换的欠驱动双足机器人控制方法的有效性.
2018, 50(4): 871-879. doi: 10.6052/0459-1879-18-049
含时滞反馈与涨落质量的记忆阻尼系统的随机共振
公徐路, 许鹏飞
针对具有记忆效应的欠阻尼系统, 存在时滞反馈与涨落质量, 本文主要研究了其输出稳态响应振幅的随机共振效应. 首先通过引入新变量和运用小时滞近似展开理论, 将具有非马尔科夫特性的原系统转化为等价的两维马尔科夫线性系统, 再利用Shapiro-Loginov公式和Laplace变换获得了系统响应的一阶稳态矩和稳态响应振幅的解析表达式. 结果表明: 当系统参数满足Routh-Hurwitz稳定条件时, 稳态响应振幅随质量涨落噪声强度、周期驱动信号频率以及时滞的变化均存在随机共振现象, 其中随机多共振现象也被观察到. 在适当范围内, 通过控制时滞反馈, 系统的随机共振效应随着时滞的增大而增强, 而较长的记忆时间及增大阻尼参数均对共振行为呈现抑制作用.有效调控时滞反馈与记忆效应的变化关系将有助于增强系统对周期驱动信号的响应强度. 最后, 通过数值模拟计算验证了理论结果的有效性.
2018, 50(4): 880-889. doi: 10.6052/0459-1879-18-051
生物、工程及交叉力学
基于孔与裂隙网络模型的平行微裂隙对驱油的影响规律研究
刘海娇, 张旭辉, 鲁晓兵
认识双重多孔介质中油水两相微观渗流机制是回答形成什么类型的裂隙网络可提高油藏采收率的关键. 微裂隙的分布可以提高多孔介质的绝对渗透率,但对于基质孔隙中的流体介质,微裂隙的存在会引起多孔介质中局部流体压力和流场的变化,导致局部流动以微裂隙流动为主,甚至出现窜流现象,降低驱油效率. 本文基于孔与裂隙双重网络模型,在网络进口设定两条平行等长且具有一定间隔的微裂隙,分析微裂隙的相对间隔(微裂隙之间距离/喉道长度)和微裂隙相对长度(微裂隙长度/喉道长度)对于微观渗流特征的影响. 结果表明:随微裂隙相对长度的增加,出现驱油效率逐渐降低,相对渗透率曲线中的油水共渗区水饱和度和等渗点增加,油水两相的共渗范围减小等现象;随着微裂隙之间相对间隔增大,周围越来越多的基质孔穴间的压力差减小,在毛管压力的限制下,驱替相绕过这些区域,而导致水窜现象.
2018, 50(4): 890-898. doi: 10.6052/0459-1879-18-025
极端工况双矩形腔静压推力轴承动态特性
于晓东, 袁腾飞, 李代阁, 曲航, 郑旭航
静压推力轴承动态特性受润滑油黏度、油膜厚度和油腔面积等因素影响, 极端工况运行过程中经常承受阶跃载荷或正弦载荷作用, 突加载荷将导致静压推力轴承动态特性改变, 表现为轴承的抗冲击能力和恢复平衡所需时间的变化. 为获得高速重载微间隙极端工况条件下双矩形腔静压推力轴承动态特性, 分别在不同油膜厚度、不同润滑油黏度以及不同油腔尺寸条件下对双矩形腔静压推力轴承的动态性能进行理论分析, 探讨了阶跃载荷作用下润滑油黏度、油膜厚度和油腔面积对轴承动态性能的影响, 揭示了油膜动态变化规律, 探究了正弦载荷作用下双矩形腔静压推力轴承的稳定性. 结果表明: 润滑油黏度、油膜厚度和油腔尺寸变化对其动态性能有很大的影响. 润滑油黏度越大、油膜厚度越小、油腔面积越大突加载荷作用下润滑油膜抵抗冲击的能力越强, 旋转工作台受到突加外力作用下恢复至平衡状态所用时间越短. 双矩形腔静压推力轴承油膜具有较大的阻尼系数, 轴承具有极强的抵抗正弦加载作用的能力
2018, 50(4): 899-907. doi: 10.6052/0459-1879-18-041
一种考虑液相惯性的一维准静力固结模型
丁洲祥
经典Terzaghi一维固结理论不考虑孔隙流体惯性影响,且该理论在不同时期模型推导和表述结果差别较大,导致当前仍存在诸多困惑甚至认识混乱的现象. 在笔者前期研究大变形动力固结理论框架内,忽略固相惯性而重点考虑液相惯性影响,经过合理简化建立反映孔隙流体惯性的一维小变形固结波动模型. 该固结波模型具有频散和耗散特性. 采用分离变量法,可得到单面排水和瞬时加载条件下无量纲形式固结波解析解答. 算例分析结果表明:固结波发展规律受无量纲数Dc变化影响而呈现不同性态;Dc数值较大时固结波响应会出现阶跃和正负波动现象;当Dc值较小时,可能出现Mandel-Cryer效应等特殊现象. 通过对早期和后期Terzaghi固结模型的分析和对比,初步探明Terzaghi固结理论模型内部的矛盾性,在普通土体坐标和固相体积坐标两种不同解读条件下,早期Terzaghi (1923,1925)模型可以分别诠释为具有小变形和大变形属性的不同固结模型. 在经典一维固结理论模型的不同诠释背景下,固结波模型也可以据此作出相应拓展和表述. 固结波理论揭示缩尺固结试验中土体物理力学参数与固结波响应两种因素之间存在一种不确定性矛盾,据此建议微观土力学研究重视尺度效应. 固结波模型的意义还在于,可为Terzaghi经典固结模型理论精度分析提供新的依据.
2018, 50(4): 908-928. doi: 10.6052/0459-1879-18-053
复杂加载条件下的砂土本构模型
万征, 孟达
试验表明,饱和砂土的应力应变关系具有显著的密度以及压力依存性,上述两点构成了描述砂土静力加载下变形特性无法忽视的因素. 此外,在循环加载等复杂加载作用下,砂土还会表现出明显的应力诱导各向异性以及相变转换特性. 基于在e--p空间中存在唯一的临界状态线这一基本假定,通过在e--p空间中引入当前状态点与临界状态线的距离R来作为反映密度与压力依存特性的状态参量, 将变相应力比以及峰值应力比表达为状态参量的指数函数,将上述应力比参量引入到统一硬化参量中可准确地反映初始状态下围压、密度 对于单调加载下应力应变关系的影响规律,能描述砂土剪缩、剪胀,应变软化、硬化等特性. 采用非相关联流动法则,p--q空间中采用水滴型屈服面,塑性势面为椭圆面,松砂在单调加载下的静态液化现象也可描述. 为反映循环加载下塑性体积应变的累积特性以及塑形偏应变的滞回特性,在循环加载下将状态参量R表达为应力比参量,并在硬化参数中引入描述应力诱导各向异性特性的旋转硬化部分,所提模型可有效地描述循环加载下剪切模量的衰减特性、刚度衰化性质、强度减小特性,在不排水约束作用下,则会产生往返活动性现象. 通过一系列的模型模拟与试验结果对比,验证了本构模型的有效性及适用性.
2018, 50(4): 929-948. doi: 10.6052/0459-1879-18-047
结构优化半解析灵敏度及误差修正改进算法
张力丹, 张盛, 陈飙松, 李云鹏
提出结构半解析灵敏度分析及其针对刚体位移的误差修正方法的改进算法, 构建灵敏度分析与误差修正项可分离形式. 该方法实现简便, 数值精度不受摄动步长与单元数目的影响. 首先从总体角度推得静力问题的误差修正半解析灵敏度分析方法, 提出了位移误差修正灵敏度列式, 并给出算法实施途径; 然后将此思路推广于自振频率、屈曲临界载荷问题, 提出了相应的计算步骤. 随后, 给出梁单元与壳单元误差修正项的具体推导方法, 并分别使用两种单元构建有限元模型进行算例测试. 结果表明, 该方法适用于多种分析类型, 数值精度不受单元数目与摄动步长的影响. 由于灵敏度分析与误差修正项可以分开计算, 该方法支持将误差修正项直接叠加于灵敏度求解结果进行误差修正, 使已有灵敏度分析程序得到充分利用. 尤其对于复杂工程结构的优化设计, 特别是形状优化设计以及尺寸、形状混合优化设计, 相比于原误差修正方法, 实现更为简便, 效率有所提升, 能为半解析灵敏度分析方法及其程序实现提供新的思路.
2018, 50(4): 949-960. doi: 10.6052/0459-1879-18-058
一种模拟大规模高频声场的双层奇异边界法
李珺璞, 陈文
大规模高频声场的数值模拟是一项非常有计算挑战性的课题. 为了解决传统边界型离散方法由于全局支撑的满阵限制, 不易应用于大规模高频声场模拟的计算瓶颈, 本文提出了一种用于模拟大规模高频声场的双层奇异边界法. 在该方法中, 通过引入双层结构, 细网格上的全局支撑的满阵被转化为局部支撑的大规模稀疏矩阵, 传统奇异边界法模拟大规模问题时所面临的高计算量以及过度存储需求遂得以解决. 其次, 双层奇异边界法仅通过粗网格评估远场作用, 且独立于特定的插值核函数. 相较于快速多级方法, 该方法具有更强的适应性和灵活性, 且多层结构使该方法具有一定的预调节作用, 非常适合求解具有大规模、高秩、高条件数特点的高频波矩阵. 在其后的散射球模型算例中, 双层奇异边界法配置10万个节点, 成功模拟了无量纲波数高达160的声散射问题. 在对于人头模型的声散射特性分析中, 双层奇异边界法比COMSOL软件计算速度快了约78.13%. 当配置8万个节点时, 双层奇异边界法成功模拟了频率高达25 kHz 的工况, 该频率已远远超出了人耳的听力极限.
2018, 50(4): 961-969. doi: 10.6052/0459-1879-18-100
科学基金
第一届全国物理力学青年学者学术研讨会报告综述
朱云飞, 韩增尧, 姜利祥, 詹世革
简要介绍了第一届全国物理力学青年学者学术研讨会的情况,对与会学者所做的报告内容进行了总结.
2018, 50(4): 970-976. doi: 10.6052/0459-1879-18-244