首页/文章/ 详情

基于 Rocky-DEM 的高压辊磨机粉碎效果影响因素分析

2月前浏览6245

机电君:最近摸索粒子仿真,故找一个案例分享,仅抛砖引玉,为大家提供多一个思路,软件。

案例做了什么:高压辊磨机的粉碎过程是通过高压慢速压缩颗粒床来实现的。案例运用离散元仿真技术来分析高压辊磨机粉碎效果的各相关因素。即对高压辊磨机的工作原理及粉碎过程进行了简要介绍,采用基于离散元仿真方法的Rocky软件,通过正交试验,定量分析辊轮转速、辊轮间隙以及入料粒径对高压辊磨机粉碎效果的影响,主要包括辊轮受力、平均破碎比、矿石与辊轮的能量耗散。


1 工作原理

高压辊磨机的结构如图1所示,主要包括机架、定辊、动辊、施压装置、传动装置和给料装置。高压辊磨机的机架是支撑整个设备的主体结构,承载着辊子、压力系统等部件。定辊和活动辊之间形成一个研磨区域,在它们的相对转动下使物料受到挤压、变形和粉碎。液压系统一方面用于控制活动辊的位置以实现物料的粉碎,另一方面能缓冲过载情况,保护设备。物料从物料仓进入高压辊磨机,在重力和辊面摩擦的作用下通过辊隙被压缩并破碎。被磨细的物料最后从高压辊磨机的两个辊子下方排出。

图 1 高压辊磨机结构


2 Rocky-DEM 仿真处理

2.1 导入仿真模型与参数设定

选用某型号高压辊磨机,辊轮直径1400mm,辊宽1100mm,静态辊隙 21mm。辊子表面装有钢钉,钢钉直径为50mm,长度为30mm。简化后的3D模型如图 2所示。

图 2 高压辊磨机简化模型图

在3D软件中将该仿真模型另存为STL文件,导入Rocky中,进行前处理,仿真流程图如图3所示。

图 3 仿真流程图

本次仿真以铁矿石作为颗粒材料,高压辊磨机选用钢材。表 1、表 2 分别为材料相关力学参数和接触参数,这些参数是基于里约热内卢联邦大学路易斯·马塞洛·塔瓦雷斯教授对中试规模高压辊磨机进行 DEM 仿真研究的标定结果。

表 1 材料相关力学参数

表 2 材料接触参数

2.2 破碎模型的选用

目前Rocky中有两种用于模拟颗粒破碎的模型:Ab-T10破碎模型和 Tavares破碎模型。Ab-T10模型将每个粒子视为单独的实体,当粒子受到定义的动能值时,会立即转化为碎片,更适合于一些需要快速破碎的场景。而Tavares 破碎能模型主要基于破碎能来描述单颗粒的破碎过程,当冲击能大于破碎能时,物料便发生破碎。该模型捕获了粒子碰撞期间发生的各种体破损机制,详细描述了脆性材料断裂机制的适应性,并解释了断裂概率的可变性和尺寸依赖性,而且能方便地对破碎后的颗粒进行粒径分布分析、计算平均破碎比等后处理,直观地展示破碎效果。Tavares破碎模型原理如图4所示。

图 4 Tavares 破碎原理图

参考罗德里格斯等校核的Tavares破碎模型参数来定义物料的破碎参数,具体破碎参数如表3所示。表中E50为破碎能量分布的中值。

表 3 Tavares 破碎模型参数

2.3 平均破碎比计算

平均破碎比可以反映物料在破碎过程中粒度变化的平均程度,是评估高压辊磨机破碎效果的关键指标之一。而中值粒径d50,即将颗粒划分为体积或质量均等的两部分时所对应的粒径,mm,常作为一个特征值表示颗粒样品平均粒径。d50可通过罗逊-拉姆勒(Rosin-Rammler)函数计算得出,该函数通常应用于描述颗粒粒径的分布情况,具体形式如下:

则高压辊磨机的理论平均破碎比 i 为:

式中,D50和d50分别表示破碎前后累计粒度分布百分数达到50%时所对应的粒径,mm。

3 正交试验模拟方案与极差分析

3.1 正交试验模拟方案的确定

本案例通过正交试验方法,考虑影响因素:辊轮转速(A)、辊轮间隙(B)以及入料粒径(C)。每个因素选取三个水平,因素水平表如表4所示,正交试验与模拟结果表如表5所示,为评估高压辊磨机的粉碎效果,主要从辊轮平均受力、平均破碎比以及矿石与辊轮平均能量耗散这三个方面进行分析。

表 4 正交试验因数水平表

表 5 正交试验表与模拟结果

以下是一组仿真试验0.5~0.8 s 的破碎过程图,直观展示了破碎效果和颗粒粒径大小。

(a)0.5 s 破碎仿真图

(b)0.6 s 破碎仿真图

(c)0.7 s 破碎仿真图

(d)0.8 s 破碎仿真图

3.2 正交试验极差分析

表6展示了基于表5中正交试验模拟结果计算得出的辊轮平均受力极差分析。其中,Ki代表i水平对应的试验指标总和,ki则为Ki的平均值。由表 6 可知,以辊轮的平均受力为试验指标,因素重要性顺序为:B>A>C,即辊轮间隙影响最大,辊轮转速次之,入料粒径影响最小。由此得出的最优解为A1B3C1,即辊轮转速30r/min(A1),辊轮间隙25mm(B3),入料粒径为 55 mm(C1)。

表 6 辊轮平均受力极差分析表

根据正交试验模拟结果,得到矿石平均破碎比极差分析表如表7所示:

表 7 破碎比极差分析表

由表7可知,以矿石平均破碎比为试验指标,因素重要性顺序为:C>A>B,即入料粒径影响最大,辊轮转速次之,辊轮间隙影响最小。由此得出的最优解为A2B1C3,即辊轮转速45 r/min(A2),辊轮间隙21 mm(B1),入料粒径为65mm(C3)。

由正交试验模拟结果计算得到矿石与辊轮平均能量耗散极差分析表如表8 所示。由表8可知,以矿石与辊轮平均能量耗散为试验指标,因素重要性顺序为:B>A>C,即辊轮间隙影响最大,辊轮转速次之,入料粒径影响最小。由此得出的最优解为A3B1C3,即辊轮转速60r/min(A3),辊轮间隙 21mm(B1),入料粒径为65 mm(C3)。

表 8 矿石与辊轮平均能量耗散极差分析表


4 单因素对矿石破碎效果的影响

4.1 辊轮转速对破碎效果的影响

为研究辊轮转速对矿石破碎效果的影响,选取30、45、60r/min 三组不同转速进行分析对比,辊轮间隙为21mm,入料粒径为65 mm。在以上参数条件下进行模拟仿真,得到模拟结果如表9所示,不同辊轮转速下矿石破碎次数随时间变化,矿石破碎后粒径分布如图6、7 所示。

由表9与图6可知,矿石破碎总数、总破碎次数、平均破碎比大小随着辊轮转速的增加而减小;辊轮平均受力随着辊轮转速的增加而增加;矿石与辊轮平均能量耗散、矿石与矿石平均能量耗散随着辊轮转速的增加而增加。这表明在辊轮间隙和入料粒度保持相同的状态下,辊轮转速提高,矿石与辊轮、矿石与矿石间产生的无用功都增加,且增量较大,对平均破碎比的下降有直接影响。由图7可以看出不同辊轮转速下矿石破碎后的粒径分布状态,辊轮转速在30r/min 时的破碎效果最好。

图 6 不同辊轮转速下矿石破碎次数随时间变化图

图 7 不同辊轮转速下矿石破碎后粒径分布

表 9 不同辊轮转速下各试验指标结果表

4.2 辊轮间隙对破碎效果的影响

为研究辊轮间隙对矿石破碎效果的影响,选取21、23、25mm三组不同辊轮间隙进行分析对比,辊轮转速为30r/min,进料粒度为65mm。在以上参数条件下进行模拟仿真,得到模拟结果如表10所示。不同辊轮间隙下矿石破碎次数随时间变化以及矿石破碎后粒径分布如图8、9所示。由表10与图8,辊轮平均受力、平均破碎比大小、总破碎次数、矿石破碎总数随辊轮间隙的增加而减小;矿石与辊轮、矿石与矿石平均能量耗散随辊轮间隙的增加而减小。这表明在辊轮转速、入料粒度保持相同状态下,辊轮间隙变大,矿石与辊轮、矿石与矿石间产生的无用功相对减少,对于破碎比的影响并不明显。由图9可以看出不同辊轮间隙下矿石破碎后的粒径分布状态,辊轮间隙在21mm时的破碎效果最好。

图 8 不同辊轮间隙下矿石破碎次数随时间变化图

图 9 不同辊轮间隙下矿石破碎后粒径分布

表 10 不同辊轮间隙下各试验指标结果

4.3 入料粒径对破碎效果的影响

为研究入料粒径对矿石破碎效果的影响,选取55、60、65mm 三组不同入料粒径进行分析对比,辊轮转速为30r/min,辊轮间隙为21mm。在以上参数条件下进行模拟仿真,得到模拟结果如表11所示。不同入料粒径下矿石破碎次数随时间变化以及矿石破碎后粒径分布如图10、11所示。由表11与图10可知,辊轮平均受力、平均破碎比大小随着入料粒径的增加而减小;总破碎次数、矿石破碎总数随着入料粒径的增加而增大;矿石与辊轮平均能量耗散随着入料粒径的增加而增加;矿石与矿石平均能量耗散随着入料粒径的增加而减小。这表明在辊轮转速、辊轮间隙保持相同的状态下,随着入料粒径的增加,矿石与辊轮之间作用产生的无用功增加,矿石与矿石之间作用产生的无用功减少,对于平均破碎比的影响并不明显。由图11可以看出不同入料粒径下矿石破碎后的粒径分布状态,入料粒径在 65 mm 时的破碎效果最好。

图 10 不同入料粒径下矿石破碎次数随时间变化图

图 11 不同入料粒径下矿石破碎后粒径分布

表 11 不同入料粒径下各试验指标结果


5 结论

1)影响辊轮平均受力的因素重要性顺序为:辊轮间隙>辊轮转速>入料粒径;影响矿石平均破碎比的因素重要性顺序为:入料粒径>辊轮转速>辊轮间隙;影响矿石与辊轮平均能量耗散的因素重要性顺序为:辊轮间隙>辊轮转速>入料粒径。即辊轮间隙对辊轮平均受力及矿石与辊轮平均能量耗散的影响最大,而入料粒径对平均破碎比的影响最大。 

2)辊轮间隙和入料粒径一定时,分析对比30、45、60r/min 三组不同转速,得出辊轮转速在30r/min 时的破碎效果最好,矿石与辊轮、矿石与矿石的平均能量耗散对于破碎效果有直接影响;辊轮转速和入料粒径一定时,分析对比21、23、25mm 三组不同间隙,得出辊轮间隙在21mm时的破碎效果最好;辊轮间隙和辊轮转速一定时,分析对比55、60、65mm三组不同粒径,得出入料粒径在65mm时的破碎效果最好。分析对比三组单因素试验,确定辊轮转速30r/min、辊轮间隙21mm、入料粒径65mm 时,破碎效果最佳,平均破碎比为2.77。 (想知道如何操作,可看下篇,有简单案例操作)


网络整理,仅限内部分享,禁止商用

公众 号:机电君

来源:机电君
Rocky断裂碰撞离散元理论材料传动控制试验
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-29
最近编辑:2月前
ErNan.Chen🍃
硕士 | CAE工程师 即物而穷其理
获赞 31粉丝 44文章 61课程 0
点赞
收藏
作者推荐

不同温度下带阻尼薄板件振动模态实验与仿真

有没有朋友做过有阻尼的模态分析,是否模态添加了阻尼之后与不添加阻尼有很大的区别?据机电君所知,阻尼比如果在10%以内,模态将会相差不大,但是阻尼将会对高频模态影响较大,并且如果关注的频率一定要精确到单位赫兹的话阻尼对模态的影响就比较大,并且与结构相关。希望有实际经验的朋友可以留言或者私信一起讨论阻尼对结构模态的影响。前言带有黏弹性阻尼材料的薄板件常用于各类高速交通运载工具的壁板结构,例如飞机舱壁、磁浮车身壁板等,用于减少车厢壁板的振动和辐射噪声。飞机和磁浮等交通运载工具会长期在相对温度较高或较低的极端温度条件下运行,温度对黏弹性阻尼材料的模量、损耗因子会有显著的影响,尤其在低温下,阻尼材料性能大幅衰减,甚至失效,因此,对带有阻尼材料薄板件在大跨度温度范围内的模态特性进行实验很有必要,分析温度对带阻尼薄板件振动特性的影响规律。识别黏弹性阻尼材料的材料参数随温度和频率的变化规律,是分析带阻尼薄板构件在极端温度下模态特性的关键。随着激光测振系统的发展,基于激光多普勒测振原理的测试得到广泛应用。对于黏弹性阻尼材料的参数识别主要有DMA(动态力学分析)方法和基于悬臂梁振动的复合损耗因子识别法。DMA 方法对阻尼材料样件加载周期变化的力载荷,将位移响应和力进行处理得到材料在特定受力模式下的储能模量和损耗因子,但要进一步推导黏弹性材料的特性总曲线,过程复杂繁琐。悬臂梁法通过测试带阻尼材料的复合样条和不带阻尼材料的基体样条的频率响应曲线,识别出各阶共振频率和复合损耗因子,再通过公式推导出阻尼材料复模量和损耗因子,误差较大且几乎无法得到连续的黏弹性材料的特性总曲线。基于以上,学者开始尝试倒推法来识别阻尼材料的参数,并取得较好的结果。但是该方法没有得到不同温度下的材料参数。案例分析方法:基于激光多普勒测振原理,在环境箱中使用非接触方法测试各个温度下带阻尼板件的振动响应,分析带阻尼薄板件在极端高、低温下的模态特性以及随温度的变化规律。尝试基于优化算法倒推识别黏弹性阻尼材料的特性总曲线,再从特性总曲线中提取不同温度下的随频率变化的频变参数。1 带阻尼薄板件在高低温下模态特性实验1 高低温下模态实验技术带阻尼薄板件的模态测试分析需要测试板件上诸多离散点的传递函数,常用方法有加速度传感器配合移动力锤的方法,也有加速度传感器配合激振器的方法。测试件需要放置于温度箱内,移动力锤法操作困难且逐个敲击测点工作量巨大。激振器法需要大量加速度传感器,带来的附加质量会严重影响测试精度。诸多离散点的传递函数,常用方法有加速度传感器配合移动力锤的方法,也有加速度传感器配合激振器的方法。测试件需要放置于温度箱内,移动力锤法操作困难且逐个敲击测点工作量巨大。激振器法需要大量加速度传感器,带来的附加质量会严重影响测试精度。激光多普勒测振原理是通过测量从振动物体表面反射回来的激光,通过信号处理解调出多普勒频移fd,从而得到物体振动速率 v,进而获得物体的振动特性。其关系式为:2 高低温下的模态实验过程测试样件为带有阻尼材料的薄板件,基底板为平直钢板,尺500mm×500mm,厚度1. 5mm;阻尼材料为高铁车厢内壁使用的水性阻尼浆,单侧涂覆,充分固化后阻尼材料平均厚度3mm。测试板件固定在刚性支座上,连同整个支座置于温度箱中,温度箱正面用透明亚克力板件密封。测试过程中,钢板面朝外,垂直面对激光测试仪。温度箱和激光测振仪的配置如图1所示。图1 温度箱与激光测振仪带阻尼薄板件的约束状态为对边固支、对边自由。激振器自由悬挂于刚性支座,一同放置在环境箱内。力传感器(图2)粘接到薄板阻尼侧,通过激振杆连接至激振器。力信号作为测试传递函数的参考量引出到数据采集系统中。图2 传感器和激振器布置测试温度为-25℃、-10℃、5℃、20℃、35℃、50℃,将环境箱温度调整至各测试温度点,并且保持30min以上再开始测试,以保证板件阻尼材料内部达到设定的温度。在每个温度下,通过激光测振仪透过透明亚克力板扫描板件上49个测点(见图1),经FFT分析得到49组速率和激振力之间的传递函数用于后续的模态参数识别。经过预测试,不带阻尼材料的薄板样件的模态参数和有限元仿真一致性较好,因此,认为实验设置是合理的。带阻尼材料的薄板样件在室温环境下前10阶的模态频率都在300Hz以下,为确保识别各个温度下的模态参数,并留有一定的余量,频率上限定为1000Hz。3 模态实验结果分析将激光扫描得到的各个温度下的49组传递函数导入Pulse Reflex core软件进行模态参数识别,使用有理分式多项式方法(rational partial polynomial)进行总体拟合,提取各个温度下前几阶清晰的模态振型和对应的模态参数。表1和表2列出了带阻尼薄板件在各个温度下的前6阶识别出的模态频率和模态阻尼比。表1 各温度下的模态频率表2 各温度下的模态阻尼在不同温度下带阻尼薄板件前几阶振型基本一致,如图3列出了50℃时的前6阶振型。图3 前6阶模态振型温度对模态频率影响很大,图4显示为各温度下前4阶振型的模态频率随温度的变化规律,温度会显著改变带阻尼薄板件的模态频率,且模态频率与温度呈负相关。图4 各温度下的模态频率根据表2数据,在同一温度下带阻尼薄板件的各阶模态阻尼比基本保持一致。将各个温度下的前几阶模态阻尼比取算术平均值,得到各温度下的平均模态阻尼比,如图5所示。平均模态阻尼比在20℃时最高,其余温度均呈下降趋势,尤其在-25℃低温时下降尤为显著。图5 各温度下的平均模态阻尼比2 阻尼材料参数识2. 1 黏弹材料定常参数识别2. 1. 1 复合损耗因子识别为识别黏弹性阻尼材料在不同温度下的材料损耗因子和模量,首先识别复合样条的复合损耗因子。图6所示是环境箱内的测试装置和样件,悬臂梁测试样件自由长度215mm,宽度10mm,基底梁厚度1mm,阻尼材料厚度3mm。激励位置位于悬臂梁末端,响应位置距固定端64mm,激励和响应信号经采集和处理,可得到振动频率响应曲线,频率上限800Hz。测试温度包括-25℃以及从-20℃到 50℃每间隔10℃一个温度点。典型的频率响应曲线包含若干个清晰的共振峰,如图7中列出的-10℃、20℃和50℃的频率响应曲线。通过第2阶和第3阶共振峰的半功率带宽,可识别出对应的频率和复合损耗因子,如表3所示。 图6 复合损耗因子测试装置图7 频率响应曲线表3 各温度下的损耗因子2. 1. 2 黏弹材料定常参数根据每一个温度下频率和复合损耗因子可得到对应的材料模量和材料损耗因子,数据如表3所示。每一个温度下,均有一组材料模量和材料损耗因子的不随频率变化的单值参数,可用于后续的仿真参数输入。2. 2 黏弹材料不同温度下频变参数识别2. 2. 1 参数识别过程为识别不同温度下阻尼材料随频率变化的模量和损耗因子参数,首先用包含待定参数的黏弹性材料分数导数模型结合温频转换系数方程,推导不同温度下带频率特性的材料参数。将材料参数代入悬臂梁复弯曲刚度表达式,结合悬臂梁强迫振动理论模型,预测不同温度下的悬臂梁的振动频率响应。将预测的频率响应和测试的频率响应构建目标函数,用Pattern Search(模式搜索)的方法寻找目标函数的最小值,即可找到阻尼材料分数导数模型和温频转换方程的参数,使得预测理论值曲线和实测曲线无限接近,这组参数就是识别出的黏弹性材料的特性参数。识别流程图如图8所示。(公式推导以及各个模型什么的不是案例重点,并且其中公式较多,此处省略,需要的可以后台私信)图8 参数识别过程2. 2. 2 黏弹性阻尼材料的特性总曲线选用了-20℃到50℃的8个温度下测试所得的频率响应曲线,以及根据待定参数建立的8组对应温度下的频响预测曲线,构建用于拟合的目标函数,并进行优化拟合。表4给出了优化过程的初始值和约束条件。表4 优化参数和优化结果以折算频率α (T) f作为横轴,以材料模量和损耗因子为纵轴,得到黏弹性材料的特性总曲线,如图9所示。储能模量随折算频率提高而提高,损耗模量和损耗因子存在明显的峰值,在峰值两边会有明显的下降。温频转换系数如图10所示。图9 黏弹性材料特性总曲线图10 温频转换系数2. 2. 3 黏弹性阻尼材料的频变参数得到温频转换系数α (T),将其代入阻尼材料分数导数模型,即可得到对应温度下自变量为f的材料复模量,根据式(4)材料损耗因子亦可得。图11给出了-25℃、-10℃、20℃和50℃温度下阻尼材料储能模量和损耗因子的频变特性曲线。图11 频变参数材料储能模量在1000Hz范围内变化相对平缓,总体趋势是频率增加、储能模量提高。损耗因子显示出更显著的频变特征,20℃时损耗因子中间高两头低,-25℃和-10℃时是由高到低的趋势,50℃时是由低到高的趋势。因此,在研究带阻尼薄板件振动响应时,需要充分考虑阻尼材料损耗因子的频变特性。3 带阻尼板件高低温下振动响应仿真分析3.1 频率响应的有限元验证使用oSH模型(偏置壳单元加六面体单元方法)建立带阻尼薄板件的有限元模型,可达到较高的有限元建模和仿真精度。具体地,薄板用四边形壳单元,厚度1.5mm,使用钢材典型材料参数,如表5所示。阻尼材料使用六面体单元,与薄钢板壳单元在交界面处共节点。薄板壳单元偏置0.5倍的钢板厚度。建立的有限元模型如图12所示,约束条件与1.2节中所述实验一致,即对边自由、对边固支;激振点和49个响应点也与实验一致。表5 有限元模型的材料参数图12 带阻尼薄板件有限元网格根据3.1.1所述的方法和表3中对应温度下的阻尼材料模量和损耗因子,在模型中设置阻尼材料的单值参数。根据3.1.2所述的方法,将图13所示的不同温度下的频变模量和损耗因子转换成对应的TR ( f ) 和 TI ( f ),在模型中设置阻尼材料的频变参数。使用直接频率响应法分析带阻尼薄板件的振动响应,分析频率上限为1000Hz,频率间隔1Hz。在模态实验己有的温度点中,选择-25℃、-10℃、20℃和50℃这几个典型且温度跨度较大的温度点,使用其对应温度下的单值和频变参数进行仿真验证。3. 3 分析结果图13所示为-25℃、-10℃、20℃和50℃下激振点和响应点之间频率响应曲线,包括测试结果、单值参数和频变参数仿真所得的结果。2种参数得到的仿真曲线与测试曲线都有较好的一致性,相比之下,频变参数得到的前几阶频率和幅值更接近测试值。将图13所示频率响应曲线的前2阶峰值的频率和幅值取出,进行误差分析,如表6所示。单值参数和频变参数仿真得到的1阶峰的频率值最大误差分别为3.0Hz和3.0Hz,2阶峰值的频率值最大误差分别为17.8Hz和7.9Hz;单值参数和频变参数仿真得到的1阶幅值最大误差分别为11.9dB和5.3dB,2阶幅值最大误差分别为-20.9dB和2.8dB。使用频变数进行仿真得到前2阶的峰值频率和幅值均比单值参数误差更小,更接近测试结果。因此,基于优化算法的倒推法识别得到的频变参数,由于考虑了参数随频率的变化,可以更准确地用于带阻尼薄板件的振动响应仿真分析。图13 频率响应曲线表6 仿真结果的误差分析4 结语对带阻尼薄板构件在大跨度的温度范围内进行模态参数识别,分析结果表明温度对带阻尼薄板件的模态特性会有显著的影响。温度与模态频率呈负相关,在极端高、低温下,模态阻尼比会有极大的衰减。材料的模量随频率的提高而提高,材料损耗因子在不同温度下都呈现明显的频变特性,且趋势各不相同。将频变参数用于有限元仿真,分析结果显示,频变参数能更准确地预测带阻尼薄板件的振动响应,其频率响应曲线前2阶的频率和幅值的精度比单值参数更高。网络整理,仅限内部分享,禁止商用公众 号:机电君来源:机电君

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈