WO2024077889A1 - Operation prompting method and apparatus for flexible instrument, device, and storage medium - Google Patents

Operation prompting method and apparatus for flexible instrument, device, and storage medium Download PDF

Info

Publication number
WO2024077889A1
WO2024077889A1 PCT/CN2023/085904 CN2023085904W WO2024077889A1 WO 2024077889 A1 WO2024077889 A1 WO 2024077889A1 CN 2023085904 W CN2023085904 W CN 2023085904W WO 2024077889 A1 WO2024077889 A1 WO 2024077889A1
Authority
WO
WIPO (PCT)
Prior art keywords
flexible
displacement
finite element
flexible instrument
instrument
Prior art date
Application number
PCT/CN2023/085904
Other languages
French (fr)
Chinese (zh)
Inventor
陈浩
刘宏斌
Original Assignee
中国科学院自动化研究所
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 中国科学院自动化研究所 filed Critical 中国科学院自动化研究所
Publication of WO2024077889A1 publication Critical patent/WO2024077889A1/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B90/00Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Definitions

  • the present application relates to the technical field of medical devices, and in particular to an operation prompt method, device, equipment and storage medium for a flexible device.
  • the embodiments of the present application provide an operation prompt method, device, equipment and storage medium for a flexible device.
  • an operation prompt method for a flexible device comprising:
  • the finite element model of the flexible device includes a plurality of finite element nodes
  • the displacement of the flexible instrument is determined;
  • the target constraint conditions are used to represent the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body;
  • the displacement of the flexible instrument includes the displacement of each finite element node;
  • operation instruction information of the flexible instrument is output; the operation instruction information of the flexible instrument is used to indicate the way of performing surgery on the natural cavity of the human body through the flexible instrument.
  • the determination of the displacement of the flexible device based on the minimum potential energy principle and the target constraint condition includes:
  • the displacement of the flexible device is determined when the potential energy change of the flexible device is minimal and the positional relationship between the flexible device and the natural cavity of the human body meets the target constraint condition;
  • the target constraint condition is that the finite element nodes in the finite element model of the flexible device and the contact point target surfaces in the three-dimensional model of the natural cavity of the human body are not penetrable.
  • the determination of the displacement of the flexible device based on the minimum potential energy principle and the target constraint condition includes:
  • determining the force result of the flexible device according to the displacement of the flexible device includes:
  • the force result of the flexible instrument is determined according to the displacement of the flexible instrument and the stiffness matrix of the flexible instrument.
  • the method further comprises:
  • the stiffness matrix of the flexible device is determined based on a rotation matrix between a local coordinate system of a unit and a global coordinate system in a finite element model of the flexible device.
  • the method further includes:
  • the natural cavity of the human body is a narrow bifurcated cavity
  • the target constraint condition is updated.
  • the embodiment of the present application further provides an operation prompt device for a flexible device, comprising:
  • An establishment module is used to establish a finite element model of a flexible device; the finite element model of the flexible device includes a plurality of finite element nodes;
  • the first determination module is used to determine the displacement of the flexible instrument based on the minimum potential energy principle and the target constraint condition;
  • the target constraint condition is used to represent the constraint of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target surface patch in the three-dimensional model of the natural cavity of the human body;
  • the displacement of the flexible instrument includes the displacement of each finite element node;
  • a second determination module is used to determine the force result of the flexible device according to the displacement of the flexible device
  • the prompt module is used to output operation instruction information of the flexible instrument according to the determined displacement of the flexible instrument and the force result of the flexible instrument; the operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
  • an embodiment of the present application further provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the operation prompt method for the flexible device as described in the first aspect is implemented.
  • an embodiment of the present application further provides a non-transitory computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the operation prompt method of the flexible instrument as described in the first aspect is implemented.
  • an embodiment of the present application further provides a computer program product, including a computer program, which, when executed by a processor, implements the operation prompt method of the flexible instrument as described in the first aspect.
  • the operation prompt method, device, equipment and storage medium of the flexible instrument are based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determine that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate; then, based on the determined displacement of the flexible instrument and the operation instruction information of the flexible instrument, the surgical personnel can accurately operate the flexible instrument, thereby improving the efficiency and safety of the operation.
  • FIG1 is a schematic flow chart of an operation prompt method for a flexible device provided in an embodiment of the present application
  • FIGS. 2a-2b are schematic diagrams of the flexible device provided in the embodiments of the present application in contact with the natural cavity of the human body;
  • FIG3 is a schematic flow chart of another operation prompt method of a flexible device provided in an embodiment of the present application.
  • FIG4 is a schematic diagram of the structure of an operation prompt device for a flexible device provided in an embodiment of the present application.
  • FIG5 is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application.
  • the method of the embodiment of the present application can be applied to surgical scenarios based on medical devices to implement operation prompts for flexible devices and improve surgical efficiency and safety.
  • the operation prompt method of the flexible instrument in the embodiment of the present application is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate; then, based on the determined displacement of the flexible instrument and the operation instruction information of the flexible instrument, the surgical personnel can accurately operate the flexible instrument, thereby improving the efficiency and safety of the operation.
  • FIG1 is a flow chart of an embodiment of an operation prompt method for a flexible device provided in an embodiment of the present application. As shown in FIG1 , the method provided in this embodiment includes:
  • Step 101 establishing a finite element model of a flexible device; the finite element model of the flexible device includes a plurality of finite element nodes;
  • the spatial information of medical devices in the patient's body plays a vital role in the accurate operation of medical devices.
  • the flexible device will inevitably change its shape, so it is necessary to accurately determine the force and displacement of the flexible device to provide doctors with accurate information to improve the efficiency and safety of the operation.
  • a finite element model of a flexible device is first established.
  • the finite element model of a flexible device is a model of a flexible device established using a finite element analysis method, and is a group of unit combinations that are connected only at nodes, transmit forces only at nodes, and are constrained only at nodes.
  • the finite element model of a flexible device includes multiple finite element units and finite element nodes.
  • the continuous geometric mechanism corresponding to the flexible device can be discretized into a finite number of flexible device finite element units, and a finite number of flexible device finite element nodes can be set in each flexible device finite element unit, so that the continuum corresponding to the flexible device is regarded as a group of flexible device finite element units that are connected only at the finite element nodes of the flexible device.
  • the flexible device After establishing the finite element model of the flexible device, when determining the displacement and force of the flexible device, the flexible device can be regarded as being composed of many flexible device finite element units and flexible device finite element nodes that are interconnected. After determining the displacement and force of each subunit of the flexible device, the displacement and force of the flexible device can be determined. The deformation of the flexible device can be simulated by a large number of flexible device finite element units and finite element nodes, thereby reducing the difficulty of determining the displacement and force results of the flexible device.
  • Step 102 Determine the displacement of the flexible device based on the minimum potential energy principle and target constraint conditions; the target constraint conditions are used to represent the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible device includes the displacement of each finite element node;
  • the displacement of the flexible device can be determined based on the minimum potential energy principle and the target constraint condition; optionally, the potential energy of the flexible device at a certain time i-1 is expressed as follows:
  • xi-1 is the generalized coordinate column vector of the node displacement at time i-1
  • F is the known total external force
  • the total external force F does not include the interaction force between the flexible device and the natural cavity of the human body
  • K( xi-1 ) is the stiffness matrix of the flexible device
  • the potential energy of the flexible device at the next time i is as follows:
  • ⁇ xi the displacement of the flexible device from time i-1 to time i
  • xi xi-1 + ⁇ xi
  • C a constant term independent of ⁇ xi .
  • the unknown potential energy change ⁇ U related to ⁇ xi is defined as follows:
  • the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement ⁇ xi of the flexible instrument; optionally, the displacement ⁇ xi of the flexible instrument includes the displacement of each finite element node, that is, the displacement ⁇ xi of the flexible instrument is composed and determined by the multi-dimensional finite element node displacement.
  • the displacement ⁇ xi of the flexible instrument is subject to environmental constraints.
  • the embodiment of the present application is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patch in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and also makes the determined displacement of the flexible instrument more accurate.
  • Step 103 determining the force result of the flexible device according to the displacement of the flexible device
  • the force result of the flexible device can be determined based on the determined displacement of the flexible device.
  • the force result of the flexible device can be determined according to Hooke's theorem.
  • the determined force result of the flexible device includes the force result of each finite element node, that is, the determined force result of the flexible device is composed and determined based on the multi-dimensional finite element node force result.
  • Step 104 outputting operation instruction information of the flexible instrument according to the determined displacement of the flexible instrument and the force result of the flexible instrument; the operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
  • the operation instruction information of the flexible instrument can be output based on the determined displacement and force results of the flexible instrument, wherein the instruction information is used to indicate how to operate the flexible instrument to perform surgery on the natural cavity of the human body in the next step.
  • the doctor can also accurately operate the flexible instrument based on the operation instruction information of the flexible instrument, thereby improving the efficiency and safety of the operation.
  • the displacement and force results of the flexible instrument can also be displayed on the screen for the doctor's reference, so as to accurately determine the next surgical action and improve the efficiency and safety of the operation.
  • the method of the embodiment of the present application will neither introduce errors of different collision detection methods, nor assume that the contact point position between the flexible instrument and the natural cavity of the human body remains unchanged during the calculation of two time frames, thereby making the displacement and force results of the flexible instrument determined in the embodiment of the present application more accurate.
  • the speed of determining the displacement and force results of the flexible instrument will be significantly reduced when the number of contact points between the flexible instrument and the natural cavity of the human body increases, and the number of contact points greatly affects the real-time and stability of determining the displacement and force results of the flexible instrument.
  • the method of the embodiment of the present application can calculate the displacement and force of the flexible instrument in real time under different contact conditions between the flexible instrument and the natural cavity of the human body, especially when there are many contact points between the flexible instrument and the natural cavity of the human body, thereby improving the stability and practicality of determining the displacement and force results of the flexible instrument, and avoiding the problem that the efficiency of determining the displacement of the flexible instrument in the related art will decrease with the increase of the collision points between the flexible instrument and the natural cavity of the human body.
  • the embodiment of the present application does not limit the relative movement speed of the contact between the flexible instrument and the natural cavity of the human body, and will not limit the relative movement speed between the two objects of the flexible instrument and the natural cavity of the human body. It can well simulate the physical movement process of an object with a relatively large relative movement speed, so that the displacement and force results of the flexible instrument determined are more accurate.
  • the method of the above embodiment is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate; then, the surgical personnel can accurately operate the flexible instrument based on the determined displacement of the flexible instrument and the operation instruction information of the flexible instrument, thereby improving the efficiency and safety of the operation.
  • determining the displacement of the flexible device includes:
  • the displacement of the flexible device is determined when the potential energy change of the flexible device is minimal and the positional relationship between the flexible device and the natural cavity of the human body meets the target constraint condition;
  • the target constraint condition is that the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the natural cavity of the human body are not penetrable.
  • the displacement of the flexible instrument is determined, which means that the displacement of the flexible instrument is determined within the natural cavity of the human body, and the displacement of the flexible instrument cannot be determined outside the natural cavity of the human body, which makes the displacement of the flexible instrument consistent with the actual surgical situation, so that the displacement of the flexible instrument is more accurate.
  • the contact point target surface is the surface closest to the three-dimensional model of the natural cavity of the human body that the finite element node of the flexible instrument contacts.
  • the method of the above embodiment when determining the displacement of the flexible instrument based on the principle of minimum potential energy, uses the impenetrable constraint condition between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body, so that the determined displacement of the flexible instrument meets the height of the natural cavity of the human body, that is, the determined displacement of the flexible instrument is the displacement within the natural cavity of the human body, and the displacement of the flexible instrument cannot be determined outside the natural cavity of the human body, which makes the determined displacement of the flexible instrument more accurate and more in line with the actual surgical situation, so that based on the determined displacement of the flexible instrument, the next surgical plan can be determined more accurately, thereby improving the efficiency and safety of the operation.
  • determining the displacement of the flexible device includes:
  • the embodiment of the present application is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate, so that the operation prompt information of the flexible instrument can be accurately provided to the surgical personnel, thereby improving the efficiency and safety of the operation.
  • one of the advantages of determining the displacement of the flexible instrument based on the minimum potential energy principle in the embodiment of the present application is the simple expression of the environmental constraints between the flexible instrument and the natural cavity of the human body, which improves the accuracy and efficiency of the determined displacement of the flexible instrument, thereby accurately providing the operating prompt information of the flexible instrument to the surgical personnel, thereby improving the efficiency and safety of the operation.
  • the displacement ⁇ xi of the flexible device is composed and determined by multi-dimensional finite element node displacements.
  • the three-dimensional model of the natural cavity of the human body reconstructed from the three-dimensional CT before surgery is composed of multiple facets in the virtual environment.
  • the inequality constraint of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target facets in the three-dimensional model of the natural cavity of the human body simplifies the collision between the flexible instrument and the natural cavity of the human body into a collision between a point and a surface.
  • the normal vector of the contact point facet of the flexible instrument and the natural cavity of the human body is For any point on the plane, position X 0 , we have or The direction of the inequality of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patch in the three-dimensional model of the natural cavity of the human body depends on the initial value, where S i is the generalized coordinate position of the flexible device in the undeformed state.
  • S i is the generalized coordinate position of the flexible device in the undeformed state.
  • the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the human body's natural cavity are not penetrated.
  • Figure 2a shows the contact point and force diagram between the finite element model of the flexible instrument and the natural cavity model of the human body between two frames in the embodiment of the present application
  • Figure 2b shows the contact point and force diagram between the finite element model of the flexible instrument and the natural cavity model of the human body between two frames in the prior art
  • n 1 to n 6 are finite element nodes of the flexible instrument
  • X is the contact point on the facet of the natural cavity model of the human body
  • F is the contact force.
  • the contact point between the flexible instrument and the natural cavity of the human body between two frames in the embodiment of the present application can be any position within the domain composed of multiple facets in the natural cavity model of the human body, while the prior art assumes that the contact point between the flexible instrument and the natural cavity of the human body between every two frames remains unchanged, and only the force direction and magnitude at the contact point change, thereby causing estimation errors.
  • the displacement of the flexible instrument is determined based on the finite element model of the flexible instrument and the principle of minimum potential energy. Compared with the prior art, the advantage is that the assumption that the contact point between the flexible instrument and the natural cavity of the human body is fixed is removed in the embodiment of the present application.
  • the contact point between the flexible instrument and the natural cavity of the human body can change between two frames, while the prior art assumes that the contact point between the flexible instrument and the natural cavity of the human body remains unchanged in two adjacent frames.
  • the flexible instrument is inserted into the natural cavity at a fast speed, and the contact point between the flexible instrument and the natural cavity of the human body moves a large distance within one frame, and the position of the contact point changes greatly, which will bring a large error in the prior art. Since the prior art assumes that the contact point between the flexible instrument and the natural cavity of the human body is fixed, when determining the displacement of the flexible instrument in the prior art, there is a large restriction on the relative movement speed between the flexible instrument and the natural cavity of the human body.
  • the embodiment of the present application establishes an inequality constraint on the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the human body's natural cavity. This only constrains the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the human body's natural cavity to be impenetrable, while the contact position and force position between the flexible device and the human body's natural cavity can be at any position within the domain composed of multiple patches of the human body's natural cavity model. This makes the displacement and force results of the flexible device determined in Figure 2a more accurate than the displacement and force results of the flexible device determined in the prior art in Figure 2b.
  • the displacement of the flexible device is determined based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patch in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy. It does not rely on the artificially given contact force of the flexible device, nor does it rely on the collision detection method to determine the contact position between the flexible device and the natural cavity of the human body.
  • the method of the embodiment of the present application will neither introduce errors from different collision detection methods, nor assume that the contact point position between the flexible device and the natural cavity of the human body remains unchanged between two frames, and will not limit the relative movement speed of the flexible device and the natural cavity of the human body when in contact (relative to the patch scale of the natural cavity model of the human body), which greatly improves the accuracy of the determined displacement and force results of the flexible device. That is, the embodiment of the present application ensures the non-inter-invasion between the finite element model of the flexible device and the natural cavity of the human body through the inequality constraint of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the natural cavity of the human body.
  • the contact and collision positions of the flexible device and the natural cavity of the human body can be directly determined by the method of the embodiment of the present application, and the contact and collision positions of the flexible device and the natural cavity of the human body are changed between different frames, so that the determined displacement and force results of the flexible device are more accurate.
  • the method of the above embodiment is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate; then, the surgical personnel can accurately operate the flexible instrument based on the determined displacement of the flexible instrument and the operation instruction information of the flexible instrument, thereby improving the efficiency and safety of the operation.
  • determining the force result of the flexible device according to the displacement of the flexible device includes:
  • the force result of the flexible device is determined according to the displacement of the flexible device and the stiffness matrix of the flexible device.
  • the force result of the flexible device can be determined according to the determined displacement ⁇ xi of the flexible device.
  • fi represents the force result of the flexible instrument.
  • the force result of the flexible instrument is composed and determined based on the multi-dimensional finite element node force result.
  • the force result of the flexible instrument includes the external force or the gravity of the flexible instrument when it contacts the natural cavity tissue of the human body;
  • fi -1 represents the force result of the flexible instrument at moment i-1, ⁇ xi represents the displacement of the flexible instrument, and
  • K( xi-1 ) represents the stiffness matrix of the flexible instrument.
  • the stiffness matrix of the flexible device is determined based on a rotation matrix between a local coordinate system of an element in a finite element model of the flexible device and a global coordinate system.
  • the stiffness matrix of the flexible device is a semi-positive definite band symmetric matrix, which is related to the physical properties of the flexible body and the deformation size under large deformation conditions.
  • the efficiency of determining the displacement and force results of the flexible device depends on the dimension and processing efficiency of the flexible device stiffness matrix K( xi-1 ).
  • the stiffness matrix of the flexible device is:
  • Te is the rotation matrix between the unit coordinate system and the global coordinate system, which is related to the displacement x of the flexible device.
  • the value of Te needs to be recalculated in each frame calculation.
  • is a stiffness matrix assembly method, which assembles unit stiffness matrices into an overall stiffness matrix.
  • the displacement of the flexible instrument is determined based on the minimum potential energy principle by determining the corresponding displacement under the condition of the minimum potential energy change of the flexible instrument, and the stiffness matrix of the finite element model of the flexible instrument is re-determined between each frame.
  • the deformation and force of the flexible instrument under large deformation conditions can be accurately simulated, and the simulation of large deformation flexible instruments can be realized, which can provide accurate prompt information of the displacement and force results of the flexible instrument for surgery, thereby improving the efficiency and safety of the surgery.
  • the method of the above embodiment determines the displacement of the flexible instrument based on the principle of minimum potential energy by determining the corresponding displacement under the condition of minimum potential energy change of the flexible instrument, and redetermines the stiffness matrix of the finite element model of the flexible instrument between each frame, thereby being able to accurately simulate the deformation and force of the flexible instrument under large deformation conditions, and can provide doctors with accurate displacement and force results of the flexible instrument, thereby improving the efficiency and safety of the operation.
  • the method further includes:
  • the natural cavity of the human body is a narrow bifurcated cavity
  • the target constraint condition is updated.
  • the natural cavity of the human body is a narrow and bifurcated cavity
  • the finite element nodes in the finite element model of the flexible instrument due to the narrowness and bifurcation of the natural cavity of the human body, there is no penetration between the finite element nodes in the finite element model of the flexible instrument and the contact point target facets in the three-dimensional model of the natural cavity of the human body.
  • the position constraint conditions between the flexible instrument and the natural cavity of the human body need to be updated; that is, according to the calculated each time
  • the node penetration surface situation determines the inequality constraints of the positional relationship between the finite element nodes in the finite element model of the new flexible instrument and the facets in the three-dimensional model of the human body's natural cavity, until the finite element nodes of the flexible instrument converge with all the facets of the human body's natural cavity model when there is no more penetration, so that the determined displacement of the flexible instrument meets the actual situation of the operation, and the determined displacement and force results of the flexible instrument are more accurate, which can provide doctors with accurate flexible instrument displacement, force results and flexible instrument operation prompt information, which can ensure that the flexible instrument can pass through the narrow cavity of the human body better, is suitable for natural cavity surgeries such as bronchial surgeries, and improve
  • the operation prompt method of the flexible instrument as shown in Figure 3 first establishes a finite element model of the flexible instrument, and then based on the principle of minimum potential energy and the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body, the displacement of the flexible instrument can be determined; after determining the displacement of the flexible instrument, the force result of the flexible instrument can be determined based on the stiffness matrix of the flexible instrument, and the determined displacement and force result of the flexible instrument can be displayed, and the operation instruction information of the flexible instrument can be output.
  • the inequality constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the facets in the three-dimensional model of the natural cavity of the human body are updated until they converge when there is no more penetration between the finite element nodes of the flexible instrument and all the facets in the three-dimensional model of the natural cavity of the human body; by updating the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the facets in the three-dimensional model of the natural cavity of the human body, the displacement and force results of the flexible instrument determined according to the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the facets in the three-dimensional model of the natural cavity of the natural
  • the inequality constraint condition for the positional relationship between the finite element nodes in the finite element model of the updated flexible instrument and the facets in the three-dimensional model of the human body's natural cavity in the black arrow dotted line in Figure 3 is not a necessary step, but an optimization step in determining the displacement and force results of the flexible instrument in a narrow and bifurcated natural cavity of the human body.
  • the method in the embodiment of the present application can more accurately determine the displacement and force results of the flexible instrument in a narrow and bifurcated natural cavity of the human body and the flexible instrument operation prompt information, which can ensure that the flexible instrument can pass through the narrow cavity of the human body better and improve the efficiency and safety of the operation.
  • the method of the above embodiment updates the constraint conditions of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the facets in the three-dimensional model of the natural cavity of the human body, thereby more accurately determining the displacement and force results of the flexible instrument in the narrow and bifurcated natural cavity of the human body and the flexible instrument operation prompt information, thereby ensuring that the flexible instrument can better pass through the narrow cavity of the human body and improving the efficiency and safety of the operation.
  • the operation prompt device of the flexible instrument provided in the present application is described below.
  • the operation prompt device of the flexible instrument described below and the operation prompt method of the flexible instrument described above can be referred to each other.
  • FIG4 is a schematic diagram of the structure of the operation prompt device of the flexible device provided in the present application.
  • the operation prompt device of the flexible device provided in this embodiment includes:
  • Establishing module 710 for establishing a finite element model of the flexible device;
  • the finite element model of the flexible device includes a plurality of finite element nodes;
  • the first determination module 720 is used to determine the displacement of the flexible device based on the minimum potential energy principle and the target constraint condition; the target constraint condition is used to represent the constraint of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target surface patch in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible device includes the displacement of each finite element node;
  • a second determination module 730 is used to determine the force result of the flexible device according to the displacement of the flexible device
  • the prompt module 740 is used to output the operation instruction information of the flexible instrument according to the determined displacement of the flexible instrument and the force result of the flexible instrument; the operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
  • the first determination module 720 is specifically used to determine the displacement of the flexible device when the potential energy change of the flexible device is minimized and the positional relationship between the flexible device and the natural cavity of the human body meets the target constraint condition;
  • the target constraint condition is that the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the natural cavity of the human body are not penetrable.
  • the second determination module 730 is specifically used to determine the force result of the flexible device according to the displacement of the flexible device, including:
  • the force result of the flexible device is determined according to the displacement of the flexible device and the stiffness matrix of the flexible device.
  • the second determination module 730 is further used to determine the stiffness matrix of the flexible device based on a rotation matrix between a unit local coordinate system and a global coordinate system in a finite element model of the flexible device.
  • the second determination module 730 is further used to: determine the positional relationship between the finite element nodes in the finite element model of the flexible device and each surface patch in the three-dimensional model of the natural cavity of the human body based on the determined displacement of the flexible device when the natural cavity of the human body is a narrow bifurcated cavity;
  • the target constraint condition is updated.
  • the device of the embodiment of the present application is used to execute the method in any of the aforementioned method embodiments. Its implementation principle and technical effects are similar and will not be repeated here.
  • FIG5 illustrates a schematic diagram of the physical structure of an electronic device, which may include: a processor 810, a communications interface 820, a memory 830 and a communications bus 840, wherein the processor 810, the communications interface 820 and the memory 830 communicate with each other through the communications bus 840.
  • the processor 810 may call the logic instructions in the memory 830 to execute the operation prompt method of the flexible instrument, which method includes: establishing a finite element model of the flexible instrument; the finite element model of the flexible instrument includes a plurality of finite element nodes; based on the minimum potential energy principle and the target constraint condition, determining the displacement of the flexible instrument; the target constraint condition is used to represent the constraint of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patch in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument includes the displacement of each finite element node; according to the displacement of the flexible instrument, determining the force result of the flexible instrument; according to the determined displacement of the flexible instrument and the force result of the flexible instrument, outputting the operation instruction information of the flexible instrument; The operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
  • the logic instructions in the above-mentioned memory 830 can be implemented in the form of software functional units and can be stored in a computer-readable storage medium when sold or used as an independent product.
  • the technical solution of the present application can be essentially or partly embodied in the form of a software product that contributes to the prior art.
  • the computer software product is stored in a storage medium, including several instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present application.
  • the aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc.
  • the present application also provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium, and the computer program includes program instructions.
  • the computer can execute the operation prompt method of the flexible instrument provided by the above-mentioned methods, and the method includes: establishing a finite element model of the flexible instrument; the finite element model of the flexible instrument includes multiple finite element nodes; based on the minimum potential energy principle and target constraint conditions, determining the displacement of the flexible instrument; the target constraint conditions are used to represent the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument includes the displacement of each finite element node; according to the displacement of the flexible instrument, the force result of the flexible instrument is determined; according to the determined displacement of the flexible instrument and the force result of the flexible instrument, the operation instruction information of the flexible instrument is output
  • the present application also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to execute the above-mentioned flexible instrument operation prompt method, the method comprising: establishing a finite element model of the flexible instrument; the finite element model of the flexible instrument comprises a plurality of finite element nodes; based on the minimum potential energy principle and target constraint conditions, determining the displacement of the flexible instrument; the target constraint conditions are used to represent the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument comprises the displacement of each finite element node; based on the displacement of the flexible instrument, determining the force result of the flexible instrument; based on the determined displacement of the flexible instrument and the force result of the flexible instrument, outputting the operation instruction information of the flexible instrument; the operation instruction information of the flexible instrument is used to indicate the method of
  • the device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Those of ordinary skill in the art may understand and implement it without creative work.
  • each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware.
  • the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM/RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Robotics (AREA)
  • Instructional Devices (AREA)

Abstract

Embodiments of the present application provide an operation prompting method and apparatus for a flexible instrument, a device, and a storage medium. The method comprises: constructing a finite element model of a flexible instrument, wherein the finite element model of the flexible instrument comprises a plurality of finite element nodes; determining the displacement of the flexible instrument on the basis of the principle of minimum potential energy and a target constraint condition, wherein the target constraint condition is used for representing the constraint of a positional relationship between the finite element nodes in the finite element model of the flexible instrument and contact target patches in a three-dimensional model of the natural orifice of the human body, and the displacement of the flexible instrument comprises the displacement of each finite element node; determining a stress result of the flexible instrument according to the displacement of the flexible instrument; and outputting operation indication information of the flexible instrument according to the determined displacement of the flexible instrument and stress result of the flexible instrument, wherein the operation indication information of the flexible instrument is used for indicating a mode of performing surgery on the natural orifice of the human body by means of the flexible instrument. The method of the embodiments of the present application can effectively improve the efficiency and safety of surgery.

Description

柔性器械的操作提示方法、装置、设备和存储介质Method, device, equipment and storage medium for operation prompting of flexible instrument
相关申请的交叉引用CROSS-REFERENCE TO RELATED APPLICATIONS
本申请要求于2022年10月13日提交的申请号为202211255852.6,发明名称为“柔性器械的操作提示方法、装置、设备和存储介质”的中国专利申请的优先权,其通过引用方式全部并入本文。This application claims priority to Chinese patent application No. 202211255852.6, filed on October 13, 2022, and entitled “Operation prompt method, device, equipment and storage medium for flexible instruments”, which is incorporated herein by reference in its entirety.
技术领域Technical Field
本申请涉及医疗器械技术领域,尤其涉及一种柔性器械的操作提示方法、装置、设备和存储介质。The present application relates to the technical field of medical devices, and in particular to an operation prompt method, device, equipment and storage medium for a flexible device.
背景技术Background technique
在微创手术过程中,患者体内医疗器械的空间位置信息对于医疗器械的准确操作起着至关重要的作用。但是很多医疗器械在进入体内后,会不可避免的产生形状变化,因此如何准确的确定出柔性器械的受力和变形情况,为手术提供提供准确的信息,具有重要的意义。During minimally invasive surgery, the spatial position information of medical devices in the patient's body plays a vital role in the accurate operation of medical devices. However, many medical devices will inevitably change their shape after entering the body, so how to accurately determine the force and deformation of flexible devices and provide accurate information for surgery is of great significance.
相关技术中,确定柔性器械的受力和变形情况时,需要依赖于人为给定柔性器械与人体自然腔道的接触受力或依赖于碰撞检测方法确定柔性器械与人体自然腔道接触位置并且在迭代过程中柔性器械与人体自然腔道的接触点不发生改变,导致确定出的柔性器械的受力和位移的误差较大,增加了病人的手术风险。In the related art, when determining the force and deformation of a flexible instrument, it is necessary to rely on artificially given contact force between the flexible instrument and the natural cavity of the human body or rely on a collision detection method to determine the contact position of the flexible instrument and the natural cavity of the human body, and the contact point between the flexible instrument and the natural cavity of the human body does not change during the iteration process, resulting in large errors in the determined force and displacement of the flexible instrument, increasing the surgical risk for the patient.
发明内容Summary of the invention
针对现有技术中的问题,本申请实施例提供一种柔性器械的操作提示方法、装置、设备和存储介质。In view of the problems in the prior art, the embodiments of the present application provide an operation prompt method, device, equipment and storage medium for a flexible device.
具体地,本申请实施例提供了以下技术方案:Specifically, the embodiments of the present application provide the following technical solutions:
第一方面,本申请实施例提供了一种柔性器械的操作提示方法,包括:In a first aspect, an embodiment of the present application provides an operation prompt method for a flexible device, comprising:
建立柔性器械的有限元模型;所述柔性器械的有限元模型包括多个有限元节点;Establishing a finite element model of a flexible device; the finite element model of the flexible device includes a plurality of finite element nodes;
基于最小势能原理和目标约束条件,确定柔性器械的位移;所述目标约束条件用于表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束;所述柔性器械的位移包括各个有限元节点位移;Based on the minimum potential energy principle and target constraint conditions, the displacement of the flexible instrument is determined; the target constraint conditions are used to represent the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument includes the displacement of each finite element node;
根据柔性器械的位移,确定柔性器械的受力结果;Determine the force result of the flexible device according to the displacement of the flexible device;
根据确定出的所述柔性器械的位移和所述柔性器械的受力结果,输出柔性器械的操作指示信息;所述柔性器械的操作指示信息用于指示通过所述柔性器械对人体自然腔道进行手术的方式。According to the determined displacement of the flexible instrument and the force result of the flexible instrument, operation instruction information of the flexible instrument is output; the operation instruction information of the flexible instrument is used to indicate the way of performing surgery on the natural cavity of the human body through the flexible instrument.
进一步地,所述基于最小势能原理和目标约束条件,确定柔性器械的位移,包括:Furthermore, the determination of the displacement of the flexible device based on the minimum potential energy principle and the target constraint condition includes:
在柔性器械的势能变化最小且所述柔性器械与人体自然腔道之间的位置关系满足目标约束条件的情况下,确定所述柔性器械的位移;所述目标约束条件为所述柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间不可穿透。The displacement of the flexible device is determined when the potential energy change of the flexible device is minimal and the positional relationship between the flexible device and the natural cavity of the human body meets the target constraint condition; the target constraint condition is that the finite element nodes in the finite element model of the flexible device and the contact point target surfaces in the three-dimensional model of the natural cavity of the human body are not penetrable.
进一步地,所述基于最小势能原理和目标约束条件,确定柔性器械的位移,包括:Furthermore, the determination of the displacement of the flexible device based on the minimum potential energy principle and the target constraint condition includes:
利用如下公式确定柔性器械的位移:

or
Axi=bi
The displacement of the flexible device is determined using the following formula:

or
Axibi
所述表示柔性器械的势能最小变化量;所述Δxi表示柔性器械i-1时刻至i时刻的位移;所述xi-1为i-1时刻有限元节点位移;所述K(xi-1)表示柔性器械在位移xi-1情况下的刚度矩阵;所述F表示所述柔性器械所受目标外力,所述目标外力不包括柔性器械与人体自然腔道之间的接触力;所述or表示目标约束条件;所述表示柔性器械与人体自然腔道的接触点面片法向量;所述xi表示i时刻有限元节点位移;所述Si表示未形变状态下柔性器械的广义坐标位置;所述X0表示人体自然腔道的三维模型接触点面片上任意位置;所述Axi=bi表示第一约束条件;所述A表示投影矩阵;所述bi表示柔性器械的目标位置;所述xi表示i时刻有限元节点位移。 Said represents the minimum change in potential energy of the flexible device; Δxi represents the displacement of the flexible device from time i-1 to time i; xi -1 is the displacement of the finite element node at time i-1; K( xi-1 ) represents the stiffness matrix of the flexible device under the displacement xi-1 ; F represents the target external force on the flexible device, and the target external force does not include the contact force between the flexible device and the natural cavity of the human body; or Represents the target constraint; represents the normal vector of the contact point patch between the flexible instrument and the natural cavity of the human body; the xi represents the finite element node displacement at time i; the Si represents the generalized coordinate position of the flexible instrument in the undeformed state; the X0 represents an arbitrary position on the contact point patch of the three-dimensional model of the natural cavity of the human body; the Axi = bi represents the first constraint condition; the A represents the projection matrix; the bi represents the target position of the flexible instrument; and the xi represents the finite element node displacement at time i.
进一步地,所述根据柔性器械的位移,确定柔性器械受力结果,包括:Furthermore, determining the force result of the flexible device according to the displacement of the flexible device includes:
根据所述柔性器械位移和所述柔性器械的刚度矩阵,确定柔性器械的受力结果。The force result of the flexible instrument is determined according to the displacement of the flexible instrument and the stiffness matrix of the flexible instrument.
进一步地,所述方法还包括:Furthermore, the method further comprises:
基于柔性器械的有限元模型中的单元局部坐标系与全局坐标系之间的旋转矩阵,确定柔性器械的刚度矩阵。The stiffness matrix of the flexible device is determined based on a rotation matrix between a local coordinate system of a unit and a global coordinate system in a finite element model of the flexible device.
进一步地,确定柔性器械的位移之后,还包括:Further, after determining the displacement of the flexible device, the method further includes:
在所述人体自然腔道为狭窄分叉腔道的情况下,基于确定出的所述柔性器械的位移,确定柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的各个面片之间的位置关系;In the case where the natural cavity of the human body is a narrow bifurcated cavity, based on the determined displacement of the flexible device, determining the positional relationship between the finite element nodes in the finite element model of the flexible device and each surface patch in the three-dimensional model of the natural cavity of the human body;
在柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的多个面片之间的至少一个面片存在穿透的情况下,对所述目标约束条件进行更新。In the case where there is penetration between a finite element node in the finite element model of the flexible device and at least one of the plurality of facets in the three-dimensional model of the natural cavity of the human body, the target constraint condition is updated.
第二方面,本申请实施例还提供了一种柔性器械的操作提示装置,包括:In a second aspect, the embodiment of the present application further provides an operation prompt device for a flexible device, comprising:
建立模块,用于建立柔性器械的有限元模型;所述柔性器械的有限元模型包括多个有限元节点;An establishment module is used to establish a finite element model of a flexible device; the finite element model of the flexible device includes a plurality of finite element nodes;
第一确定模块,用于基于最小势能原理和目标约束条件,确定柔性器械的位移;所述目标约束条件用于表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束;所述柔性器械的位移包括各个有限元节点位移;The first determination module is used to determine the displacement of the flexible instrument based on the minimum potential energy principle and the target constraint condition; the target constraint condition is used to represent the constraint of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target surface patch in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument includes the displacement of each finite element node;
第二确定模块,用于根据柔性器械的位移,确定柔性器械的受力结果;A second determination module is used to determine the force result of the flexible device according to the displacement of the flexible device;
提示模块,用于根据确定出的所述柔性器械的位移和所述柔性器械的受力结果,输出柔性器械的操作指示信息;所述柔性器械的操作指示信息用于指示通过所述柔性器械对人体自然腔道进行手术的方式。The prompt module is used to output operation instruction information of the flexible instrument according to the determined displacement of the flexible instrument and the force result of the flexible instrument; the operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
第三方面,本申请实施例还提供了一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述程序时实现如第一方面所述柔性器械的操作提示方法。In a third aspect, an embodiment of the present application further provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the operation prompt method for the flexible device as described in the first aspect is implemented.
第四方面,本申请实施例还提供了一种非暂态计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如第一方面所述柔性器械的操作提示方法。In a fourth aspect, an embodiment of the present application further provides a non-transitory computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the operation prompt method of the flexible instrument as described in the first aspect is implemented.
第五方面,本申请实施例还提供了一种计算机程序产品,包括计算机程序,所述计算机程序被处理器执行时实现如第一方面所述柔性器械的操作提示方法。In a fifth aspect, an embodiment of the present application further provides a computer program product, including a computer program, which, when executed by a processor, implements the operation prompt method of the flexible instrument as described in the first aspect.
本申请实施例提供的柔性器械的操作提示方法、装置、设备和存储介质,基于柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束及最小势能原理,确定出柔性器械最小势能变化量情况下对应的位移就是柔性器械的位移,从而使得确定出的柔性器械的位移满足人体腔道高度受限的环境,也就使得确定出的柔性器械的位移更加的准确;进而手术人员基于确定出的柔性器械的位移和柔性器械的操作指示信息也就可以准确的进行柔性器械的操作,提升手术的效率和安全性。The operation prompt method, device, equipment and storage medium of the flexible instrument provided in the embodiments of the present application are based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determine that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate; then, based on the determined displacement of the flexible instrument and the operation instruction information of the flexible instrument, the surgical personnel can accurately operate the flexible instrument, thereby improving the efficiency and safety of the operation.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
为了更清楚地说明本申请或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the present application or the prior art, a brief introduction will be given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.
图1是本申请实施例提供的柔性器械的操作提示方法的流程示意图;FIG1 is a schematic flow chart of an operation prompt method for a flexible device provided in an embodiment of the present application;
图2a-2b是本申请实施例提供的柔性器械与人体自然腔道接触的示意图;2a-2b are schematic diagrams of the flexible device provided in the embodiments of the present application in contact with the natural cavity of the human body;
图3是本申请实施例提供的另一柔性器械的操作提示方法的流程示意图;FIG3 is a schematic flow chart of another operation prompt method of a flexible device provided in an embodiment of the present application;
图4是本申请实施例提供的柔性器械的操作提示装置的结构示意图;FIG4 is a schematic diagram of the structure of an operation prompt device for a flexible device provided in an embodiment of the present application;
图5是本申请实施例提供的电子设备的结构示意图。FIG5 is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application.
具体实施方式Detailed ways
为使本申请的目的、技术方案和优点更加清楚,下面将结合本申请中的附图,对本申请中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。In order to make the purpose, technical solutions and advantages of this application clearer, the technical solutions in this application will be clearly and completely described below in conjunction with the drawings in this application. Obviously, the described embodiments are part of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.
本申请实施例的方法可以应用于基于医疗器械的手术场景中,实现柔性器械的操作提示,提升手术效率和安全性。The method of the embodiment of the present application can be applied to surgical scenarios based on medical devices to implement operation prompts for flexible devices and improve surgical efficiency and safety.
相关技术中,确定柔性器械的受力和位移情况时,需要依赖于人为给定柔性器械与人体自然腔道的接触受力或依赖于碰撞检测方法确定柔性器械与人体自然腔道接触位置并且在迭代过程中柔性器械与人体自然腔道的接触点不发生改变,导致确定出的柔性器械的受力和位移的误差较大,增加了病人的手术风险。In the related art, when determining the force and displacement of a flexible instrument, it is necessary to rely on artificially given contact force between the flexible instrument and the natural cavity of the human body or rely on a collision detection method to determine the contact position of the flexible instrument and the natural cavity of the human body, and the contact point between the flexible instrument and the natural cavity of the human body does not change during the iteration process, resulting in large errors in the determined force and displacement of the flexible instrument, increasing the surgical risk for the patient.
本申请实施例的柔性器械的操作提示方法,基于柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束及最小势能原理,确定出柔性器械最小势能变化量情况下对应的位移就是柔性器械的位移,从而使得确定出的柔性器械的位移满足人体腔道高度受限的环境,也就使得确定出的柔性器械的位移更加的准确;进而手术人员基于确定出的柔性器械的位移和柔性器械的操作指示信息也就可以准确的进行柔性器械的操作,提升手术的效率和安全性。The operation prompt method of the flexible instrument in the embodiment of the present application is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate; then, based on the determined displacement of the flexible instrument and the operation instruction information of the flexible instrument, the surgical personnel can accurately operate the flexible instrument, thereby improving the efficiency and safety of the operation.
下面结合图1-图5以具体的实施例对本申请的技术方案进行详细说明。下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例不再赘述。The technical solution of the present application is described in detail with reference to specific embodiments in conjunction with Figures 1 to 5. The following specific embodiments may be combined with each other, and the same or similar concepts or processes may not be described in detail in some embodiments.
图1是本申请实施例提供的柔性器械的操作提示方法一实施例的流程示意图。如图1所示,本实施例提供的方法,包括:FIG1 is a flow chart of an embodiment of an operation prompt method for a flexible device provided in an embodiment of the present application. As shown in FIG1 , the method provided in this embodiment includes:
步骤101、建立柔性器械的有限元模型;柔性器械的有限元模型包括多个有限元节点;Step 101, establishing a finite element model of a flexible device; the finite element model of the flexible device includes a plurality of finite element nodes;
具体地,在微创手术过程中,患者体内医疗器械的空间信息对于医疗器械的准确操作起着至关重要的作用。但是柔性器械在进入体内后,会不可避免的产生形状变化,因此需要准确的确定出柔性器械的受力和位移情况,为医生提供提供准确的信息,以提高手术的效率和手术的安全性。Specifically, during minimally invasive surgery, the spatial information of medical devices in the patient's body plays a vital role in the accurate operation of medical devices. However, after entering the body, the flexible device will inevitably change its shape, so it is necessary to accurately determine the force and displacement of the flexible device to provide doctors with accurate information to improve the efficiency and safety of the operation.
本申请实施例中,首先建立柔性器械的有限元模型,可选地,柔性器械的有限元模型是运用有限元分析方法建立的柔性器械的模型,是一组仅在节点处连接、仅靠节点传力、仅在节点处受约束的单元组合体,柔性器械的有限元模型中包括多个有限元单元和有限元节点。建立柔性器械的有限元模型就可以把柔性器械对应的连续几何机构离散成有限个柔性器械有限元单元,并在每一个柔性器械有限元单元中设定有限个柔性器械有限元节点,从而将柔性器械对应的连续体看作仅在柔性器械有限元节点处相连接的一组柔性器械有限元单元的集合体。In the embodiment of the present application, a finite element model of a flexible device is first established. Optionally, the finite element model of a flexible device is a model of a flexible device established using a finite element analysis method, and is a group of unit combinations that are connected only at nodes, transmit forces only at nodes, and are constrained only at nodes. The finite element model of a flexible device includes multiple finite element units and finite element nodes. By establishing a finite element model of a flexible device, the continuous geometric mechanism corresponding to the flexible device can be discretized into a finite number of flexible device finite element units, and a finite number of flexible device finite element nodes can be set in each flexible device finite element unit, so that the continuum corresponding to the flexible device is regarded as a group of flexible device finite element units that are connected only at the finite element nodes of the flexible device.
在建立柔性器械的有限元模型后,就可以在确定柔性器械的位移和受力时,将柔性器械看成由许多柔性器械有限元单元和柔性器械有限元节点相互连接组成,确定柔性器械的各个子单元的位移和受力情况后,进而确定出柔性器械的位移和受力,通过大量的柔性器械有限元单元和有限元节点来模拟柔性器械的变形,从而降低确定柔性器械位移和受力结果的难度。After establishing the finite element model of the flexible device, when determining the displacement and force of the flexible device, the flexible device can be regarded as being composed of many flexible device finite element units and flexible device finite element nodes that are interconnected. After determining the displacement and force of each subunit of the flexible device, the displacement and force of the flexible device can be determined. The deformation of the flexible device can be simulated by a large number of flexible device finite element units and finite element nodes, thereby reducing the difficulty of determining the displacement and force results of the flexible device.
步骤102、基于最小势能原理和目标约束条件,确定柔性器械的位移;目标约束条件用于表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束;柔性器械的位移包括各个有限元节点位移;Step 102: Determine the displacement of the flexible device based on the minimum potential energy principle and target constraint conditions; the target constraint conditions are used to represent the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible device includes the displacement of each finite element node;
具体地,在建立柔性器械的有限元模型后,就可以基于最小势能原理和目标约束条件,确定柔性器械的位移;可选地,某一时刻i-1柔性器械的势能表示如下:
Specifically, after the finite element model of the flexible device is established, the displacement of the flexible device can be determined based on the minimum potential energy principle and the target constraint condition; optionally, the potential energy of the flexible device at a certain time i-1 is expressed as follows:
其中U为总势能,xi-1为i-1时刻节点位移广义坐标列向量,F为已知总外力;总外力F不包括柔性器械与人体自然腔道之间的相互作用力;K(xi-1)为柔性器械刚度矩阵,柔性器械下一时刻i的势能如下:

Where U is the total potential energy, xi-1 is the generalized coordinate column vector of the node displacement at time i-1, and F is the known total external force; the total external force F does not include the interaction force between the flexible device and the natural cavity of the human body; K( xi-1 ) is the stiffness matrix of the flexible device, and the potential energy of the flexible device at the next time i is as follows:

其中Δxi为i-1到i时刻柔性器械的位移,xi=xi-1+Δxi,C为与Δxi无关的常数项。定义与Δxi相关的未知势能变化ΔU如下:
Where Δxi is the displacement of the flexible device from time i-1 to time i, xi = xi-1 + Δxi , and C is a constant term independent of Δxi . The unknown potential energy change ΔU related to Δxi is defined as follows:
根据最小势能原理,利用增量法,已知Ui-1需要求Ui最小值,只需要ΔU最小即可,而柔性器械最小势能变化量情况下对应的位移就是柔性器械的位移Δxi;可选地,柔性器械的位移Δxi包括各个有限元节点位移,也就是柔性器械的位移Δxi是多维的有限元节点位移组成和确定的。另外,在确定柔性器械的位移时,需要考虑到在人体腔道这种高度受限的环境下,柔性器械的位移Δxi是受到环境约束的。即本申请实施例是基于柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束及最小势能原理,确定出柔性器械最小势能变化量情况下对应的位移就是柔性器械的位移,从而使得确定出的柔性器械的位移满足人体腔道高度受限的环境,也就使得确定出的柔性器械的位移更加的准确。According to the minimum potential energy principle, using the incremental method, it is known that U i-1 needs to find the minimum value of U i , and only ΔU needs to be minimized, and the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement Δxi of the flexible instrument; optionally, the displacement Δxi of the flexible instrument includes the displacement of each finite element node, that is, the displacement Δxi of the flexible instrument is composed and determined by the multi-dimensional finite element node displacement. In addition, when determining the displacement of the flexible instrument, it is necessary to consider that in the highly restricted environment of the human cavity, the displacement Δxi of the flexible instrument is subject to environmental constraints. That is, the embodiment of the present application is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patch in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and also makes the determined displacement of the flexible instrument more accurate.
步骤103、根据柔性器械的位移,确定柔性器械的受力结果;Step 103, determining the force result of the flexible device according to the displacement of the flexible device;
具体地,在基于最小势能原理和柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束条件,确定柔性器械的位移后,进而就可以根据确定出的柔性器械的位移确定柔性器械的受力结果。可选地,可以根据胡克定理确定柔性器械的受力结果。可选地,确定出的柔性器械的受力结果包括各个有限元节点的受力结果,也就是确定出的柔性器械的受力结果是基于多维的有限元节点受力结果组成和确定的。Specifically, after determining the displacement of the flexible device based on the principle of minimum potential energy and the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patch in the three-dimensional model of the natural cavity of the human body, the force result of the flexible device can be determined based on the determined displacement of the flexible device. Optionally, the force result of the flexible device can be determined according to Hooke's theorem. Optionally, the determined force result of the flexible device includes the force result of each finite element node, that is, the determined force result of the flexible device is composed and determined based on the multi-dimensional finite element node force result.
步骤104、根据确定出的柔性器械的位移和柔性器械的受力结果,输出柔性器械的操作指示信息;柔性器械的操作指示信息用于指示通过柔性器械对人体自然腔道进行手术的方式。Step 104: outputting operation instruction information of the flexible instrument according to the determined displacement of the flexible instrument and the force result of the flexible instrument; the operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
具体地,在基于最小势能原理和柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束条件,确定柔性器械的位移和柔性器械的受力结果后,就可以基于确定出的柔性器械的位移和受力结果,输出柔性器械的操作指示信息,其中,指示信息用于指示下一步应该如何操作柔性器械对人体人体自然腔道进行手术,医生基于柔性器械的操作指示信息也就可以准确的进行柔性器械的操作,提升手术的效率和安全性。可选地,确定出柔性器械的位移和受力结果后,也可以将柔性器械的位移和受力结果在屏幕上显示出来,以供医生进行参考,从而准确确定出下一步的手术动作,提升手术的效率和安全性。Specifically, after determining the displacement of the flexible instrument and the force results of the flexible instrument based on the minimum potential energy principle and the constraint conditions of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patch in the three-dimensional model of the natural cavity of the human body, the operation instruction information of the flexible instrument can be output based on the determined displacement and force results of the flexible instrument, wherein the instruction information is used to indicate how to operate the flexible instrument to perform surgery on the natural cavity of the human body in the next step. The doctor can also accurately operate the flexible instrument based on the operation instruction information of the flexible instrument, thereby improving the efficiency and safety of the operation. Optionally, after determining the displacement and force results of the flexible instrument, the displacement and force results of the flexible instrument can also be displayed on the screen for the doctor's reference, so as to accurately determine the next surgical action and improve the efficiency and safety of the operation.
现有的技术方案通常直接给出柔性器械的受力结果或依赖于碰撞检测算法确定柔性器械与人体自然腔道的接触位置,并且一些技术方案在迭代过程中柔性器械与人体自然腔道的接触点不发生改变,导致确定出的柔性器械的位移和受力结果出现误差,准确性和实时性不足,同时导致指示信息出现误差,影响手术的效率和安全性。本申请实施例中,通过将柔性器械的位移及受力结果的确定与柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束条件相结合,不依赖于人为给定的柔性器械与人体自然腔道之间的接触受力,也不依赖于碰撞检测算法确定接触位置,因而本申请实施例的方法既不会引入不同碰撞检测方法的误差,也不必假设在两次时间帧计算的过程中柔性器械与人体自然腔道之间的接触点位置不变,从而使得本申请实施例中确定出的柔性器械的位移和受力结果更加的准确。Existing technical solutions usually directly give the force results of the flexible instrument or rely on the collision detection algorithm to determine the contact position between the flexible instrument and the natural cavity of the human body, and in some technical solutions, the contact point between the flexible instrument and the natural cavity of the human body does not change during the iteration process, resulting in errors in the displacement and force results of the determined flexible instrument, insufficient accuracy and real-time performance, and errors in the indication information, affecting the efficiency and safety of the operation. In the embodiment of the present application, by combining the determination of the displacement and force results of the flexible instrument with the constraint conditions of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patch in the three-dimensional model of the natural cavity of the human body, it does not rely on the contact force between the flexible instrument and the natural cavity of the human body given by humans, nor does it rely on the collision detection algorithm to determine the contact position. Therefore, the method of the embodiment of the present application will neither introduce errors of different collision detection methods, nor assume that the contact point position between the flexible instrument and the natural cavity of the human body remains unchanged during the calculation of two time frames, thereby making the displacement and force results of the flexible instrument determined in the embodiment of the present application more accurate.
另一方面,相关技术中,依赖于碰撞检测算法或自身在接触空间进行简化计算(如SOFA方法),在柔性器械与人体自然腔道接触点数量增加时确定柔性器械的位移和受力结果的速度会显著减少,接触点数量极大影响其确定柔性器械的位移和受力结果的实时性和稳定性。本申请实施例中确定柔性器械位移和受力的过程中,不依赖于人为给定的柔性器械与人体自然腔道之间的接触受力,也不依赖碰撞检测算法确定柔性器械与人体自然腔道的接触位置,使得本申请实施例的方法能够在柔性器械和人体自然腔道的不同接触情况下,尤其是柔性器械和人体自然腔道的接触点较多时,实时计算柔性器械的位移和受力,提高了确定柔性器械位移和受力结果的稳定性和实用性,避免了相关技术中确定柔性器械的位移会随着柔性器械和人体自然腔道碰撞点的增加而降低效率的问题。On the other hand, in the related art, the speed of determining the displacement and force results of the flexible instrument will be significantly reduced when the number of contact points between the flexible instrument and the natural cavity of the human body increases, and the number of contact points greatly affects the real-time and stability of determining the displacement and force results of the flexible instrument. In the process of determining the displacement and force of the flexible instrument in the embodiment of the present application, it does not rely on the artificially given contact force between the flexible instrument and the natural cavity of the human body, nor does it rely on the collision detection algorithm to determine the contact position of the flexible instrument and the natural cavity of the human body, so that the method of the embodiment of the present application can calculate the displacement and force of the flexible instrument in real time under different contact conditions between the flexible instrument and the natural cavity of the human body, especially when there are many contact points between the flexible instrument and the natural cavity of the human body, thereby improving the stability and practicality of determining the displacement and force results of the flexible instrument, and avoiding the problem that the efficiency of determining the displacement of the flexible instrument in the related art will decrease with the increase of the collision points between the flexible instrument and the natural cavity of the human body.
第三方面,相关技术中,在柔性器械与人体自然腔道接触过程中相对运动速度较大时(相对于人体自然腔道模型的面片尺度),模拟确定出的柔性器械的位移和受力情况的准确性严重不足。本申请实施例对柔性器械和人体自然腔道之间接触的相对运动速度没有限制,不会限制柔性器械和人体自然腔道两物体间的相对运动速度,能够很好的模拟相对运动速度较大物体的物理运动过程,使得确定出的柔性器械的位移和受力结果更加的准确。On the third aspect, in the related art, when the relative movement speed of the flexible instrument and the natural cavity of the human body is relatively large during the contact process (relative to the patch scale of the natural cavity model of the human body), the accuracy of the displacement and force of the flexible instrument determined by simulation is seriously insufficient. The embodiment of the present application does not limit the relative movement speed of the contact between the flexible instrument and the natural cavity of the human body, and will not limit the relative movement speed between the two objects of the flexible instrument and the natural cavity of the human body. It can well simulate the physical movement process of an object with a relatively large relative movement speed, so that the displacement and force results of the flexible instrument determined are more accurate.
上述实施例的方法,基于柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束及最小势能原理,确定出柔性器械最小势能变化量情况下对应的位移就是柔性器械的位移,从而使得确定出的柔性器械的位移满足人体腔道高度受限的环境,也就使得确定出的柔性器械的位移更加的准确;进而手术人员基于确定出的柔性器械的位移和柔性器械的操作指示信息也就可以准确的进行柔性器械的操作,提升手术的效率和安全性。The method of the above embodiment is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate; then, the surgical personnel can accurately operate the flexible instrument based on the determined displacement of the flexible instrument and the operation instruction information of the flexible instrument, thereby improving the efficiency and safety of the operation.
在一实施例中,基于最小势能原理和目标约束条件,确定柔性器械的位移,包括:In one embodiment, based on the minimum potential energy principle and target constraint conditions, determining the displacement of the flexible device includes:
在柔性器械的势能变化最小且柔性器械与人体自然腔道之间的位置关系满足目标约束条件的情况下,确定柔性器械的位移;目标约束条件为柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间不可穿透。The displacement of the flexible device is determined when the potential energy change of the flexible device is minimal and the positional relationship between the flexible device and the natural cavity of the human body meets the target constraint condition; the target constraint condition is that the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the natural cavity of the human body are not penetrable.
具体地,由于人体腔道高度是受限的,因而在基于最小势能原理确定柔性器械的位移时,需要使得确定出的柔性器械的位移满足人体自然腔道的高度,即在柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间不可穿透的情况下,确定出柔性器械的位移,也就使得确定出的柔性器械的位移是在人体自然腔道内的位移,而不能将柔性器械的位移确定到人体自然腔道之外,也就使得确定出的柔性器械的位移符合实际手术情况,从而确定出的柔性器械的位移也就更加的准确。可选地,接触点目标面片为柔性器械的有限元节点接触的人体自然腔道的三维模型中距离最近的面片。Specifically, since the height of the human body cavity is limited, when determining the displacement of the flexible instrument based on the principle of minimum potential energy, it is necessary to make the displacement of the flexible instrument meet the height of the natural cavity of the human body, that is, when the finite element nodes in the finite element model of the flexible instrument and the contact point target surface in the three-dimensional model of the natural cavity of the human body are impenetrable, the displacement of the flexible instrument is determined, which means that the displacement of the flexible instrument is determined within the natural cavity of the human body, and the displacement of the flexible instrument cannot be determined outside the natural cavity of the human body, which makes the displacement of the flexible instrument consistent with the actual surgical situation, so that the displacement of the flexible instrument is more accurate. Optionally, the contact point target surface is the surface closest to the three-dimensional model of the natural cavity of the human body that the finite element node of the flexible instrument contacts.
上述实施例的方法,在基于最小势能原理确定柔性器械的位移时,通过柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间不可穿透的约束条件,使得确定出的柔性器械的位移满足人体自然腔道的高度,也就使得确定出的柔性器械的位移是在人体自然腔道内的位移,而不能将柔性器械的位移确定到人体自然腔道之外,也就使得确定出的柔性器械的位移更加的准确也更加符合实际手术情况,从而基于确定出的柔性器械的位移也就可以更加准确的确定下一步的手术方案,提高手术的效率和安全性。The method of the above embodiment, when determining the displacement of the flexible instrument based on the principle of minimum potential energy, uses the impenetrable constraint condition between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body, so that the determined displacement of the flexible instrument meets the height of the natural cavity of the human body, that is, the determined displacement of the flexible instrument is the displacement within the natural cavity of the human body, and the displacement of the flexible instrument cannot be determined outside the natural cavity of the human body, which makes the determined displacement of the flexible instrument more accurate and more in line with the actual surgical situation, so that based on the determined displacement of the flexible instrument, the next surgical plan can be determined more accurately, thereby improving the efficiency and safety of the operation.
在一实施例中,基于最小势能原理和目标约束条件,确定柔性器械的位移,包括:In one embodiment, based on the minimum potential energy principle and target constraint conditions, determining the displacement of the flexible device includes:
利用如下公式确定柔性器械的位移:

or
Axi=bi
The displacement of the flexible device is determined using the following formula:

or
Axibi
表示柔性器械的势能最小变化量;Δxi表示柔性器械i-1时刻至i时刻的位移;所述xi-1为i-1时刻有限元节点位移;所述K(xi-1)表示柔性器械在位移xi-1情况下的刚度矩阵;F表示柔性器械所受目标外力,目标外力不包括柔性器械与人体自然腔道之间的接触力; or表示目标约束条件;表示柔性器械与人体自然腔道的接触点面片法向量;xi表示i时刻有限元节点位移;Si表示未形变状态下柔性器械的广义坐标位置;X0表示人体自然腔道的三维模型对应平面上任意位置;Axi=bi表示第一约束条件;A表示投影矩阵;bi表示柔性器械的已知或可测量目标位置;xi表示i时刻有限元节点位移。 represents the minimum change in potential energy of the flexible device; Δxi represents the displacement of the flexible device from time i-1 to time i; the xi -1 is the displacement of the finite element node at time i-1; the K( xi-1 ) represents the stiffness matrix of the flexible device under the displacement xi-1 ; F represents the target external force on the flexible device, and the target external force does not include the contact force between the flexible device and the natural cavity of the human body; or represents the target constraint; represents the surface normal vector of the contact point between the flexible instrument and the natural cavity of the human body; xi represents the finite element node displacement at time i; Si represents the generalized coordinate position of the flexible instrument in the undeformed state; X0 represents any position on the plane corresponding to the three-dimensional model of the natural cavity of the human body; Axi = bi represents the first constraint condition; A represents the projection matrix; bi represents the known or measurable target position of the flexible instrument; xi represents the finite element node displacement at time i.
具体地,本申请实施例基于柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束及最小势能原理,确定出柔性器械最小势能变化量情况下对应的位移就是柔性器械的位移,从而使得确定出的柔性器械的位移满足人体腔道高度受限的环境,也就使得确定出的柔性器械的位移更加的准确,从而也就可以准确的向手术人员提供柔性器械的操作提示信息,提升手术的效率和安全性。Specifically, the embodiment of the present application is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate, so that the operation prompt information of the flexible instrument can be accurately provided to the surgical personnel, thereby improving the efficiency and safety of the operation.
可选地,相较于静力平衡方法,本申请实施例中基于最小势能原理确定柔性器械位移的优势之一是柔性器械与人体自然腔道之间的环境约束的简易表达,提高了确定出的柔性器械位移的准确性和效率,从而也就可以准确的向手术人员提供柔性器械的操作提示信息,提升手术的效率和安全性。可选地,利用如下公式确定柔性器械的位移:

or
Axi=bi
Optionally, compared to the static balance method, one of the advantages of determining the displacement of the flexible instrument based on the minimum potential energy principle in the embodiment of the present application is the simple expression of the environmental constraints between the flexible instrument and the natural cavity of the human body, which improves the accuracy and efficiency of the determined displacement of the flexible instrument, thereby accurately providing the operating prompt information of the flexible instrument to the surgical personnel, thereby improving the efficiency and safety of the operation. Optionally, the displacement of the flexible instrument is determined using the following formula:

or
Axibi
表示柔性器械的势能最小变化量;Δxi表示柔性器械i-1时刻至i时刻的位移;所述xi-1为i-1时刻有限元节点位移;所述K(xi-1)表示柔性器械在位移xi-1情况下的刚度矩阵;F表示柔性器械所受目标外力,目标外力不包括柔性器械与人体自然腔道之间的接触力; or表示目标约束条件;表示柔性器械与人体自然腔道的接触点面片法向量;xi表示i时刻有限元节点位移;Si表示未形变状态下柔性器械的广义坐标位置;X0表示人体自然腔道的三维模型对应平面上任意位置;Axi=bi表示第一约束条件;A表示投影矩阵;bi表示柔性器械的目标位置;xi表示i时刻有限元节点位移。其中,Axi=bi表示将柔性器械的已知的实际的位置映射到柔性器械的有限元模型中; 用于确定柔性器械的位移,即通过确定柔性器械最小势能变化量情况下对应的位移从而确定柔性器械的位移Δxi,柔性器械的位移Δxi是多维的有限元节点位移组成和确定的。 represents the minimum change in potential energy of the flexible device; Δxi represents the displacement of the flexible device from time i-1 to time i; the xi -1 is the displacement of the finite element node at time i-1; the K( xi-1 ) represents the stiffness matrix of the flexible device under the displacement xi-1 ; F represents the target external force on the flexible device, and the target external force does not include the contact force between the flexible device and the natural cavity of the human body; or represents the target constraint; represents the surface normal vector of the contact point between the flexible device and the natural cavity of the human body; xi represents the finite element node displacement at time i; Si represents the generalized coordinate position of the flexible device in the undeformed state; X0 represents any position on the plane corresponding to the three-dimensional model of the natural cavity of the human body; Axi = bi represents the first constraint condition; A represents the projection matrix; bi represents the target position of the flexible device; xi represents the finite element node displacement at time i. Among them, Axi = bi represents mapping the known actual position of the flexible device to the finite element model of the flexible device; It is used to determine the displacement of the flexible device, that is, to determine the displacement Δxi of the flexible device by determining the displacement corresponding to the minimum potential energy change of the flexible device. The displacement Δxi of the flexible device is composed and determined by multi-dimensional finite element node displacements.
可选地,自手术前CT三维重建而来的人体自然腔道三维模型在虚拟环境中由多个面片组成,柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的不等式约束将柔性器械与人体自然腔道的碰撞简化为点与面的碰撞,柔性器械与人体自然腔道的接触点面片法向量为平面上任一点位置为X0,有 or柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的不等式方向取决于初值大小,其中Si为未形变状态下柔性器械的广义坐标位置。可选地,or表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间不可穿透,将柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的不等式作为约束加入至确定柔性器械位移的最小势能原理中,再引入代表柔性机器人通过传感器或其他算法得到的已知柔性器械位置的等式约Axi=bi,即通过如下公式就可以准确的确定出柔性器械的位移,从而使得确定出的柔性器械的位移满足人体自然腔道高度的约束条件,也就使得确定出的柔性器械的位移更加的准确:

or
Axi=bi
Optionally, the three-dimensional model of the natural cavity of the human body reconstructed from the three-dimensional CT before surgery is composed of multiple facets in the virtual environment. The inequality constraint of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target facets in the three-dimensional model of the natural cavity of the human body simplifies the collision between the flexible instrument and the natural cavity of the human body into a collision between a point and a surface. The normal vector of the contact point facet of the flexible instrument and the natural cavity of the human body is For any point on the plane, position X 0 , we have or The direction of the inequality of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patch in the three-dimensional model of the natural cavity of the human body depends on the initial value, where S i is the generalized coordinate position of the flexible device in the undeformed state. Optionally, or The finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the human body's natural cavity are not penetrated. The inequality of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the human body's natural cavity is added as a constraint to the minimum potential energy principle for determining the displacement of the flexible device. Then, the equation representing the known position of the flexible device obtained by the flexible robot through sensors or other algorithms is introduced, approximately Axi = bi . That is, the displacement of the flexible device can be accurately determined by the following formula, so that the determined displacement of the flexible device satisfies the constraint condition of the height of the human body's natural cavity, and the determined displacement of the flexible device is more accurate:

or
Axibi
例如,如图2a和2b所示,图2a表示本申请实施例中两帧之间柔性器械的有限元模型和人体自然腔道模型之间的接触点和受力示意图,图2b表示现有技术中两帧之间柔性器械的有限元模型和人体自然腔道模型之间的接触点和受力示意图,其中n1到n6为柔性器械的有限元节点,X为人体自然腔道模型面片上的接触点,F为接触受力。由图2a和2b可见,本申请实施例在两帧间柔性器械和人体自然腔道的接触点可为人体自然腔道模型中多个面片组成的域内的任意位置,而现有技术中假设每两帧间柔性器械和人体自然腔道的接触点不变,只有接触点处受力方向和大小变化,从而导致估计误差。本申请实施例中基于柔性器械的有限元模型和最小势能原理确定柔性器械的位移,相较于现有技术的优势在于,本申请实施例中去除了柔性器械和人体自然腔道之间接触点固定的假设,在两帧之间柔性器械和人体自然腔道之间接触点可以发生变化,而现有技术方案中假设在相邻两帧柔性器械和人体自然腔道之间的接触点不变。在实际应用中,如经自然腔道的过程中,柔性器械插入自然腔道速度较快,柔性器械和人体自然腔道之间接触点在一帧内移动较大距离,接触点位置变化较大,这在现有技术中会带来较大误差,由于现有技术中假设柔性器械和人体自然腔道之间接触点固定,使得现有技术中确定柔性器械的位移时,对柔性器械和人体自然腔道之间的相对运动速度有较大限制,相对运动速度太快时会发生穿模和不稳定现象。本申请实施例通过建立柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的不等式约束,仅约束了柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间不可穿透,而柔性器械和人体自然腔道之间接触位置与受力位置可以在人体自然腔道模型的多个面片组成的域内的任意位置,从而使得图2a确定出的柔性器械的位移和受力结果相较于图2b现有技术中确定出的柔性器械的位移和受力结果更加的准确。For example, as shown in Figures 2a and 2b, Figure 2a shows the contact point and force diagram between the finite element model of the flexible instrument and the natural cavity model of the human body between two frames in the embodiment of the present application, and Figure 2b shows the contact point and force diagram between the finite element model of the flexible instrument and the natural cavity model of the human body between two frames in the prior art, where n 1 to n 6 are finite element nodes of the flexible instrument, X is the contact point on the facet of the natural cavity model of the human body, and F is the contact force. As can be seen from Figures 2a and 2b, the contact point between the flexible instrument and the natural cavity of the human body between two frames in the embodiment of the present application can be any position within the domain composed of multiple facets in the natural cavity model of the human body, while the prior art assumes that the contact point between the flexible instrument and the natural cavity of the human body between every two frames remains unchanged, and only the force direction and magnitude at the contact point change, thereby causing estimation errors. In the embodiment of the present application, the displacement of the flexible instrument is determined based on the finite element model of the flexible instrument and the principle of minimum potential energy. Compared with the prior art, the advantage is that the assumption that the contact point between the flexible instrument and the natural cavity of the human body is fixed is removed in the embodiment of the present application. The contact point between the flexible instrument and the natural cavity of the human body can change between two frames, while the prior art assumes that the contact point between the flexible instrument and the natural cavity of the human body remains unchanged in two adjacent frames. In practical applications, such as in the process of passing through the natural cavity, the flexible instrument is inserted into the natural cavity at a fast speed, and the contact point between the flexible instrument and the natural cavity of the human body moves a large distance within one frame, and the position of the contact point changes greatly, which will bring a large error in the prior art. Since the prior art assumes that the contact point between the flexible instrument and the natural cavity of the human body is fixed, when determining the displacement of the flexible instrument in the prior art, there is a large restriction on the relative movement speed between the flexible instrument and the natural cavity of the human body. When the relative movement speed is too fast, mold penetration and instability will occur. The embodiment of the present application establishes an inequality constraint on the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the human body's natural cavity. This only constrains the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the human body's natural cavity to be impenetrable, while the contact position and force position between the flexible device and the human body's natural cavity can be at any position within the domain composed of multiple patches of the human body's natural cavity model. This makes the displacement and force results of the flexible device determined in Figure 2a more accurate than the displacement and force results of the flexible device determined in the prior art in Figure 2b.
本申请实施例中基于柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束及最小势能原理确定出柔性器械的位移,不依赖于人为给定的柔性器械的接触受力,也不依赖于碰撞检测方法确定柔性器械与人体自然腔道的接触位置,使得本申请实施例的方法既不会引入不同碰撞检测方法的误差,也不必假设在两帧之间柔性器械和人体自然腔道的接触点位置不变,而且不会限制柔性器械和人体自然腔道接触时的相对运动速度(相对于人体自然腔道模型的面片尺度),极大的提升了确定出的柔性器械的位移和受力结果的准确性。也就是本申请实施例通过柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的不等式约束,保证柔性器械的有限元模型和人体自然腔道之间的不可互侵性,其不依赖于碰撞检测方法,也就是不需要事先假设柔性器械与人体自然腔道之间的接触点和碰撞位置,可以通过本申请实施例的方法直接确定柔性器械与人体自然腔道的接触和碰撞位置,且确定出的不同帧之间柔性器械与人体自然腔道的接触和碰撞位置是变化的,使得确定出的柔性器械的位移和受力结果更加的准确。In the embodiment of the present application, the displacement of the flexible device is determined based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patch in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy. It does not rely on the artificially given contact force of the flexible device, nor does it rely on the collision detection method to determine the contact position between the flexible device and the natural cavity of the human body. Therefore, the method of the embodiment of the present application will neither introduce errors from different collision detection methods, nor assume that the contact point position between the flexible device and the natural cavity of the human body remains unchanged between two frames, and will not limit the relative movement speed of the flexible device and the natural cavity of the human body when in contact (relative to the patch scale of the natural cavity model of the human body), which greatly improves the accuracy of the determined displacement and force results of the flexible device. That is, the embodiment of the present application ensures the non-inter-invasion between the finite element model of the flexible device and the natural cavity of the human body through the inequality constraint of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the natural cavity of the human body. It does not rely on the collision detection method, that is, there is no need to assume the contact points and collision positions between the flexible device and the natural cavity of the human body in advance. The contact and collision positions of the flexible device and the natural cavity of the human body can be directly determined by the method of the embodiment of the present application, and the contact and collision positions of the flexible device and the natural cavity of the human body are changed between different frames, so that the determined displacement and force results of the flexible device are more accurate.
上述实施例的方法,基于柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束及最小势能原理,确定出柔性器械最小势能变化量情况下对应的位移就是柔性器械的位移,从而使得确定出的柔性器械的位移满足人体腔道高度受限的环境,也就使得确定出的柔性器械的位移更加的准确;进而手术人员基于确定出的柔性器械的位移和柔性器械的操作指示信息也就可以准确的进行柔性器械的操作,提升手术的效率和安全性。The method of the above embodiment is based on the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body and the principle of minimum potential energy, and determines that the displacement corresponding to the minimum potential energy change of the flexible instrument is the displacement of the flexible instrument, so that the determined displacement of the flexible instrument meets the highly restricted environment of the human cavity, and the determined displacement of the flexible instrument is more accurate; then, the surgical personnel can accurately operate the flexible instrument based on the determined displacement of the flexible instrument and the operation instruction information of the flexible instrument, thereby improving the efficiency and safety of the operation.
在一实施例中,根据柔性器械的位移,确定柔性器械受力结果,包括:In one embodiment, determining the force result of the flexible device according to the displacement of the flexible device includes:
根据柔性器械位移和柔性器械的刚度矩阵,确定柔性器械的受力结果。The force result of the flexible device is determined according to the displacement of the flexible device and the stiffness matrix of the flexible device.
具体地,在建立柔性器械的有限元模型,并基于最小势能原理和目标约束条件确定柔性器械的位移Δxi后,就可以根据确定出的柔性器械的位移Δxi确定柔性器械的受力结果。可选地,根据如下公式确定确定柔性器械的受力结果:
fi=fi-1+K(xi-1)Δxi
Specifically, after establishing a finite element model of the flexible device and determining the displacement Δxi of the flexible device based on the minimum potential energy principle and the target constraint condition, the force result of the flexible device can be determined according to the determined displacement Δxi of the flexible device. Optionally, the force result of the flexible device can be determined according to the following formula:
fi =fi -1 +K(xi -1 ) Δxi
其中,fi表示柔性器械的受力结果,可选地,柔性器械的受力结果是基于多维的有限元节点受力结果组成和确定的,柔性器械的受力结果包括柔性器械与人体自然腔道组织接触所受外力或本身重力;所述fi-1表示柔性器械i-1时刻的受力结果,所述Δxi柔性器械的位移,所述K(xi-1)表示柔性器械的刚度矩阵。Wherein, fi represents the force result of the flexible instrument. Optionally, the force result of the flexible instrument is composed and determined based on the multi-dimensional finite element node force result. The force result of the flexible instrument includes the external force or the gravity of the flexible instrument when it contacts the natural cavity tissue of the human body; fi -1 represents the force result of the flexible instrument at moment i-1, Δxi represents the displacement of the flexible instrument, and K( xi-1 ) represents the stiffness matrix of the flexible instrument.
在一实施例中,基于柔性器械的有限元模型中的单元局部坐标系与全局坐标系之间的旋转矩阵,确定柔性器械的刚度矩阵。In one embodiment, the stiffness matrix of the flexible device is determined based on a rotation matrix between a local coordinate system of an element in a finite element model of the flexible device and a global coordinate system.
具体地,柔性器械的刚度矩阵为一半正定带状对称矩阵,在大形变情况下与柔性体物理性质和形变大小有关,柔性器械的位移和受力结果的确定效率取决于柔性器械刚度矩阵K(xi-1)的维度和处理效率,可选地,柔性器械的刚度矩阵为:
Specifically, the stiffness matrix of the flexible device is a semi-positive definite band symmetric matrix, which is related to the physical properties of the flexible body and the deformation size under large deformation conditions. The efficiency of determining the displacement and force results of the flexible device depends on the dimension and processing efficiency of the flexible device stiffness matrix K( xi-1 ). Optionally, the stiffness matrix of the flexible device is:
其中Te是单元坐标系与全局坐标系之间的旋转矩阵,与柔性器械的位移x有关,在每一帧计算都需要重新计算Te的值, 为单元刚度,在大变形情况下每一次计算都重新计算新的刚度矩阵。∑是一种刚度矩阵组装方法,其将单元刚度矩阵组装为总体刚度矩阵。可选地,本申请实施例中基于最小势能原理通过确定柔性器械最小势能变化量情况下对应的位移来确定柔性器械的位移,并在每一帧间重新确定柔性器械的有限元模型的刚度矩阵,通过确定柔性器械的最小势能变化情况与刚度矩阵相结合,从而能够准确模拟大变形情况下柔性器械的变形和受力,实现大变形柔性器械模拟,进而也就可以为手术提供准确的柔性器械位移和受力结果的提示信息,提升手术的效率和安全性。Where Te is the rotation matrix between the unit coordinate system and the global coordinate system, which is related to the displacement x of the flexible device. The value of Te needs to be recalculated in each frame calculation. is the unit stiffness, and a new stiffness matrix is recalculated for each calculation under large deformation conditions. ∑ is a stiffness matrix assembly method, which assembles unit stiffness matrices into an overall stiffness matrix. Optionally, in an embodiment of the present application, the displacement of the flexible instrument is determined based on the minimum potential energy principle by determining the corresponding displacement under the condition of the minimum potential energy change of the flexible instrument, and the stiffness matrix of the finite element model of the flexible instrument is re-determined between each frame. By combining the minimum potential energy change of the flexible instrument with the stiffness matrix, the deformation and force of the flexible instrument under large deformation conditions can be accurately simulated, and the simulation of large deformation flexible instruments can be realized, which can provide accurate prompt information of the displacement and force results of the flexible instrument for surgery, thereby improving the efficiency and safety of the surgery.
上述实施例的方法,基于最小势能原理通过确定柔性器械最小势能变化量情况下对应的位移来确定柔性器械的位移,并在每一帧间重新确定柔性器械的有限元模型的刚度矩阵,从而能够准确模拟大变形情况下柔性器械的变形和受力,也就可以为医生提供准确的柔性器械的位移和受力结果,提升手术的效率和安全性。The method of the above embodiment determines the displacement of the flexible instrument based on the principle of minimum potential energy by determining the corresponding displacement under the condition of minimum potential energy change of the flexible instrument, and redetermines the stiffness matrix of the finite element model of the flexible instrument between each frame, thereby being able to accurately simulate the deformation and force of the flexible instrument under large deformation conditions, and can provide doctors with accurate displacement and force results of the flexible instrument, thereby improving the efficiency and safety of the operation.
在一实施例中,确定柔性器械的位移之后,还包括:In one embodiment, after determining the displacement of the flexible device, the method further includes:
在人体自然腔道为狭窄分叉腔道的情况下,基于确定出的柔性器械的位移,确定柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的各个面片之间的位置关系;In the case where the natural cavity of the human body is a narrow bifurcated cavity, based on the determined displacement of the flexible device, determining the positional relationship between the finite element nodes in the finite element model of the flexible device and each surface patch in the three-dimensional model of the natural cavity of the human body;
在柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的多个面片之间的至少一个面片存在穿透的情况下,对目标约束条件进行更新。In the case where there is penetration between a finite element node in the finite element model of the flexible device and at least one of the plurality of facets in the three-dimensional model of the natural cavity of the human body, the target constraint condition is updated.
具体地,在人体自然腔道为狭窄分叉腔道的情况下,由于人体自然腔道狭窄且分叉,导致柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间不穿透的情况下,柔性器械的有限元模型中的有限元节点仍然与人体自然腔道的三维模型中的其他面片之间可能存在穿透情况,因而在确定柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的多个面片之间的至少一个面片存在穿透的情况下,需要对柔性器械和人体自然腔道之间的位置约束条件进行更新;也就是根据每一次计算出的节点穿透表面情况确定新的柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的面片之间的位置关系的不等式约束条件,直到柔性器械的有限元节点与人体自然腔道模型的所有面片不再有穿透发生时收敛,使得确定出的柔性器械的位移满足手术的实际情况,也就使得确定出的柔性器械的位移和受力结果更加的准确,进而也就可以为医生提供准确的柔性器械位移、受力结果及柔性器械操作提示信息,能够保证柔性器械较好的通过人体狭窄腔道,适合于经支气管等人体自然腔道手术,提升手术的效率和安全性。 Specifically, in the case where the natural cavity of the human body is a narrow and bifurcated cavity, due to the narrowness and bifurcation of the natural cavity of the human body, there is no penetration between the finite element nodes in the finite element model of the flexible instrument and the contact point target facets in the three-dimensional model of the natural cavity of the human body. However, there may still be penetration between the finite element nodes in the finite element model of the flexible instrument and other facets in the three-dimensional model of the natural cavity of the human body. Therefore, when it is determined that there is penetration between the finite element nodes in the finite element model of the flexible instrument and at least one of the multiple facets in the three-dimensional model of the natural cavity of the human body, the position constraint conditions between the flexible instrument and the natural cavity of the human body need to be updated; that is, according to the calculated each time The node penetration surface situation determines the inequality constraints of the positional relationship between the finite element nodes in the finite element model of the new flexible instrument and the facets in the three-dimensional model of the human body's natural cavity, until the finite element nodes of the flexible instrument converge with all the facets of the human body's natural cavity model when there is no more penetration, so that the determined displacement of the flexible instrument meets the actual situation of the operation, and the determined displacement and force results of the flexible instrument are more accurate, which can provide doctors with accurate flexible instrument displacement, force results and flexible instrument operation prompt information, which can ensure that the flexible instrument can pass through the narrow cavity of the human body better, is suitable for natural cavity surgeries such as bronchial surgeries, and improves the efficiency and safety of the operation.
示例性的,如图3所示的柔性器械的操作提示方法,首先建立柔性器械的有限元模型,然后基于最小势能原理和柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束,就可以确定出柔性器械的位移;在确定出柔性器械的位移后,就可以基于柔性器械的刚度矩阵确定出柔性器械的受力结果,并将确定出的柔性器械的位移和受力结果进行显示,并输出柔性器械的操作指示信息。可选地,在人体自然腔道为狭窄分叉腔道的情况下,若确定柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的多个面片之间的至少一个面片存在穿透的情况下,则对柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的面片之间的位置关系的不等式约束条件进行更新,直到柔性器械的有限元节点与人体自然腔道模型的所有面片不再有穿透发生时收敛;通过更新柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的面片之间的位置关系的约束条件,使得根据更新后的柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的面片之间的位置关系的约束条件确定出的柔性器械的位移和受力结果能够保证柔性器械较好的通过狭窄腔道,适合于经支气管等空间狭小的人体自然腔道手术模拟情况。即图3黑色箭头虚线中的更新柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的面片之间的位置关系的不等式约束条件不是必须步骤,而是在确定柔性器械在狭窄分叉的人体自然腔道的位移和受力结果的情况下的优化步骤,从而使得本申请实施例中的方法可以更加准确的确定柔性器械在狭窄分叉的人体自然腔道的位移和受力结果及柔性器械操作提示信息,能够保证柔性器械较好的通过人体狭窄腔道,提升手术的效率和安全性。Exemplarily, the operation prompt method of the flexible instrument as shown in Figure 3 first establishes a finite element model of the flexible instrument, and then based on the principle of minimum potential energy and the constraints of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body, the displacement of the flexible instrument can be determined; after determining the displacement of the flexible instrument, the force result of the flexible instrument can be determined based on the stiffness matrix of the flexible instrument, and the determined displacement and force result of the flexible instrument can be displayed, and the operation instruction information of the flexible instrument can be output. Optionally, in the case where the natural cavity of the human body is a narrow and bifurcated cavity, if it is determined that there is penetration between the finite element nodes in the finite element model of the flexible instrument and at least one of the multiple facets in the three-dimensional model of the natural cavity of the human body, the inequality constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the facets in the three-dimensional model of the natural cavity of the human body are updated until they converge when there is no more penetration between the finite element nodes of the flexible instrument and all the facets in the three-dimensional model of the natural cavity of the human body; by updating the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the facets in the three-dimensional model of the natural cavity of the human body, the displacement and force results of the flexible instrument determined according to the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the facets in the three-dimensional model of the natural cavity of the human body after the update can ensure that the flexible instrument can pass through the narrow cavity better, and is suitable for simulating natural cavity surgery in narrow spaces such as bronchial tubes. That is, the inequality constraint condition for the positional relationship between the finite element nodes in the finite element model of the updated flexible instrument and the facets in the three-dimensional model of the human body's natural cavity in the black arrow dotted line in Figure 3 is not a necessary step, but an optimization step in determining the displacement and force results of the flexible instrument in a narrow and bifurcated natural cavity of the human body. As a result, the method in the embodiment of the present application can more accurately determine the displacement and force results of the flexible instrument in a narrow and bifurcated natural cavity of the human body and the flexible instrument operation prompt information, which can ensure that the flexible instrument can pass through the narrow cavity of the human body better and improve the efficiency and safety of the operation.
上述实施例的方法,通过更新柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的面片之间的位置关系的约束条件,从而更加准确的确定柔性器械在狭窄分叉的人体自然腔道的位移和受力结果及柔性器械操作提示信息,能够保证柔性器械较好的通过人体狭窄腔道,提升手术的效率和安全性。The method of the above embodiment updates the constraint conditions of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the facets in the three-dimensional model of the natural cavity of the human body, thereby more accurately determining the displacement and force results of the flexible instrument in the narrow and bifurcated natural cavity of the human body and the flexible instrument operation prompt information, thereby ensuring that the flexible instrument can better pass through the narrow cavity of the human body and improving the efficiency and safety of the operation.
下面对本申请提供的柔性器械的操作提示装置进行描述,下文描述的柔性器械的操作提示装置与上文描述的柔性器械的操作提示方法可相互对应参照。The operation prompt device of the flexible instrument provided in the present application is described below. The operation prompt device of the flexible instrument described below and the operation prompt method of the flexible instrument described above can be referred to each other.
图4是本申请提供的柔性器械的操作提示装置的结构示意图。本实施例提供的柔性器械的操作提示装置,包括:FIG4 is a schematic diagram of the structure of the operation prompt device of the flexible device provided in the present application. The operation prompt device of the flexible device provided in this embodiment includes:
建立模块710,用于建立柔性器械的有限元模型;柔性器械的有限元模型包括多个有限元节点;Establishing module 710, for establishing a finite element model of the flexible device; the finite element model of the flexible device includes a plurality of finite element nodes;
第一确定模块720,用于基于最小势能原理和目标约束条件,确定柔性器械的位移;目标约束条件用于表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束;柔性器械的位移包括各个有限元节点位移;The first determination module 720 is used to determine the displacement of the flexible device based on the minimum potential energy principle and the target constraint condition; the target constraint condition is used to represent the constraint of the positional relationship between the finite element nodes in the finite element model of the flexible device and the contact point target surface patch in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible device includes the displacement of each finite element node;
第二确定模块730,用于根据柔性器械的位移,确定柔性器械的受力结果;A second determination module 730 is used to determine the force result of the flexible device according to the displacement of the flexible device;
提示模块740,用于根据确定出的柔性器械的位移和柔性器械的受力结果,输出柔性器械的操作指示信息;柔性器械的操作指示信息用于指示通过柔性器械对人体自然腔道进行手术的方式。The prompt module 740 is used to output the operation instruction information of the flexible instrument according to the determined displacement of the flexible instrument and the force result of the flexible instrument; the operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
可选地,所述第一确定模块720,具体用于:在柔性器械的势能变化最小且柔性器械与人体自然腔道之间的位置关系满足目标约束条件的情况下,确定柔性器械的位移;目标约束条件为柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间不可穿透。Optionally, the first determination module 720 is specifically used to determine the displacement of the flexible device when the potential energy change of the flexible device is minimized and the positional relationship between the flexible device and the natural cavity of the human body meets the target constraint condition; the target constraint condition is that the finite element nodes in the finite element model of the flexible device and the contact point target patches in the three-dimensional model of the natural cavity of the human body are not penetrable.
可选地,所述第一确定模块720,具体用于:利用如下公式确定柔性器械的位移:

or
Axi=bi
Optionally, the first determination module 720 is specifically configured to determine the displacement of the flexible device using the following formula:

or
Axibi
表示柔性器械的势能最小变化量;Δxi表示柔性器械i-1时刻至i时刻的位移;K(xi-1)表示柔性器械的刚度矩阵;xi-1为i-1时刻有限元节点位移;F表示柔性器械所受目标外力,目标外力不包括柔性器械与人体自然腔道之间的接触力;or表示目标约束条件;表示柔性器械与人体自然腔道的接触点面片法向量;xi表示时刻有限元节点位移;Si表示未形变状态下柔性器械的广义坐标位置;X0表示人体自然腔道的三维模型对应平面上任意位置;Axi=bi表示第一约束条件;A表示投影矩阵;bi表示柔性器械的目标位置;xi表示i时刻有限元节点位移。 represents the minimum change in potential energy of the flexible device; Δxi represents the displacement of the flexible device from time i-1 to time i; K( xi-1 ) represents the stiffness matrix of the flexible device; xi-1 is the displacement of the finite element node at time i-1; F represents the target external force on the flexible device, and the target external force does not include the contact force between the flexible device and the natural cavity of the human body; or represents the target constraint; represents the surface normal vector of the contact point between the flexible instrument and the natural cavity of the human body; xi represents the finite element node displacement at time i; Si represents the generalized coordinate position of the flexible instrument in the undeformed state; X0 represents any position on the plane corresponding to the three-dimensional model of the natural cavity of the human body; Axi = bi represents the first constraint condition; A represents the projection matrix; bi represents the target position of the flexible instrument; xi represents the finite element node displacement at time i.
可选地,所述第二确定模块730,具体用于:根据柔性器械的位移,确定柔性器械受力结果,包括:Optionally, the second determination module 730 is specifically used to determine the force result of the flexible device according to the displacement of the flexible device, including:
根据柔性器械位移和柔性器械的刚度矩阵,确定柔性器械的受力结果。The force result of the flexible device is determined according to the displacement of the flexible device and the stiffness matrix of the flexible device.
可选地,所述第二确定模块730,还用于:基于柔性器械的有限元模型中的单元局部坐标系与全局坐标系之间的旋转矩阵,确定柔性器械的刚度矩阵。Optionally, the second determination module 730 is further used to determine the stiffness matrix of the flexible device based on a rotation matrix between a unit local coordinate system and a global coordinate system in a finite element model of the flexible device.
可选地,所述第二确定模块730,还用于:在人体自然腔道为狭窄分叉腔道的情况下,基于确定出的柔性器械的位移,确定柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的各个面片之间的位置关系;Optionally, the second determination module 730 is further used to: determine the positional relationship between the finite element nodes in the finite element model of the flexible device and each surface patch in the three-dimensional model of the natural cavity of the human body based on the determined displacement of the flexible device when the natural cavity of the human body is a narrow bifurcated cavity;
在柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的多个面片之间的至少一个面片存在穿透的情况下,对目标约束条件进行更新。In the case where there is penetration between a finite element node in the finite element model of the flexible device and at least one of the plurality of facets in the three-dimensional model of the natural cavity of the human body, the target constraint condition is updated.
本申请实施例的装置,其用于执行前述任一方法实施例中的方法,其实现原理和技术效果类似,此次不再赘述。The device of the embodiment of the present application is used to execute the method in any of the aforementioned method embodiments. Its implementation principle and technical effects are similar and will not be repeated here.
图5示例了一种电子设备的实体结构示意图,该电子设备可以包括:处理器(processor)810、通信接口(Communications Interface)820、存储器(memory)830和通信总线840,其中,处理器810,通信接口820,存储器830通过通信总线840完成相互间的通信。处理器810可以调用存储器830中的逻辑指令,以执行柔性器械的操作提示方法,该方法包括:建立柔性器械的有限元模型;柔性器械的有限元模型包括多个有限元节点;基于最小势能原理和目标约束条件,确定柔性器械的位移;目标约束条件用于表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束;柔性器械的位移包括各个有限元节点位移;根据柔性器械的位移,确定柔性器械的受力结果;根据确定出的柔性器械的位移和柔性器械的受力结果,输出柔性器械的操作指示信息; 柔性器械的操作指示信息用于指示通过柔性器械对人体自然腔道进行手术的方式。FIG5 illustrates a schematic diagram of the physical structure of an electronic device, which may include: a processor 810, a communications interface 820, a memory 830 and a communications bus 840, wherein the processor 810, the communications interface 820 and the memory 830 communicate with each other through the communications bus 840. The processor 810 may call the logic instructions in the memory 830 to execute the operation prompt method of the flexible instrument, which method includes: establishing a finite element model of the flexible instrument; the finite element model of the flexible instrument includes a plurality of finite element nodes; based on the minimum potential energy principle and the target constraint condition, determining the displacement of the flexible instrument; the target constraint condition is used to represent the constraint of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patch in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument includes the displacement of each finite element node; according to the displacement of the flexible instrument, determining the force result of the flexible instrument; according to the determined displacement of the flexible instrument and the force result of the flexible instrument, outputting the operation instruction information of the flexible instrument; The operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
此外,上述的存储器830中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。In addition, the logic instructions in the above-mentioned memory 830 can be implemented in the form of software functional units and can be stored in a computer-readable storage medium when sold or used as an independent product. Based on this understanding, the technical solution of the present application can be essentially or partly embodied in the form of a software product that contributes to the prior art. The computer software product is stored in a storage medium, including several instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present application. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc. Various media that can store program codes.
另一方面,本申请还提供一种计算机程序产品,所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机执行时,计算机能够执行上述各方法所提供的柔性器械的操作提示方法,该方法包括:建立柔性器械的有限元模型;柔性器械的有限元模型包括多个有限元节点;基于最小势能原理和目标约束条件,确定柔性器械的位移;目标约束条件用于表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束;柔性器械的位移包括各个有限元节点位移;根据柔性器械的位移,确定柔性器械的受力结果;根据确定出的柔性器械的位移和柔性器械的受力结果,输出柔性器械的操作指示信息;柔性器械的操作指示信息用于指示通过柔性器械对人体自然腔道进行手术的方式。On the other hand, the present application also provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium, and the computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the operation prompt method of the flexible instrument provided by the above-mentioned methods, and the method includes: establishing a finite element model of the flexible instrument; the finite element model of the flexible instrument includes multiple finite element nodes; based on the minimum potential energy principle and target constraint conditions, determining the displacement of the flexible instrument; the target constraint conditions are used to represent the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument includes the displacement of each finite element node; according to the displacement of the flexible instrument, the force result of the flexible instrument is determined; according to the determined displacement of the flexible instrument and the force result of the flexible instrument, the operation instruction information of the flexible instrument is output; the operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
又一方面,本申请还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各提供的柔性器械的操作提示方法,该方法包括:建立柔性器械的有限元模型;柔性器械的有限元模型包括多个有限元节点;基于最小势能原理和目标约束条件,确定柔性器械的位移;目标约束条件用于表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束;柔性器械的位移包括各个有限元节点位移;根据柔性器械的位移,确定柔性器械的受力结果;根据确定出的柔性器械的位移和柔性器械的受力结果,输出柔性器械的操作指示信息;柔性器械的操作指示信息用于指示通过柔性器械对人体自然腔道进行手术的方式。On the other hand, the present application also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to execute the above-mentioned flexible instrument operation prompt method, the method comprising: establishing a finite element model of the flexible instrument; the finite element model of the flexible instrument comprises a plurality of finite element nodes; based on the minimum potential energy principle and target constraint conditions, determining the displacement of the flexible instrument; the target constraint conditions are used to represent the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument comprises the displacement of each finite element node; based on the displacement of the flexible instrument, determining the force result of the flexible instrument; based on the determined displacement of the flexible instrument and the force result of the flexible instrument, outputting the operation instruction information of the flexible instrument; the operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Those of ordinary skill in the art may understand and implement it without creative work.
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM/RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.
最后应说明的是:以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。 Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technologies therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims (10)

  1. 一种柔性器械的操作提示方法,包括:An operation prompt method for a flexible device, comprising:
    建立柔性器械的有限元模型;所述柔性器械的有限元模型包括多个有限元节点;Establishing a finite element model of a flexible device; the finite element model of the flexible device includes a plurality of finite element nodes;
    基于最小势能原理和目标约束条件,确定柔性器械的位移;所述目标约束条件用于表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束;所述柔性器械的位移包括各个有限元节点位移;Based on the minimum potential energy principle and target constraint conditions, the displacement of the flexible instrument is determined; the target constraint conditions are used to represent the constraints on the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target patches in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument includes the displacement of each finite element node;
    根据柔性器械的位移,确定柔性器械的受力结果;Determine the force result of the flexible device according to the displacement of the flexible device;
    根据确定出的所述柔性器械的位移和所述柔性器械的受力结果,输出柔性器械的操作指示信息;所述柔性器械的操作指示信息用于指示通过所述柔性器械对人体自然腔道进行手术的方式。According to the determined displacement of the flexible instrument and the force result of the flexible instrument, operation instruction information of the flexible instrument is output; the operation instruction information of the flexible instrument is used to indicate the way of performing surgery on the natural cavity of the human body through the flexible instrument.
  2. 根据权利要求1所述的柔性器械的操作提示方法,其中,所述基于最小势能原理和目标约束条件,确定柔性器械的位移,包括:The method for prompting an operation of a flexible device according to claim 1, wherein the step of determining the displacement of the flexible device based on the minimum potential energy principle and the target constraint condition comprises:
    在柔性器械的势能变化最小且所述柔性器械与人体自然腔道之间的位置关系满足目标约束条件的情况下,确定所述柔性器械的位移;所述目标约束条件为所述柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间不可穿透。The displacement of the flexible device is determined when the potential energy change of the flexible device is minimal and the positional relationship between the flexible device and the natural cavity of the human body meets the target constraint condition; the target constraint condition is that the finite element nodes in the finite element model of the flexible device and the contact point target surfaces in the three-dimensional model of the natural cavity of the human body are not penetrable.
  3. 根据权利要求1或2所述的柔性器械的操作提示方法,其中,所述基于最小势能原理和目标约束条件,确定柔性器械的位移,包括:The operation prompt method of the flexible instrument according to claim 1 or 2, wherein the determining the displacement of the flexible instrument based on the minimum potential energy principle and the target constraint condition comprises:
    利用如下公式确定柔性器械的位移:


    Axi=bi
    The displacement of the flexible device is determined using the following formula:


    Axibi
    所述表示柔性器械的势能最小变化量;所述Δxi表示柔性器械i-1时刻至i时刻的位移;所述K(xi-1)表示柔性器械的刚度矩阵;所述xi-1为i-1时刻有限元节点位移;所述F表示所述柔性器械所受目标外力,所述目标外力不包括柔性器械与人体自然腔道之间的接触力;所述 表示目标约束条件;所述表示柔性器械与人体自然腔道的接触点面片法向量;所述xi表示i时刻柔性器械的有限元节点位移;所述Si表示未形变状态下柔性器械的广义坐标位置;所述X0表示人体自然腔道的三维模型对应平面上任意位置;所述Axi=bi表示第一约束条件;所述A表示投影矩阵;所述bi表示柔性器械的目标位置;所述xi表示i时刻柔性器械的有限元节点位移。Said represents the minimum change in potential energy of the flexible device; Δxi represents the displacement of the flexible device from time i-1 to time i; K( xi-1 ) represents the stiffness matrix of the flexible device; xi -1 is the displacement of the finite element node at time i-1; F represents the target external force on the flexible device, and the target external force does not include the contact force between the flexible device and the natural cavity of the human body; Represents the target constraint; represents the surface normal vector of the contact point between the flexible instrument and the natural cavity of the human body; the xi represents the finite element node displacement of the flexible instrument at time i; the Si represents the generalized coordinate position of the flexible instrument in the undeformed state; the X0 represents an arbitrary position on the plane corresponding to the three-dimensional model of the natural cavity of the human body; the Axi = bi represents the first constraint condition; the A represents the projection matrix; the bi represents the target position of the flexible instrument; the xi represents the finite element node displacement of the flexible instrument at time i.
  4. 根据权利要求3所述的柔性器械的操作提示方法,其中,所述根据柔性器械的位移,确定柔性器械受力结果,包括:The operation prompt method of the flexible instrument according to claim 3, wherein the step of determining the force result of the flexible instrument according to the displacement of the flexible instrument comprises:
    根据所述柔性器械位移和所述柔性器械的刚度矩阵,确定柔性器械的受力结果。The force result of the flexible instrument is determined according to the displacement of the flexible instrument and the stiffness matrix of the flexible instrument.
  5. 根据权利要求4所述的柔性器械的操作提示方法,其中,所述方法还包括:The method for prompting an operation of a flexible device according to claim 4, wherein the method further comprises:
    基于柔性器械的有限元模型中的单元局部坐标系与全局坐标系之间的旋转矩阵,确定柔性器械的刚度矩阵。The stiffness matrix of the flexible device is determined based on a rotation matrix between a local coordinate system of a unit and a global coordinate system in a finite element model of the flexible device.
  6. 根据权利要求3所述的柔性器械的操作提示方法,其中,确定柔性器械的位移之后,还包括:The method for prompting operation of a flexible device according to claim 3, wherein after determining the displacement of the flexible device, the method further comprises:
    在所述人体自然腔道为狭窄分叉腔道的情况下,基于确定出的所述柔性器械的位移,确定柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的各个面片之间的位置关系;In the case where the natural cavity of the human body is a narrow bifurcated cavity, based on the determined displacement of the flexible device, determining the positional relationship between the finite element nodes in the finite element model of the flexible device and each surface patch in the three-dimensional model of the natural cavity of the human body;
    在柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的多个面片之间的至少一个面片存在穿透的情况下,对所述目标约束条件进行更新。In the case where there is penetration between a finite element node in the finite element model of the flexible device and at least one of the plurality of facets in the three-dimensional model of the natural cavity of the human body, the target constraint condition is updated.
  7. 一种柔性器械的操作提示装置,包括:An operation prompting device for a flexible instrument, comprising:
    建立模块,用于建立柔性器械的有限元模型;所述柔性器械的有限元模型包括多个有限元节点;An establishment module is used to establish a finite element model of a flexible device; the finite element model of the flexible device includes a plurality of finite element nodes;
    第一确定模块,用于基于最小势能原理和目标约束条件,确定柔性器械的位移;所述目标约束条件用于表示柔性器械的有限元模型中的有限元节点与人体自然腔道的三维模型中的接触点目标面片之间的位置关系的约束;所述柔性器械的位移包括各个有限元节点位移;The first determination module is used to determine the displacement of the flexible instrument based on the minimum potential energy principle and the target constraint condition; the target constraint condition is used to represent the constraint of the positional relationship between the finite element nodes in the finite element model of the flexible instrument and the contact point target surface patch in the three-dimensional model of the natural cavity of the human body; the displacement of the flexible instrument includes the displacement of each finite element node;
    第二确定模块,用于根据柔性器械的位移,确定柔性器械的受力结果;A second determination module is used to determine the force result of the flexible device according to the displacement of the flexible device;
    提示模块,用于根据确定出的所述柔性器械的位移和所述柔性器械的受力结果,输出柔性器械的操作指示信息;所述柔性器械的操作指示信息用于指示通过所述柔性器械对人体自然腔道进行手术的方式。The prompt module is used to output operation instruction information of the flexible instrument according to the determined displacement of the flexible instrument and the force result of the flexible instrument; the operation instruction information of the flexible instrument is used to indicate the method of performing surgery on the natural cavity of the human body through the flexible instrument.
  8. 一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其中,所述处理器执行所述程序时实现如权利要求1至6任一项所述的柔性器械的操作提示方法的步骤。An electronic device comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the method for prompting an operation of a flexible device as claimed in any one of claims 1 to 6 are implemented.
  9. 一种非暂态计算机可读存储介质,其上存储有计算机程序,其中,该计算机程序被处理器执行时实现如权利要求1至6任一项所述的柔性器械的操作提示方法的步骤。A non-transitory computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of the operation prompt method for a flexible device as described in any one of claims 1 to 6.
  10. 一种计算机程序产品,其上存储有可执行指令,其中,该指令被处理器执行时使处理器实现如权利要求1至6中任一项所述的柔性器械的操作提示方法的步骤。 A computer program product stores executable instructions thereon, wherein when the instructions are executed by a processor, the processor implements the steps of the operation prompt method for a flexible device as claimed in any one of claims 1 to 6.
PCT/CN2023/085904 2022-10-13 2023-04-03 Operation prompting method and apparatus for flexible instrument, device, and storage medium WO2024077889A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202211255852.6A CN115906556A (en) 2022-10-13 2022-10-13 Operation prompting method, device, equipment and storage medium of flexible instrument
CN202211255852.6 2022-10-13

Publications (1)

Publication Number Publication Date
WO2024077889A1 true WO2024077889A1 (en) 2024-04-18

Family

ID=86494059

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2023/085904 WO2024077889A1 (en) 2022-10-13 2023-04-03 Operation prompting method and apparatus for flexible instrument, device, and storage medium

Country Status (2)

Country Link
CN (1) CN115906556A (en)
WO (1) WO2024077889A1 (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006020792A2 (en) * 2004-08-10 2006-02-23 The General Hospital Corporation Methods and apparatus for simulation of endovascular and endoluminal procedures
CN103961179A (en) * 2014-04-09 2014-08-06 深圳先进技术研究院 Surgical instrument moving simulation method
CN113349928A (en) * 2021-05-20 2021-09-07 清华大学 Augmented reality surgical navigation device and method for flexible instrument

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006020792A2 (en) * 2004-08-10 2006-02-23 The General Hospital Corporation Methods and apparatus for simulation of endovascular and endoluminal procedures
CN103961179A (en) * 2014-04-09 2014-08-06 深圳先进技术研究院 Surgical instrument moving simulation method
CN113349928A (en) * 2021-05-20 2021-09-07 清华大学 Augmented reality surgical navigation device and method for flexible instrument

Also Published As

Publication number Publication date
CN115906556A (en) 2023-04-04

Similar Documents

Publication Publication Date Title
KR102016959B1 (en) Method and apparatus for generating blood vessel model
US10297066B2 (en) Animating a virtual object in a virtual world
US8395626B2 (en) Method and system for interactive simulation of materials
Székely et al. Modelling of soft tissue deformation for laparoscopic surgery simulation
JP6872556B2 (en) Population-based surface mesh reconstruction
WO2019214339A1 (en) Local coordinate system setting method for three-dimensional digital model of teeth
KR20160008980A (en) System and method of recovering lagrange multipliers in modal dynamic analysis
JP6771995B2 (en) Methods and systems for obtaining vision measurements for digital human models
JP7372996B2 (en) Node grouping method, device and electronic equipment
JP4283752B2 (en) Motion of a virtual articulated object in a virtual environment while avoiding internal collisions between articulated elements
CN112862933A (en) Method, apparatus, device and storage medium for optimizing a model
EP2887241A1 (en) Cad-based initial surface geometry correction
US20230061175A1 (en) Real-Time Simulation of Elastic Body
WO2024077889A1 (en) Operation prompting method and apparatus for flexible instrument, device, and storage medium
CN111341449A (en) Simulation method for virtual blood vessel interventional operation training
Xie et al. Finite-element kalman filter with state constraint for dynamic soft tissue modelling
Marinkovic et al. Towards real-time simulation of deformable structures by means of co-rotational finite element formulation
JP2022062274A (en) Function processing method, device, and electronic apparatus
CN115249290A (en) Spatial data processing method, spatial positioning method and equipment for unilateral temporal bone
Neubauer et al. Haptic collision handling for simulation of transnasal surgery
JP2007293540A (en) Continuous elastic body deformation simulation method, program, and recording medium
Hoogendoorn et al. An atlas‐and data‐driven approach to initializing reaction‐diffusion systems in computer cardiac electrophysiology
Erkoç et al. An observation based muscle model for simulation of facial expressions
JP7457877B2 (en) Estimation device, training device, method and program
CN115937373B (en) Avatar driving method, apparatus, device and storage medium

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 23876101

Country of ref document: EP

Kind code of ref document: A1