首页/文章/ 详情

基于Edem⁃Recurdyn联合仿真的火车皮厢底冻黏浮煤清扫机构性能优化分析

1月前浏览807
以HY-0型清扫机清扫机构底面清扫滚筒为研 究对象,利用Edem离散元软件及Recurdyn进行联合仿真,建立冻煤颗粒模型和清扫机构运动学模型对清扫过程进行模拟仿真。 首先综合考虑温度、含水量煤粒度等因素对冻煤黏结参数的影响进而得到冻黏浮煤参数;再将滚筒的三维模型格式导入到Recurdyn软件中,设置其约束并根据清扫参数设置其运动参数。 仿真采用单因素到多因素的研究方法分析清扫滚筒的清扫性能受清扫转速、牵引速度和清扫深度的影响,得到拟合方程并通过Matlab目标值优化得到清扫性能最优解下的各参数值。

1 冻黏浮煤清扫滚筒仿真参数设置

1. 1 冻黏颗粒模型的建立
根据清扫机构采用铣削式清扫,针对车厢底冻黏浮煤清扫的主要执行部件为底部清扫滚筒,其清扫过程如图1所示,结构组成包括滚筒体、螺旋叶片和刀型截齿。
图1 冻煤清扫示意图
对冻黏浮煤进行模型建立的过程主要是将不规则 的具有黏性的冻黏煤颗粒,根据其自身的物理化学性质,将其简化成具有相同性能的黏结颗粒。 Edem软件能够利用离散单元算法对形状复杂的冻煤材料进行简化,将其理想化为相互独立、相互接触和相互作用的颗粒群体。
颗粒单元具有几何和物理两种类型的基本特征。这些特征受到多种因素的影响,主要包括煤粒的含水量、外界温度、下部煤受到的压力以及煤颗粒大小等, 对其进行分析如下:
1. 1. 1 含水量
煤是一种固液气并存的分散系统,在低温条件下,随着含水量的增加,煤颗粒之间会形成连续的水膜,使煤的冻黏强度增加。 一般情况下,湿煤在含水量为16%时黏附最严重,含有杂质时煤的冻黏量可达15%左右。 根据试验得出冻黏浮煤冻黏强度随含水量的变化曲线,如图2所示。
图2 冻黏强度随含水率的变化规律
1. 1. 2 温度
温度是导致车厢冻黏的重要因素,冻黏现象也只会在低温条件下才会产生。 当气温高于0℃时车厢中煤粒周围的水分子是以液态形式存在;当气温低于0℃,煤颗粒周围黏附界面间的水分开始凝结,产生冻黏现象,且与车厢直接接触的煤层会与车厢冻黏。随着温度降低,冻黏现象会逐渐增多;在温度到达一定限度时,冻黏现象基本完成,车厢内部煤粒状态趋于稳定。其变化规律如图3所示。
图3 冻黏强度随温度变化规律
1.1.3 压力
压力作为煤冻黏条件的次要因素,对煤的冻黏强度影响较大,试验得出,通常压力会使煤的冻黏强度增大,在预压力情况下法向冻黏强度是无压力下的6.2倍,切向冻黏强度是无压力下的 5.9倍。案例中的冻黏过程中,压力来源主要是上层煤炭对下层煤炭因自重产生的压力,并且压力方向具有随机性。因此,压力对冻黏的影响作用较小。
1.1.4 煤粒度
不同粒度组成的煤样黏附力不同,颗粒越小黏附力越大,颗粒越大黏附力越小,且不同大小的颗粒黏附力相差较大。随着煤粒度的增加,煤炭冻黏强度与运输设备材料的冻黏强度均逐渐降低。由图4可知,煤粒度小颗粒间隙小,其冻黏界面连续致密,颗粒间的吸附作用强,导致冻黏强度增大,在冻黏现象产生时,颗粒间空隙小容易形成连续致密的煤泥冰界面,颗粒间的接触面积较大,原因在于煤粒度小,其冻黏界面连续致密,颗粒间的吸附作用强,导致冻黏强度大。而煤粒度大的颗粒间会存在空隙缺 陷,不易形成连续的界面,颗粒间的接触面积较小。结合上述微观角度分析冻黏机制,煤粒度较小的煤颗粒容易形成物理、化学吸附作用,进而导致煤粒度小的煤炭冻黏强度大于煤粒度大的煤炭。
图4 不同粒度煤颗粒间的空隙示意图
根据以上参数及相关研究的试验结果,离散元模型建立方法在Edem软件中首先对颗粒大小进行设置,将煤颗粒简化成特定形状、特定尺寸的颗粒。煤粒特性参数如表1所示。
表1 煤粒特性参数
表2 颗粒间约束
通常情况下车厢内部冻黏总高度6-15cm,计算得出颗粒生成预设置总数约为20000个。
1.2 车厢模型建立
颗粒参数设置完成后在软件中建立车厢的物理模型,以C系列敞车为例,最大尺寸为10m×1.8m×2.6m,由于离散元软件运算复杂且耗时较长,所以 在不影响仿真真实性的前提下对车厢模型在长度和 高度方向上进行简化,选择建立box,尺寸设置为1000mm×1800mm×2000mm。添加车厢材料,敞车车厢材质一般为耐候钢,根据材料特性设置参数如表3所示。
表3 材料特性参数
1.3 Bonding模型参数设置
为模拟实际工程中的冻煤黏结情况,在Edem版本中的Phsics添加Bonding模型使颗粒与颗 粒之间具有黏结力。Bonding模型用到的是Hertz-Mindlin黏结接触模型。当需要用有限尺度的黏合剂 黏结颗粒模型时,可以应用此接触模型阻止切向和法向的相对运动;当达到最大 法向和切向应力时,这种结合就被破坏了,此后颗粒作为硬球对彼此产生作用。 这种模型适用于煤层等结构,煤颗粒间法向、切向冻黏模型都是软件自带的。想研究模型的网上找找一大把,都是利用现有软件建立的,对于我们软件使用者意义不大,我们只要知道自己使用的模型是对就好。
表4 Bonding键参数
1.4 清扫滚筒模型的建立
采用三维建模软件建立清扫滚筒模型,包括滚筒主体、螺旋叶片和刀型截齿,简化滚筒内部结构。根据相关理论研究,对滚筒进行受力分析得 到螺旋叶片最优螺旋升角为25°,螺旋叶片采用两侧 对称布置的形式,能够对物料进行集中收拢,方便后续排料;刀型截齿增加滚筒的破碎能力,布置在叶片内侧 的刀座上,布置间距为450mm,沿螺旋叶片内侧布置。 其结构参数如表5所示。
表5 滚筒结构参数
图6所示为滚筒在旋转时与煤炭颗粒接触的力学模型。确定颗粒受到滚筒叶片的切向力和摩擦力自外向内持续运动。
图6 清扫滚筒模型

2 清扫机构力学模型构建

2.1 清扫滚筒刀具三向力模型
刀具在清扫过程中承受的清扫阻力可以分解为沿 刀具轨迹切线方向的切向力Fτ 、沿刀具刀尖与清扫滚 筒刀尖点速度瞬心连线方向的径向力Fn和与清扫滚 筒轴线平行的侧向力Fz,如图7所示。
图7 刀具清扫过程清扫任一瞬时切向力、径向力和侧向力示意图
2.2 清扫滚筒单刀某时刻受力分析
冻黏浮煤清扫过程为旋转铣削,刀具的清扫厚度在不断变化,所受阻力呈动态变化。滚筒的转角和厚度在一个清扫循环中随时间变化而改变。以清扫滚筒中心点为原点建立坐标系,原点为O,原点到刀具刀头的半径为R,滚筒旋转角速度为ω,刀头线速度为vt ,关系式为vt =ωR;滚筒牵引速度为v,转角为φ,清扫深度为H,清扫厚度为h,清扫进距为fp ,关系式为fp = 2πv/ω。以刀尖点线速度切向τ和法向n建立坐标系, 沿刀具轨迹切线方向的切向力为Fτ ,径向力为Fn ,将上述力进行合成得到清扫滚筒纵向垂直平面内清扫阻力合力为Fτn ,清扫滚筒重力为G,刀具刀头速度瞬心。为O′,刀头距离速度瞬心的距离为ρ,瞬心O′到清扫滚筒轴线的距离为r,关系式为r=v/ω,刀头合速度与清扫线速度的夹角为ψ,图8中Rx 、Ry 分别为清扫滚筒轴上的水平和竖直方向上的支撑反力,驱动力矩为 MR 。清扫阻力合力可分解为Fx 和Fy 。根据上述清扫厚度计算公式,利用平衡关系可以得到刀具在清扫过程中某一时刻各力之间的关系式,如式所示。
图8 刀具清扫过程中某一时刻的受力分析

3 清扫参数对清扫阻力的影响

3.1 清扫阻力单因素仿真
为研究清扫机构在清扫过程中所受阻力的影响, 采用联合仿真的方法对其过程进行模拟,具体步骤如下:首先建立清扫机构三维模型并添加相应的材料参数,随后导入到Recurdyn中创建刚性体并添加运动副及运动参数方程;其次将Edem中的耦合文件添加到Recurdyn中,完成耦合后,将模型的零件生成可联合 仿真的wall文件导入到Edem中,如图9所示,最后对不同运动参数进行单因素变量仿真。
图9 清扫过程联合仿真模型
利用Edem和Recurdyn的后处理功能,可以分别得到煤炭颗粒的运动情况和清扫滚筒单刀受力情况, 针对不同的仿真组进行输出得到相应的仿真数据。
滚筒转速是冻黏浮煤清扫机构的重要性能参数,滚筒转速主要影响的是清扫滚筒最外圈线速度,对清扫性能有重要影响。根据其他相关铣削清扫工具的转速参数,转速通常在40-220r/min,因此以滚筒转速 为第一研究对象,控制其他因素恒定,将滚筒转速仿真 组的参数从60r/min开始,间隔20r/min,共设置 5组。
牵引速度作为清扫阻力的影响因素之一,其数值 决定了清扫耗时的长短和清扫滚筒在水平方向上的运动情况,对其进行单因素仿真试验,相关清扫设备沿清 扫方向的移动速度通常在0.04-0.1m/s。因此,为研究牵引速度对单刀清扫阻力的影响规律,设置仿真组,滚筒转速根据上述仿真结果设置为80r/min,牵引速 度以0.04m/s为起始,间隔0.01m/s,共5组。
清扫深度对单刀清扫阻力的影响规律,根据清扫机构清扫滚筒的结构参数,清扫深度范围选取60-140mm,每组仿真深度间隔20mm,共5组,滚筒转速 为80r/min、牵引速度为0.05m/s。具体参数如表6 所示。
表6 仿真组参数设置
清扫过程联合仿真结束后,通过联合仿真软件的分析处理功能得到清扫滚筒刀具的受力图和不同转速下单刀具的切向力峰值,将数据导出并绘制曲线,图10为某一转速下单刀清扫循环切向力图。由图10可以看出,切向力从切入过程中开始到切出冻黏浮煤层 后结束。
图10 清扫过程单刀切向力变化图
在冻黏浮煤清扫过程中,单把刀具的三向力包括切向力、径向力和侧向力,也就是清扫冻黏浮煤阻力在切向、径向和侧向三个方向上的分力,它们决定了清扫阻力的大小。根据仿真结果,相对应地输出了单刀在进行冻煤破碎时的三向力曲线。根据仿真结果发现,刀具在清扫过程中所收到的侧向力极小,其峰值为0-50N,大部分接近于0。其原因在于刀具所受到的侧向力两侧对称。因此,认为清扫阻力主要受切向力和径向力影响,输出其曲线如图11所示。
图11 单刀清扫切向力、径向力峰值平均值随不同因素变化曲线
通过Recurdyn和Edem的联合仿真建立仿真模型,以单因素变量对清扫滚筒刀具清扫切向力进行研究,得到了单因素对清扫刀具在清扫过程中所受阻力的影响规律。随着清扫转速的升高,单刀切向力在一 定范围内波动;随着牵引速度的提高,单刀切向力在0.05m/s达到最低随后一直上升;随清扫深度的增加,单刀切向力逐渐增大。
根据刀具三向受力情况,确定了此清扫过程中,决定清扫阻力的主要力为切向力和径向力,且径向力变化趋势与切向力相同,因此确定切向力为清扫阻力的主要因素。
滚筒转速、牵引速度和清扫深度对刀具所受到的 切向力和径向力变化规律趋势相同且数值相近,与理论计算得到的切向力与径向力峰值的比值相同,说明切向力和径向力决定了清扫机构清扫阻力的大小。因此减小切向力大小以降低清扫机构清扫过程中清扫阻力,提高清扫性能。
3.2 清扫滚筒的性能优化分析
根据对单因素运动参数研究,得到了滚筒转速n、牵引速度v、清扫深度H分别对清扫刀具阻力的影响规律,并得到了清扫阻力最小时的三种参数的取值范围,现利用响应面法对这三种影响因素进行正交仿真试验,得出三因素对清扫滚筒清扫比能耗的综合影响,并利用Matlab多项式拟合的方法得出三因素的匹配模型,得到最优清扫效率下各参数的最优取值。
根据公式推导出清扫过程清扫滚筒的整体清扫比能耗,它是评价清扫滚筒清扫性能的重要指标,代表了清扫滚筒破碎单位体积冻黏浮煤所消耗的能量,花费的能量越少,清扫效率越高,清扫滚筒螺旋叶片和刀具的损耗越小,其计算式为
利用Recurdyn后处理功能导出清扫滚筒平均转矩,计算冻黏浮煤所破碎的体积,最后对17组仿真数据进行计算得到各自的清扫比能耗,如表7所示。
表7 清扫比能耗仿真试验
对清扫比能耗仿真数据进行分析,得到结果,利用Matlab对数据进行2阶拟合,得到清扫比能耗与运动参数的函数关系式为
拟合后匹配得到数学模型,对其可靠性进行分析,以清扫比能耗的置信度和方差为评价指标,判断拟合是否可行。
方差分析如表8所示。模型F值为30.3%意味 着模型显著,只有0.01%的可能性才会出现如此大的F值。对其P值和失拟项分析后,P值小于0.05表示模型项显著,由于噪声,有49.71%的可能性拟合不足F值,是不显著的。因此,模型符合要求,得到的拟合的数学模型能够合理地描述清扫比能耗与三个运动参数之间的影响变化规律,模型准确度较高。
表8 清扫比能耗方差分析
对结果进行拟合如表9所示。由表9可知,预测 R2为0.8107,与调整R2为0.9428合理一致;信噪比最好大于4,实际为16.952,表示信号充足,模型精度较高。
表9 拟合统计
由图12可以看出,残差各点合理的分布在拟合直线附近,说明残差值小,数学模型较合理准确。
图12 回归拟合的残差分布
根据各因素对清扫比能耗的影响规律,导出单因素对清扫比能耗的影响规律,如图13所示。图13(a) 所示为滚筒转速对清扫比能耗的影响,当牵引速度为0.05m/s,清扫深度为100mm时,可以看出随着滚筒 转速的增加,清扫比能耗先缓慢降低后快速升高。 图13(b)所示为牵引速度对清扫比能耗的影响,当清扫转速为80r/min,清扫深度为100mm,可以看出随着牵引速度的增加,清扫比能耗在逐渐降低。图13 (c)所示为清扫深度对清扫比能耗的影响,滚筒转速为80r/min,牵引速度为0.05m/s时,随着清扫深度的增加,清扫比能耗逐渐降低。说明清扫滚筒的牵引速度和清扫深度的增加有助于减少能耗,成反比关系。
图13 单因素对清扫比能耗的影响规律
得到目标值的二次拟合多项式后,对此类无约束优化问题采用多目标粒子群算法寻求拟合函数保证清扫能力的情况下能耗最低的性能参数匹配方案。根据单因素的仿真结果范围设置变量的上下界,迭代次数共50次,迭代结果如图14所示。
图14 粒子群算法优化结果
由图14可知,方程在迭代12次后稳定,最终得到清扫转速为78r/min、牵引速度为0.05m/s、清扫深度为100mm时,清扫滚筒比能耗最低,数值为 0.22kWh/m^3,在此参数下清扫机构性能达到最优。

网络整理,仅限内部分享,禁止商用
公 众号:机电君


来源:机电君
RecurDyn化学MATLAB煤炭离散元理论EDEM材料控制试验
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-02-25
最近编辑:1月前
ErNan.Chen🍃
硕士 | CAE工程师 即物而穷其理
获赞 34粉丝 58文章 73课程 0
点赞
收藏
作者推荐

案例-不同铺层顺序对GO-CF/BF混杂增强复合 材料微观组织及弯曲性能的影响

案例给到了一种研可行的研究方法,思路:采用仿真分析与试验验证对真空浸渗热压成型工艺(VIHPS)制备的GO CF/BF混杂增强复合材料的微观组织与弯曲性能进行了探究。首先用扫描电镜对GO-CF/BF混杂增强复合材料微观组织进行表征,并用COMSOL软件对其浸渗效果进行仿真。然后将DIGIMAT软件建立的代表性体积单元模型导入ANSYS进行三点弯曲仿真。研究机构、学校学生发文章可以适用试试文章的思路。铺层顺序对GO-CF/BF浸渗效果的影响COMSOL仿真仿真流程 利用多孔介质理论研究不同铺层顺序对GO-CF/BF混杂增强复合材料浸渗效果的影响规律。在模拟计算之前,通过扫描提取的方法构建不同的模型。首先,对复合材料进行微观组织形貌表征,从图像中提取纤维的分布。其次,在COMSOL软件中,选择二维模型中的两相流,将“层流”与“相 场”接口相结合。根据提取的复合材料微观组织形貌建立湿相与非湿相的两相流模型,其中湿相是纤维表面涂抹的基体溶液,而非湿相是纤维与纤维间的孔隙。将两相流的属性映射到复合材料相 应的区域上,通过联合体构建和复合边形成,运用 求和与差集操作精确获取目标模型部位。然后选择材料属性并设置流体参数:密度为1.36 g·mm-3, 动力粘度为0.562 Pa·s等,指定模型的进出口,设置入口压力为0.5 MPa,选择其边界、初始值域等, 设置两种纤维边界的浸润角。在多物理场接口中, 设定流体的相态属性,如密度和粘度。随后,系统会自动创建有限元网格。最后计算仿真结果。仿真流程如图所示。通过该方法分别对不同的纤维分布的模型进行模拟,并与微观组织样貌图进行对比分析,从而得到浸渗效果最优的铺层顺序。COMSOL仿真结果基于复合材料的微观组织样貌图,利用COMSOL软件对GO-CF/BF混杂增强复合材料进行仿真模拟。仿真结果如图5所示,其中白色部分为纤维分布,浸渗效果右侧表示固化混合溶液的体积分数, 颜色越深表明固化溶液的体积分数越高,浸渗效果越好。从图(a)中可以看出当铺层顺序为CCCBBB时,大部分区域基体体积分数较高,表明浸渗效果较为理想,但仍有小部分区域浸渗效果较差,使得材料内部产生孔洞缺陷,这主要是因为该处附近的纤维较为密集,出现了团聚现象,基体溶液难以进入该区域。铺层顺序为CBBBCC如图(b)所示,此时上方与下方碳纤维区域固化溶液体积分数相对中间玄武岩纤维区域较小,不同纤维区域间浸渗效果相差较大,从而削弱了复合材料内部的热量与力的传递,使得复合材料整体性能下降。铺层顺序为BBBCCC如图(c)所示,可以观察到碳纤维区域相较玄武岩纤维区域被更多的基体所包裹,表明碳纤维区域的浸渗效果更为充分。铺层顺序为BCBCBC如图(d)所示,从中可见基体溶液在材料内部有较高的体积分数,并且分布均匀,这表明在此铺层顺序下浸渗效果最优。不同铺层顺序GO-CF/BF混杂增强复合材料的COMSOL浸渗仿真图GO-CF/BF混杂增强复合材料微观表征分析由于碳纤维与玄武岩纤维的密度存在差异,在纤维层数相同条件下,玄武岩区域的孔隙率相比碳纤维区域高,因此GO-CF/BF混杂增强复合材料的孔隙率将随着铺层顺序的改变而发生变化,不同铺层顺序的复合材料微观形貌如图所示。当铺层顺序为CCCBBB时,复合材料的微观浸渗形貌如图(a)所示,表面的碳纤维被基体紧密的包裹,从侧面观察,纤维束的端口被基体牢固地粘结,这种结构没有出现明显的孔洞或缺陷。图(b)为GO-CF/BF混杂增强复合材料在CBBBCC铺层顺序下的微观组织形貌,此时在上方的碳纤维表面观察到有少量的基体附着,纤维之间存在明显的缝隙。玄武岩纤维位于中间区域,呈现亮白色。下方的碳纤维被基体紧密包裹,纤维间的缝隙已被基体填充。这是因为在挤压力和重力的作用下,基体自上而下浸渗。并且玄武岩纤维的单丝间缝隙较碳纤维的单丝间缝隙大,导致底部碳纤维的浸渗效果优于顶部。铺层顺序为BBBCCC时,GO CF/BF混杂增强复合材料的微观形貌如图(c)所示。玄武岩纤维主要以簇状结构粘结,少量的玄武岩纤维之间存在缝隙。在虚线下方的碳纤维区域,其端口表面被大量的基体包裹,浸渗效果相比玄武岩纤维区域要更充分,这是由于碳纤维相比玄武岩纤维有更好的吸附能力,从而碳纤维表面粘附大量基体。然而,仍有少量的碳纤维之间存在缝隙。当铺层顺序为BCBCBC时,复合材料的微观浸渗形貌如图(d)所示,玄武岩纤维与碳纤维均被基体紧密包裹,且相邻纤维之间未观察到缝隙缺陷,与 COMSOL仿真结果一致。这主要是因为碳纤维与玄武岩纤维交替排列使得纤维间的孔隙分布更加均匀,且这种纤维交替排列的铺层顺序可能导致基体在浸渗过程中受到不同的压力与阻力,使其内部产生较大的压力梯度,从而使复合材料浸渗充分。不同铺层顺序GO-CF/BF混杂增强复合材料的微观形貌图铺层顺序对GO-CF/BF混杂增强复合材料弯曲性能的影响复合材料失效准则复合材料层合板是由若干层单层板按特定的铺层顺序与铺层角度粘合而成。相比于单层板,层合板的力学性能取决于单层板性能与铺层方式,具有良好的可设计性。在复合材料结构设计的经典层合板理论框架下,首先基于所给定的单层板弹性性能参数及其铺层方式,通过计算得出各层合板的刚度与柔度模量,并进一步确定每一单层板中的应力与应变分布。随后,根据选定的失效准则来评估单层板的破坏载荷,从而确定极限强度。针对GO-CF/BF混杂增强复合材料层合板,当承受的载荷超过了材料的强度时会发生失效破坏。 蔡-吴(Tsai-Wu)准则综合考虑了复合材料在纤维方向的拉伸与压缩强度、垂直于纤维方向的强度,以及材料的纵向剪切强度。该准则在使用时考虑了纤维、基体的拉伸与压缩性能。根据这一准则, 在应力空间中表面破坏存在下列形式,如式所示。蔡-吴强度理论的失效准则函数为复合材料细观力学理论复合材料的宏观力学性能,如刚度与强度,以及其失效机制,均深受其细观结构特征与组分材料的属性所影响。细观力学的建立需具备以下假设条件:(1)单向复合材料具有宏观均匀性、线弹性,无初应力;(2)纤维为横观各向同性、线弹性,基体为均匀性、各向同性、线弹性;(3)纤维与基体间 的粘结效果理想,复合材料所含孔隙与初始微裂纹总体积之和小,可忽略不计[25]。假设复合材料内有一微小的代表性体积单元(Representative Volume Element,RVE),在该单元体积内含有纤维和基体,对于GO-CF/BF混杂增强复合材料,将碳纤维与玄武岩纤维看成圆柱状均匀地排列在 基体材料中,如图所示。根据Eshelby 等效夹杂理论(Eshelby Equivalent Inclusion),现将GO看作夹杂,EP视为基体的一种模型。 当设定GO-CF/BF混杂增强复合材料的边界上受到远场均匀应力为σ0时,另有一种形状与弹性性质与上述复合材料相同的均质材料,在同等外力条件下,基体材料中的平均应力σm为:GO-CF/BF混杂增强复合材料微观三点弯曲力学模型在DIGIMAT软件环境中,为了构建GO-CF/BF混杂增强复合材料的数值模型,首先对复合材料的微观结构进行表征,即分别定义CF、BF、EP和氧化GO不同的微观相。CF与BF设置为Continuous fiber 相、EP设置为Matrix相、GO设置为Inclusion相,并在Continuous fiber相、Matrix相、Inclusion相中设 置相应的参数,具体如下:在Type中将Interface behavior设置为Perfectly bonded,Phase material设置为CF;在 Parameters中将CF的质量分数设置成10%,Diameter定义为7×10-3 mm,Orientation选择Fixed项,其中θ=90°,Φ=0°。同理对BF建立的过程和参数也一样,其中BF的Diameter定义为1.2×10-2 mm。 GO质量分数设置为0.05%,Size设置为0.335 nm,Orientation设置为Random 3D。根据上述建立过程, 在RVE选项中设置6层结构,每一层都重复上述的相定义和参数设置过程,以构建一个完整的RVE模 型。CF、BF、GO和EP在DIGIMAT软件中的相关性能参数如表所示,GO-CF/BF混杂增强复合材料微观力学模型选用Mori-Tanaka模型。GO、CF、BF与EP的相关性能参数为研究不同铺层顺序对GO-CF/BF混杂复合材料弯曲性能的影响,需要构建4种不同铺层顺序的细观力学模型,如图(a)~(d)所示。在这些模型中,CF、BF与GO呈现随机分布,图(e)为局部放大视图,用以展示这一微观结构的特征。不同铺层顺序的GO-CF/BF混杂增强复合材料的RVE模型GO-CF/BF混杂增强复合材料仿真与试验结果分析仿真流程ANSYS中的三点弯曲仿真 主要用于模拟和分析材料在弯曲载荷作用下的行为。将ANSYS软件的前/后处理模块 (ANSYS Composite PrePost,ACP)、Static Structural模块关联共享,根据Tsai-Wu强度准则判定GO-CF/BF混杂增强复合材料层合板试样是否弯曲失效。通过这些综合分析,揭示 了GO-CF/BF混杂增强复合材料在受到弯曲载荷时的力学响应特征,相关流程如下: (1)导入模型: 在ANSYS的预处理模块中导入由DIGIMAT-FE软件生成的层合板实体单元RVE模型。 随后,建立与实验设置相对应的压头和支座的模型,如图所示。GO-CF/BF混杂增强复合材料的三点弯曲仿真模型(2)设置材料属性: 在ANSYS软件中,对RVE模型各部位进行选择并赋予相应的材料属性,其中圆柱体为纤维,球体代表GO,而立方体则代表EP。此步骤的目的是将DIGIMAT软件中各材料的属性与 ANSYS软件相关联。(3)设置接触条件:仿真模型使用"面-面"接触方式来模拟压头与层合板之间的接触。其中压头被定义为接触元件,层合板的上表面则被定义为目标元件。接触属性设置为“摩擦”,摩擦系数值设为0.2。 对于支座与层合板间的接触,同样使用了“面-面”接触模拟,其中支座作为接触体,层合板的下表面作为目标体,接触类型则设定为“无摩擦”,以模拟理想光滑支撑条件。(4)划分网格:鉴于压头与支座被建模为刚体,它们无需进行网格细化。为了确保GO-CF/BF混杂增强 复合材料仿真结果的可靠性,研究采用了四面体网格划分方法,并自动调整网格单元大小。四面体网格能够更好地简化模型,提高计算效率,同时增强数值解的稳定性和收敛性。(5)加载求解:为了模拟准静态加载条件下的固定支撑环境,将支座的6个自由度全部约束。在压头处施加沿层合板法线方向的5mm强制位移,并在预定时间内完成该位移,以复现实验条件下的准静态加载场景。结果分析在层合板弯曲受载过程中,不同铺层顺序的GO-CF/BF混杂增强复合材料抵御外载荷的能力将存在差异。当外部施加的载荷超出了最弱层的临界承载极限时,将引发该层材料断裂失效。随着外载荷的持续增加,这种局部的失效会扩展至整个层合板,进而引发整体结构的破坏。在GO-CF/BF混杂增强复合材料中,碳纤维和玄武岩纤维起承载作用,基体起传递载荷的作用,复合材料的弯曲性能与纤维铺层顺序紧密相关。不同铺层顺序的GO-CF/BF混杂增强复合材料的弯曲应力云图如图所示。 通过分析4组应力图发现,复合材料在压头与支撑座附近均出现了应力集中现象,并向周围递减,这表明试样的失效是从加载位置的拉伸与压缩表面开始的。由图(c)可以看出,当铺层顺序为BBBCCC时,随着载荷增加,碳纤维周围的基体受到很高的应力集中,基体出现裂纹,并逐渐出现裂纹扩展。此时复合材料变形后的弯曲应力峰值最大,为442.89 MPa。CCCBBB、CBBBCC和BCBCBC铺层顺序的复合材料弯曲强度分别为325.18、296.17、336.02 MPa,与BBBCCC铺层相比,分别下降了26.58%、33.13%和24.13%。特别是铺层顺序选择CBBBCC时,如图(b)所示,根据前文对该铺层 的浸渗仿真与微观形貌的研究可知,此时复合材料内部基体含量较低,纤维之间存在明显缝隙缺陷,因此其弯曲强度和弯曲模量相比其余铺层下降幅度极大,变形后的弯曲应力峰值最小,为296.17MPa。不同铺层顺序GO-CF/BF混杂增强复合材料的等效应力云图结合上图可知,在对复合材料施加载荷过程中,试样的中性层上侧受挤压力,下侧受拉伸力, 且距离中性层越远,所受应力越大,中性层处应力最小,接近于0,如图(a)所示。因此,本文主要 分析复合材料受压部分的上、下表面。中性层上侧区域的应力分布情况大致呈现下凸的弧形,下侧 区域与之相反,如图(b)和图(c)所示。铺层顺序为CBBBCC时,在三点弯曲载荷作用下,层合板 内的纤维单元所承受的等效应力较小,所能承受的载荷也随之较小,裂纹扩展的速度也较快。当铺 层顺序为BBBCCC时,压头处的弯曲应力达到峰值,在变形过程中,层合板层内纤维单元等效应力 随着压头加载而增大,但裂纹扩展速度较缓慢,因此失效变形区域呈现越来越弱的趋势。从图(d) 可以看出GO附近的应力较小,这表明了GO的加入可以改善其界面性能,并在其内部生成网状结构, 从而分散应力,使复合材料强度增强。不同铺层顺序GO-CF/BF混杂增强复合材料的局部放大应力云图试验验证三点弯曲试验材料的弯曲性能按照试验标准(GB/T1449—2005) ,采用三点弯曲方法进行测试,跨距L和厚度h比值为32∶1,将试验材料居中放在试验台支座上,保证材料两端与两侧试验台支座中心等距离。 在实验室温条件相同的自然环境下,给预制的试验材料以5mm/min的速度持续加载,直到试验材料 发生明显的脆性断裂为止,记录此时的极限载荷为弯曲强度值。其中加载压头与支座半径均为5mm, 弯曲试样尺寸为80.0mm×12.5mm×2.0mm,跨距取64mm,加载示意图如图所示。铺层顺序为BBBCCC时GO-CF/BF混杂增强复合材料的三点弯曲模型试验验证与结果分析在复合材料层合板的弯曲破坏过程中,其外层受压缩和拉伸的交互作用,这会导致外表面逐渐产生裂纹,且随着应力的增加,裂纹逐渐向内层扩展。在三点弯曲试验中存在着多种失效形式,主要包括:基体开裂、纤维断裂、界面脱粘和层间剥离。由于CF与BF在强度上的差异,不同的铺层顺 序将影响层合板在弯曲试验中的失效和破坏模式。试验结果表明,铺层顺序为BBBCCC时,复合材 料的弯曲强度最大,高达410.59 MPa。CCCBBB、CBBBCC和BCBCBC铺层顺序的复合材料弯曲强度 分别为305.78、277.52、326.31 MPa,与BBBCCC铺层相比,分别下降了25.53%、32.40%和20.53%。 特别是铺层顺序为CBBBCC时,复合材的弯曲强度 最小,低至277.52 MPa。这是因为弯曲受载时,在采 用CBBBCC铺层顺序的复合材料中,压缩侧的主要 承载部分是碳纤维,其具备优越的强度和刚度。然 而,当这些碳纤维发生脆性断裂失效时,载荷逐渐转 移到具有较低模量与刚度的玄武岩纤维上,载荷下 降至玄武岩纤维承载能力水平,其弯曲性能下降。 试验结果与仿真结果的弯曲强度对比图如图13所 示,比较4组模拟数据弯曲强度,当铺层顺序为BCBCBC时,其试验值与仿真值误差最小,此时复合 材料的试验弯曲强度为326.31 MPa;而仿真结果表 明,复合材料的弯曲强度为336.02 MPa。试验与仿真 结果误差小于10%,两者较为符合。比较其余3组模拟弯曲强度,实验结果全部略小于模拟结果,主要原因是仿真过程中材料属性为理想或平均情况下 的数值,与试验材料略有差异,而且在GO-CF/BF混杂复合材料的制备过程中,材料内部可能存在如 孔隙缺陷或浸透不充分等问题,使材料的力学性能受到影响。但模拟结果中弯曲强度的变化趋势与实验结果仍然吻合,所以有限元分析方法对GO-CF/BF混杂增强复合材料力学性能的预测具有参考价值。结论基于VIHPS制备了GO-CF/BF混杂增强复合材料,并通过COMSOL与ANSYS有限元仿真软 件,分别建立了GO-CF/BF混杂增强复合材料浸渗模型与三点弯曲模型,探究了不同铺层顺序对该复合材料微观组织与弯曲性能的影响,主要结论如下: (1)通过COMSOL软件对不同铺层顺序的GO-CF/BF混杂增强复合材料进行浸渗效果仿真,并将仿真 结果与其微观形貌对比。结果表明,相比其余3种铺层顺序,CF与BF交替排列(BCBCBC)时,纤维间 的孔隙分布更加均匀,使得纤维被基体包裹,相邻纤维间不存在缝隙等缺陷,浸渗与界面结合效果 好。因此铺层顺序对复合材料微观组织的浸渗效果有显著影响。 (2)通过DIGIMAT软件建立不同铺层顺序的RVE模型,再将RVE模型导入ANSYS软件构建GO-CF/BF 混杂增强复合材料有限元模型,并进行三点弯曲仿真。结果表明,BBBCCC的铺层顺序复合材料弯 曲强度最大,为442.89 MPa。其原因是在此铺层下,复合材料浸渗效果较好,相同纤维间的结合使复 合材料内部更加紧密,减少层间界面处的缺陷和应力集中,且外侧受压的玄武岩纤维具有较好的韧 性,可以更好的分散应力。CCCBBB、CBBBCC和BCBCBC铺层顺序的复合材料与BBBCCC相比弯曲 强度分别下降了26.58%、33.13%和24.13%。因此GO-CF/BF复合材料的铺层顺序可以显著的影响其宏 观力学性能。 (3)将三点弯曲仿真结果与试验结果进行对比发现,当铺层顺序为BBBCCC时,此时复合材料弯曲强 度的仿真结果与试验结果均达到最高值,分别为442.89 MPa与410.59 MPa,试验与仿真结果误差小于 10%,两者较为符合。因此有限元分析方法在预测GO-CF/BF复合材料的力学性能上有一定参考价值。网络整理,禁止商用,仅限内部分享公众 号:机电君来源:机电君

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