留言板

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

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

倒刺型接触超材料力学性能研究

田耕鑫 曹升虎 张健

田耕鑫, 曹升虎, 张健. 倒刺型接触超材料力学性能研究[J]. 应用数学和力学, 2024, 45(9): 1172-1181. doi: 10.21656/1000-0887.440285
引用本文: 田耕鑫, 曹升虎, 张健. 倒刺型接触超材料力学性能研究[J]. 应用数学和力学, 2024, 45(9): 1172-1181. doi: 10.21656/1000-0887.440285
TIAN Gengxin, CAO Shenghu, ZHANG Jian. Study on Mechanical Properties of Barbed Contact Metamaterials[J]. Applied Mathematics and Mechanics, 2024, 45(9): 1172-1181. doi: 10.21656/1000-0887.440285
Citation: TIAN Gengxin, CAO Shenghu, ZHANG Jian. Study on Mechanical Properties of Barbed Contact Metamaterials[J]. Applied Mathematics and Mechanics, 2024, 45(9): 1172-1181. doi: 10.21656/1000-0887.440285

倒刺型接触超材料力学性能研究

doi: 10.21656/1000-0887.440285
详细信息
    作者简介:

    田耕鑫(1999—),男,硕士生(E-mail: 1300907140@qq.com)

    通讯作者:

    曹升虎(1972—),男,副教授,博士,硕士生导师(通讯作者. E-mail: caosh@xaut.edu.cn)

  • 中图分类号: O34

Study on Mechanical Properties of Barbed Contact Metamaterials

  • 摘要: 受到一些植物茎部和部分猫科动物舌头上的倒刺结构在不同方向上表现出差异性的启发,该研究设计了一种可重复使用且易恢复的倒刺型超材料,并对其力学性能进行了理论分析和数值分析. 研究结果表明,在倒刺的往复运动中,当选取矩形截面高×宽为1 mm×1 mm,长为20 mm,且与竖直方向呈60°角的倒刺时,与阻挡凸台正向接触过程中的最大反力约为逆向接触过程中最大反力的20倍左右,而所消耗能量则相差约200倍. 当倒刺与竖直方向的角度减小时,倒刺结构吸收能量更多,恢复所需能量更少;当倒刺长度增大时,倒刺结构吸收能量减小,恢复时所需能量也减小. 这表明该结构具有出色的抗冲击性能和能量吸收能力,正逆向所需能量差异较大,结构更容易恢复,并且可以通过对倒刺角度和长度设计提升吸能效率.
  • 超材料不同于传统意义上的普通材料,它们是人工制造的一类材料,无法在自然环境中形成. 力学超材料是一种人工结构体,其可以具有与常规力学性质相反的特点,这是由其微观结构单元的几何形状所决定[1],而非其组成部分的性质. 超材料具备多种独特的力学性能,例如可编程性[2]、可重构性[3]、负刚度特性、复杂的双稳定性以及可调刚度[4-5]. 这些独特的力学性能使超材料在吸能领域具有广泛的应用潜力,并为解决工程中复杂的吸能问题提供了新的可能性[6-7].

    优秀的吸能材料应具备轻质、高能量密度和可重复使用的特点. 然而,传统吸能材料往往无法满足这些要求,如金属的塑性变形和陶瓷的破裂等[8]. 主要原因在于它们过于依赖材料本身的特性和破坏变形来实现吸能任务[9-10],无法实现材料的可重复使用,并且在生产和日常生活中面临不同的吸能需求. 因此,针对依赖于材料属性的吸能材料(如金属材料)[11],需要进行重新设计,但这通常需要耗费大量时间. 为了解决每次面对不同吸能任务时需要重新设计材料的问题,我们需要开发一种新型的吸能材料,以能够快速适应不同吸能任务的需求.

    超材料的广泛应用为解决这个问题提供了解决方案. 超材料由人工设计和组装的微观结构单元构成,可以通过排列这些微观结构单元来实现各种吸能任务[12-14]. 进一步设计微观结构单元也可以解决塑性材料吸能材料无法重复使用的问题,例如卡扣型超材料和多稳态吸能材料. 然而,对这些材料分析后可知,它们在吸能和恢复过程中的能量差异不大. 换句话说,在恢复过程中需要消耗相对较多的能量. 为了解决恢复过程中能量消耗较多的问题,需要找到新的方法以节省成本. 通过观察自然界中的动植物,我们发现了一种特殊的倒刺结构,该结构在不同方向上施加相同载荷时产生不同的力学响应. 图 1展示了自然界中三种常见的倒刺结构.

    图  1  自然界中具有倒刺结构的生物
    Figure  1.  Organisms with barbed structures in the nature

    受到这些结构的启发,本文提出了一种新型材料,即倒刺型摩擦吸能超材料,通过结构变形来实现能量吸收. 这种超材料不仅具有轻质、高能量密度和可重复使用的优点,还具备超材料的可编程性和可恢复性. 与吸能过程相比,它在恢复过程中所需能量较少,同时具有更高的吸能性能. 设计这种材料的灵感来源于某些植物茎部和部分猫科动物舌头上的倒刺结构. 某些猫科动物为了吃到骨头上的肉,会竖起舌头上的倒刺增加接触力,从而舔舐下肉. 然而,倒刺不总是竖起的,当它们不需要捕猎或进食时会收起倒刺,以减少接触力,避免伤害到自己. 因此,在梳理毛发时,它们会收起倒刺. 这种结构正是模拟了这种舔舐过程.

    笔者受生物体倒刺结构的启发,设计了如图 2所示的倒刺吸能结构. 该结构由内筒和外筒组成,其内径和外径不同. 内筒上均匀排列着倒刺,并在外筒内侧边缘设计了凸台,用于阻止倒刺向下运动. 为了提高内筒在向下运动时的稳定性,在内筒和外筒之间设计了滑轨和滑槽.

    图  2  倒刺吸能超材料示意图
    Figure  2.  Schematic diagram of the barb-absorbing metamaterial

    当系统的上部受到冲击或负载时,内筒会向下运动,倒刺受到凸台的阻挡,从而发生变形以吸收外部加载的能量. 随着内筒的向下运动,倒刺逐个与外筒内的凸台接触,重复这一接触过程,直到外部加载无法完成倒刺的接触变形为止. 经过变形的倒刺通过往复振动将储存的能量耗散掉. 当内筒运动到底端时,吸能过程完成. 完成吸能过程后,当施加向上的外力于内筒时,由于倒刺结构的特殊性,只需施加很小的力即可使内筒恢复到原始状态,完成恢复过程.

    图 3所示,为了简化分析,我们选择如图 2中的一个倒刺与凸台进行接触. 由于这两个构件是相对运动的,可以固定倒刺的一端,让凸台向上(或向下)运动并与倒刺接触变形. 进一步简化,可以将倒刺和凸台视为平面结构. 因此,可以考虑一个悬臂梁,其左端固定,长度为L,与竖直方向夹角为φ,与一个厚度可以忽略不计的矩形凸台进行接触分析. 在这种情况下,X表示凸台与倒刺末端的距离.

    图  3  倒刺吸能结构简化平面图
    Figure  3.  The simplified plan of the barbed energy absorption structure

    在Euler-Bernoulli梁理论下(即只考虑弯曲形变,忽略横截面的切应变)进行理论分析. 在此假设中,认为凸台是一个刚性部件,在受负载和恢复过程中与其他部件接触时不发生变形. 结构的能量吸收能力主要受到梁的弯曲影响[15-21]. 因此,在分析倒刺结构的弯曲行为时,只需要考虑梁的弯曲部分.

    在凸台向上运动与倒刺接触过程中建立如下模型,如图 4所示,由于凸台向上运动极为缓慢,可以将其看作静力问题. 为了分析此问题,以梁的固定端为原点杆方向建立坐标系xOy,由于杆受到力F作用,杆件自由端在变形过程是不确定的. 为了使计算简便,建立一个恒定不变的坐标系,以A为原点沿着梁方向建立曲线坐标系AO,使用s的值可以表示在任意变形状态下杆件同一位置,在As=0,Os=L. 在梁变形过程中,坐标系AO随着梁一同变形,如图 4所示,图中悬臂梁末端在刚体作用下受到力F的作用出现图示变形. Fx轴夹角为γF指向x负半轴时γ取零,顺时针增大,逆时针减小,取左侧截面,上侧受拉弯矩为正. α是梁上一点s变形后,切线与梁末端受力载荷方向的夹角.

    图  4  倒刺接触变形图
    Figure  4.  The barb contact deformation diagram

    图 4xy分别是在跟随着梁的坐标系中,坐标值s处对应的点在xOy坐标系中的横纵坐标.

    力学方程  梁上任意一点s所受弯矩为

    $$ \begin{equation*} M(x, y)=F(a-x) \sin \gamma+F(b-y) \cos \gamma , \end{equation*} $$ (1)

    ab是梁自由端A点的坐标.

    几何方程   梁自由端沿梁建立坐标系s该坐标系随着梁一同变形,这个坐标系可以保证梁上每一个点对应每一个s,若取梁上任一点s处的微段ds,有以下方程:

    $$ \begin{equation*} \frac{\mathrm{d} x}{\mathrm{~d} s}=-\cos \beta, \frac{\mathrm{d} y}{\mathrm{~d} s}=-\sin \beta, -\kappa=\frac{\mathrm{d} \beta}{\mathrm{~d} s}, \end{equation*} $$ (2)

    κ(s)是该点处的曲率,β(s)是该点处梁的切线与x轴的夹角.

    对式(1)两边对弧长s求导可得

    $$ \begin{equation*} \frac{\mathrm{d} M(x, y)}{\mathrm{d} s}=F \sin \gamma \cdot \cos \beta+F \cos \gamma \cdot \sin \beta=F \sin (\gamma+\beta) . \end{equation*} $$ (3)

    由Euler-Bernoulli梁理论,梁上任意一点处的弯矩和该点处的曲率成正比:

    $$ \begin{equation*} M=E I \kappa(s), \end{equation*} $$ (4)

    E是弹性模量,I是惯性矩.

    为了得到κ(s)和α(s), 设α0为已知, 即α(0)=β(0)+γ=α0,引入参数$\omega=\sqrt{\frac{F}{E I}}, q=\sin \left(\frac{\alpha_{0}}{2}\right)$, 代入边界条件:

    $$ \alpha(L)=\gamma, \kappa(0)=0 . $$

    在此基础上通过Jacobi椭圆函数进行求解得到κ(s), α(s):

    $$ \alpha(s)=2 \arcsin (q \cdot \operatorname{sn}(\omega s+c, q)) , $$ (5)
    $$ \kappa(s)=-2 q \omega \cdot \operatorname{cn}(\omega s+c, q), $$ (6)

    其中

    $$ c=\mathrm{sn}^{-1}\left(\frac{\sin \left(\frac{\alpha_{0}}{2}\right)}{q}, q\right). $$

    α=β+γ,代入$\frac{\mathrm{d} x}{\mathrm{~d} s}=-\cos \beta, \frac{\mathrm{d} y}{\mathrm{~d} s}=-\sin \beta$可得

    $$ \frac{\mathrm{d} x}{\mathrm{~d} s}=-\cos \beta=-\cos \alpha \cdot \cos \gamma-\sin \alpha \cdot \sin \gamma , $$ (7)
    $$ \frac{\mathrm{d} y}{\mathrm{~d} s}=-\sin \beta=-\sin \alpha \cdot \cos \gamma+\cos \alpha \cdot \sin \gamma . $$ (8)

    令$\xi(s)=-\int_{L}^{s} \cos \alpha \cdot \mathrm{~d} s, \eta(s)=-\int_{L}^{s} \sin \alpha \cdot \mathrm{~d} s$,通过计算整理,则有

    $$ x(s)=\xi(s) \cdot \cos \gamma+\eta(s) \cdot \sin \gamma, $$ (9)
    $$ y(s)=\eta(s) \cdot \cos \gamma-\xi(s) \cdot \sin \gamma. $$ (10)

    因此求解x(s), y(s)转为求解ξ(s), η(s). 将α(s)代入,可得

    $$ \begin{align*} & \xi(s)=\left[\frac{2 \mathrm{E}(q)}{\mathrm{K}(q)}-1\right](L-s)+\frac{2}{\omega}[Z(\omega L+c, q)-Z(\omega s+c, q)] , \end{align*} $$ (11)
    $$ \eta(s)=-\frac{2 q}{\omega} \cdot[\operatorname{cn}(\omega L+c, q)-\operatorname{cn}(\omega s+c, q)] , $$ (12)

    K是第一类完全椭圆积分函数(Elliptic K函数的简写),E是第二类完全椭圆积分函数.

    想要求解在凸台向上运动过程中所受反力F,通过上述理论分析只需得到α0, γab中至少两个参数,便可求解F.

    应变能

    $$ \begin{equation*} U=\int_{0}^{L} \mathrm{~d} U=\int_{0}^{L} \frac{M^{2}}{2 E I} \mathrm{~d} s=\frac{1}{2 E I} \int_{0}^{L}\left[\frac{2 F k}{\omega}(\operatorname{cn}(c, k)-\operatorname{cn}(\omega s+c, k))\right]^{2} \mathrm{~d} s . \end{equation*} $$ (13)

    根据本研究设计的模型进行理论推导和数值模拟,以确保结果的准确性. 本文所设计的模型应用机理是通过结构的大变形来完成吸能任务. 因此在选择材料时,必须考虑到其能否完成大变形,并且在变形后能否恢复原状以实现可重复使用.

    图 5所示,为了确定倒刺所需的材料,通过文献调研选择了TPU材料进行仿真模拟. TPU具有较大的变形伸长率,不容易破坏,并且在本研究中,随着应变的变化,其弹性模量变化较小,可以作为线弹性材料进行计算. 通过对实验数据进行拟合[22],得到了TPU材料的弹性模量为E=104.3 MPa. 观察实验得到的应力-应变数据和数值模拟过程中材料的应力-应变曲线,可以看出仿真所使用的材料与实验数据拟合较好,说明了仿真模拟的可靠性.

    图  5  实验数据拟合图
    Figure  5.  The experimental data fitting curve

    接下来,选取一根长为L=60 mm,截面为1 mm×1 mm的矩形倒刺,并使倒刺与固定端夹角为60°,同时使倒刺与凸台恰好接触即X=0时,进行理论求解和数值分析. 将求解在凸台向上(或向下)运动到不同位置时力F的大小,以及整个过程中所需要的能量.

    表 1表 2中是选取在凸台向上或向下移动时不同位移下得到接触力的理论解和数值解,通过数据对比可以看出,在整个接触过程中理论解与数值解之间虽然很接近但是还是有一些误差,发生这种误差可能是因为剪切效应以及弹性模量变化的影响,通过计算可以得出相对误差小于5%,说明椭圆函数法和数值方法分析梁大变形是相对准确的.

    表  1  吸能过程理论解与数值模拟凸台在不同位置时力的大小和凸台与倒刺脱离时的能量
    Table  1.  Theoretical solution of the energy absorption process and numerical simulation of the force of the rib in different positions and the energy of the rib at the moment of sepration from the barb
    displacement D/mm 1 3 5 6 7 10 15 17.84(separation) Eint/mJ
    FFEM/mN 4.04 12.35 20.01 23.58 33.38 51.97 71.32 72.91 0.644 2
    Ftheory/mN 4.01 12.18 19.71 23.14 32.52 51.06 69.95 70.22 0.627 9
    relative error δ 0.007 0.013 0.015 0.019 0.026 0.027 0.029 0.037 0.026
    下载: 导出CSV 
    | 显示表格
    表  2  恢复过程理论解与数值模拟凸台在不同位置时力的大小和凸台与倒刺脱离时的能量
    Table  2.  The theoretical solution of the recovery process and the numerical simulation of the force of the rib in different positions and the energy of the rib at the moment of sepration from the barb
    displacement D/mm 1 1.5 2 2.14(separation) Eres/mJ
    FFEM/mN 1.53 2.56 4.03 4.52 0.003 30
    Ftheory/mN 1.52 2.54 3.99 4.45 0.003 26
    relative error δ 0.006 0.008 0.010 0.015 0.011
    下载: 导出CSV 
    | 显示表格

    图 6(a)所示的力-位移曲线可以观察到,随着凸台的向上运动倒刺发生弯曲变形,接触力F在逐渐增大到最大值时倒刺与凸台脱离,此后F瞬间变为0,在接触变形过程中会分为图所示3个阶段:在ⅰ阶段倒刺末端受到凸台施加的竖直向上力F作用,可直接进行计算,倒刺末端会随着凸台右端滑移. 当凸台运动到一定位置,倒刺末端转角为30°时进入ⅱ阶段,此时力F方向不变但是接触由点接触变为面接触,因此力-位移曲线会出现一小段阶跃,在计算时由于过程较短使用点接触代替面接触进行计算,凸台继续运动由于凸台尖端与倒刺固定端的水平距离是固定的,通过计算可以得到在倒刺上一点施加竖直向上的力F,该点变形后距固定端的水平距离为凸台尖端与倒刺固定端的水平距离,且转角为30°. 接下来到ⅲ阶段,面接触转为点接触,从曲线上可以看出,力F增大幅度逐渐减小主要是由于力F方向的改变使得力F增大幅度变小. 凸台上端与倒刺接触点向倒刺末端进行滑移直到脱离为止.

    图  6  吸能与恢复过程理论分析与数值模拟力-位移对比图
      为了解释图中的颜色,读者可以参考本文的电子网页版本,后同.
    Figure  6.  Theoretical analysis of energy absorption and recovery processes and numerical simulation of force-displacement curves

    图 6(b)所示的恢复过程,由于倒刺结构角度差异所以只出现一种情况,随着凸台的向下运动,力F在逐渐增大,由于角度的原因. 倒刺只需要一点变形,即可完成恢复过程,所以在接触过程的反力和位移都较小.

    结合上述图表可以看出,在凸台向上向下运动的两个过程中,接触力F最大值出现在倒刺与凸台脱离接触时,它们之间相差16倍. 因此可以得出,在两个运动过程中倒刺最大的变形程度相差较大,即使是接触位移相同时,接触力F也是相差较大的,通过分析可知,两个过程差距主要在于力对整个倒刺结构的作用方向.

    图 6中,Eint是力-位移曲线与x轴所围成的面积,Eres是在恢复过程中力-位移曲线与x轴所围成的面积,容易得出Eint>Eres. 由表 1表 2中数据可以得出Eint的数值是Eres的195倍,结合力-位移曲线可以看出两个过程能量之间差距如此之大的原因:凸台向上运动过程与倒刺接触位移相较与凸台向下运动过程与倒刺接触位移更大,并且在接触位移相同时凸台向上运动时的接触力大于凸台向下运动时的接触力.

    综上所述,由于倒刺结构的设计角度,倒刺型摩擦超材料在吸能过程中表现出更大的变形以及与凸台更长的接触位移,因此会消耗更多的能量;而在凸台恢复过程中,变形相对较小,倒刺与凸台的接触位移更短,因此消耗的能量相对较少. 由此可知,倒刺型摩擦超材料是一种相对理想的吸能超材料,可应用于各种吸能任务的不同场景中.

    通过上述内容,可以得知倒刺型超材料在吸能和恢复过程中消耗能量的差异很大. 接下来,将进一步讨论几何尺寸变化对吸能和恢复过程的影响.

    ① 当倒刺长度不变且与凸台恰好接触时,使倒刺与竖直方向的角度改变时吸能和恢复过程中力-位移曲线以及消耗能量的变化.

    图 7中可以观察到,角度的增大和减小对吸能和恢复过程产生了重大影响. 在倒刺与竖直方向夹角为60°时,出现两种不同的变化. 当夹角大于60°时,在图 7(a)中关于倒刺与凸台接触的吸能过程的叙述中,力-位移曲线的3个阶段随着角度增大变得越来越不明显,并且接触力和接触位移随角度的增大而减小,在恢复过程中力-位移曲线形状没有明显变化,但接触力和接触位移随角度的增大而减小. 当夹角小于60°时,力-位移曲线3个阶段更为明显,并且接触力和接触位移随角度的减小而增大,恢复过程趋势与之相同. 当角度小于60°时,在吸能过程中,力-位移曲线在倒刺与凸台接触即将脱离的阶段显示下降趋势,即倒刺与凸台脱离时的接触力不是最大接触力. 这种情况主要是由于夹角减小和接触力从竖直方向转向水平方向的共同作用所导致. 恢复过程的力-位移曲线趋势没有发生很大变化,但是随着角度增加,恢复所需的最大接触力和接触位移逐渐增大.

    图  7  倒刺与竖直方向不同角度下吸能和恢复方向的力-位移曲线
    Figure  7.  Force-displacement curves of energy absorption and recovery processes under different angles of barbed and vertical directions

    表 3为倒刺与竖直方向不同角度下吸能和恢复所用能量及比能耗. 从表 3中可以看出,角度的变化对吸能和恢复过程中消耗的能量有很大影响. 随着角度增加,吸能过程消耗的能量减少,恢复过程消耗的能量增大,即比能耗减小. 即角度越小材料吸能更多并且更容易回复. 然而,需要注意的是,角度的增加会导致接触位移的增加,即变形程度增大. 本研究使用的材料为TPU,可以在较大夹角下实现变形. 如果使用其他材料,则需要谨慎考虑材料是否会发生破坏的问题.

    表  3  倒刺与竖直方向不同角度下吸能和恢复所用能量及比能耗
    Table  3.  The absorption and recovery of energy and specific energy consumption in barbed and vertical directions at different angles
    angle φ/(°) 40 50 60 70 80
    energy-absorption Eint/mJ 1.592 1 1.059 55 0.644 2 0.372 7 0.093
    energy-recovery Eres/mJ 0.001 9 0.003 0.003 3 0.005 0.007 3
    specific energy consumption Eint/Eres 837.95 353.18 195.21 74.54 12.74
    下载: 导出CSV 
    | 显示表格

    ② 当倒刺角度为60°,且与凸台恰好接触时,改变倒刺长度吸能和恢复过程中力-位移曲线和所消耗能量的变化.

    图 8可以看出,随着倒刺长度的增加,在接触过程的力-位移曲线的整体趋势没有明显变化. 然而当其他条件不变时,倒刺长度的增加会使结构在吸能过程的接触力减小,倒刺与凸台的接触位移增大;在恢复过程中倒刺长度的增加会使的接触力减小,倒刺与凸台的接触位移同时也会减小.

    图  8  不同倒刺长度的吸能和恢复方向的力-位移曲线
    Figure  8.  Force-displacement curves of energy absorption and recovery processes with of different barb lengths

    表 4为不同倒刺长度的吸能和恢复所用能量及比能耗. 根据表 4中数据可以看出,长度变化会让吸能和恢复过程所需能量减小但比能耗增大,说明长度的改变会使结构在吸能过程减小一些能量但恢复过程中会更容易恢复. 需要注意的是,吸能过程能量虽然相对减少程度不是很大,但接触距离增长在应对不同问题时也是需要考虑的因素.

    表  4  不同倒刺长度的吸能和恢复所用能量及比能耗
    Table  4.  Energy absorption and recovery energy and specific energy consumption with different barb lengths
    length L/mm 20 22 24 26 28
    energy-absorption Eint/mJ 0.644 2 0.551 5 0.530 6 0.453 8 0.437 3
    energy-recovery Eres/mJ 0.003 3 0.002 5 0.002 3 0.001 6 0.001 3
    specific energy consumption Eint/Eres 195.21 220.60 230.70 283.63 336.38
    下载: 导出CSV 
    | 显示表格

    经过对植物茎部和猫科动物舌头上的倒刺结构进行研究,本文设计了一种倒刺型超材料. 通过对该超材料进行模型简化的理论分析和数值分析,得出了以下结论:

    1) 使用椭圆函数理论和有限元方法对倒刺的大变形过程进行了比较,验证了这两种方法的可行性.

    2) 倒刺型超材料表现出明显的各向异性,吸能过程中消耗的能量远大于恢复过程所需的能量.

    3) 通过数值模拟不同几何尺寸下的倒刺型超材料发现,倒刺与竖直方向的角度和倒刺长度对吸能恢复过程的能量有重大影响,为倒刺型超材料在不同应用情况下的能量吸收提供了重要参考.

    4) 倒刺型吸能超材料具备可重复使用性,吸能过程完成后更容易恢复到初始状态,可以多次使用,降低了材料更换成本.

    因此,本文设计的倒刺型超材料展现出高能量吸收能力和可重复使用性,为解决各种复杂吸能问题提供了潜在的解决方案. 这项研究对吸能超材料的进一步发展和应用提供了有益的参考.

  • 图  1  自然界中具有倒刺结构的生物

    Figure  1.  Organisms with barbed structures in the nature

    图  2  倒刺吸能超材料示意图

    Figure  2.  Schematic diagram of the barb-absorbing metamaterial

    图  3  倒刺吸能结构简化平面图

    Figure  3.  The simplified plan of the barbed energy absorption structure

    图  4  倒刺接触变形图

    Figure  4.  The barb contact deformation diagram

    图  5  实验数据拟合图

    Figure  5.  The experimental data fitting curve

    图  6  吸能与恢复过程理论分析与数值模拟力-位移对比图

      为了解释图中的颜色,读者可以参考本文的电子网页版本,后同.

    Figure  6.  Theoretical analysis of energy absorption and recovery processes and numerical simulation of force-displacement curves

    图  7  倒刺与竖直方向不同角度下吸能和恢复方向的力-位移曲线

    Figure  7.  Force-displacement curves of energy absorption and recovery processes under different angles of barbed and vertical directions

    图  8  不同倒刺长度的吸能和恢复方向的力-位移曲线

    Figure  8.  Force-displacement curves of energy absorption and recovery processes with of different barb lengths

    表  1  吸能过程理论解与数值模拟凸台在不同位置时力的大小和凸台与倒刺脱离时的能量

    Table  1.   Theoretical solution of the energy absorption process and numerical simulation of the force of the rib in different positions and the energy of the rib at the moment of sepration from the barb

    displacement D/mm 1 3 5 6 7 10 15 17.84(separation) Eint/mJ
    FFEM/mN 4.04 12.35 20.01 23.58 33.38 51.97 71.32 72.91 0.644 2
    Ftheory/mN 4.01 12.18 19.71 23.14 32.52 51.06 69.95 70.22 0.627 9
    relative error δ 0.007 0.013 0.015 0.019 0.026 0.027 0.029 0.037 0.026
    下载: 导出CSV

    表  2  恢复过程理论解与数值模拟凸台在不同位置时力的大小和凸台与倒刺脱离时的能量

    Table  2.   The theoretical solution of the recovery process and the numerical simulation of the force of the rib in different positions and the energy of the rib at the moment of sepration from the barb

    displacement D/mm 1 1.5 2 2.14(separation) Eres/mJ
    FFEM/mN 1.53 2.56 4.03 4.52 0.003 30
    Ftheory/mN 1.52 2.54 3.99 4.45 0.003 26
    relative error δ 0.006 0.008 0.010 0.015 0.011
    下载: 导出CSV

    表  3  倒刺与竖直方向不同角度下吸能和恢复所用能量及比能耗

    Table  3.   The absorption and recovery of energy and specific energy consumption in barbed and vertical directions at different angles

    angle φ/(°) 40 50 60 70 80
    energy-absorption Eint/mJ 1.592 1 1.059 55 0.644 2 0.372 7 0.093
    energy-recovery Eres/mJ 0.001 9 0.003 0.003 3 0.005 0.007 3
    specific energy consumption Eint/Eres 837.95 353.18 195.21 74.54 12.74
    下载: 导出CSV

    表  4  不同倒刺长度的吸能和恢复所用能量及比能耗

    Table  4.   Energy absorption and recovery energy and specific energy consumption with different barb lengths

    length L/mm 20 22 24 26 28
    energy-absorption Eint/mJ 0.644 2 0.551 5 0.530 6 0.453 8 0.437 3
    energy-recovery Eres/mJ 0.003 3 0.002 5 0.002 3 0.001 6 0.001 3
    specific energy consumption Eint/Eres 195.21 220.60 230.70 283.63 336.38
    下载: 导出CSV
  • [1] XU R, HE Y, LI X, et al. Snap-fit mechanical metamaterials[J]. Applied Materials Today, 2023, 30: 101714. doi: 10.1016/j.apmt.2022.101714
    [2] PAN F, LI Y, LI Z, et al. 3D pixel mechanical metamaterials[J]. Advanced Materials, 2019, 31(25): 1900548. doi: 10.1002/adma.201900548
    [3] 王竞哲, 陈保才, 朱绍伟, 等. 圆锥形负刚度超材料吸能性能研究[J]. 应用数学和力学, 2023, 44(10): 1172-1179. doi: 10.21656/1000-0887.440055

    WANG Jingzhe, CHEN Baocai, ZHU Shaowei, et al. Study on energy absorption performances of conical negative stiffness metamaterials[J]. Applied Mathematics and Mechanics, 2023, 44(10): 1172-1179. (in Chinese) doi: 10.21656/1000-0887.440055
    [4] 姬忠莹, 闫昌友, 张晓琴, 等. 仿生取向结构表界面及其摩擦各向异性研究进展[J]. 表面技术, 2018, 47(6): 112-121.

    JI Zhongying, YAN Changyou, ZHANG Xiaoqin, et al. Research advances in biomimetic surfaces with oriented structures and its frictional anisotropy[J]. Surface Technology, 2018, 47(6): 112-121. (in Chinese)
    [5] MA S, SCARAGGI M, YAN C, et al. Bioinspired 3D printed locomotion devices based on anisotropic friction[J]. Small, 2018, 15(1): 1802931.
    [6] KLATT T, HABERMAN M R. A nonlinear negative stiffness metamaterial unit cell and small-on-large multiscale material model[J]. Journal of Applied Physics, 2013, 114(3): 033503. doi: 10.1063/1.4813233
    [7] FANG N, XI D, XU J, et al. Ultrasonic metamaterials with negative modulus[J]. Nature Materials, 2006, 5(6): 452-456. doi: 10.1038/nmat1644
    [8] DING Y, LIU Z, QIU C, et al. Metamaterial with simultaneously negative bulk modulus and mass density[J]. Physical Review Letters, 2007, 99(9): 093904. doi: 10.1103/PhysRevLett.99.093904
    [9] LAKES R. Foam structures with a negative Poisson's ratio[J]. Science, 1987, 235(4792): 1038-1040. doi: 10.1126/science.235.4792.1038
    [10] ABRAMOWICZ W, JONES N. Dynamic axial crushing of circular tubes[J]. International Journal of Impact Engineering, 1984, 2(3): 263-281. doi: 10.1016/0734-743X(84)90010-1
    [11] REID S R, REDDY T Y, GRAY M D. Static and dynamic axial crushing of foam-filled sheet metal tubes[J]. International Journal of Mechanical Sciences, 2011, 28(5): 295-322.
    [12] REDDY T Y, REID S R. Axial splitting of circular metal tubes[J]. International Journal of Mechanical Sciences, 1986, 28(2): 111-131. doi: 10.1016/0020-7403(86)90018-4
    [13] 侯海量, 朱锡, 李伟. 轻型陶瓷/金属复合装甲抗弹机理研究[J]. 兵工学报, 2013, 34(1): 105-114.

    HOU Hailiang, ZHU Xi, LI Wei. Investigation on bullet proof mechanism of light ceramic/steel composite armor[J]. Acta Armamentarii, 2013, 34(1): 105-114. (in Chinese)
    [14] 侯海量, 朱锡, 阚于龙. 陶瓷材料抗冲击响应特性研究进展[J]. 兵工学报, 2008, 29(1): 94-99.

    HOU Hailiang, ZHU Xi, KAN Yulong. Advance ot dynamic behavior of ceramic material under the impact of projectile[J]. Acta Armamentarii, 2008, 29(1): 94-99. (in Chinese)
    [15] CONWAY H D. The nonlinear bending of thin circular rods[J]. Journal of Applied Mechanics, 1956, 23(1): 7-10.
    [16] 裴晓辉. 悬臂梁平面大变形的椭圆函数解[D]. 西安: 西安电子科技大学, 2014.

    PEI Xiaohui. Elliptic function solution to large-deflection problems of cantilever beams[D]. Xi'an: Xidian University, 2014. (in Chinese)
    [17] 刘鸿文. 材料力学[M]. 北京: 高等教育出版社, 2011: 175-189.

    LIU Hongwen. Mechanics of Materials[M]. Beijing: Higher Education Press, 2011: 175-189. (in Chinese)
    [18] SHVARTSMAN B S. Large deflections of a cantilever beam subjected to a follower force[J]. Journal of Sound and Vibration, 2007, 304(3): 969-973.
    [19] GREENHILL G. The Applications of Elliptic Functions[M]. London: Macmillan, 1892: 277-305.
    [20] ARMITAGE J V, EBERLEIN W F. Elliptic Functions[M]. Cambridge: Cambridge University Press, 2006: 324-366.
    [21] HOWELL L L, MIDHA A. A method for the design of compliant mechanisms with small-length flexural pivots[J]. Journal of Mechanical Design, 1994, 116(1): 280-290.
    [22] JONSSON A. Some guidelines for implicit analyses using LS-DYNA[Z]. 4th revision. 2014.
  • 加载中
图(8) / 表(4)
计量
  • 文章访问数:  275
  • HTML全文浏览量:  112
  • PDF下载量:  42
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-09-20
  • 修回日期:  2023-11-25
  • 刊出日期:  2024-09-01

目录

/

返回文章
返回