CN115438533A - 基于反节点振动位移的超声骨切削力的间接测量方法 - Google Patents

基于反节点振动位移的超声骨切削力的间接测量方法 Download PDF

Info

Publication number
CN115438533A
CN115438533A CN202210936098.6A CN202210936098A CN115438533A CN 115438533 A CN115438533 A CN 115438533A CN 202210936098 A CN202210936098 A CN 202210936098A CN 115438533 A CN115438533 A CN 115438533A
Authority
CN
China
Prior art keywords
model
cutting force
vibration displacement
ultrasonic
nodal
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.)
Pending
Application number
CN202210936098.6A
Other languages
English (en)
Inventor
王林伟
刘宇
邹启磊
舒彬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN202210936098.6A priority Critical patent/CN115438533A/zh
Publication of CN115438533A publication Critical patent/CN115438533A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

基于反节点振动位移的超声骨切削力的间接测量方法，包括：确定超声骨刀节点位置，得到其简化模型；建立压电换能器模型，计算模型参数，获得其振动位移；建立阶梯变幅杆模型，根据特征值、特征向量和能量守恒原则计算动力学参数；建立骨组织的力负载模型，辨识模型参数；基于上述模型，搭建Matlab‑Simulink仿真程序；在此基础上，利用数据驱动建模方法，以反节点振动位移x1的有效值x1rms、振动频率fx1为输入，输出切削力的有效值Frms，获得反节点振动位移与切削力的关系；利用传感器测量反节点振动位移的有效值x1rms及振动频率fx1，根据上述步骤获得的关系计算切削力Frms。本发明提供了超声骨刀‑骨组织相互作用的一维数值计算模型，基于此提供超声骨切削力的间接测量方法。

Description

基于反节点振动位移的超声骨切削力的间接测量方法
发明名称
基于反节点振动位移的超声骨切削力的间接测量方法。
技术领域
本发明属于机器人超声骨刀加工应用领域，具体涉及基于反节点振动位移的超声骨切削力的间接测量方法。
背景技术
超声骨刀作为一种新型手术器械，高频电信号通过磁致伸缩换能器或压电陶瓷换能器转化为高频机械振动，激起刀具的高频纵向振动模态，实现上百微米级振幅的超声振动进行切削。相对于传统的骨切削工具，超声骨刀具有低切削力、低切削温度和组织选择性等优点，近年来一直作为手持工具被外科医生广泛应用于脊柱外科、颌面外科、整形外科等领域，然而超声骨刀在骨科机器人中的应用尚未见报道。超声骨刀对末端切削力敏感，切削力的大小会影响切削稳定性(改变骨刀谐振频率、振动幅值)、切削质量，进而影响手术效果。因此，超声骨刀切削力的检测模型，是超声骨刀应用在机器人辅助手术中的一个亟待解决的问题，对在手术过程中实现骨质层感知、切削力控制具有重要意义。
发明内容
针对现有技术存在的问题，本发明提供基于反节点振动位移的超声骨切削力的间接测量方法。
本发明的技术方案是：
基于反节点振动位移的超声骨切削力的间接测量方法，包括以下步骤：
步骤1：确定超声骨刀节点位置，得到超声骨刀的简化模型。
步骤2：建立压电换能器模型，计算模型参数，获得压电换能器的振动位移。
步骤3：建立阶梯变幅杆模型，根据特征值、特征向量和能量守恒原则计算动力学参数。
步骤4：建立骨组织的力负载模型，辨识模型参数。
步骤5：建立超声骨刀-骨组织相互作用模型，搭建Matlab-Simulink仿真程序。
步骤6：基于Matlab-Simulink仿真程序，利用数据驱动建模方法，以反节点振动位移x1的有效值x1rms、振动频率fx1为输入，输出切削力的有效值Frms，获得反节点振动位移x1与切削力F的关系。
步骤7：利用电涡流位移传感器测量反节点振动位移的有效值x1rms及振动频率fx1，根据步骤6获得的关系计算切削力的有效值Frms。
有益效果
本发明将超声骨刀这一个复杂的连续系统(具有无限多个自由度)简化为一个2自由度的离散系统，该模型具有原始系统的属性，同时可以用于仿真。开发的一维数值计算模型处理所需信息的速度比有限元、实验快得多，易于建立数据集。超声骨刀切削力的间接检测方法将促进超声骨刀在骨科机器人中的应用，为骨科机器人在手术过程中实现骨质层感知、切削力控制提供理论基础。
附图说明
图1是本发明具体实施方式的超声骨刀简化模型图；
图2是本发明具体实施方式的特征值特征向量获得示意图；
图3是本发明具体实施方式的骨组织的力负载模型图；
图4是本发明具体实施方式的Matlab-Simulink仿真模型示意图；
图5是本发明具体实施方式的反节点振动位移测量示意图；
图6是本发明具体实施方式的基于反节点振动位移的超声骨切削力的间接测量方法流程图。
具体实施方式
下面结合附图对本发明的具体实施方式做详细说明。
基于反节点振动位移的超声骨切削力的间接测量方法，如图6所示，包括以下步骤：
步骤1：确定超声骨刀节点位置，得到超声骨刀的简化模型。
超声骨刀的简化模型如图1，超声骨刀在纵向模态振动工作时，整个超声骨刀(不考虑刀头)有两个节点：一个位于两片压电陶瓷中间，另一个位于变幅杆台阶处；由于两片压电陶瓷中间存在节点，可防止超声波向后传播，因此忽略了节点左侧部分，超声骨刀的简化模型由两部分组成：压电换能器模型、阶梯变幅杆模型。
根据反节点个数将变幅杆简化为2自由度质量-弹簧-阻尼(MSD)系统，它具有两个谐振频率，第一谐振频率对应于m1，m2同相振动。第二谐振频率对应于m1，m2反相振动。
步骤2：建立压电换能器模型，计算模型参数，获得压电换能器的振动位移。
压电换能器的工作原理是通过在压电陶瓷两端施加交变电压激励，利用压电材料的逆压电特性，使其产生纵向振动的压缩或膨胀变形。在不考虑压电材料的非线性的假设下，压电陶瓷的电气和机械行为之间的相互作用可以用以下线性关系来描述：
Figure BDA0003783453420000031
式中：S是应变，T是应力，SE是恒定电场下的弹性柔度，d是压电电荷常数，D是介电位移，E是电场强度，εT是恒定应力下的介电常数。
Figure BDA0003783453420000032
式中：S0是压电陶瓷截面积，l0是压电陶瓷厚度，x0是压电陶瓷的振动位移，F0是变幅杆施加到压电换能器的力，u是激励电压。
将式(2)带入式(1)，可得压电换能器的振动位移为
Figure BDA0003783453420000033
式中：
Figure BDA0003783453420000034
为弹性柔度，d33为电荷常数，
Figure BDA0003783453420000035
为介电常数。
步骤3：建立阶梯变幅杆模型，根据特征值、特征向量和能量守恒原则计算动力学参数。
步骤3-1：动力学建模
根据Euler-Lagrange方程，无负载超声骨刀的动力学方程可以写成：
Figure BDA0003783453420000036
压电陶瓷与2-DOF模型之间的相互作用力F0可表示为：
Figure BDA0003783453420000037
式中：m1，m2，k1，k2，c1，c2分别为等效质量、等效刚度和等效阻尼。x0为压电陶瓷的振动位移，x1，x2分别为m1，m2的振动位移。
根据边界条件，超声骨刀的纵向振动位移可表示为：
Figure BDA0003783453420000038
式中：V1，V2为m1，m2的振动幅值，κ为圆波数，l0为变幅杆大圆柱长度，l1为变幅杆全长。
步骤3-2：动力学参数辨识
系统的质量矩阵、刚度矩阵、特征值、特征向量可以表示为：
Figure BDA0003783453420000039
式中：ω1为1阶固有频率，ω2为2阶固有频率，ψ1为与ω1对应的振幅比，ψ2为与ω2对应的振幅比。
由特征方程MλΨ＝KΨ，可得
Figure BDA0003783453420000041
然而，对于式(8)有四个未知参数(m1，m2，k1，k2)只有三个独立的方程，因此解是无限组合(无限个解)。根据能量守恒原则可以补充方程
Figure BDA0003783453420000042
联立方程式(8)、式(9)可得
Figure BDA0003783453420000043
由此，可以求解系统动力学参数(m1，m2，k1，k2)。
系统特征值测量：利用ANSYS Workbench 20.0软件中模态分析模块Modal对超声骨刀进行模态分析，以获得2个纵向振动模态对应的振动频率ω1、ω2。
系统特征向量测量：如图2特征值特征向量获得示意图，在超声骨刀模型前后两端创建两个参考点，连接参考点形成路径，通过ANSYSWorkbench 20.0软件中的MEMS ACT插件和谐响应分析模块Harmonic Response对超声骨刀进行谐响应分析，对超声骨刀分别施加ω1，ω2频率的电压，得到超声骨刀沿路径的轴向振动位移分布曲线，计算对应频率的振幅比ψ1、ψ2后即可求得特征向量。
通过调整系数c1，c2以调整模拟的超声波振动幅值x1，x2与实验测量值V1，V2相匹配。
步骤4：建立骨组织的力负载模型，辨识模型参数。
步骤4-1：力负载模型建立。
超声骨刀切削过程中，刀具与骨组织之间的相互作用会对振动系统产生力负载，可用Kelvin-Voigt模型表示如图3，负载可以等效成一个弹性系数为k的线性弹簧和一个阻尼系数为c的阻尼器的模型。负载模型的动态响应可以由下式表示：
Figure BDA0003783453420000051
式中：k为刚度系数；c为阻尼系数；x2为刀具的振动位移；Δ为初始干涉。
步骤4-2：力负载模型参数辨识
通过调整初始干涉、刚度系数和阻尼系数的值以匹配在超声骨切削有限元仿真中获得的切削力。
步骤5：建立超声骨刀-骨组织相互作用模型，搭建Matlab-Simulink仿真程序。
有负载超声骨刀的动力学方程可以写成：
Figure BDA0003783453420000052
式中：ψ2(ωF)为与ωF对应的振幅比。
由式(12)可以看出，切削力F变化将导致超声骨刀谐振频率、振型发生变化，进而使振动位移x1或x2的变化，切削力F和振动位移x1或x2存在某种对应关系。因此，可以通过测量振动位移x1或x2的变化来反映切削力F。然而，一方面狭窄的切削区，不易固定安装传感器；另一方面冷却水雾、血液等恶劣的工作环境。这两方面导致振动位移x2无法测量，因此本文提出一种基于反节点振动位移x1的超声骨切削力的间接测量方法。
根据式(3)、(5)、(11)建立MATLAB-Simulink仿真模型如图4。切削力将导致超声骨刀谐振频率发生漂移，为保证骨刀以谐振频率工作，需要进行频率追踪。通过给系统输入50KHz脉冲激励，测量系统输出，时间步长为1e-6，采样频率分辨率1Hz，根据频响函数定义，将输入、输出数据进行傅里叶变换再相除，得到超声骨刀的频响函数，由此确定谐振频率。在仿真过程中超声骨刀采用恒压驱动策略，且使驱动电压频率＝谐振频率，以模拟骨刀的追频功能。
步骤6：基于Matlab-Simulink仿真程序，利用数据驱动建模方法，以反节点振动位移x1的有效值x1rms、振动频率fx1为输入，输出切削力的有效值Frms，获得反节点振动位移x1与切削力F的关系。
首先，将刚度系数设置为m水平，将阻尼系数设置为n水平，利用全因子实验设计m*n组仿真实验。接着，根据设计的实验在MATLAB-Simulink仿真模型改变刚度系数和阻尼系数，并重置驱动电压频率，使其与刚度系数和阻尼系数对应的谐振频率相等进行仿真。然后记录反节点振动位移x1的有效值x1rms、振动频率fx1(与谐振频率相等)和切削力F的有效值Frms。最后，利用数据驱动建模方法获得反节点振动位移x1与切削力F的关系。
步骤7：利用电涡流位移传感器测量反节点振动位移的有效值x1rms及振动频率fx1，根据步骤6获得的关系计算切削力的有效值Frms。
如图5反节点振动位移测量示意图，包括外壳1、电涡流位移传感器2、标记物3、变幅杆4、安装有分析软件的计算机5、信号发生器6、功率放大器7。在变幅杆反节点x1位置安装标记物，标记物随变幅杆一起轴向方向往复振动。在骨刀外壳安装电涡流位移传感器，外壳与变幅杆节点连接，无运动。
电涡流位移传感器2采用基恩士EX-305V，采样频率40KHz，用于测量反节点振动位移x1信号。
信号采集系统5采用NI 9234采集卡，用于采集来自位移位传感器的信号；
安装有分析软件的计算机1用于对信号采集系统获得的位移信号进行分析，获得反节点振动位移x1的有效值x1rms、振动频率fx，进而计算切削力F的有效值Frms。

Claims (8)

1.基于反节点振动位移的超声骨切削力的间接测量方法，其特征在于，包括以下步骤：
步骤1，确定超声骨刀节点位置，得到超声骨刀的简化模型；
步骤2，建立压电换能器模型，计算模型参数，获得压电换能器的振动位移；
步骤3，建立阶梯变幅杆模型，根据特征值、特征向量和能量守恒原则计算动力学参数；
步骤4，建立骨组织的力负载模型，辨识模型参数；
步骤5，建立超声骨刀-骨组织相互作用模型，搭建Matlab-Simulink仿真程序；
步骤6，基于Matlab-Simulink仿真程序，利用数据驱动建模方法，以反节点振动位移x1的有效值x1rms、振动频率fx1为输入，输出切削力的有效值Frms，获得反节点振动位移x1与切削力F的关系；
步骤7，利用电涡流位移传感器测量反节点振动位移的有效值x1rms及振动频率fx1，根据步骤6获得的关系计算切削力的有效值Frms。
2.根据权利要求1所述的基于反节点振动位移的超声骨切削力的间接测量方法，其特征在于，步骤1所述根据超声骨刀节点位置、个数，将超声骨刀简化为压电换能器模型、阶梯变幅杆模型两部分。
3.根据权利要求1所述的基于反节点振动位移的超声骨切削力的间接测量方法，其特征在于，步骤2所述在不考虑压电材料的非线性的假设下，利用线性关系来描述压电陶瓷的电气和机械行为之间的相互作用，建立压电换能器模型。
4.根据权利要求1所述的基于反节点振动位移的超声骨切削力的间接测量方法，其特征在于，步骤3所述根据Euler-Lagrange方程建立阶梯变幅杆模型，根据特征值、特征向量和能量守恒原则计算其动力学参数。
5.根据权利要求1所述的基于反节点振动位移的超声骨切削力的间接测量方法，其特征在于，步骤4所述在超声骨刀切削过程中，刀具与骨组织之间的相互作用会对振动系统产生力负载，其可以等效成一个弹性系数为k的线性弹簧和一个阻尼系数为c的阻尼器的模型。
6.根据权利要求1所述的基于反节点振动位移的超声骨切削力的间接测量方法，其特征在于，步骤5所述根据压电换能器模型、阶梯变幅杆模型和力负载模型建立超声骨刀-骨组织相互作用模型，搭建Matlab-Simulink仿真程序。
7.根据权利要求1所述的基于反节点振动位移的超声骨切削力的间接测量方法，其特征在于，步骤6所述基于MATLAB-Simulink仿真程序，建立以反节点振动位移x1的有效值x1rms、振动频率fx1为输入，以切削力的有效值Frms为输出的数据集，利用数据驱动建模方法，获得反节点振动位移x1与切削力F的关系。
8.根权利要求1所述的基于反节点振动位移的超声骨切削力的间接测量方法，其特征在于，步骤7所述反节点振动位移测量平台组成，包括外壳，电涡流位移传感器，标记物，变幅杆，安装有分析软件的计算机，信号发生器，功率放大器，在变幅杆反节点x1位置安装标记物，标记物随变幅杆一起轴向方向往复振动，在骨刀外壳安装电涡流位移传感器，外壳与变幅杆节点连接，无运动。
CN202210936098.6A 2022-08-05 2022-08-05 基于反节点振动位移的超声骨切削力的间接测量方法 Pending CN115438533A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210936098.6A CN115438533A (zh) 2022-08-05 2022-08-05 基于反节点振动位移的超声骨切削力的间接测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210936098.6A CN115438533A (zh) 2022-08-05 2022-08-05 基于反节点振动位移的超声骨切削力的间接测量方法

Publications (1)

Publication Number Publication Date
CN115438533A true CN115438533A (zh) 2022-12-06

Family

ID=84243193

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210936098.6A Pending CN115438533A (zh) 2022-08-05 2022-08-05 基于反节点振动位移的超声骨切削力的间接测量方法

Country Status (1)

Country Link
CN (1) CN115438533A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117574741A (zh) * 2024-01-17 2024-02-20 西安稀有金属材料研究院有限公司 一种有限元分析下二次结构调整超声骨刀频率的方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117574741A (zh) * 2024-01-17 2024-02-20 西安稀有金属材料研究院有限公司 一种有限元分析下二次结构调整超声骨刀频率的方法
CN117574741B (zh) * 2024-01-17 2024-04-30 西安稀有金属材料研究院有限公司 一种有限元分析下二次结构调整超声骨刀频率的方法

Similar Documents

Publication Publication Date Title
Rani et al. Computational modeling and experimental studies of the dynamic performance of ultrasonic horn profiles used in plastic welding
Rupitsch et al. Inverse method to estimate material parameters for piezoceramic disc actuators
Ghenna et al. Vector control of piezoelectric transducers and ultrasonic actuators
Li et al. A parametric study for the design of an optimized ultrasonic percussive planetary drill tool
CN115438533A (zh) 基于反节点振动位移的超声骨切削力的间接测量方法
Leniowska et al. MFC sensors and actuators in active vibration control of the circular plate
Pérez-Sánchez et al. Numerical design and analysis of a langevin power ultrasonic transducer for acoustic cavitation generation
Giraud et al. Vector control method applied to a traveling wave in a finite beam
Shiyang et al. Analysis of the temperature field distribution for piezoelectric plate-type ultrasonic motor
JP2001350741A (ja) 振動解析の方法および装置ならびにコンピュータ読み取り可能な記録媒体
Zhang et al. A novel Z-shaped elastic flange structure for increasing the amplitude output of a piezoelectric ultrasonic transducer
Jalili et al. Modeling a multilayer piezo-electric transducer by equivalent electro-mechanical admittance matrix
Zhang et al. Numerical analysis and experimental verification of synchronized switching damping systems
Reichel et al. Efficient numerical modeling of oscillatory fluid-structure interaction
Al Sarraf et al. Design optimization of two diagonal slotted block horn for ultrasonic plastic welding
Klinge Determination of the geometry of the RVE for cancellous bone by using the effective complex shear modulus
Nader et al. Determination of piezoelectric transducer damping by using experimental and finite element simulations
JP2002148169A (ja) 定断面作動体の出力圧の測定方法及び測定装置
Bellizzi et al. Response spectral densities and identification of a randomly excited non-linear squeeze film oscillator
Qian et al. Design of an ultrasonic concentrator for vibro-tactile sensors using electro-mechanical analogy
Du et al. A 2D Ultrasonic-Assisted Turning Tool With Function of In Situ Vibration Amplitude Self-Sensing
Simkovics et al. Finite element analysis of hysteresis effects in piezoelectric transducers
Campos‐Pozuelo et al. Finite amplitude standing waves in metallic rods
Jiang et al. Design of Sandwich Transducer Based on the Equivalent Length Algorithm.
Pan et al. A comparison of modal expansion and travelling wave methods for predicting energy flow in beam structures

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