删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

一种滞弹簧耗能的新型离散元滚动阻力模型研究

本站小编 Free考研考试/2022-01-01


引 言


滚动阻力问题涉及车辆工程[1], 土木工程[2]和农业[3-4]等许多领域. 与滑动阻力相比滚动阻力通常较小[5], 如硬质圆形或球型颗粒相对滚动时滚动摩擦系数一般为10?5 ~ 10?3, 滑动摩擦系数通常为0.1 ~ 1. 从数值来看, 滚动摩擦系数比滑动摩擦系数要小得多, 但是滚动摩擦对颗粒的宏观力学特性有着非常重要的影响. Bardet等[6]在离散元模型中最早考虑了滚动约束的作用, 他们发现不考虑滚动约束作用, 离散元模拟结果得到的力学参数在理论值范围之外. 他们认为接触力偶矩可能起着重要作用, 他们指出当约束颗粒滚动自由度后颗粒体系的内摩擦角增大. Morgan等[7]直接将转动阻尼引入离散元模型进行断层泥的模拟, 他们得到了与实验室评估接近的结果. Sakaguchi等[8]将“滚动摩擦”的概念引入离散元模型, 进行谷仓清空过程颗粒流阻塞试验与数值模拟研究. 他们在离散元计算程序中引入了一个滚动摩擦力矩, 发现常规离散元数值模拟中圆盘颗粒成拱不稳定, 易于破坏, 考虑滚动摩擦后能有效的模拟出物理试验中得到的阻塞成拱现象.

目前在滚动阻力理论及工程应用方面许多****开展了相关研究[9-16]. 然而, 颗粒体系稳定过程中滚动阻力作用机制仍不十分清楚, 与材料[17-18]及滚动阻力有关的许多力学模型和方法仍有待研究.

滚动阻力通常被****们称为“滚动摩擦”. 根据摩擦学理论, 滚动阻力的主要来源是接触表面上的微小滑移、塑性变形、材料的黏滞性、表面附着力和形状效应[19].

1875年, Reynolds[20]发现, 当金属圆柱体在橡胶表面上滚动时, 在竖向压力作用下接触面的切向位移会发生微小差异, 称为微滑移. 类似地, Tabor等[21]研究了弹性范围内, 硬质球体和硬质圆柱体在橡胶软基上滚动时的摩擦作用. 他们发现接触面上的切向力始终小于在滚动过程中存在润滑剂的情况下产生的界面滑移值. 因此, 他们认为微滑移并非弹性范围内滚动摩擦的主要原因. 为了解释这种情况下的滚动特性, 他们提出了弹性滞后是滚动摩擦的主要原因.

Flom和Bueche[22]采用具有黏滞特性的Voigt模型对弹性范围内材料的弹性滞后进行了研究. 根据该模型, 当坚硬球体或圆柱体在软基上以一定速度滚动时, 接触面的后缘将与软基脱离形成接触面压力的不均匀分布. 而滚动阻力力矩正是源自这种不对称分布的接触压力. 他们认为, 相对坚硬的球体在较软的基材上的滚动摩擦力会随滚动速度的变化而改变. 对于黏弹性材料的滚动颗粒, 滚动阻力主要是由接触面上体积变形产生的能量耗散引起的, 而与表面附着力的关系较小[23]. 在该研究中, 弹性滞后模型被当作具有黏弹特性的力学模型, 是速度相关型的. 当滚动速度非常小时, 滚动阻力接近于零. 同样的, 对于在硬质基体材料上滚动的软黏弹性球, 滚动摩擦阻力取决于滚动速度, 当滚动速度为零时, 滚动阻力也将为零[24]. 但是, 当采用这种速度相关型黏弹性模型研究颗粒堆积稳定时, 若颗粒处于稳定状态, 滚动阻力将会消失. 显然, 这种稳定机制是有问题的. 当滚动速度为零时, 由于阻力的消失颗粒体系稳定性会降低, 较颗粒滚动速度不为零时更容易崩塌是不合理的.

Greenwood等[23]对橡胶薄壁管进行了扭转和拉伸共同作用试验, 来研究材料的弹性滞后现象. 试验结果显示, 橡胶具有与加载速度无关的弹性滞后. 当在弹性范围内进行加载和卸载时, 由于应变变化落后于应力, 使应力应变曲线的加载和卸载过程不一致, 形成了闭环, 进而产生了能量耗散. 因此, 对于弹性材料, 弹性滞后引起的滚动阻力包括速度相关型的和速度无关型的两部分.

1979年, Cundall和Strack[25]提出了离散元方法(DEM)并迅速引起关注[26-27], 它可以方便地再现物理试验过程细节, 并且能有效降低试验成本, 目前已成为一种被广泛接受的数值计算方法. 离散元模型中滚动阻力力学模型的建立在颗粒体系力学行为模拟中至关重要.

为了研究在剪切带试验中观察到的颗粒间的巨大空隙和高旋转梯度, Iwashita和Oda[28]在DEM模型中考虑了滚动阻力, 建立了改进的离散元模型(MDEM). 该模型滚动阻力由转动方向的弹簧、黏壶、非承拉节点和摩擦型滑阻器元件表达, 其中由弹簧和黏壶组成的Voigt模型, 是一种速度相关型的阻力元件. 而速度无关型的滚动阻力通常由摩擦型滑阻器元件表达. 研究结果表明该模型可以有效预测剪切带的剪胀行为. MDEM是目前应用比较成功的离散元模型, 在此基础上, 许多****开展了相关的计算应用与模型改进工作[9,29-30]. 在本文称MDEM模型为常规DEM模型.

目前对于滚动阻力与滚动速度之间的相关性仍没有统一的认识, 在离散元模拟中通常采用两种典型的滚动阻力公式[31-33]: 第一种滚动力矩与滚动速度无关, 力矩大小正比于法向接触力, 与滚动方向相反. 当颗粒间接触力恒定时, 滚动力矩是一个恒定值; 第二种滚动力矩与相对角速度成正比, 表现为一种黏滞力特征. 离散元模拟结果显示按照速度相关型黏滞滚动阻力公式不能很好地模拟颗粒的稳定堆积. 而单独按照速度无关型恒定的滚动阻力公式虽然模拟颗粒堆积稳定性有所提高, 但颗粒临近稳定静止时会在平衡位置往复振动, 此时滚动阻力大小不变, 方向正负变化, 微观上无法达到静止状态[19]. 因此, 在离散元算法中, 摩擦型滑阻器元件表达的速度无关型的滚动阻力当颗粒临近静止时要退出工作, 只有速度相关型黏滞力元件工作. 所以, 在颗粒堆积问题离散元模拟中, 弹性滞后引起的速度无关型滚动阻力是建立颗粒临近静止状态颗粒接触力学模型的一个重要思路.

直接采用滚动阻力公式表达的力学元件建立离散元模型应用方便, 但由于机理上缺乏深入认识, 一些滚动阻力参数确定往往通过经验或试算得到. 滚动阻力通常较小, 通过试验方法直接识别滚动阻力参数的难度较大[34].

笔者提出了通过圆形颗粒在刚性平面上滚动停止前的往复摆动现象测量滚动阻力参数的方法,研发了一个颗粒微动力光学试验检测系统[35], 可通过测量颗粒的往复摆动曲线识别Voigt模型的转动刚度与黏滞阻尼参数. 试验研究发现通过识别参数的常规DEM模型计算得到的颗粒摆动位移曲线与试验曲线在临近静止时刻吻合程度剧烈下降, 计算得到的颗粒稳定时间较试验长, 且试验曲线在临近静止时刻出现摆动频率的改变, 常规DEM模型不能从机理上解释这一现象.

为此, 本文基于弹性滞后理论研究建立了一种滞弹簧力学元件, 将与速度无关的材料弹性滞后特性引入, 提出一种滞弹性滚动阻力模型, 以此建立对试验中颗粒临近静止状态滞弹性耗能机理的解释. 与传统DEM模型相比, 改进后的滞弹性滚动阻力模型计算结果与试验结果更为符合, 验证了滞弹簧滚动阻力模型的有效性.


1.
弹性滞后


弹性范围内材料加载卸载时, 应变往往落后于应力, 使得应力?应变加载线与卸载线不重合围成一封闭回线, 形成弹性滞后现象. 与速度相关的弹性滞后现象通常表达为弹性材料的黏性行为, 而与速度无关的弹性滞后源于材料加载与卸载过程应力应变不能同步, 这一现象的主要原因在于材料分子间的相互作用和弛豫时间, 与加载速度无关. 通过薄壁橡胶管的纯拉伸试验, Tabor[21]获得应力?应变滞回曲线. 当外力没有达到极限应变卸载时, 应力?应变曲线上将形成一个转折点(εrev1, σrev1), 如图1所示. 当卸载不完全, 再次加载时, 应变将在最后一次卸载的终点(εrev2, σrev2)继续加载. 定义弹性滞后上升曲线为加载曲线, 下降曲线定义为卸载曲线, 并以(εe, σe)表示弹性极限点. 因此, 可以通过如下幂指数表达式定义应力和应变之间的关系



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-1.jpg'"
class="figure_img
figure_type1 bbb " id="Figure1" />



图
1
弹性滞后示意图



Figure
1.
Schematic diagram of elastic hysteresis



下载:
全尺寸图片
幻灯片


加载







$$ {sigma }_{
m{load}}(varepsilon,{varepsilon }_{
m{rev}})=({sigma }_{
m{e}}-{sigma }_{
m{rev}})frac{varepsilon -{varepsilon }_{
m{rev}}}{{varepsilon }_{
m{e}}-{varepsilon }_{
m{rev}}}+{sigma }_{
m{rev}}$$
(1)

卸载







$${sigma _{{
m{unload}}}}(varepsilon,{varepsilon _{{
m{rev}}}}) = {sigma _{{
m{rev}}}}left[ {1 - {{left(1 - frac{varepsilon }{{{varepsilon _{{
m{rev}}}}}}
ight)}^beta }}
ight]$$
(2)

若加载起始点是(0, 0), 加载公式可简化为







$${sigma _{{
m{load}}}}(varepsilon,{varepsilon _{{
m{rev}}}}) = {sigma _{
m{e}}}{left(frac{varepsilon }{{{varepsilon _{
m{e}}}}}
ight)^beta }$$
(3)

当应变ε ∈ (0, εe), 表示材料处在弹性范围. 定义参数β ∈ (0, 1), 表示应变滞后于应力的程度. β值越接近0, 加载曲线和卸载曲线所包围的面积越大, 弹性滞后所引起的能量耗散也就越大.

为定义在复杂应力下的弹性滞后的能量耗散过程, 可根据加载和卸载构造相应的能量耗散过程. 以单个完整加载或单个完整卸载定义为一个子过程, 将整个荷载过程分为多个子过程组合, 可表达为







$$varepsilon {
m{ = }}Delta {varepsilon _{
m{1}}}{
m{ + }}Delta {varepsilon _{
m{2}}}{
m{ + }} cdots {
m{ + }}Delta {varepsilon _n}{
m{ + }}Delta varepsilon {
m{ = }}sumlimits_{k = 1}^n {Delta {varepsilon _{{k}}} + Delta varepsilon } $$
(4)

其中







$$Delta {varepsilon _{{k}}} = {varepsilon _{{{{
m{rev}}(k)}}}} - {varepsilon _{{{{
m{rev}}(k - 1)}}}}$$
(5)

Δε是不足一个完整子过程的多余应变. 子过程的能量密度表达式可写为

加载







$${e_{{k}}} = int_{{varepsilon _{{
m{rev}}({{k - 1}})}}}^{{varepsilon _{{
m{rev}}({{k - 1}})}} + Delta {varepsilon _{{k}}}} {{sigma _{{
m{load}}}}(} varepsilon ,{varepsilon _{{{{
m{rev}}}}({{k}})}}){
m{d}}varepsilon $$
(6)

卸载







$${e_{{{k + 1}}}} = int_{{varepsilon _{{{{
m{rev}}(k)}}}}}^{{varepsilon _{{{{
m{rev}}(k)}}}} + Delta {varepsilon _{{{(k + 1)}}}}} {{sigma _{{
m{unload}}}}(varepsilon,{varepsilon _{{{{
m{rev}}(k + 1)}}}}){
m{d}}varepsilon } $$
(7)

而整个过程的能量密度可通过每个子过程累加得到







$$begin{split}&{e_{{
m{sum}}}} = int_0^{Delta {varepsilon _1}} {{sigma _{{
m{load}}}}{
m{d}}varepsilon + int_{{varepsilon _{{
m{rev(1)}}}} + Delta {varepsilon _2}}^{Delta {varepsilon _2}} {{sigma _{{
m{unload}}}}{
m{d}}} } varepsilon + cdots +& int_{{varepsilon _{{{{
m{rev}}(n - 2)}}}}}^{{varepsilon _{{{{
m{rev}}(n - 2)}}}} + Delta {varepsilon _{{{n - 1}}}}} {{sigma _{{
m{load}}}}{
m{d}}varepsilon + int_{{varepsilon _{{{{
m{rev}}(n - 1)}}}}}^{{varepsilon _{{{{
m{rev}}(n - 1)}}}} + Delta {varepsilon _2}} {{sigma _{{
m{unload}}}}{
m{d}}} } varepsilon+ int_{{varepsilon _{{{{
m{rev}}(n)}}}}}^{Delta varepsilon } {sigma {
m{d}}} varepsilon end{split}$$
(8)

进一步可简化为







$$begin{split}{e_{{
m{sum}}}} = &sumlimits_{k = 1}^n {left[ {int_{{varepsilon _{{{{
m{rev}}(k - 1)}}}}}^{{varepsilon _{{{{
m{rev}}(k - 1)}}}} + Delta {varepsilon _{{k}}}} {{sigma _{{
m{load}}}}(varepsilon,{varepsilon _{{{{
m{rev}}}}}}){
m{d}}varepsilon } }
ight.}+&left. {{
m{ }}int_{{varepsilon _{{{{
m{rev}}(k)}}}}}^{{varepsilon _{{{{
m{rev}}(k)}}}} + Delta {varepsilon _{{{k + 1}}}}} {{sigma _{{
m{unload}}}}(varepsilon,{varepsilon _{{{{
m{rev}}}}}}){
m{d}}varepsilon } }
ight] + int_{{varepsilon _{{{{
m{rev}}(n)}}}}}^{Delta varepsilon } {sigma {
m{d}}varepsilon } end{split}$$
(9)

$e$表示每个荷载过程的能量密度, $displaystyleint_{{varepsilon _{{{{
m{rev}}(n)}}}}}^{Delta varepsilon } {sigma {
m{d}}varepsilon }$表示不能构成一个子过程的多余应变的能量密度.


2.
滞弹簧与HDEM模型


根据接触方向, 常规DEM接触模型[28]可分为法向接触模型, 切向接触模型和滚动阻力模型, 如图2所示. 滚动方向阻力模型由滚动弹簧、滚动黏壶、摩擦器和非承拉节点组成. 滚动模型所提供的滚动阻力可表示为







$$ M=mathrm{min}({K}_{
m{r}}theta +{c}_{
m{r}}dot{theta },;; {mu }_{
m{r}}{F}_{
m{n}})$$
(10)

式中Kr是弹簧刚度系数, cr是滚动阻尼系数, μr是滚动摩擦系数.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-2.jpg'"
class="figure_img
figure_type1 bbb " id="Figure2" />



图
2
常规DEM模型



Figure
2.
DEM model



下载:
全尺寸图片
幻灯片


从上式可看出, 颗粒发生持续滚动时, 滚动力矩总是等于最大摩擦力矩μrFn, 当颗粒滚动不能持续时, 滚动阻力小于最大摩擦力矩μrFn, 常规DEM模型所提供的阻力值与滚动角速度有关. 滚动角速度越大, 滚动模型所提供的阻力值越大. 滚动弹簧是线弹性的, 不能耗散能量, 动能的耗散仅依靠黏壶黏滞力做功实现. 为表征滚动摩擦中速度无关的摩擦力, 根据上述建立弹性滞后的应力?应变关系, 提出图3所示的滞弹簧元件.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-3.jpg'"
class="figure_img
figure_type1 bbb " id="Figure3" />



图
3
滞弹簧



Figure
3.
Hysteresis spring



下载:
全尺寸图片
幻灯片


滞弹簧的角位移与弹性力不同于常规DEM滚动模型中滚动弹簧的线弹性关系. 参照弹性滞后应力应变表达, 滞弹簧的加载和卸载遵循以下关系

加载







$$ {F}_{
m{load}}(varDelta,{varDelta }_{
m{rev}})=({F}_{
m{e}}-{F}_{
m{rev}})left(frac{varDelta -{varDelta }_{
m{rev}}}{{varDelta }_{
m{e}}-{varDelta }_{
m{rev}}}
ight)+{F}_{
m{rev}}$$
(11)

卸载







$${F_{{
m{unload}}}}(varDelta,{varDelta _{{
m{rev}}}}) = {F_{{
m{rev}}}}left[ {1 - {{left(1 - frac{varDelta }{{{varDelta _{{
m{rev}}}}}}
ight)}^beta }}
ight]$$
(12)

加载时, 若加载起点为原点(0, 0), 加载方程可简化为







$${F_{{
m{load}}}}(varDelta,{varDelta _{{
m{rev}}}}) = {F_{
m{e}}}{left(frac{varDelta }{{{varDelta _{
m{e}}}}}
ight)^beta }$$
(13)

式中Δ表示滞弹簧的位移变形, Δe是滞弹簧弹性变形量的极限值, F表示滞弹簧的恢复力值, Fe是弹簧弹性力的极限值.

滞弹簧耗散的能量可以表示为







$$begin{split}{E_{{
m{sum}}}} = &sumlimits_{k = 1}^n {left[ {int_{{varDelta _{{{{
m{rev}}(k - 1)}}}}}^{{varDelta _{{{{
m{rev}}(k - 1)}}}} + {varDelta _{{k}}}} {{F_{{
m{load}}}}(delta,{varDelta _{{{{
m{rev}}}}}}){
m{d}}delta } }
ight.}+&left. {{
m{ }}int_{{varDelta _{{{{
m{rev}}(k)}}}}}^{{varDelta _{{{{
m{rev}}(k)}}}} + {varDelta _{{{k + 1}}}}} {{F_{{
m{unload}}}}(delta,{varDelta _{{{{
m{rev}}}}}}){
m{d}}delta } }
ight] + int_{{varDelta _{{{{
m{rev}}(n)}}}}}^varDelta {F{
m{d}}delta } end{split}$$
(14)

其中







$$varDelta = {varDelta _{{{{
m{rev}}(k)}}}} - {varDelta _{{{{
m{rev}}(k - 1)}}}}$$
(15)

Esum表示整个滞弹簧运动过程中由于弹性滞后效应而耗散的能量, Δrev表示加载和卸载过程中滞弹簧转折点的位移值, Frev表示加载和卸载过程中滞弹簧转折点的力值.

当滚动颗粒在平衡位置往复摆动时, 滚动过程可分为正向加载、正向卸载、负向加载、负向卸载四个阶段. 根据滞弹簧变形与滚动角速度, 式(11)和式(12)计算卸载与卸载时滞弹簧元件的弹性恢复力, 负向时弹性恢复力取负.

将滞弹簧与黏壶、摩擦器及非承拉节点等元件进行组合, 提出一种新的滞弹性滚动阻力离散元模型(HDEM), 如图4所示. 滞弹簧可以表征颗粒堆积稳定过程中与运动速度无关的滚动阻力耗能.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-4.jpg'"
class="figure_img
figure_type1 bbb " id="Figure4" />



图
4
HDEM模型



Figure
4.
HDEM model



下载:
全尺寸图片
幻灯片



3.
验证



3.1
滞弹簧有效性验证


当一个在刚性平面上的滚动圆柱试件速度减小到一定程度后, 会在一个平衡位置往复摆动. 由于动能耗散, 摆动幅度逐渐减小直至为零, 如图5所示. 常规DEM滚动阻力模型参数可以通过测量试件的摆动来识别[36].



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-5.jpg'"
class="figure_img
figure_type1 bbb " id="Figure5" />



图
5
颗粒自由滚动示意图



Figure
5.
Particle free-rolling on a flat surface



下载:
全尺寸图片
幻灯片


基于常规DEM模型, 自由滚动的圆柱试件的运动平衡方程为






$$sum {{F_{{x}}} = 0} qquadqquadqquad$$







$$mddot x + {c_{
m{s}}}dot x{
m{ + }}{K_{
m{s}}}x + frac{{{M_{
m{r}}}}}{R} = 0$$
(16)







$${M_{
m{r}}} = {J_{
m{z}}}ddot theta + {c_{
m{r}}}dot theta + {K_{
m{r}}}theta qquad $$
(17)

或






$$sum {{M_{
m{z}}} = 0}qquadqquadqquadqquad $$







$$({J_{
m{c}}} + m{R^2})ddot theta + {c_{
m{r}}}dot theta + {K_{
m{r}}}theta + {F_{{x}}}R = 0$$
(18)







$${J_{
m{z}}} = {J_{
m{c}}} + m{R^2}qquadqquadqquadqquad$$
(19)







$${F_{{x}}} = mddot xqquadqquadqquadqquadqquad$$
(20)

其中, θ为滚动角位移, Ks和Kr分别表示切向和滚动弹簧刚度系数, cs和cr分别表示切向和滚动方向阻尼系数, Jz是对接触点的转动惯量, Jc是对颗粒形心处的转动惯量, Fx是水平惯性力, Mr是惯性力矩, R为圆柱试件半径.

方程(18)为有阻尼振动方程, 可得到滚动角位移表达式为







$$theta = {{
m{e}}^{ - xi omega t}}Asin left( {omega t + alpha }
ight)$$
(21)

其中ω为振动圆频率, ξ为阻尼比, A为最大幅值, α为相位角. 而阻尼系数和滚动刚度可表达为







$${c_{
m{r}}} = 2{J_{
m{z}}}omega xi $$
(22)







$${K_{
m{r}}} = {omega ^2}{J_{
m{z}}}$$
(23)

因此, 只需测得颗粒的往复摆动曲线, 按照式(21) ~ 式(23)可识别出振动圆频率ω和阻尼比ξ, 进而识别出阻尼系数cr和滚动刚度Kr.

采用激光位移传感器微振动位移检测试验装置, 可试验测得圆柱体试件往复摆动曲线, 试验装置如图6.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-6.jpg'"
class="figure_img
figure_type1 bbb " id="Figure6" />



图
6
试验装置图



Figure
6.
Experimental device diagram



下载:
全尺寸图片
幻灯片


通过检测试验, 我们可以得到圆柱滚动的时间历程曲线, 如图7所示.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-7.jpg'"
class="figure_img
figure_type1 bbb " id="Figure7" />



图
7
圆柱滚动角位移时程曲线



Figure
7.
Angular displacement variation versus time for a rolling cylinder



下载:
全尺寸图片
幻灯片


为验证滞弹簧元件有效性, 分别使用常规DEM滚动阻力模型与HDEM模型对颗粒的纯滚动过程进行了离散元模拟, 如图8所示. 图9是圆柱体在平衡位置摆动过程的数值模拟结果. 图中常规DEM模型中弹簧角位移与弹性力的线性变化关系, 而HDEM模型中滞弹簧的角位移与弹性力形成滞回环, 与所建立的位移?荷载公式一致.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-8.jpg'"
class="figure_img
figure_type1 bbb " id="Figure8" />



图
8
颗粒转动模型示意图



Figure
8.
Particle rotation model



下载:
全尺寸图片
幻灯片




onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-9.jpg'"
class="figure_img
figure_type1 bbb " id="Figure9" />



图
9
滚动过程



Figure
9.
Rolling process



下载:
全尺寸图片
幻灯片


图10中黑色曲线是采用激光位移传感试验测得的滚动弹簧刚度值, 并使用常规DEM模型进行离散元数值模拟得到的动能衰减包络线. 紫色曲线是将常规DEM模型中阻尼系数放大1.5倍后的动能衰减包络线, 蓝色曲线是HDEM模型模拟出的动能衰减包络线. 在摆动起始时刻将黏壶的阻尼系数设置为0, 可以看出将阻尼系数放大后, 动能衰减曲线近似向下平移. 而HDEM模拟的动能衰减曲线在起始时刻7 s时与蓝色曲线交汇, 交汇前位于蓝色曲线上方, 交汇后位于下方, 说明在临近静止过程中滞弹簧的耗能大于黏壶.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-10.jpg'"
class="figure_img
figure_type1 bbb " id="Figure10" />



图
10
HDEM与DEM模拟动能变化



Figure
10.
Kinetic energy evolution with HDEM and DEM simulation



下载:
全尺寸图片
幻灯片


参数β能反映滞弹簧的耗能能力, 其值越小耗能能力越强. 通过与试验数据进行对比, 经过试算拟合可得β值. 对10个不同材质圆柱形试件测量结果拟合得到的β值如图11所示. 聚氨酯圆柱β平均值为0.844, 铝圆柱β平均值为0.874. 聚氨酯的平均值要小于铝的平均值, 分析原因为聚氨酯材料质地较软, 在弹性滞后过程中会产生更多的能耗.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-11.jpg'"
class="figure_img
figure_type1 bbb " id="Figure11" />



图
11
橡胶圆柱与铝圆柱的β值



Figure
11.
β values for the rubber and aluminum cylinder



下载:
全尺寸图片
幻灯片


虽然材料的弹性滞后耗能同样存在法向与切向. 但由于离散元接触模型中法向与切向刚度比转动方向刚度大得多, 法向与切向振动频率高, 速度相关型的黏壶元件耗能要比滞弹簧耗能大得多, 因此HDEM模型中虑法向、切向滞弹簧与转动方向滞弹簧相比作用不显著, 在模型中可不考虑使用.


3.2
耗能分析


图12为橡胶圆柱自由滚动激光位移传感器试验、常规DEM模型和HDEM模型离散元数值模拟得到的时间?相对位移图. 从图中可以看出, 常规DEM和HDEM在振荡早期的动能衰减差异不大. 但从整体上, 常规DEM模型计算得到的滚动到静止时间较试验结果要长, 显示常规DEM模型不易达到静止稳定; HDEM模型与试验结果更为接近. 因此可以得出HDEM模型在接近静止时具有更强的能量耗散能力, 与试验结果吻合更好. 试验曲线不光滑是由于仪器采集误差造成的.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-12.jpg'"
class="figure_img
figure_type1 bbb " id="Figure12" />



图
12
DEM滚动模型与HDEM模型模拟结果



Figure
12.
Rolling angle versus time for DEM and HDEM model



下载:
全尺寸图片
幻灯片


HDEM模型中有滞弹簧和黏壶两个耗能元件. 为分析速度无关滚动阻力与速度相关滚动阻力关系, 提取图12中数值模拟结果,得到体系总动能变化如图13(a)与黏壶作用耗能变化如图13(b). 可以看出, 体系动能不断衰减, 黏壶能量耗散逐渐减小.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-13.jpg'"
class="figure_img
figure_type1 bbb " id="Figure13" />



图
13
体系动能与黏壶耗能



Figure
13.
Kinetic energy and energy dissipated by the damper



下载:
全尺寸图片
幻灯片


图14(a)是HDEM模型中滞弹簧耗能与黏壶耗能变化. 图14(b)是他们的比值随时间变化趋势. 可以看出随着滚动速度降低, HDEM模型中滞弹簧元件所代表的与速度无关的能量耗散比例越来越大.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-14.jpg'"
class="figure_img
figure_type1 bbb " id="Figure14" />



图
14
滞弹簧与黏壶耗能



Figure
14.
Energy consumption ratio of hysteresis spring and damper



下载:
全尺寸图片
幻灯片



3.3
频率拟合结果


选取4组橡胶圆柱和铝圆柱试件的试验检测数据和常规DEM模型、HDEM模型计算结果数据进行对比(如图15和图16). 由于常规DEM模型表达的振动具有频率不变特性, 常规DEM模型与观测试验结果对比可以看出, 摆动速度接近0时, 试验中圆柱体摆动频率有增大现象; HDEM模型的数值模拟结果与试验结果一致, 在摆动速度接近0时同样表现出频率增大的现象. HDEM模型能够模拟滚动试件临近静止时刻摆动频率变高的现象.



onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-15.jpg'"
class="figure_img
figure_type1 bbb " id="Figure15" />



图
15
橡胶圆柱数值模拟与试验对比结果



Figure
15.
Comparison of the numerical simulation and experimental results for the rubber cylinder



下载:
全尺寸图片
幻灯片




onerror="this.onerror=null;this.src='https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/9//lxxb2021-236-16.jpg'"
class="figure_img
figure_type1 bbb " id="Figure16" />



图
16
铝圆柱数值模拟与试验对比结果



Figure
16.
Comparison of the numerical simulation and experimental results for the aluminum cylinder



下载:
全尺寸图片
幻灯片


从模型计算结果与试验结果对比可以看出, 本文建立的滞弹簧滚动阻力模型能够很好地模拟颗粒滚动速度接近于零状态时的能量耗散过程, 能够对颗粒滚动阻力现象进行合理的滚动阻力机理解释. 建立的滞弹性滚动阻力可为颗粒材料堆积稳定问题研究提供方法.


4.
结 论


本文研究建立了滚动阻力滞弹性表达的HDEM模型, 与常规DEM模型的数值模拟结果进行了比较. 通过圆柱试件在平台上的自由滚动试验验证了该模型的有效性, 并得出以下结论:

(1) HDEM滚动阻力模型模拟结果与试验现象吻合, 能较好地解释颗粒在临近静止阶段的能量耗散特性;

(2) 弹性滞后引起的滚动阻力包括速度相关型的和速度无关型的两部分, 滚动试件临近静止时刻, 与速度无关的阻力成分占比越来越大;

(3) HDEM滚动阻力模型能较好地拟合橡胶材料与铝材料圆柱试件的摆动频率, 且能很好地模拟试验中临近静止时刻频率变高的现象.

相关话题/过程 图片 材料 计算 力学
  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效，10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材，产品类型包含电子书、题库、全套资料以及视频，无论您是考研复习、考证刷题，还是考前冲刺等，不同类型的产品可满足您学习上的不同需求。 ...
    考试优惠券 本站小编 Free壹佰分学习网 2022-09-19
  • 基于深度学习和细观力学的颗粒材料本构关系研究
    引言颗粒材料广泛存在于各类自然环境和工程活动中,人类几乎所有的基础设施都建造于岩土材料这类典型的颗粒材料上,理解和预测此类材料在外载作用下的力学响应具有重要的意义.岩土材料通常由粒径不一,形状各异的矿物颗粒以及颗粒之间的空隙组成.颗粒间的互相滑动会引起不可恢复的塑性变形,因而,颗粒材料在剪切作用下极 ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01
  • 冲击载荷下颗粒材料临边界区域的波动行为及变形特征分析
    引言作为离散颗粒材料的土壤(或岩体),是日常生活中最常见的随机多孔介质,冲击载荷作用下材料的动力学响应对于工程实践和地质灾害防控有着重要的指导意义.注意到应力波的传播行为与介质性质密切相关,基于应力波的行为特征分析可提供一种非侵入式研究颗粒材料性质的途径.因此有非常多的****致力于研究冲击载荷作用 ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01
  • 一种新的玻尔兹曼方程可计算模型构造与分析
    引言近年来,随着对临近空间战略意义认识的加深,各航空航天大国发展了众多临近空间飞行器.为满足长航时、大载荷、高超声速巡航等需求,这类临近空间飞行器通常整体尺寸较大,各部件间差异显著,飞行过程会出现多物理多场耦合、多流区共存的多尺度复杂流动现象,传统的数值模拟方法如基于宏观连续介质假设的纳维?斯托克斯 ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01
  • 数据驱动印度洋海域全局动力学研究
    引言海洋系统天然受到多种环境因素(风速、温度、海水含盐量等)、人类活动和天体运动等不确定扰动的影响,导致其洋流运动复杂而多样[1-2].然而,一些观测和研究发现,其内部蕴含着能够对洋流的运动方向和趋势起限制和导向作用的复杂动力结构和特征,影响甚至主导了表面漂浮物的流动形式和演化路径.例如,经志友等[ ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01
  • 非常规油气资源渗流力学进展专题序
    石油和天然气在我国经济发展中起到了积极支撑作用.随着经济发展的后石油时代到来,我国油气资源对外依存度已超出了能源安全的警戒线.为保障国家能源安全,我国对非常规油气资源的开发力度大幅增加.自1996年以来,我国在页岩油气、煤层气、致密油气以及天然气水合物等非常规油气资源勘探开发领域取得了显著进展.作为 ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01
  • 压电超构材料及其波动控制研究: 现状与展望
    引言超构材料或超材料(metamaterials)是将精心设计的基本单元通过一定的空间排列来实现普通材料所不具有的奇异或反常性能,如带隙、波导、负折射、负模量、负密度、超透镜、声学聚焦、声学隐身和拓扑态等[1-2],已成为一个多学科交叉的前沿研究领域.过去几十多年来陆续出现的左手材料[3]、光子晶体 ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01
  • 含天然气水合物土微观力学特性研究进展
    引言天然气水合物(简称水合物)是由天然气分子与水分子在高压、低温环境下形成的非化学计量笼型结晶物质[1],广泛分布于全球近海大陆架及陆地永久冻土地带[2],资源储量巨大,被认为是21世纪最具潜力的替代能源之一[3].自20世纪开始,世界各国陆续成立了专门的水合物研究机构并开展了相关研究[4],在此基 ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01
  • 页岩凝析气藏相平衡的快速准确计算方法
    引言随着以煤炭、石油和木柴为代表的传统能源消耗带来的巨量的空气污染日益引起全社会的关注,天然气作为一种相对较为清洁的能源得到了广泛的重视[1].近年来在非常规油气藏勘探和开发技术上的迅猛发展,使得具有巨大潜在储量但一直得不到有效开发的以页岩气为代表的非常规油气资源成为了投资热点和研究焦点[2].中国 ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01
  • 挠性航天器太阳翼全局模态动力学建模与实验研究
    引言随着航天领域的发展,现代大型柔性航天器往往安装有诸如大型太阳翼等柔性结构,为了得到可持续的能源,其尺寸日益增大、结构重量越来越轻,太阳翼的弹性振动不可避免地与航天器主体平台的运动相互耦合,这种耦合效应随着太阳翼尺寸的增大显著增强.此时,单个太阳能帆板的模态(假设模态)并不能准确反映整个太阳翼在系 ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01
  • 隧道开挖影响下地层-基础体系的接触力学响应分析
    引言在地层条件软弱、周围建筑物繁多的城市环境中进行隧道的施工,不可避免需要穿越和邻近施工,隧道施工扰动将对既有建(构)筑物的力学状态产生较大影响,进而导致建筑物的沉降、变形,甚至发生破坏[1-3].在城市运行高可靠度的要求下,隧道施工引起既有结构的工程响应预测与安全控制问题受到众多专家****与工程 ...
    中科院力学研究所 本站小编 Free考研考试 2022-01-01

Free考研考试FreeKaoYan.Com
欢迎来到Free考研考试，"为实现人生的Free而奋斗"

© 2020 FreeKaoYan! . All rights reserved.