Please wait a minute...

当期目录

    2021年 第53卷 第4期    刊出日期:2021-04-08
    上一期   
    综述
    原子力显微镜在二维材料力学性能测试中的应用综述
    高扬
    2021, 53(4):  929-943.  DOI: 10.6052/0459-1879-20-354
    摘要 ( 451 )   PDF (21904KB) ( 590 )  
    相关文章 | 计量指标
    以石墨稀为代表,二维材料有着诸多优异的性质,在下一代电子器件等领域拥有广阔的应用前景。目前绝大多数关于二维材料的研究都集中在其电子学和光学性质和应用,对于其力学性质的研究则相对欠缺。原子力显微镜是低维材料力学性质表征的主要手段,例如基于原子力显微镜的纳米压痕技术。在这篇综述中,作者首先简要介绍了二维材料的基本背景以及原子力显微镜的工作原理。作者进一步展示了纳米压痕技术的工作原理和理论背景,并回顾了利用纳米压痕技术研究二维材料面内力学性质的相关实验和理论工作,同时探讨了原子力显微镜在表征二维材料力学性能中存在的测量误差及来源。由于二维材料展现出强烈的各向异性,纳米压痕技术在能够很好地测量二维材料面内/横向的力学性质的同时,在二维材料层间/纵向的力学性质表征等方面也有明显的局限性。作者在第三部分中介绍了一种全新的基于原子力显微镜的埃(?)压痕技术,该技术能够将纵向位移尺度控制在0.1纳米以内,从而精确的表征和调控二维材料的层间范德华作用力,即层间/纵向的力学性质。最后作者简要介绍了范德华异质结材料的基本性质,探讨了埃压痕技术在该材料力学性质研究中的潜在应用。
    流体力学
    水中高压脉动气泡与浮体流固耦合特性研究
    胡振宇 曹卓尔 李帅 张阿漫
    2021, 53(4):  944-961.  DOI: 10.6052/0459-1879-20-357
    摘要 ( 288 )   PDF (19341KB) ( 227 )  
    参考文献 | 相关文章 | 计量指标
    本文针对水中电火花气泡与水面浮体流固耦合作用开展实验和数值研究,采用边界积分法对气泡运动进行数值模拟,利用辅助函数法提高非线性流固耦合问题的计算精度,同时运用双节点法保证气-液-固三相交界线的计算稳定性。在实验中,采用水下放电技术生成气泡,使用高速摄影捕捉气泡动力学行为与浮体运动响应。首先对比数值与实验结果,取得了良好的一致性。然后通过对气泡与浮体的无量纲距离 (气泡最大半径为特征长度)系统研究发现:(i)\gamma_s从0.2增大至2时,气泡在坍塌阶段分别形成了颈缩型环状射流(0.2<=\gamma_s<=0.3)、接触射流(0.4<=\gamma_s<=0.6)、非接触射流(0.7<=\gamma_s<=1)、对射流(1.1<=\gamma_s<=1.3)和反射流(1.4<=\gamma_s<=2)等五种典型射流模式;(ii)数值研究表明,当0.7<=\gamma_s<=0.9时,射流砰击速度可达~1000m/s;(iii)在本文实验条件下,\gamma_s<1.5时浮体对气泡的Bjerknes吸引力强于自由液面的Bjerknes排斥力导致气泡在坍塌阶段向浮体迁移;当\gamma_s>1.5时自由液面对气泡的排斥作用更强,使气泡在坍塌阶段远离自由液面。
    不同壁面取向下超疏水平面直轨道上的气泡滑移
    叶煜航 凃程旭 包福兵 汪钰锟 杨森森
    2021, 53(4):  962-972.  DOI: 10.6052/0459-1879-20-405
    摘要 ( 75 )   PDF (4828KB) ( 46 )  
    参考文献 | 相关文章 | 计量指标
    利用特定几何的超疏水表面实现气泡定向输运在矿物浮选和生物孵化等领域具有广阔的应用前景,对平面直线超疏水轨道而言,其壁面取向是相关工程结构的关键参数,但超疏水壁面取向对倾斜壁面气泡滑移的影响尚不清楚。本文采用高速阴影成像系统研究了不同壁面取向(-90° ≤ β ≤ 90°)及轨道倾角(45° ≤ α ≤ 75°)下,气泡(Deq=2.4mm, Re=500~700,We=7~13)在轨道宽度为2mm的超疏水直线轨道上的运动特性。气泡在轨道上的滑移近似为匀速,形状为具有多脊的半子弹型。根据气液界面波动程度的不同滑移气泡可分为波动型和稳定型,稳定型气泡只在较小倾角且较大方位角时出现(α<70°, 45°≤∣β∣)。根据倾角不同滑移速度关于β有2种变化规律:当α≤65°,气泡滑移速度近似为关于β=0°的单峰分布(β=0°时,气泡滑移速度最大);当α≥70°,气泡滑移速度在不同的方位角下基本保持稳定。轨道取向(方位角β)及轨道倾角(α)通过改变气泡沿轨道方向的驱动力和气泡迎风面积影响气泡的滑移速度和气液界面稳定性。
    基于WENO-THINC/WLIC模型的水气二相流数值模拟
    韦志龙 蒋勤
    2021, 53(4):  973-985.  DOI: 10.6052/0459-1879-20-430
    摘要 ( 86 )   PDF (3347KB) ( 37 )  
    相关文章 | 计量指标
    水气二相流与诸多领域的实际工程问题密切相关. 对二相流运动进行高精度的数值模拟是计算流体力学研究的难点和热点.针对开敞水域的自由表面流运动问题, 将水和空气均视为不可压缩流体, 采用五阶加权基本无震荡(weighted essentially non-oscillatory, WENO)格式求解描述流体运动的纳维斯托克斯(Navier-Stokes, NS)方程, 利用以加权线性界面算法改进的多维双曲正切函数界面捕捉法(tangent of hyperbola for interface capturing with weighed line interface calculation, THINC/WLIC)追踪水气界面, 建立WENO-THINC/WLIC水气二相流运动数值模型. 模型采用分步计算法离散求解控制方程, 通过压力投影法求解压强场, 并利用三阶总变差递减(total variation diminishing, TVD)龙格库塔(Runge-Kutta, RK)法对时间项进行离散求解. 通过对环境速度场下Zalesak's disk和shearing vortex界面运动问题, 线性液舱晃荡问题以及溃坝问题的模拟结果与理论分析或试验结果的比较, 对所建立的水气二相流数值模型的适用性及模拟精度进行了验证. 结果表明, 本模型的模拟结果与物理模型或理论分析结果吻合良好, 能较为准确地再现不可压缩水气二相流运动现象. 鉴于WENO格式和THINC法本身在算法及应用等方面仍在不断改进, 本研究提出的WENO-THINC耦合模型为后续更高精度的二相流计算模型开发提供了一种研究思路.
    基于特征分裂有限元准隐格式的共轭传热整体耦合数值模拟方法
    刘瑜 邓家钰 王成恩 苏红星
    2021, 53(4):  986-997.  DOI: 10.6052/0459-1879-20-299
    摘要 ( 55 )   PDF (3578KB) ( 36 )  
    相关文章 | 计量指标
    共轭传热现象在科学和工程领域中大量存在. 随着计算能力的发展,对共轭传热现象进行准确有效的数值模拟,成为科学研究和工程设计上的重要挑战.共轭传热数值模拟的方法可以分为两大类:分区耦合和整体耦合.本文采用有限元法对共轭传热问题进行整体耦合模拟. 固体传热求解采用标准的伽辽金有限元方法.流动求解采用基于特征分裂的有限元方法. 该方法是一种重要的求解流动问题的有限元方法,可以使用等阶有限元. 该方法的准隐格式与其他格式相比,具有时间步长大的特点. 将稳定项中的时间步长与全局时间步长分开,改进了准隐格式的稳定性. 基于改进的特征分裂有限元方法的准隐格式,发展了一种层流共轭传热数值模拟的整体耦合方法. 采用这种方法可以将流体计算域和固体计算域作为一个整体划分有限元网格,并且所有变量都可以采用相同的插值函数,从而有利于程序的实现. 通过对典型问题的模拟,验证了这种方法的准确性. 我们还研究了固体区域时间步长对定常共轭传热问题数值模拟收敛性的影响.
    向列相溶致液晶旋转黏度研究
    陈凌峰 于佳佳 李友荣 黄映洲 李谷元
    2021, 53(4):  998-1007.  DOI: 10.6052/0459-1879-20-272
    摘要 ( 146 )   PDF (1705KB) ( 27 )  
    参考文献 | 相关文章 | 计量指标
    溶致液晶具有良好的生物相容性、无毒性和生物降解性等,在细胞相互作用、神经刺激传递、脂肪吸收、药物智能输运等生命活动研究、医药工程和液晶显示等领域有广泛应用。本文采用旋转磁场法测得溶致液晶Sunset Yellow在不同温度和溶液浓度下向列相时旋转粘度,并结合溶致液晶分子自组装过程,理论分析了向列相溶致液晶旋转粘度随温度和溶液浓度变化规律。研究结果表明:溶致液晶旋转粘度随溶液浓度的增大而增大,但随温度的增加逐渐减小。构建了与温度和浓度相关的旋转粘度经验表达式,经验表达式计算结果与实验值吻合较好,最大误差仅为18.56%。提出利用旋转流变仪间接获得溶致液晶剪切能的新方法,溶致液晶剪切能随溶液浓度的增大而增大,但在实验温度范围内,溶致液晶剪切能几乎不随温度而改变。成功地利用了液晶分子自组装能力随温度的变化规律来研究柱状体长径比对旋转粘度的影响规律,创新地提出了“一步法”测量研究,大大减少相关实验研究的成本和复杂性。
    湿润性对孔隙介质两相渗流驱替效率的影响
    魏鹳举 胡冉 廖震 陈益峰
    2021, 53(4):  1008-1017.  DOI: 10.6052/0459-1879-20-403
    摘要 ( 158 )   PDF (14597KB) ( 136 )  
    相关文章 | 计量指标
    孔隙介质中多相渗流的驱替效率对二氧化碳封存效率和石油采收率具有决定性影响,是实际工程调控中的一个关键指标. 湿润性是影响多相渗流驱替模式及其效率的一个重要因素. 本文通过微流体模型-显微镜-高速相机实验平台,对基于砂岩孔隙结构的微流体模型进行湿润性修饰,开展了5种流量和2种湿润性的两相驱替可视化实验,研究了湿润性对砂岩孔隙结构中两相渗流驱替模式及其效率的重要影响. 实验结果表明:随着流速的增加,两相渗流驱替模式由毛细指流向稳定流发生转变;在低流速条件下,由于毛细力的主导效应,亲水性介质中指进的宽度和被驱替流体团簇的数目均小于疏水性介质. 实验结果还证实了亲水性介质中由于单支优势通道和 “绕流”现象的发生,驱替效率显著小于疏水性介质. 最后,通过对毛细管数进行修正,建立了考虑湿润性影响的驱替效率和毛细管数之间的统一关系式,为不同湿润性条件下驱替效率的预测提供了一种潜在方法.
    波浪作用下柔性草本植物受力特性研究
    张小霞 林鹏智
    2021, 53(4):  1018-1027.  DOI: 10.6052/0459-1879-20-429
    摘要 ( 59 )   PDF (7422KB) ( 77 )  
    相关文章 | 计量指标
    草本盐沼植物广泛分布于海岸带形成滨海湿地系统,这些滨海湿地具有生态、旅游、养殖及减灾等多方面价值。草本植物的消浪作用在海岸稳固和沿海防灾减灾等方面发挥了重要作用。已有研究对草本植物消浪特性的预测主要依赖于调节经验拖曳力系数,而对柔性植物与波浪的相互作用机理研究还很欠缺。草本植物通常由多片柔性叶片和一根柔性较小的杆茎组成。在波浪作用下,植物的叶片与杆茎产生不同程度的动态变形以及相互作用,从而改变波浪与植物间的相对运动速度,使得柔性植物受力特性十分复杂。本文采用简单柔性结构物受力尺度定律计算单片叶片和单根杆茎受力,并用遮蔽系数量化叶片与杆茎间相互作用引起的叶片受力降低程度,提出了波浪作用下柔性草本植物动态受力计算公式。与单片叶片、单根杆茎和既有叶片又有杆茎的模型植物及真实植物样本测量受力比较表明本文提出的植物受力模型可较好地预测植物最大受力及植物受力随波浪周期的变化情况。
    固体力学
    考虑缠结效应的超弹性本构模型
    肖锐 向玉海 钟旦明 曲绍兴
    2021, 53(4):  1028-1037.  DOI: 10.6052/0459-1879-21-008
    摘要 ( 108 )   PDF (339KB) ( 118 )  
    相关文章 | 计量指标
    经典熵弹性模型,如 Neo-Hookean模型和Arruda-Boyce八链模型,被广泛的应用于预测橡胶等软材料的超弹性力学行为,然而,大量实验结果也显示仅采用一套模型参数这类模型不能同时准确地描述橡胶在多种加载模式下的应力响应。为了克服上述模型的不足,本文在熵弹性的模型基础上引入缠结约束效应。微观上,采用Langevin统计模型来表征熵弹性变形自由能,通过管模型(Tube model)引入缠结约束自由能,并基于仿射假设,建立微观变形与宏观变形之间的映射关系。在宏观上,所建立的超弹性模型的Helmholtz自由能同时包含熵弹性和缠结约束两部分,其中熵弹性自由能与经典的Arruda-Boyce八链模型一致,依赖于柯西-格林应变张量的第一不变量,而缠结约束自由能依赖于柯西-格林应变张量的第二不变量。与文献中的实验结果对比发现,该三参数的模型能准确地预测实验中所测得的橡胶材料在单轴拉伸、纯剪切和等双轴拉伸变形条件下的应力响应,也能较好地描述不同预拉伸比条件下双轴拉伸实验结果。最后,本文比较了所建立的基于应变不变量的缠结约束模型与文献中相关的缠结约束模型在多种加载模式下自由能的异同。总的来说,本文所建立的本构理论能准确的模拟橡胶等软材料的大变形力学行为,对其工程应用有促进作用。
    二维正交各向异性位势问题的高阶单元快速多极边界元法
    李聪
    2021, 53(4):  1038-1048.  DOI: 10.6052/0459-1879-20-455
    摘要 ( 58 )   PDF (2159KB) ( 42 )  
    参考文献 | 相关文章 | 计量指标
    本文研制了一种适用于二维正交各向异性位势问题的高阶单元(线性单元和二次单元)快速多极边界元法。在快速多极边界元法中,源点对于远场区域的积分采用多极展开式计算,而对于近场区域的积分则直接进行计算。高阶单元的使用使得近场积分,尤其是奇异和几乎奇异积分的计算更加复杂,本文通过引入复数表达对其进行简化。若边界采用线性单元插值,近场积分可直接解析计算;若采用二次单元插值,则给出一个半解析算法计算近场积分。高阶单元奇异积分和几乎奇异积分的有效解决,使得高阶单元快速多极边界元法不仅能够计算一般的结构,也能被应用于薄体结构,具有更加广泛的应用前景。数值算例表明,高阶单元快速多极边界元法有效提高了计算精度,且适用于求解大规模问题。
    双材料界面裂纹复应力强度因子的正则化边界元法
    谷岩 张耀明
    2021, 53(4):  1049-1058.  DOI: 10.6052/0459-1879-20-440
    摘要 ( 88 )   PDF (1712KB) ( 77 )  
    参考文献 | 相关文章 | 计量指标
    双材料界面裂纹渐近位移和应力场表现出剧烈的振荡特性,许多用于表征经典平方根和负平方根渐近物理场的传统数值方法失效,给界面裂纹复应力强度因子的精确求解增加了难度。本文引入一种含有复振荡因子的新型“特殊裂尖单元”,可精确表征裂纹尖端渐近位移和应力场的振荡特性,在避免裂尖区域高密度网格剖分的情况下,可实现双材料界面裂纹复应力强度因子的精确求解。此外,结合边界元法中计算近奇异积分的正则化算法,成功求解了大尺寸比(超薄)双材料界面裂纹的断裂力学参数。数值算例表面,所提算法稳定,效率高,在不增加计算量的前提下,显著提高了裂尖近场力学参量和断裂力学参数的求解精度和数值稳定性。
    脆性岩石蠕变裂纹成核宏-细观力学机理研究
    李晓照
    2021, 53(4):  1059-1069.  DOI: 10.6052/0459-1879-20-400
    摘要 ( 57 )   PDF (1274KB) ( 65 )  
    参考文献 | 相关文章 | 计量指标
    脆性岩石内部细观裂纹扩展、贯通及成核影响下的脆性蠕变行为,对深部地下工程围岩微震及岩爆事件评价有着重要意义。然而,目前能够解释裂纹成核损伤突变影响下,脆性岩石完整蠕变宏细观力学机理模型的研究很少。本文基于脆性岩石亚临界裂纹扩展模型、裂纹-应变-声发射事件相关的损伤模型及裂纹成核损伤时间演化路径函数,提出了一种脆性岩石裂纹成核损伤突变影响下的蠕变宏细观力学模型。裂纹成核损伤时间演化路径函数通过岩石内部裂纹成核损伤突变大小ΔDCN及相邻裂纹成核损伤时间差Δt进行定义,该函数可以结合岩石声发射监测试验数据定义的岩石损伤数据确定。通过与试验结果对比分析验证模型的合理性。并讨论了裂纹成核损伤大小、相邻裂纹成核损伤时间间隔、及裂纹成核数量对脆性岩蠕变裂纹长度、裂纹速率、轴向应变及应变率的影响。该结果对于更加合理、经济、高效的深部地下工程施工及设计提供了一定的理论支持。
    面向增材制造的应力最小化连通性拓扑优化
    王超 徐斌 段尊义 荣见华
    2021, 53(4):  1070-1080.  DOI: 10.6052/0459-1879-20-389
    摘要 ( 135 )   PDF (7957KB) ( 328 )  
    相关文章 | 计量指标
    增材制造与拓扑优化的有机结合将极大促进高性能产品的研发,但现有基于拓扑优化的设计性能和可制造性研究多是独立开展,或常局限于传统的刚度问题,缺乏对工程中至关重要的强度问题的考虑. 本文面向增材制造,针对协同考虑强度和可制造连通性的结构优化问题,建立了材料体积和连通性标量场约束下的结构应力最小化拓扑优化模型. 针对求解过程中的不同数值困难问题,提出了有效的优化求解策略. 引入基于P范数的全局标量场约束度量,并结合稳定转换误差修正技术来实现对局部标量场的有效控制. 详细推导了相关灵敏度,然后通过典型数值算例论证了文中模型及方法的合理有效性. 结果表明,仅考虑连通性约束的刚度最大化设计不一定能避免局部高应力集中,而该设计也不一定等同于应力最小化连通性设计;充足的材料许用量和恰当的连通性约束边界条件对提高所研究设计的性能至关重要,而应力凝聚参数取值并非越大越好,合理取值才能有助于获取高性能设计. 此外,优化结果也在一定程度上论证了可制造性拓扑优化中考虑强度问题的必要性和可行性.
    饱和多孔粘弹地基热-水-力耦合动力响应分析
    郭颖 李文杰 马建军 梁斌 熊春宝
    2021, 53(4):  1081-1092.  DOI: 10.6052/0459-1879-20-385
    摘要 ( 44 )   PDF (543KB) ( 29 )  
    参考文献 | 相关文章 | 计量指标
    天然土体由于沉积条件和应力状态不同,往往会表现出一定的流变性。本文重点研究了具有流变性特征的土体中渗透系数和孔隙率变化对外荷载作用下热-水-力多场耦合的饱和地基的影响。基于Biot波动方程、达西定律和Lord-Shulman广义热弹性理论,并引入了考虑粘弹性松弛时间因子的Kelvin-Voigt粘弹性模型研究地基上表面受热/力源作用时孔隙率和渗透系数变化对均质各向同性饱和多孔粘弹性地基中所考虑的各无量纲量的影响。根据不同的边界条件采用正则模态法推导出无量纲竖向位移、超孔隙水压力、竖向应力和温度的解析表达式,结合算例分析了不同变量对各物理量的影响并加以图示。正则模态法是一种利用加权残差法,可不经正、反积分变换将方程快速解耦并消除数值反变换的局限性。结果表明:无论何种荷载作用时,荷载频率变化对所有考虑的物理量均有明显的影响;孔隙率和渗透系数均对无量纲超孔隙水压力有明显的影响,当仅考虑热荷载作用时,孔隙率和渗透系数变化对无量纲温度均没有什么影响。正则模态法可广泛应用于岩土工程领域,尤其适用于商业建筑、高速铁路和公路能源基础的热、力学特性研究中。该研究结果可为工程施工奠定一定的理论基础,具有一定的指导性意义。
    动力学与控制
    基于压电纤维复合材料的旋转叶片主动控制
    张博 丁虎 陈立群
    2021, 53(4):  1093-1102.  DOI: 10.6052/0459-1879-20-448
    摘要 ( 68 )   PDF (1213KB) ( 93 )  
    参考文献 | 相关文章 | 计量指标
    长期以来,旋转叶片结构的振动失效占据了航空发动机整机故障的相当比重。发展针对旋转叶片结构的减振技术对于减轻叶片重量,提升叶片性能,延长叶片寿命具有重要意义。通过引入压电纤维复合材料(Macro Fiber Composite, MFC)传感器和作动器,研究预变形旋转叶片2:1内共振的主动控制。建立考虑时滞效应的旋转叶片比例微分闭环控制系统运动方程。通过摄动分析推导出受控叶片的演化方程,并结合延拓法揭示速度增益、位移增益、时滞量等系统参数对受控系统稳态响应及稳定性的影响规律。理论研究结果与数值结果得到相互验证。研究发现时滞量对系统稳定性影响显著,位移增益存在一个范围使得系统出现多值稳态响应,进而破坏了增益平面内系统稳定区和非稳定区域的直线边界。研究结果为叶片结构的减振提供了理论基础。
    基于流形切空间插值的折叠翼参数化气动弹性建模
    詹玖榆 周兴华 黄锐
    2021, 53(4):  1103-1113.  DOI: 10.6052/0459-1879-20-376
    摘要 ( 52 )   PDF (2945KB) ( 43 )  
    参考文献 | 相关文章 | 计量指标
    变体飞行器的气动弹性力学建模是当前先进飞行器设计的研究热点和难点. 然而传统的气动弹性力学分析方法对于变体飞行器问题存在建模效率低、计算复杂等问题. 本研究提出了一种基于流形切空间插值的折叠翼参数化气动弹性建模方法. 首先, 该方法建立数个折叠角的折叠翼有限元模型, 通过流形切空间插值方法建立折叠翼参数化结构动力学模型. 其次, 由偶极子网格法得到参数化非定常气动力模型, 进而得到结构和气动耦合的折叠翼参数化气动弹性模型. 为了验证该参数化模型在气动弹性计算中的准确性,本文以折叠翼为研究对象, 从折叠翼自由振动、颤振边界预测两方面进行了算例验证, 并与直接计算方法进行了对比, 进一步验证了参数化气动弹性模型的有效性. 研究结果表明, 该参数化气动弹性模型对上述两类问题计算精度与直接计算方法一致, 并且有着计算效率更高的优势.
    局域共振型声子晶体板缺陷态带隙及其俘能特性研究
    卢一铭 曹东兴 申永军 陈许敏
    2021, 53(4):  1114-1123.  DOI: 10.6052/0459-1879-20-436
    摘要 ( 93 )   PDF (9151KB) ( 89 )  
    参考文献 | 相关文章 | 计量指标
    设计了一种由圆柱形散射体嵌入环氧树脂基体而组成的周期阵列局域共振型声子晶体板结构,分析了其平直带区域以及缺陷态的能量集中特性,并研究了其振动能量采集特性。首先基于超元胞法结合有限元方法分析了5×5完美声子晶体结构和缺陷态声子晶体结构的能带曲线和能量传输特性;考虑点缺陷局域共振声子晶体结构的能量集中特性,利用压电材料代替超元胞中某点的散射体材料引入点缺陷,分析其振动能量采集特性,结果表明单个5×5点缺陷超胞结构共振频带较窄;为提升俘能效率,提出两种由3个具有不同缺陷态数量和构型的5×5超元胞结构并联而成的5×15声子晶体板结构,机电耦合特性分析结果表明:所提出的局域共振型声子晶体板结构克服了单个点缺陷超胞结构缺陷模过少、共振频带过窄的局限性,拓宽了俘能器的工作频带,提高了输出电压。
    含电学边界的压电层合梁的非线性弯曲波
    赵希宁 杨晓东 张伟
    2021, 53(4):  1124-1137.  DOI: 10.6052/0459-1879-20-409
    摘要 ( 49 )   PDF (4267KB) ( 27 )  
    相关文章 | 计量指标
    本文研究了含电学边界条件的压电层合梁的非线性弯曲波传播特性。首先,考虑几何非线性效应和压电耦合效应,利用哈密顿原理建立了一维无限长矩形压电层合梁弯曲波的非线性方程。其次,采用Jacobi椭圆函数展开法对非线性弯曲波方程进行求解,得到了非线性弯曲波动方程在近似情况下对应的冲击波解和孤波解。最后,利用约化摄动法得到了非线性薛定谔方程,进一步得到了亮孤子和暗孤子解。基于两种方法具体研究了外加电压、压电层厚度等参数对冲击波和孤立波以及亮孤子和暗孤子特性的影响。研究结果表明,在波速较小时,外加电压对冲击波的影响较大,波速较大时,外加电压对孤立波影响减弱。通过调整作用在压电层合梁上的电压发现了存在亮孤子和暗孤子,分析结果表明随着外加电压值的增加,亮孤子和暗孤子的振幅都增加。
    空间双臂机器人抓捕翻滚目标后的鲁棒稳定控制
    夏鹏程 罗建军 王明明
    2021, 53(4):  1138-1155.  DOI: 10.6052/0459-1879-20-449
    摘要 ( 44 )   PDF (1565KB) ( 35 )  
    相关文章 | 计量指标
    针对空间双臂机器人抓捕翻滚目标的稳定控制问题,由于目标惯性参数的不确定性以及双臂同时作用于目标存在内力挤压,已有稳定控制方法无法有效地约束机械臂末端与目标的接触力与力矩,无法保证控制过程中抓捕点处的接触安全。为此,本文考虑被抓捕目标惯性参数不确定性与双臂内力挤压对抓捕后阶段组合体稳定控制的影响,提出了一种保证接触安全的鲁棒稳定控制方法。首先,根据目标的有界不确定性构造鲁棒正不变集,同时考虑双臂与目标接触内力的影响,利用预测控制方法在该不变集内规划控制目标运动的虚拟接触力与力矩。然后,根据目标动力学模型以及空间双臂机器人的运动约束,规划空间机器人的鲁棒安全期望运动轨迹。最后,利用障碍李雅普诺夫函数构造跟踪控制器,使空间双臂机器人有效地跟踪期望轨迹,从而使得目标运动受到等效于虚拟控制律的作用,在保证接触安全的同时实现鲁棒的稳定控制。通过空间双臂机器人稳定翻滚目标的仿真算例验证了所提方法的有效性。
    基于LuGre摩擦模型的接触约束法旋转柔性梁斜碰撞研究
    范纪华 谷通顺 王明强 谌宏 陈立威
    2021, 53(4):  1156-1169.  DOI: 10.6052/0459-1879-20-350
    摘要 ( 43 )   PDF (1003KB) ( 66 )  
    参考文献 | 相关文章 | 计量指标
    基于接触约束法和LuGre摩擦模型对在重力场作用下作大范围旋转运动的柔性梁系统和斜坡发生含摩擦斜碰撞的动力学问题进行研究。首先运用刚柔耦合的多体系统动力学理论对大范围运动的柔性梁进行离散化和动力学建模,在碰撞时采用冲量-动量法求出跳跃速度,其次在法向上引入接触约束求解出碰撞力,在切向上采用LuGre摩擦模型分两种方式求解摩擦力,第一种是在滑动时摩擦力由摩擦系数和碰撞力计算得出,黏滞状态下引入切向约束计算拉格朗日乘子反应实际摩擦力,根据黏滞/滑动切换判断计算出碰撞过程摩擦力(与Coulomb摩擦模型计算摩擦力一致);第二种根据LuGre摩擦模型摩擦系数和法向碰撞力计算其摩擦力,从而在碰撞时无需黏滞/滑动切换,采用相同的摩擦力计算公式。通过与Coulomb摩擦模型对比发现,LuGre摩擦模型描述碰撞切向摩擦过程更精确,LuGre摩擦模型黏滞时建立约束方程和碰撞采用统一的摩擦力公式这两种建模方式描述的斜碰撞动力学特性没有区别,进而说明采用法向接触约束和LuGre摩擦模型具有满足碰撞非嵌入情况、避免黏滞/滑动切换、描述摩擦力相对准确的优势。
    生物、工程及交叉力学
    零重力条件下低温射流抑制大尺寸液氢储罐热分层的数值研究
    郭斌 赵建福 李凯 胡文瑞
    2021, 53(4):  1170-1182.  DOI: 10.6052/0459-1879-20-343
    摘要 ( 68 )   PDF (11136KB) ( 250 )  
    相关文章 | 计量指标
    液氢是一种常沸点低、易蒸发的空间低温推进剂。空间微重力环境中浮力对流被极大减弱,当推进剂储罐壁面存在局部漏热时,储罐内部气液两相流体体系会出现环绕漏热源的热分层现象,引起局部过热沸腾,导致储罐内部压力急剧增大,危害系统结构安全。利用低温射流抑制储罐热分层现象是一种有效手段。低温流体通过设置在储罐内部的射流喷嘴与储罐内部的流体混合,消减局部高温,实现温度的均匀化。本文采用全充满的二维大尺寸储罐模型,对零重力条件下液氢储罐内局部漏热引起的热分层现象开展了数值模拟,分析了不同低温射流条件对于消除零重力条件下液氢储罐内部热分层效果的影响。研究结果表明对于大尺寸储罐,当采用圆形射流喷嘴且低温射流条件相同时,射流喷嘴的位置对罐体内部热分层消除效果影响不是很明显。当射流喷嘴位于储罐内部同一相对位置且入射流量相同时,圆形射流喷嘴因出流方向更集中,罐内流场演变更快,消除热分层比半球形射流喷嘴更有效。
    广义有限差分法在含阻抗边界空腔声学分析中的应用
    陈增涛 王发杰 王超
    2021, 53(4):  1183-1195.  DOI: 10.6052/0459-1879-20-311
    摘要 ( 140 )   PDF (7341KB) ( 54 )  
    参考文献 | 相关文章 | 计量指标
    声学分析在噪声控制、室内隔音等工程计算中有着重要的作用。由于现实生活中的声学模型往往伴随着吸声材料,因此分析含阻抗边界条件的声学问题显得十分必要。广义有限差分法是一种新型区域型无网格数值离散方法,该方法基于多元函数泰勒级数展开式和加权最小二乘拟合,将控制方程中未知参量的各阶偏导数表示为相邻节点函数值的线性组合。本文首次将广义有限差分法应用于含阻抗边界条件空腔声学问题的分析中,建立了空腔声场问题的广义有限差分法数值离散格式。与传统算法相比,所建立的数值模型具有无需网格剖分和数值积分、计算精度高、适用于大规模声学分析等优点。通过具有解析解的经典算例,研究了总节点数目和局部支撑点数目对数值结果的影响,得到了最大计算频率与节点间距之间关系的经验公式。此外,将广义有限差分法应用于无解析解的二维和三维复杂声学模型,并与COMSOL Multiphysics软件所得的有限元结果进行了比较分析。数值实验表明,该算法是一种高效、精确、稳定、收敛的数值模拟方法,在含阻抗边界空腔声学分析中具有广阔的应用前景。
    基于一类非局部宏-微观损伤模型的混凝土典型试件力学行为模拟
    任宇东 陈建兵
    2021, 53(4):  1196-1211.  DOI: 10.6052/0459-1879-20-427
    摘要 ( 58 )   PDF (9926KB) ( 50 )  
    参考文献 | 相关文章 | 计量指标
    混凝土是一类典型的准脆性材料,其受力过程中的非线性分析与裂纹模拟依然是具有挑战性的问题。经典的断裂力学与损伤力学分别从间断与连续的视角对裂纹拓扑进行了描述,是早期人们研究固体破坏问题的有力工具。本世纪以来,相场理论和近场动力学在预测裂纹的萌生、扩展与非线性分析方面取得了重要的进展。最近,结合统一相场理论与近场动力学的基本思想,发展了一类非局部宏-微观损伤模型。该模型引入物质点对的概念来刻画由于变形引起的微细观损伤,对微细观损伤在作用域中进行加权平均得到定量描述物质不连续程度的拓扑损伤。通过具有物理机制的能量退化函数,将拓扑损伤嵌入到连续介质-损伤力学的框架中,这使得该模型在进行非线性分析的同时可以自然地进行裂纹模拟,而毋须预设初始裂纹与裂纹扩展路径。本文考虑细观物理参数的空间变异性,采用非局部宏-微观损伤模型进行单轴受拉混凝土板式试件受力全过程的精细化模拟。通过一维建模标定模型细观参数,并探讨了细观参数与混凝土材料细观物理-几何特性之间的内在关联,在此基础上采用二维模型进行精细化分析。进而,考察了材料参数空间变异性对混凝土单轴受拉试件和带缺口三点弯曲试件力学行为的重要影响。本文的研究工作为非局部宏-微观损伤模型细观参数的试验标定与复杂应力状态下混凝土等准脆性材料的非线性力学行为研究提供了有意义的参考。
    危岩体崩塌灾害监测预警试验研究
    杜岩
    2021, 53(4):  1212-1221.  DOI: 10.6052/0459-1879-20-441
    摘要 ( 58 )   PDF (1504KB) ( 46 )  
    相关文章 | 计量指标
    岩体崩塌破坏的突发性使其成为最难预防的地质灾害之一,严重威胁人类的生命财产安全。边坡岩块体崩塌破坏多是系统不稳定导致的动力破坏,因此应用动力学指标进行监测预警更为有效。本研究通过引入多种时域动力学监测指标,开展了岩体崩塌破坏全过程的监测预警实验研究。通过振动幅值、峭度指标等时域动力学指标监测,可有效识别岩体分离破坏前兆现象,提前55s实现岩块体崩塌的早期预警。多个时域动力学指标均可识别岩体破坏前的非协调性动力特征,其中变异系数在识别这一振荡特征上优势明显。此外,破坏前岩体的振动速度是稳定岩体的2.1倍,岩体破坏发生时刻赋存较大的冲击能量,是岩体启程剧动的主要原因之一。本研究不仅为岩体崩塌灾害早期预警提供了新的数据与技术方法支持,也为崩塌岩体启程剧动机制与破坏后运动特性研究提供了新的启示。