https://lxxb.www.cn100led.com/ Chinese Journal of Theoretical and Applied Mechanics day 15 2024-02-07 14:14:58 lxxb@www.cn100led.com lxxb@www.cn100led.com 2024-02-07 14:14:58 zh lxxb@www.cn100led.com 0459-1879 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-18-321 <br/><p>青光眼是世界上第一大不可逆致盲眼病.其病变与眼内压直接相关,控制眼内压是目前控制青光眼发展的唯一有效途径,但发病的确切机制尚未明确.现已证实,青光眼的原发部位是巩膜筛板：由筛板前后分别承受的眼内压与颅内压产生的压力差会导致筛板结构与形态发生变化,进而挤压穿过筛板的视觉神经,造成视觉神经损伤,产生不可逆的视觉损失.因此,青光眼的发病机理与筛板的力学特性及其周围的力学环境密切相关.自从筛板被确定为青光眼视神经损害的原发部位,筛板便成为该领域的研究热点.其中,通过建立筛板力学模型,研究眼内压与颅内压作用下筛板的受力变形,进而分析筛板变形对视神经的损伤,有助于揭示青光眼视神经损伤机制及青光眼的发病机理.本文将从相关实验、理论和计算以及临床等方面介绍青光眼发病机理中筛板变形的研究进展以及目前存在的问题.</p> <br/>力学学报. 2019 51(5): 1273-1284. <br/><p>青光眼是世界上第一大不可逆致盲眼病.其病变与眼内压直接相关,控制眼内压是目前控制青光眼发展的唯一有效途径,但发病的确切机制尚未明确.现已证实,青光眼的原发部位是巩膜筛板：由筛板前后分别承受的眼内压与颅内压产生的压力差会导致筛板结构与形态发生变化,进而挤压穿过筛板的视觉神经,造成视觉神经损伤,产生不可逆的视觉损失.因此,青光眼的发病机理与筛板的力学特性及其周围的力学环境密切相关.自从筛板被确定为青光眼视神经损害的原发部位,筛板便成为该领域的研究热点.其中,通过建立筛板力学模型,研究眼内压与颅内压作用下筛板的受力变形,进而分析筛板变形对视神经的损伤,有助于揭示青光眼视神经损伤机制及青光眼的发病机理.本文将从相关实验、理论和计算以及临床等方面介绍青光眼发病机理中筛板变形的研究进展以及目前存在的问题.</p> <br/>力学学报. 2019 51(5): 1273-1284. 青光眼发病机理--筛板变形研究进展 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1273-1284. article doi:10.6052/0459-1879-18-321 10.6052/0459-1879-18-321 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-18-321 1273-1284 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-071 <br/><p>黏性液体中的气泡浮升运动有趣而又复杂,而气泡与固壁边界的相互作用更是广泛存在于实际工程中.基于轴对称数值计算,模拟了浮力驱动下气泡在黏性液体中上升并与顶部水平固壁面碰撞、回弹的过程.采用考虑表面张力的不可压、变密度Navier-Stokes方程来描述气液两相流流动,并通过基于分级八叉树的有限体积法进行数值求解.为准确捕捉气泡在回弹过程中局部而迅速的拓扑变化,采用了动态自适应网格技术耦合流体体积法(volume of fluid,VOF)来重构气泡的形状. 从气泡对壁面的碰撞和回弹的基本现象入手,研究了伽利略数 Ga和接触速度$U_{a}$对气泡回弹动力学特性的影响, 分析了气泡碰撞过程中涡结构的变化.用回弹高度$H$、回弹周期$T$、长宽比{$A_{r}$}、浮升速度$U$、轴向位置$z$和回复系数$C_{r}$等参数来表征不同条件时气泡的运动和形状特性. 研究结果表明,气泡的回弹运动特性对 Ga十分敏感. Ga的增大可加剧气泡形变, 促进气泡的回弹运动, 增多回弹次数,增大回弹参数($T$和$H)$, 提升回复系数. 然而,接触速度并非决定气泡回弹动力学的控制参数, $U_{a}$的改变并不会改变回复系数.</p> <br/>力学学报. 2019 51(5): 1285-1295. <br/><p>黏性液体中的气泡浮升运动有趣而又复杂,而气泡与固壁边界的相互作用更是广泛存在于实际工程中.基于轴对称数值计算,模拟了浮力驱动下气泡在黏性液体中上升并与顶部水平固壁面碰撞、回弹的过程.采用考虑表面张力的不可压、变密度Navier-Stokes方程来描述气液两相流流动,并通过基于分级八叉树的有限体积法进行数值求解.为准确捕捉气泡在回弹过程中局部而迅速的拓扑变化,采用了动态自适应网格技术耦合流体体积法(volume of fluid,VOF)来重构气泡的形状. 从气泡对壁面的碰撞和回弹的基本现象入手,研究了伽利略数 Ga和接触速度$U_{a}$对气泡回弹动力学特性的影响, 分析了气泡碰撞过程中涡结构的变化.用回弹高度$H$、回弹周期$T$、长宽比{$A_{r}$}、浮升速度$U$、轴向位置$z$和回复系数$C_{r}$等参数来表征不同条件时气泡的运动和形状特性. 研究结果表明,气泡的回弹运动特性对 Ga十分敏感. Ga的增大可加剧气泡形变, 促进气泡的回弹运动, 增多回弹次数,增大回弹参数($T$和$H)$, 提升回复系数. 然而,接触速度并非决定气泡回弹动力学的控制参数, $U_{a}$的改变并不会改变回复系数.</p> <br/>力学学报. 2019 51(5): 1285-1295. 浮力气泡对水平壁面的回弹动力学特性 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1285-1295. article doi:10.6052/0459-1879-19-071 10.6052/0459-1879-19-071 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-071 1285-1295 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-128 <br/><p>空化流动具有高度的压缩性,空化流动非定常特性及其流体动力与压缩性密切相关.为研究可压缩空化流动空泡脱落的回射流和激波机制下周期性空穴结构演化及其诱导流体动力特性,本文采用多场同步测试方法对典型云状空化流动进行了实验研究,获得了文丘里管扩张段内部云状空化空穴形态演化及其诱导同步壁面压力脉动信号.并基于数字图像处理技术,对附着型片状空穴和脱落型云状空穴结构演化进行了精细化定量分析.结果表明：可压缩空化流动回射流机制下,空穴演化呈现附着型空穴生长$\!$-$\!$-$\!$回射流产生及发展$\!$-$\!$-$\!$附着型空穴失稳断裂及大尺度空泡云团产生脱落的非定常过程,激波机制下空穴演化具有附着型空穴生长$\!$-$\!$-$\!$激波产生及传播$\!$-$\!$-$\!$附着型空穴失稳断裂及大尺度空泡云团脱落的非定常特征,激波传播时间占空穴脱落周期小于回射流推进.激波与空穴相互作用导致空穴内部含气率瞬间大范围大幅度下降,诱导复杂流体动力.激波传播过程中,空泡内部压力脉动大幅增加,激波前缘诱导压力脉冲,而回射流推进过程中,壁面压力脉动相对平稳,回射流头部存在小幅增加. 不同机制下空穴结构存在显著差异,具有不同的相间质量传输过程.</p> <br/>力学学报. 2019 51(5): 1296-1309. <br/><p>空化流动具有高度的压缩性,空化流动非定常特性及其流体动力与压缩性密切相关.为研究可压缩空化流动空泡脱落的回射流和激波机制下周期性空穴结构演化及其诱导流体动力特性,本文采用多场同步测试方法对典型云状空化流动进行了实验研究,获得了文丘里管扩张段内部云状空化空穴形态演化及其诱导同步壁面压力脉动信号.并基于数字图像处理技术,对附着型片状空穴和脱落型云状空穴结构演化进行了精细化定量分析.结果表明：可压缩空化流动回射流机制下,空穴演化呈现附着型空穴生长$\!$-$\!$-$\!$回射流产生及发展$\!$-$\!$-$\!$附着型空穴失稳断裂及大尺度空泡云团产生脱落的非定常过程,激波机制下空穴演化具有附着型空穴生长$\!$-$\!$-$\!$激波产生及传播$\!$-$\!$-$\!$附着型空穴失稳断裂及大尺度空泡云团脱落的非定常特征,激波传播时间占空穴脱落周期小于回射流推进.激波与空穴相互作用导致空穴内部含气率瞬间大范围大幅度下降,诱导复杂流体动力.激波传播过程中,空泡内部压力脉动大幅增加,激波前缘诱导压力脉冲,而回射流推进过程中,壁面压力脉动相对平稳,回射流头部存在小幅增加. 不同机制下空穴结构存在显著差异,具有不同的相间质量传输过程.</p> <br/>力学学报. 2019 51(5): 1296-1309. 可压缩空化流动空穴演化及压力脉动特性实验研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1296-1309. article doi:10.6052/0459-1879-19-128 10.6052/0459-1879-19-128 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-128 1296-1309 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-079 <br/><p>高速列车具有细长形状, 数值评估气动噪声往往需要巨大的计算量.目前对高速列车气动噪声的数值模拟大多基于对简化短编组列车的评估,而实际列车通常具有较长的8$\sim$16节编组.如何基于现有条件合理评价真实长度列车的气动噪声,是一个急需探讨的问题. 本文应用非线性声学求解器(NLAS)和FW--H声学比拟法的混合算法, 先求解噪声积分面上的声场脉动,再进行远场积分, 引入多噪声面积分技术,通过对三种不同长度(3节、4节、6节)列车模型的气动性能和噪声数值模拟,分析了车体长度对列车气动噪声的影响. 结果表明,同一列车模型的各节车厢具有相似的沿线噪声分布,其噪声曲线在量值上十分接近,只是主峰位置会随着车厢空间位置的不同而相应地发生偏移;不同长度编组列车对应部位之间的远场噪声特性具有较强的关联性,它们的远场噪声具有接近的总声压级和噪声频谱.通过利用短编组计算数据进行分解、平移和叠加,成功重构了4编组和6编组列车远场噪声特性,与直接计算结果相比误差在可接受范围内.由此发展了基于短编组列车噪声的数值结果,重构长编组列车沿线噪声的近似评估方法.</p> <br/>力学学报. 2019 51(5): 1310-1320. <br/><p>高速列车具有细长形状, 数值评估气动噪声往往需要巨大的计算量.目前对高速列车气动噪声的数值模拟大多基于对简化短编组列车的评估,而实际列车通常具有较长的8$\sim$16节编组.如何基于现有条件合理评价真实长度列车的气动噪声,是一个急需探讨的问题. 本文应用非线性声学求解器(NLAS)和FW--H声学比拟法的混合算法, 先求解噪声积分面上的声场脉动,再进行远场积分, 引入多噪声面积分技术,通过对三种不同长度(3节、4节、6节)列车模型的气动性能和噪声数值模拟,分析了车体长度对列车气动噪声的影响. 结果表明,同一列车模型的各节车厢具有相似的沿线噪声分布,其噪声曲线在量值上十分接近,只是主峰位置会随着车厢空间位置的不同而相应地发生偏移;不同长度编组列车对应部位之间的远场噪声特性具有较强的关联性,它们的远场噪声具有接近的总声压级和噪声频谱.通过利用短编组计算数据进行分解、平移和叠加,成功重构了4编组和6编组列车远场噪声特性,与直接计算结果相比误差在可接受范围内.由此发展了基于短编组列车噪声的数值结果,重构长编组列车沿线噪声的近似评估方法.</p> <br/>力学学报. 2019 51(5): 1310-1320. 高速列车车体长度对气动噪声影响的数值研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1310-1320. article doi:10.6052/0459-1879-19-079 10.6052/0459-1879-19-079 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-079 1310-1320 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-051 <br/><p>基于半隐式特征线分裂算子有限元法,对低雷诺数下串列布置上游静止方柱--下游双自由度运动方柱体结构的尾激振动问题进行了研究.首先与现有文献结果进行对比验证该方法的正确性.然后着重分析了雷诺数($Re$)与折减速度$(U_{\rm r})$两个关键参数对下游方柱尾激振动响应的影响,同时将计算结果与单方柱工况进行了对比. 数值计算结果表明,雷诺数和折减速度对下游方柱的振幅、振动频率和运动轨迹等动力响应特性的影响较大.随着雷诺数的增大,双柱系统的互扰效应从以涡激效应为主逐渐转变为尾激效应发挥主导作用,从而导致下游方柱的振动响应增强.单方柱工况结构运动轨迹均呈"8"字形. 然而,下游方柱的运动轨迹会随着雷诺数的增加而变得复杂.雷诺数较小时($Re\!=\!40$, 80),下游方柱的运动轨迹基本为"8"字形. 雷诺数较大时($Re\!=\!120$, 160,200), 下游方柱的运动轨迹会出现双"8"字形. 同时,下游方柱的尾流场特性主要呈现2S, 2S*, 2P, 2T, P+S和稳态6种模式.最后, 通过对流场特性进行分析,揭示了串列双方柱系统尾激振动效应的作用机理.</p> <br/>力学学报. 2019 51(5): 1321-1335. <br/><p>基于半隐式特征线分裂算子有限元法,对低雷诺数下串列布置上游静止方柱--下游双自由度运动方柱体结构的尾激振动问题进行了研究.首先与现有文献结果进行对比验证该方法的正确性.然后着重分析了雷诺数($Re$)与折减速度$(U_{\rm r})$两个关键参数对下游方柱尾激振动响应的影响,同时将计算结果与单方柱工况进行了对比. 数值计算结果表明,雷诺数和折减速度对下游方柱的振幅、振动频率和运动轨迹等动力响应特性的影响较大.随着雷诺数的增大,双柱系统的互扰效应从以涡激效应为主逐渐转变为尾激效应发挥主导作用,从而导致下游方柱的振动响应增强.单方柱工况结构运动轨迹均呈"8"字形. 然而,下游方柱的运动轨迹会随着雷诺数的增加而变得复杂.雷诺数较小时($Re\!=\!40$, 80),下游方柱的运动轨迹基本为"8"字形. 雷诺数较大时($Re\!=\!120$, 160,200), 下游方柱的运动轨迹会出现双"8"字形. 同时,下游方柱的尾流场特性主要呈现2S, 2S*, 2P, 2T, P+S和稳态6种模式.最后, 通过对流场特性进行分析,揭示了串列双方柱系统尾激振动效应的作用机理.</p> <br/>力学学报. 2019 51(5): 1321-1335. 上游静止方柱尾流对下游方柱体尾激振动效应影响 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1321-1335. article doi:10.6052/0459-1879-19-051 10.6052/0459-1879-19-051 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-051 1321-1335 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-149 <br/><p>利用自主设计的等离子喷注器采用介质阻挡放电方式产生非平衡等离子体,首先利用纹影技术、热电偶、单点红外测温等多种诊断方法实验研究了纯氧放电等离子体的电学特性、热效应及气动效应,然后通过可见光和化学自发辐射成像技术获得了火焰形态及特征参数,详细分析了等离子体对甲烷--纯氧扩散火焰形态和释热的影响,并计算了放电功率及费效比. 结果表明, 燃烧导致放电电流显著增大,其中电压幅值与氧气流速对放电电流大小的影响规律正好相反;与空气等离子体相比, 相同流量与电压条件下氧等离子体放电功率较高,但其发光强度明显较弱; 氧等离子体热效应微弱, 对燃烧的影响可以忽略,放电反应中释热过程主要由含氧组分决定;放电产生了具有3个速度分量的诱导射流, 增大了氧射流角,且电压越大越显著.等离子体主要通过气动效应改变了燃料与氧化剂的掺混,使得一定条件下火焰变得更稳定、释热更强.在所研究的范围内等离子体作用的费效比最低仅为2.2％,大流量、小混合比更有利.</p> <br/>力学学报. 2019 51(5): 1336-1349. <br/><p>利用自主设计的等离子喷注器采用介质阻挡放电方式产生非平衡等离子体,首先利用纹影技术、热电偶、单点红外测温等多种诊断方法实验研究了纯氧放电等离子体的电学特性、热效应及气动效应,然后通过可见光和化学自发辐射成像技术获得了火焰形态及特征参数,详细分析了等离子体对甲烷--纯氧扩散火焰形态和释热的影响,并计算了放电功率及费效比. 结果表明, 燃烧导致放电电流显著增大,其中电压幅值与氧气流速对放电电流大小的影响规律正好相反;与空气等离子体相比, 相同流量与电压条件下氧等离子体放电功率较高,但其发光强度明显较弱; 氧等离子体热效应微弱, 对燃烧的影响可以忽略,放电反应中释热过程主要由含氧组分决定;放电产生了具有3个速度分量的诱导射流, 增大了氧射流角,且电压越大越显著.等离子体主要通过气动效应改变了燃料与氧化剂的掺混,使得一定条件下火焰变得更稳定、释热更强.在所研究的范围内等离子体作用的费效比最低仅为2.2％,大流量、小混合比更有利.</p> <br/>力学学报. 2019 51(5): 1336-1349. 非平衡等离子体对甲烷--氧扩散火焰影响的实验研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1336-1349. article doi:10.6052/0459-1879-19-149 10.6052/0459-1879-19-149 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-149 1336-1349 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-107 <br/><p>针对复合材料水翼存在的流固耦合求解问题,结合其自身特有属性,对复合材料水翼结构变形特性进行了数值仿真计算研究.研究建立了复合材料水翼流固耦合数值计算模型,并将数值计算结果与Zarruk等的实验结果进行对比,验证模型的正确性,得出复合材料水翼尖端扭转角随雷诺数的增加而增加的研究结论.基于数值计算模型,系统地研究了不同铺层角对复合材料水翼水动力特性及强度特性的影响,结果表明：不同铺层角复合材料水翼的尖端扭转角,随铺层角的增大而减小,而其尖端位移量随铺层角的增大先减小后增大.为了削弱工程常数的影响对复合材料水翼变形的影响,研究提出了无量纲扭转角和无量纲位移量,进一步探究复合材料水翼结构的弯扭耦合作用对其变形特性的影响.最后利用蔡$\!$-$\!$-$\!$吴失效准则进行复合材料水翼强度特性的判断和分析,结果表明：不同铺层角复合材料水翼的蔡$\!$-$\!$-$\!$吴系数,随铺层角的增大呈现先减小后增大的趋势,其中0$^\circ$铺层时的复合材料水翼蔡$\!$-$\!$-$\!$吴系数最小,50$^\circ$铺层时的复合材料水翼的蔡$\!$-$\!$-$\!$吴系数最大.</p> <br/>力学学报. 2019 51(5): 1350-1362. <br/><p>针对复合材料水翼存在的流固耦合求解问题,结合其自身特有属性,对复合材料水翼结构变形特性进行了数值仿真计算研究.研究建立了复合材料水翼流固耦合数值计算模型,并将数值计算结果与Zarruk等的实验结果进行对比,验证模型的正确性,得出复合材料水翼尖端扭转角随雷诺数的增加而增加的研究结论.基于数值计算模型,系统地研究了不同铺层角对复合材料水翼水动力特性及强度特性的影响,结果表明：不同铺层角复合材料水翼的尖端扭转角,随铺层角的增大而减小,而其尖端位移量随铺层角的增大先减小后增大.为了削弱工程常数的影响对复合材料水翼变形的影响,研究提出了无量纲扭转角和无量纲位移量,进一步探究复合材料水翼结构的弯扭耦合作用对其变形特性的影响.最后利用蔡$\!$-$\!$-$\!$吴失效准则进行复合材料水翼强度特性的判断和分析,结果表明：不同铺层角复合材料水翼的蔡$\!$-$\!$-$\!$吴系数,随铺层角的增大呈现先减小后增大的趋势,其中0$^\circ$铺层时的复合材料水翼蔡$\!$-$\!$-$\!$吴系数最小,50$^\circ$铺层时的复合材料水翼的蔡$\!$-$\!$-$\!$吴系数最大.</p> <br/>力学学报. 2019 51(5): 1350-1362. 复合材料水翼水动力与结构强度特性数值研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1350-1362. article doi:10.6052/0459-1879-19-107 10.6052/0459-1879-19-107 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-107 1350-1362 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-118 <br/><p>试样尺度、缺口和加载方式通常对材料的疲劳性能具有重要影响. 因此,发展关联试样尺度、缺口和加载方式对疲劳强度影响的方法对于从材料疲劳性能到结构件疲劳性能的预测具有重要意义.首先,采用旋转弯曲加载和轴向加载方式对不同几何形状EA4T车轴钢试样进行了疲劳实验.实验结果表明, 由于试样尺度的增加,轴向加载下狗骨形试样的疲劳强度明显低于沙漏形试样; 相同寿命下,缺口显著降低试样的疲劳强度. 疲劳断口扫描电镜观测结果表明,疲劳裂纹均起源于试样表面.沙漏形试样和狗骨形试样疲劳断口大多只有一个裂纹源,而缺口试样疲劳断口均具有多裂纹源特征. 然后,采用概率控制体积方法研究了试样尺度、缺口和加载方式对疲劳强度的影响,并与临界距离和应变能密度方法进行了比较. 结果表明,概率控制体积方法能够更好地关联试样尺度、缺口和加载方式对EA4T车轴钢疲劳强度的影响.最后, 提出一种基于控制体积的结构件疲劳强度预测方法,并用于具有不连续高应力区域车轴钢试样的疲劳强度预测,预测结果与实验结果 吻合.</p> <br/>力学学报. 2019 51(5): 1363-1371. <br/><p>试样尺度、缺口和加载方式通常对材料的疲劳性能具有重要影响. 因此,发展关联试样尺度、缺口和加载方式对疲劳强度影响的方法对于从材料疲劳性能到结构件疲劳性能的预测具有重要意义.首先,采用旋转弯曲加载和轴向加载方式对不同几何形状EA4T车轴钢试样进行了疲劳实验.实验结果表明, 由于试样尺度的增加,轴向加载下狗骨形试样的疲劳强度明显低于沙漏形试样; 相同寿命下,缺口显著降低试样的疲劳强度. 疲劳断口扫描电镜观测结果表明,疲劳裂纹均起源于试样表面.沙漏形试样和狗骨形试样疲劳断口大多只有一个裂纹源,而缺口试样疲劳断口均具有多裂纹源特征. 然后,采用概率控制体积方法研究了试样尺度、缺口和加载方式对疲劳强度的影响,并与临界距离和应变能密度方法进行了比较. 结果表明,概率控制体积方法能够更好地关联试样尺度、缺口和加载方式对EA4T车轴钢疲劳强度的影响.最后, 提出一种基于控制体积的结构件疲劳强度预测方法,并用于具有不连续高应力区域车轴钢试样的疲劳强度预测,预测结果与实验结果 吻合.</p> <br/>力学学报. 2019 51(5): 1363-1371. 试样疲劳性能尺度效应的概率控制体积方法 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1363-1371. article doi:10.6052/0459-1879-19-118 10.6052/0459-1879-19-118 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-118 1363-1371 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-120 <br/><p>响应量在临近破坏时呈现出临界幂律奇异性加速特征,是一种被广泛证实的灾变破坏前兆,并被火山、滑坡和岩石破坏实验等后验预测结果证实为一种对破坏时间进行短临期预测的可行方法.但是,奇异性指数测量值的较大分散性导致了对其具体取值的争议和预测效果的不确定性.因此,理解奇异性指数取值特征及其内在物理控制因素,成为了一个核心问题.本文基于连续介质损伤力学和材料时间相关失效特征,构建了刻画损伤加速发展通向破坏过程的力学模型.导出了恒名义应力蠕变加载和控制名义应力随时间线性增大两种典型加载方式下,损伤和应变率加速发展通向破坏的临界幂律奇异性前兆特征.阐明了临界幂律奇异性指数取值依赖于材料损伤与承受真应力之间的非线性关系这一内在物理根源,表明了实际测量中奇异性指数的分散性不完全归结于测量数据误差,而是有着内在物理控制因素.针对破坏前奇异性指数的不确定性,建议了在未知奇异性指数条件下预测破坏时间的方法,并基于花岗岩脆性蠕变破坏实验进行了验证和说明.</p> <br/>力学学报. 2019 51(5): 1372-1380. <br/><p>响应量在临近破坏时呈现出临界幂律奇异性加速特征,是一种被广泛证实的灾变破坏前兆,并被火山、滑坡和岩石破坏实验等后验预测结果证实为一种对破坏时间进行短临期预测的可行方法.但是,奇异性指数测量值的较大分散性导致了对其具体取值的争议和预测效果的不确定性.因此,理解奇异性指数取值特征及其内在物理控制因素,成为了一个核心问题.本文基于连续介质损伤力学和材料时间相关失效特征,构建了刻画损伤加速发展通向破坏过程的力学模型.导出了恒名义应力蠕变加载和控制名义应力随时间线性增大两种典型加载方式下,损伤和应变率加速发展通向破坏的临界幂律奇异性前兆特征.阐明了临界幂律奇异性指数取值依赖于材料损伤与承受真应力之间的非线性关系这一内在物理根源,表明了实际测量中奇异性指数的分散性不完全归结于测量数据误差,而是有着内在物理控制因素.针对破坏前奇异性指数的不确定性,建议了在未知奇异性指数条件下预测破坏时间的方法,并基于花岗岩脆性蠕变破坏实验进行了验证和说明.</p> <br/>力学学报. 2019 51(5): 1372-1380. 连续损伤力学基临界奇异指数与破坏时间预测 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1372-1380. article doi:10.6052/0459-1879-19-120 10.6052/0459-1879-19-120 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-120 1372-1380 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-181 <br/><p>由CVD方法制备的石墨烯含有大量的晶界,通常还带有许多褶皱,本文通过分子动力学方法研究了具有褶皱和晶界的石墨烯 平面拉伸断裂行为,结果显示,在垂直晶界方向,褶皱能够显著提高小角度双晶石墨烯的断裂应力,断裂应力增幅最 大约为50%,褶皱对断裂应力的影响随晶界角的增大减弱,导致双晶石墨烯断裂应力对晶界角不敏感,只略低于单晶石墨 烯,和实验结果完全吻合;在沿晶界方向,褶皱对双晶石墨烯断裂应力影响不明显. 另外,褶皱可以显著提高双晶石墨烯的断 裂应变,增幅最大约为100%. 增强机制归纳主要如下：通过面外变形,褶皱可以部分释放晶界5-7环中C---C键的预拉伸变形, 提高双晶石墨烯的断裂应力;褶皱可以降低相邻5-7环之间相互作用,导致断裂应力对晶界角不敏感;在拉伸作用下,褶皱被部分 拉平,这可以显著降低C---C键面内拉伸变形,导致断裂应变显著增大. 本研究为准确理解多晶石墨烯断裂行为提供重要帮助.</p> <br/>力学学报. 2019 51(5): 1381-1392. <br/><p>由CVD方法制备的石墨烯含有大量的晶界,通常还带有许多褶皱,本文通过分子动力学方法研究了具有褶皱和晶界的石墨烯 平面拉伸断裂行为,结果显示,在垂直晶界方向,褶皱能够显著提高小角度双晶石墨烯的断裂应力,断裂应力增幅最 大约为50%,褶皱对断裂应力的影响随晶界角的增大减弱,导致双晶石墨烯断裂应力对晶界角不敏感,只略低于单晶石墨 烯,和实验结果完全吻合;在沿晶界方向,褶皱对双晶石墨烯断裂应力影响不明显. 另外,褶皱可以显著提高双晶石墨烯的断 裂应变,增幅最大约为100%. 增强机制归纳主要如下：通过面外变形,褶皱可以部分释放晶界5-7环中C---C键的预拉伸变形, 提高双晶石墨烯的断裂应力;褶皱可以降低相邻5-7环之间相互作用,导致断裂应力对晶界角不敏感;在拉伸作用下,褶皱被部分 拉平,这可以显著降低C---C键面内拉伸变形,导致断裂应变显著增大. 本研究为准确理解多晶石墨烯断裂行为提供重要帮助.</p> <br/>力学学报. 2019 51(5): 1381-1392. 褶皱与晶界耦合作用对石墨烯断裂行为的影响 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1381-1392. article doi:10.6052/0459-1879-19-181 10.6052/0459-1879-19-181 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-181 1381-1392 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-175 <br/><p>与石墨烯相比,氧化石墨烯(graphene oxide, GO)的亲水性、分散性和反应活性更好,更易于作为增强材料而研发生成性能超常的复合材料,但另一方面,由于其电子结构较为复杂,致使目前有关力学方面的研究存在一定差异.本文利用分子动力学方法,建立了羟基、羧基和环氧基等官能团随机分布的GO原子模型;通过单向拉伸模拟,分析了其断裂行为,结果表明,远离羟基和羧基的环氧基对断裂具有"诱导"作用,并从化学成键、体系能量和应力分布三个角度对其机理进行了阐释;此外,进一步研究了拉伸应力$\!$-$\!$-$\!$应变曲线、极限强度、极限应变等力学性能与含氧官能团覆盖度间的关系,结果表明,极限强度、极限应变均随含氧官能团覆盖度的增大而呈减小趋势.分析认为,主要原因是官能团的出现对石墨烯面内的sp$^{2}$杂化形式造成了破坏,进而使得原子间键合能弱化,随着含氧官能团的覆盖度的增大,被弱化的键合能的数量和程度将越大,从而使得GO的极限强度、极限应变等越低. 研究结果可为GO的基础研究和工程应用提供参考.</p> <br/>力学学报. 2019 51(5): 1393-1402. <br/><p>与石墨烯相比,氧化石墨烯(graphene oxide, GO)的亲水性、分散性和反应活性更好,更易于作为增强材料而研发生成性能超常的复合材料,但另一方面,由于其电子结构较为复杂,致使目前有关力学方面的研究存在一定差异.本文利用分子动力学方法,建立了羟基、羧基和环氧基等官能团随机分布的GO原子模型;通过单向拉伸模拟,分析了其断裂行为,结果表明,远离羟基和羧基的环氧基对断裂具有"诱导"作用,并从化学成键、体系能量和应力分布三个角度对其机理进行了阐释;此外,进一步研究了拉伸应力$\!$-$\!$-$\!$应变曲线、极限强度、极限应变等力学性能与含氧官能团覆盖度间的关系,结果表明,极限强度、极限应变均随含氧官能团覆盖度的增大而呈减小趋势.分析认为,主要原因是官能团的出现对石墨烯面内的sp$^{2}$杂化形式造成了破坏,进而使得原子间键合能弱化,随着含氧官能团的覆盖度的增大,被弱化的键合能的数量和程度将越大,从而使得GO的极限强度、极限应变等越低. 研究结果可为GO的基础研究和工程应用提供参考.</p> <br/>力学学报. 2019 51(5): 1393-1402. 基于分子动力学的氧化石墨烯拉伸断裂行为与力学性能研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1393-1402. article doi:10.6052/0459-1879-19-175 10.6052/0459-1879-19-175 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-175 1393-1402 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-109 <br/><p>柔性薄膜结构广泛应用于航天飞行器的关键部件.形面平整度是影响薄膜结构性能的主要因素之一.褶皱幅值是评价薄膜反射面天线形面平整度的重要指标.褶皱幅值的大小与垂直于褶皱方向的横向应变密切相关.本文基于薄板稳定性理论,针对对角受拉方形薄膜建立了一个能够准确预测其褶皱变形幅值的理论模型.该模型考虑了横向拉伸力对薄膜变形的影响,将垂直于褶皱方向的位移分解为由泊松效应造成的横向收缩位移、由面外变形造成的褶皱位移以及由横向拉伸力造成的拉伸位移三个部分,重新推导了褶皱幅值的理论公式. 基于数字图像相关技术,对受拉方形薄膜进行了散斑实验测试. 利用双目视觉三维测量系统,测量了方形薄膜的三维位移场, 获得了薄膜的三维变形形貌和褶皱波形图,研究了褶皱幅值与拉伸载荷之间的非线性关系. 与已有理论模型相比,该模型进一步提高了褶皱幅值的计算精度, 与实验结果吻合良好.本文呈现的理论研究可为精确数值模型的建立及算法实现提供有意义的指导.</p> <br/>力学学报. 2019 51(5): 1403-1410. <br/><p>柔性薄膜结构广泛应用于航天飞行器的关键部件.形面平整度是影响薄膜结构性能的主要因素之一.褶皱幅值是评价薄膜反射面天线形面平整度的重要指标.褶皱幅值的大小与垂直于褶皱方向的横向应变密切相关.本文基于薄板稳定性理论,针对对角受拉方形薄膜建立了一个能够准确预测其褶皱变形幅值的理论模型.该模型考虑了横向拉伸力对薄膜变形的影响,将垂直于褶皱方向的位移分解为由泊松效应造成的横向收缩位移、由面外变形造成的褶皱位移以及由横向拉伸力造成的拉伸位移三个部分,重新推导了褶皱幅值的理论公式. 基于数字图像相关技术,对受拉方形薄膜进行了散斑实验测试. 利用双目视觉三维测量系统,测量了方形薄膜的三维位移场, 获得了薄膜的三维变形形貌和褶皱波形图,研究了褶皱幅值与拉伸载荷之间的非线性关系. 与已有理论模型相比,该模型进一步提高了褶皱幅值的计算精度, 与实验结果吻合良好.本文呈现的理论研究可为精确数值模型的建立及算法实现提供有意义的指导.</p> <br/>力学学报. 2019 51(5): 1403-1410. 对角受拉方膜褶皱变形幅值的理论预测及实验验证 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1403-1410. article doi:10.6052/0459-1879-19-109 10.6052/0459-1879-19-109 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-109 1403-1410 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-133 <br/><p>针对平纹编织复合材料低速冲击响应和损伤问题,提出了一种多尺度分析方法. 首先, 建立微观尺度单胞模型,引入周期性边界条件,采用最大主应力失效准则和直接刚度退化模型表征纤维丝和基体的损伤起始与演化,预测了纤维束的弹性性能和强度性能. 其次,将这些性能参数代入介观尺度单胞模型,基于Hashin和Hou的混合失效准则以及连续介质损伤模型对介观尺度单胞进行6种边界条件下的渐进损伤模拟.然后采用渐进均匀化方法,以介观尺度单胞为媒介预测了0$^\circ$和90$^\circ$子胞的性能参数,并建立平纹编织复合材料的子胞模型,进而扩展成为材料的宏观尺度低速冲击模型. 在此基础上,研究了平纹编织复合材料低速冲击下的力学响应与损伤特征.结果表明：宏观冲击仿真和试验吻合较好, 验证了多尺度方法的正确性;最大接触力、材料吸能和分层面积均随冲击能量的增大而增大,分层损伤轮廓逐渐从椭圆形向圆形转化;基体拉伸和压缩损伤的长轴方向分别与子胞材料主方向正交和一致,损伤面积前者远大于后者.</p> <br/>力学学报. 2019 51(5): 1411-1423. <br/><p>针对平纹编织复合材料低速冲击响应和损伤问题,提出了一种多尺度分析方法. 首先, 建立微观尺度单胞模型,引入周期性边界条件,采用最大主应力失效准则和直接刚度退化模型表征纤维丝和基体的损伤起始与演化,预测了纤维束的弹性性能和强度性能. 其次,将这些性能参数代入介观尺度单胞模型,基于Hashin和Hou的混合失效准则以及连续介质损伤模型对介观尺度单胞进行6种边界条件下的渐进损伤模拟.然后采用渐进均匀化方法,以介观尺度单胞为媒介预测了0$^\circ$和90$^\circ$子胞的性能参数,并建立平纹编织复合材料的子胞模型,进而扩展成为材料的宏观尺度低速冲击模型. 在此基础上,研究了平纹编织复合材料低速冲击下的力学响应与损伤特征.结果表明：宏观冲击仿真和试验吻合较好, 验证了多尺度方法的正确性;最大接触力、材料吸能和分层面积均随冲击能量的增大而增大,分层损伤轮廓逐渐从椭圆形向圆形转化;基体拉伸和压缩损伤的长轴方向分别与子胞材料主方向正交和一致,损伤面积前者远大于后者.</p> <br/>力学学报. 2019 51(5): 1411-1423. 基于渐进均匀化的平纹编织复合材料低速冲击多尺度方法 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1411-1423. article doi:10.6052/0459-1879-19-133 10.6052/0459-1879-19-133 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-133 1411-1423 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-099 <br/><p>提出应用混合摄动$\!$-$\!$-$\!$伽辽金法求解随机桁架结构的几何非线性问题.将含位移项的随机割线弹性模量以及随机响应表示为幂多项式展开,利用高阶摄动方法确定随机结构几何非线性响应的幂多项式展开的各项系数.将随机响应的各阶摄动项假定为伽辽金试函数,运用伽辽金投影对试函数系数进行求解,从而得到随机桁架结构几何非线性响应的显式表达式.同已有的随机伽辽金法相比,本文所给的试函数由摄动解的线性组合而成,在求解非线性问题时,试函数的获取具有自适应性.数值算例结果表明,对于具有不同概率分布的多随机变量问题,本文方法无需对随机变量的概率分布形式进行转换,避免了转换误差,因而比同阶的广义正交多项式方法(generalizedpolynomial chaos, GPC)计算精度高.同时,在结果精度相当时,和GPC方法相比,本文方法得到的试函数系数的非线性方程维度不大,方程的求解工作量小且更易求解.当随机量涨落较大时,混合摄动$\!$-$\!$-$\!$伽辽金法计算所得的结构响应的各阶统计矩比高阶摄动法所得结果更逼近于蒙特卡洛模拟结果,显示了该方法对几何非线性随机问题求解的有效性.</p> <br/>力学学报. 2019 51(5): 1424-1436. <br/><p>提出应用混合摄动$\!$-$\!$-$\!$伽辽金法求解随机桁架结构的几何非线性问题.将含位移项的随机割线弹性模量以及随机响应表示为幂多项式展开,利用高阶摄动方法确定随机结构几何非线性响应的幂多项式展开的各项系数.将随机响应的各阶摄动项假定为伽辽金试函数,运用伽辽金投影对试函数系数进行求解,从而得到随机桁架结构几何非线性响应的显式表达式.同已有的随机伽辽金法相比,本文所给的试函数由摄动解的线性组合而成,在求解非线性问题时,试函数的获取具有自适应性.数值算例结果表明,对于具有不同概率分布的多随机变量问题,本文方法无需对随机变量的概率分布形式进行转换,避免了转换误差,因而比同阶的广义正交多项式方法(generalizedpolynomial chaos, GPC)计算精度高.同时,在结果精度相当时,和GPC方法相比,本文方法得到的试函数系数的非线性方程维度不大,方程的求解工作量小且更易求解.当随机量涨落较大时,混合摄动$\!$-$\!$-$\!$伽辽金法计算所得的结构响应的各阶统计矩比高阶摄动法所得结果更逼近于蒙特卡洛模拟结果,显示了该方法对几何非线性随机问题求解的有效性.</p> <br/>力学学报. 2019 51(5): 1424-1436. 随机桁架结构几何非线性问题的混合摄动-伽辽金法求解 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1424-1436. article doi:10.6052/0459-1879-19-099 10.6052/0459-1879-19-099 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-099 1424-1436 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-104 <br/><p>随机过程或随机系统响应的最大绝对值概率分布往往是科学与工程中关心的重要挑战性问题.本文从理论与数值上进行了Markov过程的时变最大绝对值过程及其概率分布研究.文中,通过引入扩展状态向量,构造了最大绝对值$\!$-$\!$-$\!$状态量联合向量过程,由此将不具有Markov性的最大值过程转化为具有Markov性的向量随机过程.在此基础上,通过最大绝对值$\!$-$\!$-$\!$状态量之间的关系,建立了联合向量过程的转移概率密度函数.进而,结合Chapman-Kolmogorov方程和路径积分方法,提出了最大绝对值概率密度函数求解的数值方法.由此,可以得到Markov过程最大绝对值过程的时变概率密度函数,可进一步用于结构动力可靠度分析等.通过数值算例,验证了本文所提方法的有效性. 该方法有望推广到更一般随机系统的极值分布估计之中.</p> <br/>力学学报. 2019 51(5): 1437-1447. <br/><p>随机过程或随机系统响应的最大绝对值概率分布往往是科学与工程中关心的重要挑战性问题.本文从理论与数值上进行了Markov过程的时变最大绝对值过程及其概率分布研究.文中,通过引入扩展状态向量,构造了最大绝对值$\!$-$\!$-$\!$状态量联合向量过程,由此将不具有Markov性的最大值过程转化为具有Markov性的向量随机过程.在此基础上,通过最大绝对值$\!$-$\!$-$\!$状态量之间的关系,建立了联合向量过程的转移概率密度函数.进而,结合Chapman-Kolmogorov方程和路径积分方法,提出了最大绝对值概率密度函数求解的数值方法.由此,可以得到Markov过程最大绝对值过程的时变概率密度函数,可进一步用于结构动力可靠度分析等.通过数值算例,验证了本文所提方法的有效性. 该方法有望推广到更一般随机系统的极值分布估计之中.</p> <br/>力学学报. 2019 51(5): 1437-1447. 一类Markov过程的最大绝对值过程概率密度求解的新方法 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1437-1447. article doi:10.6052/0459-1879-19-104 10.6052/0459-1879-19-104 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-104 1437-1447 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-111 <br/><p>以充液航天器为工程背景,借助多尺度方法研究刚--液耦合动力学系统非线性动力学特性.利用多维模态方法,将描述横向外激励下圆柱贮箱中液体非线性晃动的自由边界问题转换为液体模态系数相互耦合的有限维非线性常微分方程组.推导液体晃动产生的作用于贮箱壁的晃动力和晃动力矩的解析表达式,进而建立航天器刚体部分平动和液体晃动耦合的非线性动力学方程组.应用多尺度方法对刚--液耦合系统的动力学特性进行解析分析,通过固有频率的特征方程求解耦合系统固有频率,推导外激励频率接近耦合系统第一阶固有频率时液体晃动稳态解的幅值频率响应方程.结合数值方法,研究了液体晃动稳态解的幅值频率响应曲线和激励--幅值响应曲线.结果表明, 随充液比变化,液体晃动稳态解的幅值频率响应曲线会发生软、硬弹簧特性转换现象和"跳跃"现象;幅值频率响应曲线的软、硬弹簧特性转换点受重力加速度和弹簧刚度系数影响;以上所得研究结果表明,考虑非线性效应时的刚--液耦合系统动力学特性与传统的线性系统模型所显示的动力学特性具有本质区别.本文的研究工作对进一步分析充液航天器刚--液耦合非线性动力学特性具有重要参考价值.</p> <br/>力学学报. 2019 51(5): 1448-1454. <br/><p>以充液航天器为工程背景,借助多尺度方法研究刚--液耦合动力学系统非线性动力学特性.利用多维模态方法,将描述横向外激励下圆柱贮箱中液体非线性晃动的自由边界问题转换为液体模态系数相互耦合的有限维非线性常微分方程组.推导液体晃动产生的作用于贮箱壁的晃动力和晃动力矩的解析表达式,进而建立航天器刚体部分平动和液体晃动耦合的非线性动力学方程组.应用多尺度方法对刚--液耦合系统的动力学特性进行解析分析,通过固有频率的特征方程求解耦合系统固有频率,推导外激励频率接近耦合系统第一阶固有频率时液体晃动稳态解的幅值频率响应方程.结合数值方法,研究了液体晃动稳态解的幅值频率响应曲线和激励--幅值响应曲线.结果表明, 随充液比变化,液体晃动稳态解的幅值频率响应曲线会发生软、硬弹簧特性转换现象和"跳跃"现象;幅值频率响应曲线的软、硬弹簧特性转换点受重力加速度和弹簧刚度系数影响;以上所得研究结果表明,考虑非线性效应时的刚--液耦合系统动力学特性与传统的线性系统模型所显示的动力学特性具有本质区别.本文的研究工作对进一步分析充液航天器刚--液耦合非线性动力学特性具有重要参考价值.</p> <br/>力学学报. 2019 51(5): 1448-1454. 基于多尺度方法的平动圆柱贮箱航天器刚--液耦合动力学研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1448-1454. article doi:10.6052/0459-1879-19-111 10.6052/0459-1879-19-111 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-111 1448-1454 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-076 <br/><p>相比于浮动坐标系法, 绝对节点坐标法(absolute nodal coordinateformulation, ANCF)在处理柔性体非线性大变形问题上具有显著优势,ANCF将单元节点坐标定义在全局坐标系下,采用斜率矢量代替节点转角坐标, 具有常数质量阵,不存在科氏离心力等优点, 然而弹性力阵为非线性项,其求解将比较耗时且占用资源. 据此, 在弹性力求解方法中,引入弹性线方法(elastic line method, ELM),该方法将格林--拉格朗日应变张量定义在中心线上,采用曲率公式来定义弯曲应变, 转角公式来定义扭转应变.同时采用有限元法对三维柔性梁位移场进行离散,求解梁单元常数质量阵、广义刚度阵、广义力阵,进而得到单元的动力学方程, 通过转换矩阵得到三维梁的动力学方程.接着从理论上指出连续介质力学方法(continuum mechanics method,CMM)和弹性线方法在求解弹性力上的不同点, 并编制动力学仿真软件.最后分别采用连续介质力学方法和弹性线方法对柔性单摆以及履带式车辆的动力学问题进行仿真分析,结果表明：弹性线方法能在保证精度的前提下有效提高计算效率.</p> <br/>力学学报. 2019 51(5): 1455-1465. <br/><p>相比于浮动坐标系法, 绝对节点坐标法(absolute nodal coordinateformulation, ANCF)在处理柔性体非线性大变形问题上具有显著优势,ANCF将单元节点坐标定义在全局坐标系下,采用斜率矢量代替节点转角坐标, 具有常数质量阵,不存在科氏离心力等优点, 然而弹性力阵为非线性项,其求解将比较耗时且占用资源. 据此, 在弹性力求解方法中,引入弹性线方法(elastic line method, ELM),该方法将格林--拉格朗日应变张量定义在中心线上,采用曲率公式来定义弯曲应变, 转角公式来定义扭转应变.同时采用有限元法对三维柔性梁位移场进行离散,求解梁单元常数质量阵、广义刚度阵、广义力阵,进而得到单元的动力学方程, 通过转换矩阵得到三维梁的动力学方程.接着从理论上指出连续介质力学方法(continuum mechanics method,CMM)和弹性线方法在求解弹性力上的不同点, 并编制动力学仿真软件.最后分别采用连续介质力学方法和弹性线方法对柔性单摆以及履带式车辆的动力学问题进行仿真分析,结果表明：弹性线方法能在保证精度的前提下有效提高计算效率.</p> <br/>力学学报. 2019 51(5): 1455-1465. 基于绝对节点坐标法的弹性线方法研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1455-1465. article doi:10.6052/0459-1879-19-076 10.6052/0459-1879-19-076 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-076 1455-1465 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-154 <br/><p>利用固定点理论优化接地类型的动力吸振器得到的结果可能不是全局最优参数,在选择其他参数时主系统可以获得更小的振幅, 接地类型动力吸振器的优化问题值得进一步研究. 因此,以一种接地式三要素型动力吸振器为对象,通过研究系统参数变化对固 定点位置与主系统最大振幅的影响,得到了此吸振器的局部最优参数并分析了它的性能. 首先建立了此系统模型的运动微分方程, 得到了主系统振幅放大因子,发现系统存在3个与阻尼无关的固定点. 固定点中幅值较大点随系统参数变化的趋势可以代表最大振 幅随系统参数变化的趋势,因此利用盛金公式得到了固定点幅值的表达式. 为了更加精确,进一步使用数值算法得到了最大振幅与 系统参数的关系图,发现系统中存在局部最优参数. 通过对比接地式吸振器与接地三要素型吸振器的最大振幅随系统参数变化的趋 势,得到了接地式三要素型吸振器的局部最优参数,并发现当固有频率比小于局部最优频率比时,接地式三要素型吸振器模型主系 统的最大振幅要远小于接地式动力吸振器模型.</p> <br/>力学学报. 2019 51(5): 1466-1475. <br/><p>利用固定点理论优化接地类型的动力吸振器得到的结果可能不是全局最优参数,在选择其他参数时主系统可以获得更小的振幅, 接地类型动力吸振器的优化问题值得进一步研究. 因此,以一种接地式三要素型动力吸振器为对象,通过研究系统参数变化对固 定点位置与主系统最大振幅的影响,得到了此吸振器的局部最优参数并分析了它的性能. 首先建立了此系统模型的运动微分方程, 得到了主系统振幅放大因子,发现系统存在3个与阻尼无关的固定点. 固定点中幅值较大点随系统参数变化的趋势可以代表最大振 幅随系统参数变化的趋势,因此利用盛金公式得到了固定点幅值的表达式. 为了更加精确,进一步使用数值算法得到了最大振幅与 系统参数的关系图,发现系统中存在局部最优参数. 通过对比接地式吸振器与接地三要素型吸振器的最大振幅随系统参数变化的趋 势,得到了接地式三要素型吸振器的局部最优参数,并发现当固有频率比小于局部最优频率比时,接地式三要素型吸振器模型主系 统的最大振幅要远小于接地式动力吸振器模型.</p> <br/>力学学报. 2019 51(5): 1466-1475. 接地式三要素型动力吸振器性能分析 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1466-1475. article doi:10.6052/0459-1879-19-154 10.6052/0459-1879-19-154 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-154 1466-1475 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-119 <br/><p>针对折叠舵面内、外舵铰接处存在的间隙对地面振动响应的影响及间隙处的非线性建模方法展开研究.消除间隙,利用锤击法对线性折叠舵面进行模态实验,得到了前五阶模态参数;打开间隙,进行振动台扫频基础激励,实验结果表明间隙的存在会使结构的动力学响应产生非线性现象,如正反向扫描差异、跳跃、多谐波及频率漂移.非线性的影响主要体现在一阶弯曲模态上,激励量级的增大和间隙的减小均会使基频增大,且逐渐趋向于无间隙的结果,但对第二阶扭转模态的影响与第一阶相比较小.建立了折叠舵面的有限元模型. 提出了一种适用于具有集中非线性的折叠机构的模型缩减方法,并对舵面进行了模态缩减.根据Hertz接触理论,用具有线性和3/2次刚度组合形式的非线性扭转弹簧来模拟铰接处的间隙和接触.通过比较锤击实验与数值计算得到的前四阶频率和振型对模型的线性部分进行验证.通过Bathe两子步隐式复合算法计算基础激励下非线性结构的动力学响应,得到的传递函数可以模拟实验中出现的频率变化特征,验证了连接处非线性建模方法的合理性.</p> <br/>力学学报. 2019 51(5): 1476-1488. <br/><p>针对折叠舵面内、外舵铰接处存在的间隙对地面振动响应的影响及间隙处的非线性建模方法展开研究.消除间隙,利用锤击法对线性折叠舵面进行模态实验,得到了前五阶模态参数;打开间隙,进行振动台扫频基础激励,实验结果表明间隙的存在会使结构的动力学响应产生非线性现象,如正反向扫描差异、跳跃、多谐波及频率漂移.非线性的影响主要体现在一阶弯曲模态上,激励量级的增大和间隙的减小均会使基频增大,且逐渐趋向于无间隙的结果,但对第二阶扭转模态的影响与第一阶相比较小.建立了折叠舵面的有限元模型. 提出了一种适用于具有集中非线性的折叠机构的模型缩减方法,并对舵面进行了模态缩减.根据Hertz接触理论,用具有线性和3/2次刚度组合形式的非线性扭转弹簧来模拟铰接处的间隙和接触.通过比较锤击实验与数值计算得到的前四阶频率和振型对模型的线性部分进行验证.通过Bathe两子步隐式复合算法计算基础激励下非线性结构的动力学响应,得到的传递函数可以模拟实验中出现的频率变化特征,验证了连接处非线性建模方法的合理性.</p> <br/>力学学报. 2019 51(5): 1476-1488. 有间隙折叠舵面的振动实验与非线性建模研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1476-1488. article doi:10.6052/0459-1879-19-119 10.6052/0459-1879-19-119 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-119 1476-1488 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-093 <br/><p>通过将系统参数定义为参数变量, 构成参数空间,研究齿轮传动系统在参数空间和状态空间耦合下的非线性全局动力学特性,以及多参数、多初值和多稳态行为之间的关联特性.首先设计了一个两空间耦合下非线性系统多稳态行为的计算和辨识方法.其次,基于该方法并结合相图、Poincaré映射图、分岔图、最大Lyapunov指数、吸引域等,研究齿轮传动系统在不同参数平面上多稳态行为的存在区域和分布特性,以及多稳态行为在状态平面上的分布特性,揭示了参数平面和状态平面上系统可能隐藏的多稳态行为和分岔,并分析了多稳态行为的形成机理. 结果发现,两空间耦合下系统在参数平面上存在大量多稳态行为并呈"带状"分布, 状态平面上多稳态行为出现两种不同的侵蚀现象, 即内部侵蚀和边界侵蚀.分岔点或分岔曲线对初值的敏感性导致多稳态行为的出现.当齿侧间隙和误差波动在较小的范围内变化时,系统全局动力学特性受间隙和误差扰动的影响较小,受啮合频率的影响较大.两空间耦合下系统全局动力学特性变得丰富和复杂.</p> <br/>力学学报. 2019 51(5): 1489-1499. <br/><p>通过将系统参数定义为参数变量, 构成参数空间,研究齿轮传动系统在参数空间和状态空间耦合下的非线性全局动力学特性,以及多参数、多初值和多稳态行为之间的关联特性.首先设计了一个两空间耦合下非线性系统多稳态行为的计算和辨识方法.其次,基于该方法并结合相图、Poincaré映射图、分岔图、最大Lyapunov指数、吸引域等,研究齿轮传动系统在不同参数平面上多稳态行为的存在区域和分布特性,以及多稳态行为在状态平面上的分布特性,揭示了参数平面和状态平面上系统可能隐藏的多稳态行为和分岔,并分析了多稳态行为的形成机理. 结果发现,两空间耦合下系统在参数平面上存在大量多稳态行为并呈"带状"分布, 状态平面上多稳态行为出现两种不同的侵蚀现象, 即内部侵蚀和边界侵蚀.分岔点或分岔曲线对初值的敏感性导致多稳态行为的出现.当齿侧间隙和误差波动在较小的范围内变化时,系统全局动力学特性受间隙和误差扰动的影响较小,受啮合频率的影响较大.两空间耦合下系统全局动力学特性变得丰富和复杂.</p> <br/>力学学报. 2019 51(5): 1489-1499. 两空间耦合下齿轮传动系统多稳态特性研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1489-1499. article doi:10.6052/0459-1879-19-093 10.6052/0459-1879-19-093 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-093 1489-1499 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-049 <br/><p>超声共振谱技术通过测量样本在超声激励下产生的固有共振频率来计算弹性参数,而共振频率的提取是整个测量过程的关键.低$Q$值(品质因数)材料由于其衰减特性,导致共振谱平缓并无法直观地从谱图上观察得到共振频率,为从中提取更为有效的共振频率, 本文提出了一种新的共振频率提取方法.采用经验模态分解法将材料频率响应自适应分解为有限个具有特殊振荡特性的固有模态函数分量,根据材料的超声共振谱先验信息选择具有共振频率特性的固有模态函数分量,并从中提取共振频率. 以短切纤维环氧树脂材料(仿骨材料, $Q \approx$25)为例, 通过实验与传统线性预测方法进行对比,计算弹性系数和工程模量. 实验结果表明新方法的计算效率高,对弱激发模态更为敏感,共振频率的匹配数量(26)多于传统方法(21)且满足5倍于弹性系数的估计要求,优化后的弹性模量更接近标准值.新方法可从低$Q$值材料平缓的频谱中提取数量足够且有效的共振频率,不仅有效提升了力学参数估计的可靠性,而且拓展了超声共振谱技术的应用范围.</p> <br/>力学学报. 2019 51(5): 1500-1506. <br/><p>超声共振谱技术通过测量样本在超声激励下产生的固有共振频率来计算弹性参数,而共振频率的提取是整个测量过程的关键.低$Q$值(品质因数)材料由于其衰减特性,导致共振谱平缓并无法直观地从谱图上观察得到共振频率,为从中提取更为有效的共振频率, 本文提出了一种新的共振频率提取方法.采用经验模态分解法将材料频率响应自适应分解为有限个具有特殊振荡特性的固有模态函数分量,根据材料的超声共振谱先验信息选择具有共振频率特性的固有模态函数分量,并从中提取共振频率. 以短切纤维环氧树脂材料(仿骨材料, $Q \approx$25)为例, 通过实验与传统线性预测方法进行对比,计算弹性系数和工程模量. 实验结果表明新方法的计算效率高,对弱激发模态更为敏感,共振频率的匹配数量(26)多于传统方法(21)且满足5倍于弹性系数的估计要求,优化后的弹性模量更接近标准值.新方法可从低$Q$值材料平缓的频谱中提取数量足够且有效的共振频率,不仅有效提升了力学参数估计的可靠性,而且拓展了超声共振谱技术的应用范围.</p> <br/>力学学报. 2019 51(5): 1500-1506. 一种基于超声共振谱的低&lt;em&lt;Q&lt;/em&lt;值材料共振频率提取方法 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1500-1506. article doi:10.6052/0459-1879-19-049 10.6052/0459-1879-19-049 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-049 1500-1506 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-105 <br/><p>传统采用微分求积(differential quadrature,DQ)法求解动力问题时都是以位移响应作为基本未知量,而将速度响应和加速度响应表示为位移响应的加权和的形式.如此做法需要处理线性方程组或者矩阵方程(Sylvester方程)才能求得动力响应,导出的算法一般为有条件稳定算法.本文利用动力响应的Duhamel积分解,逆用DQ原理,提出了一种计算卷积的高精度显式算法.该算法可以逐时段地求解出动力时程响应,当各时段内DQ节点分布完全一致时,仅须进行一次Vandermonde矩阵求逆计算即可应用于各个时段,一次性获得时段内多个时刻的位移响应值,因而具有计算效率高的优点.通过分析动力方程积分格式,证明本文动力算法传递矩阵的谱半径恒等于1,因而该算法具有无条件稳定特性,且计算过程中不会产生数值耗散. 本文算法的数值精度取决于分析时段内布置的DQ节点数量$N$,具有$N-1$阶代数精度.实际操作时可以取10个甚至更多的DQ节点数,从而获得比较高的数值精度.</p> <br/>力学学报. 2019 51(5): 1507-1516. <br/><p>传统采用微分求积(differential quadrature,DQ)法求解动力问题时都是以位移响应作为基本未知量,而将速度响应和加速度响应表示为位移响应的加权和的形式.如此做法需要处理线性方程组或者矩阵方程(Sylvester方程)才能求得动力响应,导出的算法一般为有条件稳定算法.本文利用动力响应的Duhamel积分解,逆用DQ原理,提出了一种计算卷积的高精度显式算法.该算法可以逐时段地求解出动力时程响应,当各时段内DQ节点分布完全一致时,仅须进行一次Vandermonde矩阵求逆计算即可应用于各个时段,一次性获得时段内多个时刻的位移响应值,因而具有计算效率高的优点.通过分析动力方程积分格式,证明本文动力算法传递矩阵的谱半径恒等于1,因而该算法具有无条件稳定特性,且计算过程中不会产生数值耗散. 本文算法的数值精度取决于分析时段内布置的DQ节点数量$N$,具有$N-1$阶代数精度.实际操作时可以取10个甚至更多的DQ节点数,从而获得比较高的数值精度.</p> <br/>力学学报. 2019 51(5): 1507-1516. 一种结构动力时程分析的积分求微方法 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1507-1516. article doi:10.6052/0459-1879-19-105 10.6052/0459-1879-19-105 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-105 1507-1516 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-156 <br/><p>海底地震动场及海洋声场的模拟中,需要考虑复杂海床介质及海底地形的影响,涉及到海水、饱和海床、弹性基岩之间的相互耦合.传统的方法分别采用声波方程描述理想流体、Biot方程描述饱和海床、弹性波方程描述基岩,分别进行空间离散和界面耦合, 十分不便.本文基于理想流体、固体分别为饱和多孔介质的特殊情形(孔隙率分别为1和0),由饱和多孔介质的Biot方程可退化得到理想流体的声波方程和固体的弹性波方程.然后, 以饱和多孔介质方程为基础, 经集中质量有限元离散,严格考虑不同孔隙率的饱和多孔介质在不规则界面的耦合条件,通过求解法向和切向界面力的途径,建立了不同孔隙率的饱和多孔介质耦合情形的求解方法,将流体、固体、饱和多孔介质间的耦合问题纳入到统一计算框架,并编制了相应的三维并行分析程序.考虑海水--弹性基岩、海水--饱和海床--弹性基岩体系中凹陷地形情形,采用本文提出的统一计算框架, 结合透射边界条件,分析了P波入射时的动力反应, 并通过结果是否满足界面条件,验证了该统一计算框架的有效性以及并行计算的可行性.</p> <br/>力学学报. 2019 51(5): 1517-1529. <br/><p>海底地震动场及海洋声场的模拟中,需要考虑复杂海床介质及海底地形的影响,涉及到海水、饱和海床、弹性基岩之间的相互耦合.传统的方法分别采用声波方程描述理想流体、Biot方程描述饱和海床、弹性波方程描述基岩,分别进行空间离散和界面耦合, 十分不便.本文基于理想流体、固体分别为饱和多孔介质的特殊情形(孔隙率分别为1和0),由饱和多孔介质的Biot方程可退化得到理想流体的声波方程和固体的弹性波方程.然后, 以饱和多孔介质方程为基础, 经集中质量有限元离散,严格考虑不同孔隙率的饱和多孔介质在不规则界面的耦合条件,通过求解法向和切向界面力的途径,建立了不同孔隙率的饱和多孔介质耦合情形的求解方法,将流体、固体、饱和多孔介质间的耦合问题纳入到统一计算框架,并编制了相应的三维并行分析程序.考虑海水--弹性基岩、海水--饱和海床--弹性基岩体系中凹陷地形情形,采用本文提出的统一计算框架, 结合透射边界条件,分析了P波入射时的动力反应, 并通过结果是否满足界面条件,验证了该统一计算框架的有效性以及并行计算的可行性.</p> <br/>力学学报. 2019 51(5): 1517-1529. 海洋地震工程流固耦合问题统一计算框架------不规则界面情形 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1517-1529. article doi:10.6052/0459-1879-19-156 10.6052/0459-1879-19-156 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-156 1517-1529 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-18-440 <br/><p>涡核模型中的涡核尺寸对自由涡尾迹(free vortex wake,FVW)方法准确预估风力机气动特性至关重要,涡核尺寸包括初始涡核半径和由于耗散效应涡核半径在尾迹中的增长.FVW方法中涡线控制方程离散采用三步三阶预估校正格式,涡核模型采用经典Lamb-Oseen模型,并考虑了涡耗散效应和拉伸效应.首先,通过气动载荷和叶尖涡涡量平均值的分析得到初始涡核半径的取值范围;然后,根据叶尖涡耗散特性的分析,确定体现涡黏性耗散效应涡核半径增长的经验常数的取值;最后,分析了涡核尺寸对叶尖涡结构的影响,进一步验证初始涡核半径和涡黏性耗散经验常数的取值对风力机气动计算的影响.结果表明：当初始涡核半径大于50%弦长时,FVW方法收敛稳定且能准确预估风轮气动载荷;综合风轮气动载荷和叶尖涡耗散特性,初始涡核半径取60%到70%弦长为宜,且对应的涡黏性耗散经验常数取值也不同;风轮气动载荷和叶尖涡结构的准确预估主要受初始涡核半径影响,经验常数对其影响不大,而经验常数主要影响风轮下游尾流场叶尖涡的耗散特性.</p> <br/>力学学报. 2019 51(5): 1530-1537. <br/><p>涡核模型中的涡核尺寸对自由涡尾迹(free vortex wake,FVW)方法准确预估风力机气动特性至关重要,涡核尺寸包括初始涡核半径和由于耗散效应涡核半径在尾迹中的增长.FVW方法中涡线控制方程离散采用三步三阶预估校正格式,涡核模型采用经典Lamb-Oseen模型,并考虑了涡耗散效应和拉伸效应.首先,通过气动载荷和叶尖涡涡量平均值的分析得到初始涡核半径的取值范围;然后,根据叶尖涡耗散特性的分析,确定体现涡黏性耗散效应涡核半径增长的经验常数的取值;最后,分析了涡核尺寸对叶尖涡结构的影响,进一步验证初始涡核半径和涡黏性耗散经验常数的取值对风力机气动计算的影响.结果表明：当初始涡核半径大于50%弦长时,FVW方法收敛稳定且能准确预估风轮气动载荷;综合风轮气动载荷和叶尖涡耗散特性,初始涡核半径取60%到70%弦长为宜,且对应的涡黏性耗散经验常数取值也不同;风轮气动载荷和叶尖涡结构的准确预估主要受初始涡核半径影响,经验常数对其影响不大,而经验常数主要影响风轮下游尾流场叶尖涡的耗散特性.</p> <br/>力学学报. 2019 51(5): 1530-1537. 自由涡尾迹方法中涡核尺寸对风力机气动计算的影响 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1530-1537. article doi:10.6052/0459-1879-18-440 10.6052/0459-1879-18-440 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-18-440 1530-1537 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-025 <br/><p>金属材料的辐照脆化问题一直以来都是核能安全领域亟待解决的关键问题之一.为了更准确地预测金属材料的辐照脆化行为,基于Johnson-Cook本构模型,将未辐照金属材料的断裂真应力取作辐照材料的断裂真应力,建立了通过辐照退火态金属材料屈服强度就能够预测其整个真应力$\!$-$\!$-$\!$应变曲线,以及断裂真应变的辐照脆化模型.实验研究了不同中子剂量辐照退火态高纯铝的准静态拉伸真应力$\!$-$\!$-$\!$应变曲线、断裂真应力和断裂真应变随辐照剂量的变化规律.结果表明,辐照剂量越高,高纯铝的屈服强度越高,断裂真应变越低,但断裂真应力几乎不变.通过TEM显微分析获得了高纯铝内部辐照缺陷的尺寸和数密度随辐照剂量的变化规律,结果表明,辐照剂量越高,孔洞的尺寸和数密度越高,但位错环尺寸和数密度始终很小,难以准确统计.由辐照高纯铝实验数据拟合得到了辐照脆化模型所需参数,并检验了该模型的预测效果.结果表明,无论是通过实验还是显微分析得到辐照高纯铝的屈服强度,模型的预测结果均能够与实验结果较好地吻合,且模型对退火态高纯铝临界中子剂量的预测值也与文献结果一致.</p> <br/>力学学报. 2019 51(5): 1538-1544. <br/><p>金属材料的辐照脆化问题一直以来都是核能安全领域亟待解决的关键问题之一.为了更准确地预测金属材料的辐照脆化行为,基于Johnson-Cook本构模型,将未辐照金属材料的断裂真应力取作辐照材料的断裂真应力,建立了通过辐照退火态金属材料屈服强度就能够预测其整个真应力$\!$-$\!$-$\!$应变曲线,以及断裂真应变的辐照脆化模型.实验研究了不同中子剂量辐照退火态高纯铝的准静态拉伸真应力$\!$-$\!$-$\!$应变曲线、断裂真应力和断裂真应变随辐照剂量的变化规律.结果表明,辐照剂量越高,高纯铝的屈服强度越高,断裂真应变越低,但断裂真应力几乎不变.通过TEM显微分析获得了高纯铝内部辐照缺陷的尺寸和数密度随辐照剂量的变化规律,结果表明,辐照剂量越高,孔洞的尺寸和数密度越高,但位错环尺寸和数密度始终很小,难以准确统计.由辐照高纯铝实验数据拟合得到了辐照脆化模型所需参数,并检验了该模型的预测效果.结果表明,无论是通过实验还是显微分析得到辐照高纯铝的屈服强度,模型的预测结果均能够与实验结果较好地吻合,且模型对退火态高纯铝临界中子剂量的预测值也与文献结果一致.</p> <br/>力学学报. 2019 51(5): 1538-1544. 中子辐照金属材料的脆化模型研究 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1538-1544. article doi:10.6052/0459-1879-19-025 10.6052/0459-1879-19-025 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-025 1538-1544 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-074 <br/><p>土壤材料是一种典型的摩擦型材料,然而天然岩石却具有一定的凝聚力,而金属材料则完全是凝聚型材料. 在分析三种典型的材料强度准则表达式基础上,即SMP,Lade-Duncan以及广义Von-Mises准则,通过利用应力张量的不变量表达形式,提出了一种扩展准则即VML准则,该准则能够分别退化为上述3种典型准则. 在偏平面上,新准则能够描述从曲边三角形到圆形在内的多种开口形态;在子午面上,采用幂函数作为破坏准则公式,能够描述静水压力对于强度特性影响的非线性性质. 而对于土壤的屈服性质,岩土材料具有典型的压剪耦合特性,因此,为了描述剪切与等方向压缩两种路径下的体积耦合现象,采用水滴型屈服面作为屈服准则. 对于偏平面上的截面形状,讨论了给定球应力下偏应力强度值的分布形式及特点,讨论了应力罗德角对于偏平面上强度曲线的凹凸性的影响. 最后,通过多种材料的破坏与屈服试验成果,用所提新准则进行了验证. 通过强度以及屈服特性测试对比,验证了所提VML准则的合理性.</p> <br/>力学学报. 2019 51(5): 1545-1556. <br/><p>土壤材料是一种典型的摩擦型材料,然而天然岩石却具有一定的凝聚力,而金属材料则完全是凝聚型材料. 在分析三种典型的材料强度准则表达式基础上,即SMP,Lade-Duncan以及广义Von-Mises准则,通过利用应力张量的不变量表达形式,提出了一种扩展准则即VML准则,该准则能够分别退化为上述3种典型准则. 在偏平面上,新准则能够描述从曲边三角形到圆形在内的多种开口形态;在子午面上,采用幂函数作为破坏准则公式,能够描述静水压力对于强度特性影响的非线性性质. 而对于土壤的屈服性质,岩土材料具有典型的压剪耦合特性,因此,为了描述剪切与等方向压缩两种路径下的体积耦合现象,采用水滴型屈服面作为屈服准则. 对于偏平面上的截面形状,讨论了给定球应力下偏应力强度值的分布形式及特点,讨论了应力罗德角对于偏平面上强度曲线的凹凸性的影响. 最后,通过多种材料的破坏与屈服试验成果,用所提新准则进行了验证. 通过强度以及屈服特性测试对比,验证了所提VML准则的合理性.</p> <br/>力学学报. 2019 51(5): 1545-1556. 一种适用于岩土的扩展强度及屈服准则 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1545-1556. article doi:10.6052/0459-1879-19-074 10.6052/0459-1879-19-074 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-074 1545-1556 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-252 <br/><p>在钱学森先生领导下,于1957年创刊的《力学学报》作为力学学科历史最悠久的综合性期刊,她见证了20世纪50年代以来中国力学学科的发展轨迹,也见证了中国几代力学人才的成长和发展.本文采用文献计量指标,对《力学学报》从创刊至2018年的产出及学科影响力情况进行分析,用数据说明她对推动我国力学学科的发展、学术成果交流所做出的贡献.本文首先梳理了《力学学报》历年载文量、发文机构、发文作者等情况;然后分析了该刊在力学类期刊中的排名位置;最后,对载文内容特点进行分析,包括主要学科分布及近年刊文的研究热点.</p> <br/>力学学报. 2019 51(5): 1557-1564. <br/><p>在钱学森先生领导下,于1957年创刊的《力学学报》作为力学学科历史最悠久的综合性期刊,她见证了20世纪50年代以来中国力学学科的发展轨迹,也见证了中国几代力学人才的成长和发展.本文采用文献计量指标,对《力学学报》从创刊至2018年的产出及学科影响力情况进行分析,用数据说明她对推动我国力学学科的发展、学术成果交流所做出的贡献.本文首先梳理了《力学学报》历年载文量、发文机构、发文作者等情况;然后分析了该刊在力学类期刊中的排名位置;最后,对载文内容特点进行分析,包括主要学科分布及近年刊文的研究热点.</p> <br/>力学学报. 2019 51(5): 1557-1564. 《力学学报》的产出及学术影响力分析 Personal use only, all commercial or other reuse prohibited . 2019 51(5): 1557-1564. article doi:10.6052/0459-1879-19-252 10.6052/0459-1879-19-252 力学学报 51 5 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-19-252 1557-1564

map