CN112451089A - 一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法 - Google Patents

一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法 Download PDF

Info

Publication number
CN112451089A
CN112451089A CN202011092538.1A CN202011092538A CN112451089A CN 112451089 A CN112451089 A CN 112451089A CN 202011092538 A CN202011092538 A CN 202011092538A CN 112451089 A CN112451089 A CN 112451089A
Authority
CN
China
Prior art keywords
bone
axis
coordinate system
distal
xyz
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.)
Granted
Application number
CN202011092538.1A
Other languages
English (en)
Other versions
CN112451089B (zh
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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202011092538.1A priority Critical patent/CN112451089B/zh
Publication of CN112451089A publication Critical patent/CN112451089A/zh
Application granted granted Critical
Publication of CN112451089B publication Critical patent/CN112451089B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/101—Computer-aided simulation of surgical operations
    • A61B2034/102—Modelling of surgical devices, implants or prosthesis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/101—Computer-aided simulation of surgical operations
    • A61B2034/105—Modelling of the patient, e.g. for ligaments or bones
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/107—Visualisation of planned trajectories or target regions
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/108—Computer aided selection or customisation of medical implants or cutting guides

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Primary Health Care (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • Epidemiology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Robotics (AREA)
  • Urology & Nephrology (AREA)
  • Surgical Instruments (AREA)

Abstract

本发明涉及一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法，首先建立骨‑固定器系统模型，进而建立基准坐标系，根据远端骨末端中心点从初始畸形位置到目标位置的线性轨迹约束方程，分别建立远端骨末端中心点的位置等间距离散约束方程和远端骨的姿态等间距离散约束方程，进而完成胫骨畸形矫正的位姿空间等间距轨迹规划，最终得到光顺平滑且稳定的矫正轨迹和新生骨形状。该方法可以为临床医生在术前轨迹规划、术前矫正策略评估、框架配置提供一定的指导。

Description

一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法
技术领域
本发明属于骨科畸形矫正技术领域，具体涉及一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法。
背景技术
胫骨畸形严重影响着人们的生活和健康，临床上常见的胫骨畸形主要表现为：内/外旋、侧方成角及短缩畸形。重塑骨轴线和骨形态可以稳定患者的行走姿态，提高骨的质量和力学强度。并联外固定器允许对骨畸形进行逐渐矫正，通过调节各支链的长度，以1mm/天的速度将骨段以离散的步骤拉开，进而刺激骨生长，直到新的骨组织填充骨间隙。
临床上，医生通常手动调整每个支链以完成畸形矫正，这是一个复杂的过程，如果没有准确的术前轨迹规划，通常会产生大的矫正误差和不良的临床结果。在进行骨畸形矫正的过程中，期望的骨矫正轨迹被认为是远端骨末端从初始畸形位置到目标位置的空间等距离散轨迹。临床医生通常将外固定器将初始畸形位置到目标位置时每个支链的总变化量进行均等分割，进而对每个支链各自以相等的调节量调节支链的长度，以获得远端骨末端位置和姿态的等距离散变化。然而，由于并联外固定器在关节空间和操作空间具有强非线性映射关系，因此，通过该矫正方法将导致矫正轨迹偏离期望的轨迹。
因此，为了避免这种由于并联外固定器在关节空间和操作空间具有强非线性映射关系而导致的矫正轨迹偏离期望的轨迹，本发明提出了一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法，以实现对胫骨畸形矫正的位姿空间等间距轨迹规划，进而得到从初始畸形位置到目标位置的光顺平滑且稳定的矫正轨迹和新生骨形状。该方法可以为临床医生在术前轨迹规划、术前矫正策略评估、框架配置提供一定的指导。
发明内容
本发明的目的是为了避免由于并联外固定器在关节空间和操作空间具有强非线性映射关系而导致的矫正轨迹偏离期望的轨迹的问题，以实现对胫骨畸形矫正的位姿空间等间距轨迹规划，进而得到从初始畸形位置到目标位置的光顺平滑且稳定的矫正轨迹和新生骨形状。
为了实现以上目的，本发明所采用的技术方案是一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法，包括如下步骤：
(1)确定胫骨的畸形参数和外固定器参数；
(2)骨-固定器系统模型建立：近端骨垂直且刚性地固定于近端环的中心，远端骨垂直且刚性地固定器远端环中心，进而建立骨-固定器系统模型；
(3)基准坐标系建立：在近端环的中心建立固定坐标系O-xyz，在近端骨末端中心点位置建立坐标系O1-x1y1z1，在远端环中心建立动坐标系G-uvw，在远端骨末端中心点位置建立坐标系G1-u1v1w1；
(4)求解外固定器运动学位置逆解；
(5)位姿空间等间距轨迹规划方法建立：首先建立畸形矫正轨迹约束方程，随后分别对远端骨末端的位置和姿态进行等间距离散规划；
(6)获取矫正处方；
(7)生成骨的矫正轨迹及新生骨的形状；
所述步骤(1)中，设定胫骨的畸形参数：假设畸形的胫骨绕固定坐标系O-xyz的x轴的角度畸形为θ1，绕固定坐标系O-xyz的y轴的角度畸形为θ2，绕固定坐标系O-xyz的z轴的旋转畸形为θ3，沿固定坐标系O-xyz的x轴的平移畸形为ex，沿固定坐标系O-xyz的y轴的平移畸形为ey，沿固定坐标系O-xyz的z轴的缩短畸形为ez，近端骨段OO1和远端骨段GG1的长度均为HL，外固定器的近端环和远端环直径均为dR。假设总的畸形矫正步数为n，因此，远端骨末端的初始畸形位置和姿态可以表示为Γ0(XG10,YG10,ZG10,θ1,θ2,θ3)，远端骨末端在步骤h时的位置和姿态可以表示为Γh(XG1h,YG1h,ZG1h,θ1h,θ2h,θ3h)，远端骨末端的目标位置和姿态表示为Γn(0,0,HL+ez,0,0,0)；
所述步骤(2)中，在角度旋转中心处将畸形骨截骨为近端骨和远端骨两段，在初始畸形位置，畸形骨的近端骨通过克氏针垂直且刚性地固定在近端环中心，畸形骨的远端骨通过克氏针垂直且刚性地固定在远端环中心，近端骨作为参考骨，当畸形被完全矫正时，外固定器到达中性位置，远端骨与近端骨轴线对齐，进而完成骨-固定器系统建立；
所述步骤(3)中，固定坐标系O-xyz建立在近端环的中心，z轴垂直于近端环向下，坐标系O1-x1y1z1建立在近端骨末端中心点位置并且与固定坐标系O-xyz平行且方向一致；
相应地，动坐标系G-uvw建立在远端环中心，w轴垂直于远端环向下，坐标系G1-u1v1w1建立在远端骨末端中心点位置并且与动定坐标系G-uvw平行且方向一致，外固定器在中性位置时，固定坐标系O-xyz和动坐标系G-uvw平行且方向一致，进而完成基准坐标系建立；
所述步骤(4)中，根据封闭矢量环法求解外固定器运动学位置逆解；
所述步骤(5)中，设总的畸形矫正步数为n，远端骨末端中心点从初始畸形位置到目标位置的均匀间隔的线性轨迹就是目标矫正轨迹，该轨迹的约束方程用式(1)表示：
Figure BDA0002722626310000021
其中G10(XG10,YG10,ZG10)和G1n(XG1n,YG1n,ZG1n)分别是固定坐标系O-xyz中步骤0和步骤n处远端骨末端中心点G1的位置矢量，其中步骤0表示初始畸形位置，步骤n表示目标位置；
将远端骨末端中心点G1分别沿固定坐标系O-xyz的x轴、y轴、z轴的位置总调节量进行等距离散分割为n个部分，进而均等规划远端骨末端中心点的G1位置轨迹，可用式(2)表示：
Figure BDA0002722626310000031
进一步，远端骨末端中心点G1在固定坐标系O-xyz中步骤h(h＝1、2,...,n)处的位置G1h(XG1h,YG1h,ZG1h)的约束方程用式(3)表示，且该方程为远端骨末端中心点G1的位置等间距离散约束方程：
Figure BDA0002722626310000032
将远端骨分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态总调节量进行等距离散分割为n个部分，进而均等规划远端骨的姿态变化轨迹，可用式(4)表示：
Figure BDA0002722626310000033
其中θ1,θ2,θ3分别表示远端骨在初始畸形位置时，分别绕固定坐标系O-xyz的x轴、y轴、z轴的畸形角度，也是远端骨在初始畸形位置时的姿态，为了公式的统一性，令θ10,θ20,θ30表示在固定坐标系O-xyz中步骤0处分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态，θ1n,θ2n,θ3n表示在固定坐标系O-xyz中步骤n处分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态；
远端骨在固定坐标系O-xyz中步骤h(h＝1、2,...,n)处的姿态θh(θ1h,θ2h,θ3h)的约束方程用式(5)表示，且该方程为远端骨的姿态等间距离散约束方程：
Figure BDA0002722626310000041
其中θ1h,θ2h,θ3h(h＝1,2,...,n)表示在固定坐标系O-xyz中步骤h处分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态；
通过分别建立远端骨末端中心点G1的位置等间距离散约束方程和远端骨的姿态等间距离散约束方程，进而完成胫骨畸形矫正的位姿空间等间距轨迹规划。
所述步骤(6)中，基于位姿空间等间距轨迹规划方法，确定矫正过程中每个矫正步骤时远端骨末端的位置和姿态Γh(XG1h,YG1h,ZG1h,θ1h,θ2h,θ3h)，通过运动学位置逆解求得在矫正步h处的支链长度lh(l1h，l2h，l3h，l4h，l5h，l6h)，在矫正过程中得到的每个支链的长度就是所求得畸形矫正处方，其中lih(i＝1,2,...,6)表示在矫正步h处并联外固定器标号为i的支链的长度。
所述步骤(7)中，将步骤(6)中通过位姿空间等间距轨迹规划方法得到的一些列远端骨末端的位置和姿态Γh(XG1h,YG1h,ZG1h,θ1h,θ2h,θ3h)通过Matlab数值软件拟合生成骨的矫正轨迹，采用以位置为中心点的圆面在空间的姿态表示远端骨末端面的姿态，进而生成从初始畸形位置到目标位置的均匀间隔的线性轨迹，新生骨的形状规则且光顺。
本发明与现有技术相比，具有以下有益效果：
本发明的一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法，通过分别对远端骨末端的位置和姿态进行空间等间距离散，得到了从初始畸形位置到目标位置的光顺平滑且稳定骨的矫正轨迹和新生骨形状，避免了由于并联外固定器在关节空间和操作空间具有强非线性映射关系而导致的矫正轨迹偏离期望的轨迹的问题，该方法可以为临床医生在术前轨迹规划、术前矫正策略评估、框架配置提供一定的指导。
附图说明
图1是本发明的一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法流程图；
图2是畸形骨的初始畸形位置和目标位置示意图；
图3是骨-固定器系统初始畸形位置；
图4是骨-固定器系统目标位置；
图5是矫正过程中每个支链长度变化图；
图6是畸形矫正轨迹及新生骨的形状。
具体实施方式
下面结合附图及具体的实施例对本发明进一步介绍，所述是对本发明的解释而不是限定。
参图1所示为本发明的一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法流程图，由图可知，包括以下步骤：
(1)确定胫骨的畸形参数和外固定器参数
所述步骤(1)中，参图2、图3和图4所示，设定胫骨的畸形参数和外固定器参数：假设畸形的胫骨绕固定坐标系O-xyz的x轴的角度畸形为14°，绕固定坐标系O-xyz的y轴的角度畸形为22°，绕固定坐标系O-xyz的z轴的旋转畸形为20°，沿固定坐标系O-xyz的x轴的平移畸形为3mm，沿固定坐标系O-xyz的y轴的平移畸形为2mm，沿固定坐标系O-xyz的z轴的缩短畸形为20mm，近端骨段OO1和远端骨段GG1的长度均为70mm，外固定器的近端环和远端环直径均为162mm，假设总的畸形矫正步数为20，以保证沿固定坐标系O-xyz的z轴为1mm/day的最佳牵引速率，因此，远端骨末端的初始畸形位置和姿态可以表示为Γ0(3,2,70,14,22,20)，远端骨末端在步骤h时的位置和姿态可以表示为Γh(XG1h,YG1h,ZG1h,θ1h,θ2h,θ3h)，远端骨末端的目标位置和姿态表示为Γn(0,0,90,0,0,0)；
(2)骨-固定器系统模型建立
所述步骤(2)中，参图3和图4所示，在角度旋转中心处将畸形骨截骨为近端骨和远端骨两段，在初始畸形位置，畸形骨的近端骨通过克氏针垂直且刚性地固定在近端环中心，畸形骨的远端骨通过克氏针垂直且刚性地固定在远端环中心，近端骨作为参考骨，当畸形被完全矫正时，外固定器到达中性位置，也是目标位置，远端骨轴线与近端骨轴线对齐，进而完成骨-固定器系统建立；
(3)基准坐标系建立
所述步骤(3)中，参图3和图4所示，固定坐标系O-xyz建立在近端环的中心，z轴垂直于近端环向下，坐标系O1-x1y1z1建立在近端骨末端中心点位置并且与固定坐标系O-xyz平行且方向一致；
相应地，动坐标系G-uvw建立在远端环中心，w轴垂直于远端环向下，坐标系G1-u1v1w1建立在远端骨末端中心点位置并且与动定坐标系G-uvw平行且方向一致，外固定器在中性位置时，固定坐标系O-xyz和动坐标系G-uvw平行且方向一致，进而完成基准坐标系建立；
(4)求解外固定器运动学位置逆解
所述步骤(4)中，根据封闭矢量环法求解外固定器运动学位置逆解；
(5)位姿空间等间距轨迹规划方法建立
所述步骤(5)中，设总的畸形矫正步数为n，远端骨末端中心点G1从初始畸形位置到目标位置的均匀间隔的线性轨迹就是目标矫正轨迹，该轨迹的约束方程用式(1)表示：
Figure BDA0002722626310000061
其中G10(XG10,YG10,ZG10)和G1n(XG1n,YG1n,ZG1n)分别是固定坐标系O-xyz中步骤0和步骤n处远端骨末端中心点G1的位置矢量，其中步骤0表示初始畸形位置，步骤n表示目标位置；
将远端骨末端中心点G1分别沿固定坐标系O-xyz的x轴、y轴、z轴的位置总调节量进行等距离散分割为n个部分，进而均等规划远端骨末端中心点的G1位置轨迹，可用式(2)表示：
Figure BDA0002722626310000062
进一步，远端骨末端中心点G1在固定坐标系O-xyz中步骤h(h＝1、2,...,n)处的位置G1h(XG1h,YG1h,ZG1h)的约束方程用式(3)表示，且该方程为远端骨末端中心点G1的位置等间距离散约束方程：
Figure BDA0002722626310000063
将远端骨分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态总调节量进行等距离散分割为n个部分，进而均等规划远端骨的姿态变化轨迹，可用式(4)表示：
Figure BDA0002722626310000064
其中θ1,θ2,θ3分别表示远端骨在初始畸形位置时，分别绕固定坐标系O-xyz的x轴、y轴、z轴的畸形角度，也是远端骨在初始畸形位置时的姿态，为了公式的统一性，令θ10,θ20,θ30表示在固定坐标系O-xyz中步骤0处分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态，θ1n,θ2n,θ3n表示在固定坐标系O-xyz中步骤n处分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态；
进一步，远端骨在固定坐标系O-xyz中步骤h(h＝1、2,...,n)处的姿态θh(θ1h,θ2h,θ3h)的约束方程用式(5)表示，且该方程为远端骨的姿态等间距离散约束方程：
Figure BDA0002722626310000071
其中θ1h,θ2h,θ3h(h＝1,2,...,n)表示在固定坐标系O-xyz中步骤h处分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态；
通过分别建立远端骨末端中心点G1的位置等间距离散约束方程和远端骨的姿态等间距离散约束方程，进而完成胫骨畸形矫正的位姿空间等间距轨迹规划。
(6)获取矫正处方；
所述步骤(6)中，参图(5)所示，将步骤(1)中胫骨的畸形参数和外固定器参数代入位姿空间等间距轨迹规划方法，通过Matlab数值软件求解矫正过程中每个矫正步骤的远端骨末端的位置和姿态Γh(XG1h,YG1h,ZG1h,θ1h,θ2h,θ3h)，通过运动学位置逆解求得每个矫正步骤处的支链长度lh(l1h，l2h，l3h，l4h，l5h，l6h)，在矫正过程中得到的每个支链的长度就是所求得的畸形矫正处方，其中lih(i＝1,2,...,6)表示在矫正步h处并联外固定器标号为i的支链的长度。
(7)生成骨的矫正轨迹及新生骨的形状；
所述步骤(7)中，参图(6)所示，将步骤(6)中通过位姿空间等间距轨迹规划方法得到的一系列远端骨末端的位置和姿态Γh(XG1h,YG1h,ZG1h,θ1h,θ2h,θ3h)，代入Matlab数值软件拟合生成骨的矫正轨迹，以远端骨末端中心点G1的位置为中心点的圆面在空间的姿态表示远端骨末端面的姿态，进而生成位姿空间等间距轨迹规划下的骨的矫正轨迹及新生骨的形状，参图(6)所示，生成了从初始畸形位置到目标位置的均匀间隔的线性轨迹，新生骨的形状规则且光顺。
上述实施例仅是针对本发明的可行性实施方式的具体说明，它们并非用以限制本发明的保护范围，对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。

Claims (1)

1.一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法，其特征在于，包括以下三个步骤：
(1)骨-固定器系统模型建立：近端骨垂直且刚性地固定于近端环的中心，远端骨垂直且刚性地固定器远端环中心，进而建立骨-固定器系统模型；
(2)基准坐标系建立；
(3)位姿空间等间距轨迹规划方法建立：首先建立畸形矫正轨迹约束方程，随后分别对远端骨末端的位置和姿态进行等间距离散规划；
步骤(1)中，在角度旋转中心处将畸形骨截骨为近端骨和远端骨两段，在初始畸形位置，畸形骨的近端骨通过克氏针垂直且刚性地固定在近端环中心，畸形骨的远端骨通过克氏针垂直且刚性地固定在远端环中心，近端骨作为参考骨，当畸形被完全矫正时，外固定器到达中性位置，远端骨与近端骨轴线对齐；
步骤(2)中，固定坐标系O-xyz建立在近端环的中心，z轴垂直于近端环向下，坐标系O1-x1y1z1建立在近端骨末端中心点位置并且与固定坐标系O-xyz平行且方向一致；
相应地，动坐标系G-uvw建立在远端环中心，w轴垂直于远端环向下，坐标系G1-u1v1w1建立在远端骨末端中心点位置并且与动定坐标系G-uvw平行且方向一致，外固定器在中性位置时，固定坐标系O-xyz和动坐标系G-uvw平行且方向一致；
步骤(3)中，设总的畸形矫正步数为n，远端骨末端中心点从初始畸形位置到目标位置的均匀间隔的线性轨迹就是目标矫正轨迹，该轨迹的约束方程用式(1)表示：
Figure FDA0002722626300000011
其中G10(XG10,YG10,ZG10)和G1n(XG1n,YG1n,ZG1n)分别是固定坐标系O-xyz中在步骤0和步骤n处远端骨末端中心点G1的位置矢量，其中步骤0表示位于初始畸形位置，步骤n表示到达目标位置；
将远端骨末端中心点G1分别沿固定坐标系O-xyz的x轴、y轴、z轴的位置总调节量进行等距离散分割为n个部分，进而均等规划远端骨末端中心点的G1位置轨迹，可用式(2)表示：
Figure FDA0002722626300000012
进一步，远端骨末端中心点G1在固定坐标系O-xyz中步骤h(h＝1、2,...,n)处的位置G1h(XG1h,YG1h,ZG1h)的约束方程用式(3)表示，且该方程为远端骨末端中心点G1的位置等间距离散约束方程：
Figure FDA0002722626300000021
将远端骨分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态总调节量进行等距离散分割为n个部分，进而均等规划远端骨的姿态变化轨迹，用式(4)表示：
Figure FDA0002722626300000022
其中θ1,θ2,θ3分别表示远端骨在初始畸形位置时，分别绕固定坐标系O-xyz的x轴、y轴、z轴的畸形角度，也是远端骨在初始畸形位置时的姿态，令θ10,θ20,θ30表示在固定坐标系O-xyz中步骤0处分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态，θ1n,θ2n,θ3n表示在固定坐标系O-xyz中步骤n处分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态；
远端骨在固定坐标系O-xyz中步骤h处的姿态θh(θ1h,θ2h,θ3h)的约束方程用式(5)表示，其中h＝1、2,...,n；且该方程为远端骨的姿态等间距离散约束方程：
Figure FDA0002722626300000023
其中θ1h,θ2h,θ3h(h＝1,2,...,n)表示在固定坐标系O-xyz中步骤h处分别绕固定坐标系O-xyz的x轴、y轴、z轴的姿态；
通过分别建立远端骨末端中心点G1的位置等间距离散约束方程和远端骨的姿态等间距离散约束方程，进而完成胫骨畸形矫正的位姿空间等间距轨迹规划。
CN202011092538.1A 2020-10-13 2020-10-13 一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法 Active CN112451089B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011092538.1A CN112451089B (zh) 2020-10-13 2020-10-13 一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011092538.1A CN112451089B (zh) 2020-10-13 2020-10-13 一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法

Publications (2)

Publication Number Publication Date
CN112451089A true CN112451089A (zh) 2021-03-09
CN112451089B CN112451089B (zh) 2023-01-20

Family

ID=74833784

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011092538.1A Active CN112451089B (zh) 2020-10-13 2020-10-13 一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法

Country Status (1)

Country Link
CN (1) CN112451089B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05297916A (ja) * 1992-04-16 1993-11-12 Nachi Fujikoshi Corp ロボットの軌跡制御方法
RU2013154004A (ru) * 2013-12-04 2015-07-20 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования Нижегородский государственный технический университет им. Р.Е. Алексеева (НГТУ) Способ восстановления целостности и функциональности пяточной кости человека и аппарат для его осуществления
CN108245122A (zh) * 2018-01-12 2018-07-06 北京理工大学 磁引导式胶囊内窥镜系统及轨迹规划方法
CN110009724A (zh) * 2019-02-28 2019-07-12 成都真实维度科技有限公司 一种根据关节非骨质建模分区体积比的矫正分析方法
CN110104216A (zh) * 2019-01-28 2019-08-09 西北工业大学深圳研究院 一种针对运动学冗余双臂空间机器人的协同路径规划方法
WO2020029378A1 (zh) * 2018-08-10 2020-02-13 天津大学 骨折并联外固定支架空间位姿的自动识别方法
CN110974414A (zh) * 2019-12-11 2020-04-10 北京工业大学 一种用于胫骨畸形矫正的空间锥-柱干涉检验方法
US20200316779A1 (en) * 2019-04-08 2020-10-08 Teradyne, Inc. System and method for constraint management of one or more robots

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05297916A (ja) * 1992-04-16 1993-11-12 Nachi Fujikoshi Corp ロボットの軌跡制御方法
RU2013154004A (ru) * 2013-12-04 2015-07-20 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования Нижегородский государственный технический университет им. Р.Е. Алексеева (НГТУ) Способ восстановления целостности и функциональности пяточной кости человека и аппарат для его осуществления
CN108245122A (zh) * 2018-01-12 2018-07-06 北京理工大学 磁引导式胶囊内窥镜系统及轨迹规划方法
WO2020029378A1 (zh) * 2018-08-10 2020-02-13 天津大学 骨折并联外固定支架空间位姿的自动识别方法
CN110104216A (zh) * 2019-01-28 2019-08-09 西北工业大学深圳研究院 一种针对运动学冗余双臂空间机器人的协同路径规划方法
CN110009724A (zh) * 2019-02-28 2019-07-12 成都真实维度科技有限公司 一种根据关节非骨质建模分区体积比的矫正分析方法
US20200316779A1 (en) * 2019-04-08 2020-10-08 Teradyne, Inc. System and method for constraint management of one or more robots
CN110974414A (zh) * 2019-12-11 2020-04-10 北京工业大学 一种用于胫骨畸形矫正的空间锥-柱干涉检验方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
孙昊等: "基于蚁群算法的骨外固定器处方优化", 《仪器仪表学报》 *
廖筑秀等: "新型微创技术在胎儿外科中的应用与展望", 《医疗卫生装备》 *
李剑锋等: "Tricept型并联机床曲面加工的刀轨规划", 《北京工业大学学报》 *
熊涛等: "下肢康复机器人屈伸训练轨迹规划", 《机械设计与制造》 *

Also Published As

Publication number Publication date
CN112451089B (zh) 2023-01-20

Similar Documents

Publication Publication Date Title
US11653978B2 (en) External fixator deformity correction systems and methods
US20210216671A1 (en) Method making it possible to produce the ideal curvature of a rod of vertebral osteosynthesis material designed to support a patient's vertebral column
Yang et al. Three-dimensional printing of patient-specific surgical plates in head and neck reconstruction: a prospective pilot study
CN108742804B (zh) 一种用于骨折复位的并联外固定支架
EP2767252B1 (en) Computer implemented method for planning deformity correction
US10321959B2 (en) Method of surgical planning
JP6486913B2 (ja) ヒトまたは動物の骨組織の切断および切断の計画
Majdani et al. A robot-guided minimally invasive approach for cochlear implant surgery: preliminary results of a temporal bone study
US10555815B2 (en) Surgical kit for cartilage repair comprising implant and a set of tools
JP6521707B2 (ja) 穿刺プランニング装置及び穿刺システム
US10952820B2 (en) System of manufacturing orthodontic wire, method for manufacturing the orthodontic wire using the same, and orthodontic wire bending machine for performing the same
EP2680233A2 (en) Method and system for providing a face adjustment image
CN112451089B (zh) 一种用于胫骨畸形矫正的位姿空间等间距轨迹规划方法
CN110141363B (zh) 一种基于结构光扫描的脊柱多阶段配准系统
US9028252B1 (en) Device, system and method for in-situ drill guide sleeve orientation
CN109106457A (zh) 一种基于精准化定制的牙齿活动矫治器用卡环的制造方法
CN110974414B (zh) 一种用于胫骨畸形矫正的空间锥-柱干涉检验方法
US20060239577A1 (en) Process of using computer modeling, reconstructive modeling and simulation modeling for image guided reconstructive surgery
Servi et al. CAD-based automatic modelling of customized cutting templates for pectus arcuatum surgical correction
KR100572259B1 (ko) 원형 체외 고정기의 스트럿 길이 결정 방법 및 기록매체
Zhou et al. The accuracy of computer-aided simulation system protocol for positioning the maxilla with a intermediate splint in orthognathic surgery
CN117133398A (zh) 基于六轴空间外固定支架的骨折骨畸治疗用电子处方生成方法
WO2018191903A1 (zh) 植入牙机器人系统及其操作方法
CN117174224A (zh) 基于Ortho-SUV外固定支架的骨折骨畸治疗电子处方生成方法
CN118000761A (zh) 一种种植体稳定性相关参数的测量方法和系统

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