CN115008263B - 一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法 - Google Patents
一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法 Download PDFInfo
- Publication number
- CN115008263B CN115008263B CN202210739186.7A CN202210739186A CN115008263B CN 115008263 B CN115008263 B CN 115008263B CN 202210739186 A CN202210739186 A CN 202210739186A CN 115008263 B CN115008263 B CN 115008263B
- Authority
- CN
- China
- Prior art keywords
- grinding
- force
- vibration
- single abrasive
- friction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002245 particle Substances 0.000 title claims abstract description 108
- 238000000034 method Methods 0.000 title claims abstract description 57
- 239000006061 abrasive grain Substances 0.000 claims abstract description 40
- 238000013178 mathematical model Methods 0.000 claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 32
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000003786 synthesis reaction Methods 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 230000000704 physical effect Effects 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000002604 ultrasonography Methods 0.000 claims 1
- 239000000463 material Substances 0.000 abstract description 10
- 238000003754 machining Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 239000002131 composite material Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005299 abrasion Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 229910001651 emery Inorganic materials 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000005304 optical glass Substances 0.000 description 1
- 238000006748 scratching Methods 0.000 description 1
- 230000002393 scratching effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B24—GRINDING; POLISHING
- B24B—MACHINES, DEVICES, OR PROCESSES FOR GRINDING OR POLISHING; DRESSING OR CONDITIONING OF ABRADING SURFACES; FEEDING OF GRINDING, POLISHING, OR LAPPING AGENTS
- B24B1/00—Processes of grinding or polishing; Use of auxiliary equipment in connection with such processes
- B24B1/04—Processes of grinding or polishing; Use of auxiliary equipment in connection with such processes subjecting the grinding or polishing tools, the abrading or polishing medium or work to vibration, e.g. grinding with ultrasonic frequency
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Abstract
本发明涉及一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，包括：获取单颗磨粒在超声椭圆振动辅助条件下的运动学特征；建立各磨削分力的数学模型，磨削分力包括：切屑变形力、磨粒切刃与工件滑擦的第一摩擦力以及切屑流出工件表面的第二摩擦力；每种磨削分力包括法向力、切向力和轴向力；基于各个磨削分力的数学模型，建立单颗磨粒磨削力初始模型，初始模型适用于无超声振动辅助的普通磨削、超声振动辅助磨削及超声椭圆振动磨削三种加工条件；基于磨削试验获取各加工条件下初始模型的模型参数，获得各加工条件下的单颗磨粒磨削力预测模型；采用预测模型预测磨削力。本发明可高效准确地对单颗磨粒磨削力进行预测，为材料磨削加工提供参考。
Description
技术领域
本发明涉及超声振动辅助磨削加工领域，尤其涉及一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法。
背景技术
以先进陶瓷、光学玻璃为代表的硬脆材料因具有高强硬度、低密度、耐高温、耐腐蚀、比刚度大等优良性能得到了广泛应用。一直以来，磨削是对硬脆材料进行机械加工的主要方式，但由于硬脆材料具有高硬脆性、低韧性的特点，属于难加工材料，所以常规的磨削方式在加工硬脆材料时往往会出现表面裂纹、亚表面损伤等现象，且砂轮磨损较严重、加工效率难以提高。超声振动辅助磨削在硬脆材料精密加工方面具有显著优势，如降低磨削力和磨削温度、提高材料去除率、减少砂轮磨损、提高加工表面完整性等。
磨削力是磨削加工中的一个重要参数，是加工过程的集中体现，也与磨粒磨损、工件加工质量存在密切的规律关系。研究磨削力对揭示加工机理并了解加工机理与加工结果的规律具有重要作用。从单颗磨粒角度研究磨削加工机理是一种简化而有效的手段，尽管与实际的磨削加工工艺有很大差别，但复杂的磨削过程是由若干个单颗磨粒在磨削区域的集成结果，且单颗磨粒不会受到其他磨粒的影响，很大程度上降低了研究难度。因此，利用单颗磨粒运动或加工机理建立磨削力预测模型是探究总磨削力的基础，对了解磨削过程具有很大帮助。尽管目前研究从不同角度构建了单颗磨粒磨削力预测模型，但并未形成完善体系，且预测的准确度和普适性都不高。
发明内容
本发明提供一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，以解决上述技术问题。
为解决上述技术问题，本发明提供一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，包括如下步骤：
S100：基于超声振动原理，获取单颗磨粒在超声椭圆振动辅助条件下的运动学特征；
S200：基于磨削原理与超声椭圆振动对磨削力的影响规律，分别建立各磨削分力的数学模型，所述磨削分力至少包括：切屑变形力、磨粒切刃与工件滑擦的第一摩擦力以及切屑流出工件表面的第二摩擦力；每种所述磨削分力包括法向力、切向力和轴向力；
S300：基于各个所述磨削分力的数学模型，建立单颗磨粒磨削力初始模型，所述初始模型适用于无超声振动辅助的普通磨削、超声振动辅助磨削以及超声椭圆振动磨削三种加工条件；
S400：基于磨削试验获取各加工条件下所述初始模型的模型参数，获得各加工条件下的单颗磨粒磨削力预测模型；
S500：采用所述预测模型预测磨削力。
较佳地，步骤S100包括：
S110：根据砂轮与工件的位置建立三维坐标系M-XYZ，获取单颗磨粒在坐标系M-XYZ中的运动学方程，并基于所述运动学方程计算单颗磨粒在各坐标方向上的运动分速度vx、vy、vz；
S120：分别获取椭圆振动信号和组成所述椭圆振动信号的两个方向垂直的原始振动信号的表达方程，并求取各振动信号的振动速度vu、vy1、vz1；
S130：获取单颗磨粒在磨削弧区的合成速度ve；
S140：分别计算椭圆振动信号的振动速度vu和合成速度ve与磨削方向的夹角β、α，其中，β、α随时间t周期性变化。
较佳地，所述椭圆振动信号的振动速度vu和合成速度ve与磨削方向的夹角β、α的计算公式分别为：
其中，vs为砂轮线速度；t为运动时间；vw为工件进给速度；Ay、Az分别为y、z方向的超声振幅，且Az/Ay＝λ；w为超声振动角频率，w＝2πf，f为超声振动频率；分别为y、z方向振动的初始相位角，/>
较佳地，步骤S200中，建立所述切屑变形力的数学模型的方法包括：
S211：基于实际工作磨粒的几何模型，建立超声椭圆振动辅助条件下单颗磨粒的受力模型，进行力学分析；
S212：获取所述磨粒前刀面上的微小单元面积的切屑变形力dFgc、磨粒与工件接触的微小单元面积dA；
S213：对微小单元面积的所述切屑变形力dFgc作积分运算，获取切屑变形力Fgc；
S214：基于步骤S100获取的运动学特征以及所述切屑变形力的法向力、切向力和轴向力的关系，获得所述切屑变形力的数学模型。
较佳地，步骤S200中，建立所述第一摩擦力的数学模型的方法包括：
S221：借鉴于通过磨削原理获取普通磨削下单颗磨粒由摩擦引起的法向力、切向力，通过摩擦二项式定理获取普通磨削的摩擦系数，可得具体加工条件下的摩擦系数μ：
其中，ds为砂轮直径，为磨粒实际磨损平面与工件间的平均接触压强，p0为系数，能够通过试验得到，m、n为系数，取决于接触面的物理性能；
S222：计算振动条件和无振动条件的摩擦力之间的比值r，
其中，ζ为物体宏观速度与振动速度的幅值之比；β为物体宏观速度与振动速度的夹角；τ为规范化时间，τ＝ωt；
S223：基于所述摩擦系数μ和比值r，获得所述第一摩擦力的数学模型。
较佳地，步骤S200中，建立所述第二摩擦力的数学模型的方法包括：
S231：基于所述第二摩擦力与所述切屑变形力存在摩擦系数的关系，求得所述第二摩擦力的数学模型。
较佳地，所述单颗磨粒的形状为截角八面体，所述单颗磨粒磨削力的初始模型是基于磨粒的几何形貌建立的。
较佳地，所述单颗磨粒磨削力初始模型为：
其中，γ为磨粒侧面角；a为截角八面体的边长；θ为磨粒的负前角；ap为磨削深度；为系数。
较佳地，当r＝1、α＝0时，所述单颗磨粒磨削力初始模型为普通磨削加工条件下的单颗磨粒磨削力初始模型。
较佳地，步骤S400中，基于拟合法、代入法或遗传算法计算所述模型参数m、n、P0、Fp、
与现有技术相比，本发明提供的超声振动辅助磨削下基于单颗磨粒的磨削力预测方法具有如下优点：
1、本发明基于单颗磨粒运动学规律建立磨削力预测模型，使单颗磨粒磨削力预测模型不局限于具体加工材料，更具有普适性；
2、本发明将磨削力分为切屑变形力、磨粒切刃与工件滑擦的第一摩擦力以及切屑流出工件表面的第二摩擦力，分别分析超声振动对各磨削分力的影响，并在预测模型中有具体参数体现，使预测模型考虑更全面，也更符合实际加工过程；
3、该预测模型适用于普通磨削、超声振动辅助磨削、超声椭圆振动磨削多个加工背景；
4、本发明基于实际磨粒的几何形貌特征进行受力分析，在很大程度上使磨削力的预测结果更符合实际。
附图说明
图1为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法的流程示意图；
图2为单颗磨粒超声椭圆振动磨削示意图；
图3为超声振动原理图；
图4为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法中椭圆振动信号合成示意图；
图5为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法中β求解分析示意图；
图6为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法中α求解分析示意图；
图7为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法中磨粒的几何形貌示意图；
图8为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法中单颗磨粒受力分析图；
图9为图8的Y-Y剖视图；
图10为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法中切屑变形力受力分解图；
图11a和11b分别为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法中不同超声振幅下cosα的变化规律图；
图12a和12b分别为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法中不同超声振幅下r的变化规律图；
图13a和13b分别为本发明一具体实施方式中超声振动辅助磨削下基于单颗磨粒的磨削力预测方法中通过拟合法计算模型参数值的拟合结果仿真示意图。
图中：10-砂轮、20-工件、30-磨粒。
具体实施方式
为了更详尽的表述上述发明的技术方案，以下列举出具体的实施例来证明技术效果；需要强调的是，这些实施例用于说明本发明而不限于限制本发明的范围。
本发明提供的超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，如图1所示，包括如下步骤：
S100：基于超声振动原理，获取单颗磨粒在超声椭圆振动辅助条件下的运动学特征；
S200：基于磨削原理与超声椭圆振动对磨削力的影响规律，分别建立各磨削分力的数学模型，所述磨削分力至少包括：切屑变形力、磨粒切刃与工件滑擦的第一摩擦力以及切屑流出工件表面的第二摩擦力；每种所述磨削分力包括法向力、切向力和轴向力；
S300：基于各个所述磨削分力的数学模型，建立单颗磨粒磨削力初始模型，所述初始模型适用于无超声振动辅助的普通磨削、超声振动辅助磨削以及超声椭圆振动磨削三种加工条件；
S400：基于磨削试验获取各加工条件下所述初始模型的模型参数，获得各加工条件下的单颗磨粒磨削力预测模型；
S500：采用所述预测模型预测磨削力。
本发明基于单颗磨粒运动学规律建立磨削力预测模型，使单颗磨粒磨削力预测模型不局限于具体加工材料，更具有普适性；将磨削力分为切屑变形力、磨粒切刃与工件滑擦的第一摩擦力以及切屑流出工件表面的第二摩擦力，分别分析超声振动对各磨削分力的影响，并在预测模型中有具体参数体现，使预测模型考虑更全面，也更符合实际加工过程；该预测模型适用于普通磨削、超声振动辅助磨削、超声椭圆振动磨削多个加工背景。
请重点参考图2，在一些实施例中，步骤S100包括：
S110：根据砂轮10与工件20的位置建立三维坐标系M-XYZ，获取单颗磨粒在坐标系M-XYZ中的运动学方程，并基于所述运动学方程计算单颗磨粒在各坐标方向上的运动分速度vx、vy、vz；
S120：椭圆振动信号是一种合成振动信号，它由两个方向垂直的振动信号形成。分别获取椭圆振动信号和组成所述椭圆振动信号的两个方向垂直的原始振动信号的表达方程，并求取各振动信号的振动速度vu、vy1、vz1；
S130：获取单颗磨粒在磨削弧区的合成速度ve；
S140：分别计算椭圆振动信号的振动速度vu和合成速度ve与磨削方向的夹角β、α，其中，β、α随时间t周期性变化。
在一些实施例中，如图2和图3所示，单颗磨粒在坐标系M-XYZ中的运动学方程为：
式中，ds为砂轮直径；vs为砂轮线速度；t为运动时间；vw为工件进给速度；Ay、Az分别为y、z方向的超声振幅(Az/Ay＝λ)；w为超声振动角频率(w＝2πf，f为超声振动频率)；分别为y、z方向振动的初始相位角/>
单颗磨粒在各坐标方向上的运动分速度vx、vy、vz为：
如图4所示，椭圆振动信号(合成振动信号)由两个方向垂直的振动信号(原始振动信号)形成，合成振动信号、原始振动信号的振动速度(vu、v1、v2)为：
当φz-φy＝π/2时，Φ取负号；当φz-φy＝-π/2时，Φ取正号。在磨削过程中，由于磨削深度相比砂轮10周长非常之小，以致磨削角θg很小(cosθg≈1)，所以可近似地将砂轮10在微观磨削弧区内的运动看作直线运动。单颗磨粒沿砂轮运动方向(Y方向)的运动速度简化为：
vy(t)＝vw+vs+vy1(t) (5)
进一步，单颗磨粒在磨削弧区的合成速度为：
请重点参考图5和图6，在一些实施例中，椭圆振动信号的振动速度vu和合成速度ve与磨削方向(磨粒运动方向)的夹角β、α的计算公式分别为：
在一些实施例中，步骤S200中，建立所述切屑变形力的数学模型的方法包括：
S211：基于实际工作磨粒的几何模型，建立超声椭圆振动辅助条件下单颗磨粒的受力模型，进行力学分析；
S212：获取所述磨粒前刀面上的微小单元面积的切屑变形力dFgc、磨粒与工件接触的微小单元面积dA；
S213：对微小单元面积的所述切屑变形力dFgc作积分运算，获取切屑变形力Fgc；
S214：基于步骤S100获取的运动学特征以及所述切屑变形力的法向力、切向力和轴向力的关系，获得所述切屑变形力的数学模型。
在一些实施例中，磨粒30的总体几何形状近似于截角八面体，如图7所示。图8为单颗磨粒的受力分解图，在不考虑摩擦作用的情况下，单颗磨粒的切屑变形力将垂直作用于磨粒30的前刀面上。在Y-Y截面内，磨粒前刀面上的微小单元面积的切屑变形力dFgc为：
dFgc＝FpcosθdA (9)
式中，Fp为单位磨削力；θ为磨粒的负前角；dA为磨粒30与工件20接触的微小单元面积；
结合磨粒30的形状特征，磨粒与工件接触的微小单元面积dA可表示为：
式中，dA1为前刀面正面的微小单元面积；dA2为前刀面侧面的微小单元面积；dL为前刀面正面的微小单元长度；a为截角八面体的边长。
进一步，对微小单元面积的切屑变形力dFgc作积分运算，获取切屑变形力Fgc数值大小，求解得：
根据对超声椭圆振动下单颗磨粒的运动学分析，切屑变形力的法向分力、切向分力和轴向分力的关系如图9所示，单颗磨粒的切屑变形力为：
在一些实施例中，步骤S200中，建立所述第一摩擦力的数学模型的方法包括：
S221：通过磨削原理获取普通磨削下单颗磨粒由摩擦引起的法向力、切向力，在一些实施例中，单颗磨粒由摩擦引起的法向力、切向力为：
式中，δ为磨粒与工件的实际接触面积；p为磨粒实际磨损平面与工件间的平均接触压强，μ为磨粒与工件之间的摩擦系数。
借鉴于上述原理，可以通过摩擦二项式定理，得到具体加工条件下的摩擦系数μ为：
其中，ds为砂轮直径，为磨粒实际磨损平面与工件间的平均接触压强，p0为系数，能够通过试验得到，m、n为系数，取决于接触面的物理性能。
S222：研究表明，振动具有减小摩擦的效果，在超声振动下，单颗磨粒的摩擦力得到降低，有无振动的摩擦力存在比值关系，比值r的计算公式如下：
其中，ζ为物体宏观速度与振动速度的幅值之比；β为物体宏观速度与振动速度的夹角；τ为规范化时间，τ＝ωt；
近似地将砂轮10在微观磨削弧区内的运动看作直线运动，得ζ的计算公式如下：
β随t周期性变化，可通过式(7)求得；
S223：基于所述摩擦系数μ和比值r，获得所述第一摩擦力的数学模型为：
在一些实施例中，步骤S200中，建立所述第二摩擦力的数学模型的方法包括：
S231：基于所述第二摩擦力与所述切屑变形力存在摩擦系数的关系，求得所述第二摩擦力的数学模型。
在一些实施例中，所述单颗磨粒磨削力的初始模型是基于磨粒的几何形貌建立的。
在一些实施例中，切屑在流出过程中与磨粒前刀面发生摩擦，该处的摩擦力为：
dFgf＝μdFgc (18)
在超声椭圆振动下，切屑流出的摩擦力也可分解为法向力、切向力和轴向力，因此，可将所述第二摩擦力的数学模型表示为：
在一些实施例中，基于建立的切屑变形力的数学模型、第一摩擦力的数学模型以及第二摩擦力的数学模型，获取单颗磨粒磨削力初始模型的计算公为：
其中，γ为磨粒侧面角，如图8所示，γ决定了磨粒两侧的磨粒-工件实际接触面积，当磨粒为如图7所示的截角八面体时，γ＝30°；a为截角八面体的边长；θ为磨粒的负前角；ap为磨削深度；为系数。
特别地，当r＝1、α＝0时，所述单颗磨粒磨削力初始模型为普通磨削加工条件下的单颗磨粒磨削力初始模型。
在一些实施例中，通过如图2所示的超声椭圆振动磨削方式，开展SiC陶瓷的单颗磨粒划擦试验，利用Kistler三向测力仪获取实际磨削力，本实施例采用单因素试验，自变量参数为磨削速度vs、单颗磨粒切厚agmax，试验磨削工艺参数如表1所示，其中，vw的试验参数值通过公式(21)求解，具体数值如
表2所示，
表1工艺参数
表2 vw的试验参数值
在一些实施例中，获取磨削加工条件下cosα与摩擦力比r，具体计算公式如下：
在超声椭圆振动下，合成速度ve随振动时间呈现周期性变化，以致角度α呈周期性变化，因此可通过式(22)求解一个周期内cosα的平均值同样，施加超声椭圆振动与无振动的摩擦力比值r也随时间周期性变化，可通过式(23)求解一个周期内的平均值。
如图11a和11b所示，在加工条件下，cosα随单颗磨粒切厚几乎没有变化，因此可使用各单颗磨粒切厚(0.05～0.4μm)下cosα均值表示某磨削速度下的cosα值。Ay＝0μm表示普通磨削下的cosα值(α＝π/2)。当Ay＝0.3μm时，cosα值基本与普通磨削下的情况相同，可近似认为cosα＝1，这表明此时轴向磨削力很小，可以忽略。但Ay较大时，cosα随磨削速度的增大而增加，逐渐趋近于1。相同磨削速度下，cosα随Ay的增大而减小。
如图12a和12b所示，摩擦力比值r与cosα具有相似的变化规律。但当振幅较小，如Ay＝0.3μm时，超声椭圆振动磨削与普通磨削下的r存在较明显的差别，这表明此时超声椭圆振动对摩擦力依然有较显著的影响。
在一些实施例中，通过开展单颗磨粒超声椭圆振动磨削单因素试验，通过测力仪测量各加工条件下的法向力、切向力与轴向力，代入单颗磨粒磨削力预测模型的计算公式，分别得到各加工条件对应的模型参数下m、n、P0、Fp、的值，进而计算得到多组模型参数值，计算各参数平均值；基于不同加工参数下cosα、r的计算值以及模型参数平均值，获取单颗磨粒磨削力计算公式。
为获取磨削加工条件下m、n、P0、Fp、的值，所述的计算模型参数值的方法可以为拟合法、代入法或遗传算法，图13a和13b即为通过拟合法获得的拟合结果仿真示意图。
表3、4分别为普通磨削、超声椭圆振动磨削下单颗磨粒磨削力的预测结果。
表3普通磨削、vs＝8m/s时的切向力、法向力
表4超声椭圆振动磨削、不同磨削速度下的切向力、法向力
/>
综上所述，本发明提供的超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，包括如下步骤：S100：基于超声振动原理，获取单颗磨粒在超声椭圆振动辅助条件下的运动学特征；S200：基于磨削原理与超声椭圆振动对磨削力的影响规律，分别建立各磨削分力的数学模型，所述磨削分力至少包括：切屑变形力、磨粒切刃与工件滑擦的第一摩擦力以及切屑流出工件表面的第二摩擦力；每种所述磨削分力包括法向力、切向力和轴向力；S300：基于各个所述磨削分力的数学模型，建立单颗磨粒磨削力初始模型，所述初始模型适用于无超声振动辅助的普通磨削、超声振动辅助磨削以及超声椭圆振动磨削三种加工条件；S400：基于磨削试验获取各加工条件下所述初始模型的模型参数，获得各加工条件下的单颗磨粒磨削力预测模型；S500：采用所述预测模型预测磨削力。本发明通过分析单颗磨粒在超声椭圆振动辅助条件下的运动学特征，再将磨削力分为切屑变形力、磨刃与工件滑擦的第一摩擦力、切屑流出工件表面的第二摩擦力，并基于磨粒实际几何形貌、通过力学分析分别探究超声椭圆振动对各组成力的影响规律，与普通磨削不同的变化参数包括合成速度与运动方向夹角α、超声-普通摩擦力比r，进而构建超声椭圆振动下单颗磨粒磨削力模型；通过开展单颗磨粒普通磨削和超声椭圆振动磨削SiC试验，确定两种磨削方式下的预测模型，对单颗磨粒磨削力进行了有效预测，更符合实际加工过程。
显然，本领域的技术人员可以对发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包括这些改动和变型在内。
Claims (5)
1.一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，其特征在于，包括如下步骤：
S100：基于超声振动原理，获取单颗磨粒在超声椭圆振动辅助条件下的运动学特征，包括：
S110：根据砂轮与工件的位置建立三维坐标系M-XYZ，获取单颗磨粒在坐标系M-XYZ中的运动学方程，并基于所述运动学方程计算单颗磨粒在各坐标方向上的运动分速度v x 、v y 、v z ；
S120：分别获取椭圆振动信号和组成所述椭圆振动信号的两个方向垂直的原始振动信号的表达方程，并求取各振动信号的振动速度v u 、v y1 、v z1 ；
S130：获取单颗磨粒在磨削弧区的合成速度v e ；
S140：分别计算椭圆振动信号的振动速度v u 和合成速度v e 与磨削方向的夹角β、α，其中，β、α随时间t周期性变化，β、α的计算公式分别为：
，
，
其中，为砂轮线速度；/>为运动时间；/>为工件进给速度；/>、/>分别为/>、/>方向的超声振幅，且/>；/>为超声振动角频率，/>，/>为超声振动频率；/>、/>分别为/>、方向振动的初始相位角，/>；
S200：基于磨削原理与超声椭圆振动对磨削力的影响规律，分别建立各磨削分力的数学模型，所述磨削分力至少包括：切屑变形力、磨粒切刃与工件滑擦的第一摩擦力以及切屑流出工件表面的第二摩擦力；每种所述磨削分力包括法向力、切向力和轴向力；
所述步骤S200中，建立所述第一摩擦力的数学模型的方法包括：
S221：借鉴于通过磨削原理获取普通磨削下单颗磨粒由摩擦引起的法向力、切向力，通过摩擦二项式定理获取普通磨削的摩擦系数，可得具体加工条件下摩擦系数μ：
，
其中，为砂轮直径，/>为磨粒实际磨损平面与工件间的平均接触压强，p 0 为系数，能够通过试验得到，m、n为系数，取决于接触面的物理性能；
S222：计算振动条件和无振动条件的摩擦力之间的比值r，
，
其中，ζ为物体宏观速度与振动速度的幅值之比；β为物体宏观速度与振动速度的夹角；τ为规范化时间，τ=ωt；
S223：基于所述摩擦系数μ和比值r，获得所述第一摩擦力的数学模型；
S300：所述单颗磨粒的形状为截角八面体，基于磨粒的几何形貌和各个所述磨削分力的数学模型，建立单颗磨粒磨削力初始模型如下：
，
其中，为磨粒侧面角；/>为截角八面体的边长；/>为磨粒的负前角；a p 为磨削深度；/>为系数；
所述初始模型适用于无超声振动辅助的普通磨削、超声振动辅助磨削以及超声椭圆振动磨削三种加工条件；
S400：基于磨削试验获取各加工条件下所述初始模型的模型参数，获得各加工条件下的单颗磨粒磨削力预测模型；
S500：采用所述预测模型预测磨削力。
2.权利要求1所述的超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，其特征在于，步骤S200中，建立所述切屑变形力的数学模型的方法包括：
S211：基于实际工作磨粒的几何模型，建立超声椭圆振动辅助条件下单颗磨粒的受力模型，进行力学分析；
S212：获取所述磨粒前刀面上的微小单元面积的切屑变形力dF gc 、磨粒与工件接触的微小单元面积d A ；
S213：对微小单元面积的所述切屑变形力dF gc 作积分运算，获取切屑变形力F gc ；
S214：基于步骤S100获取的运动学特征以及所述切屑变形力的法向力、切向力和轴向力的关系，获得所述切屑变形力的数学模型。
3.如权利要求1所述的超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，其特征在于，步骤S200中，建立所述第二摩擦力的数学模型的方法包括：
S231：基于所述第二摩擦力与所述切屑变形力存在摩擦系数的关系，求得所述第二摩擦力的数学模型。
4.如权利要求1所述的超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，其特征在于，当r=1、α=0时，所述单颗磨粒磨削力初始模型为普通磨削加工条件下的单颗磨粒磨削力初始模型。
5.如权利要求1所述的超声振动辅助磨削下基于单颗磨粒的磨削力预测方法，其特征在于，步骤S400中，基于拟合法、代入法或遗传算法计算模型参数m、n、P 0 、F p 、。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210739186.7A CN115008263B (zh) | 2022-06-14 | 2022-06-14 | 一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210739186.7A CN115008263B (zh) | 2022-06-14 | 2022-06-14 | 一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115008263A CN115008263A (zh) | 2022-09-06 |
| CN115008263B true CN115008263B (zh) | 2023-12-22 |
Family
ID=83076289
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210739186.7A Active CN115008263B (zh) | 2022-06-14 | 2022-06-14 | 一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115008263B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6239174A (ja) * | 1985-07-29 | 1987-02-20 | エナジ− − アダプテイブ グラインデイング インコ−ポレイテイツド | 研削方法及びその実施システム |
| CN207923641U (zh) * | 2017-12-27 | 2018-09-28 | 青岛理工大学 | 支持不同润滑工况的单颗磨粒速度及尺寸效应试验系统 |
| CN108723900A (zh) * | 2018-04-25 | 2018-11-02 | 北京航空航天大学 | 一种旋转超声振动铣磨的全域切削力建模方法与装置 |
| CN108897955A (zh) * | 2018-07-01 | 2018-11-27 | 北京工业大学 | 一种考虑螺纹升角和磨粒重叠系数的高速内螺纹磨削力计算方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019109920A1 (zh) * | 2017-12-06 | 2019-06-13 | 青岛理工大学 | 超声振动辅助磨削液微通道浸润的纳米流体微量润滑磨削装置 |
-
2022
- 2022-06-14 CN CN202210739186.7A patent/CN115008263B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6239174A (ja) * | 1985-07-29 | 1987-02-20 | エナジ− − アダプテイブ グラインデイング インコ−ポレイテイツド | 研削方法及びその実施システム |
| CN207923641U (zh) * | 2017-12-27 | 2018-09-28 | 青岛理工大学 | 支持不同润滑工况的单颗磨粒速度及尺寸效应试验系统 |
| CN108723900A (zh) * | 2018-04-25 | 2018-11-02 | 北京航空航天大学 | 一种旋转超声振动铣磨的全域切削力建模方法与装置 |
| CN108897955A (zh) * | 2018-07-01 | 2018-11-27 | 北京工业大学 | 一种考虑螺纹升角和磨粒重叠系数的高速内螺纹磨削力计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 轴向超声振动辅助磨削的磨削力建模;唐进元;周伟华;黄于林;;机械工程学报;52(第15期);第184-189页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115008263A (zh) | 2022-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shen et al. | Ultrasonic vibration-assisted milling of aluminum alloy | |
| Pei et al. | Modeling of ductile-mode material removal in rotary ultrasonic machining | |
| Cao et al. | Effect of intermittent cutting behavior on the ultrasonic vibration-assisted grinding performance of Inconel718 nickel-based superalloy | |
| Qingyu et al. | Modelling of grinding mechanics: A review | |
| Gopal et al. | A new chip-thickness model for performance assessment of silicon carbide grinding | |
| Zhou et al. | Towards understanding the ploughing friction mechanism in ultrasonic assisted grinding with single grain | |
| Salvatore et al. | Experimental and numerical study of media action during tribofinishing in the case of SLM titanium parts | |
| Li et al. | Experimental study on vibration-assisted grinding | |
| Liu et al. | Experimental study of tool wear and its effects on cutting process of ultrasonic-assisted milling of Ti6Al4V | |
| Walia et al. | Determining dynamically active abrasive particles in the media used in centrifugal force assisted abrasive flow machining process | |
| Mali et al. | Current status and application of abrasive flow finishing processes: a review | |
| Haidong et al. | A study on ultrasonic elliptical vibration cutting of Inconel 718 | |
| CN115008263B (zh) | 一种超声振动辅助磨削下基于单颗磨粒的磨削力预测方法 | |
| Li et al. | Improving the working surface condition of electroplated cubic boron nitride grinding quill in surface grinding of Inconel 718 by the assistance of ultrasonic vibration | |
| Lu et al. | Friction modeling of tool-chip interface based on shear-slip theory for vibration assisted swing cutting | |
| Liu et al. | Experimental and numerical study on surface generated mechanism of robotic belt grinding process considering the dynamic deformation of elastic contact wheel | |
| Misra et al. | 4 Magnetic abrasive finishing process | |
| Ghule et al. | Ultrasonic Vibrations Assisted Turning (UAT): A Review | |
| Jiao et al. | Research on ultrasonic-assisted fixed-abrasive lapping technology for engineering ceramics cylindrical part | |
| Hassan et al. | A review of the milling process to fabricate a dimple structure | |
| Zhan et al. | Study on the specific grinding energy of cemented carbide (YG8) grinding with a vitrified diamond wheel in high speed regime | |
| Ma et al. | Investigation of intermittent cutting mechanism of abrasive particles in ultrasonic-assisted external cylindrical honing | |
| Lam et al. | Experimental Study on Cutting Force in Shear Thickening Polishing for Spherical Steel | |
| Rai et al. | Role of sacrificing layer on glasses during ultrasonic machining (USM)-a review | |
| Senthilkumar et al. | State of the art of micro turning process |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |