US20080051704A1 - Catheter and system for using same - Google Patents

Abstract

The present invention describes a system for controlling a thin flexible thermoplastic catheter. The system includes a plurality of shape memory alloy filaments attached to the distal end of the catheter, each filament having a phase and a temperature; a means for receiving a strain value for at least one of the filaments; a means for determining the phase change that will results in the strain value, whereby the phase change is dependent on the temperature and the temperature is dependent on a voltage; and a means for setting the voltage in each filament thereby resulting in movement of the catheter. In addition, a three-degrees of freedom force sensor measures the magnitude and direction of the force exerted on the tip of the catheter. The catheter can be autonomously guided or it can be interfaced to a haptic device. The catheter can also be fitted with microtools.

Description

FIELD OF THE INVENTION
• This invention relates to catheters and in particular a haptically controlled active catheter instrumented with shape memory alloy actuators and a three degrees of freedom force sensor.
BACKGROUND OF THE INVENTION
• Angioplasty is a minimally invasive procedure that involves the insertion of a catheter into the blood vessel for removal of blockages in blood flow. In the conventional approach, the catheter is inserted into the body through the femoral artery and is guided through the lumen of the blood vessel till it reaches a blockage. A stent (superelastic Shape Memory Alloy) is then deployed to open the blood vessel at the blockage and let normal blood flow resume. The surgeon is provided with images that are obtained either by X-ray imaging or by Magnetic Resonance Imaging. These images enable the healthcare provider to track the end point of the catheter (position in absolute coordinates) in real time and determine the future course of insertion.
• There are a number of problems associated with the conventional way of performing angioplasty. Specifically, the catheter insertion completely depends on the expertise and dexterity of the surgeon. In the case of intravascular neurosurgery, the catheter is pushed through extremely delicate and complex cranial blood vessels to treat aneurisms in the brain. The repeated insertion of the catheter through several trials could tear the blood vessel at the junction and cause heavy bleeding. This could also result in prolonged operating times and fatigue to clinicians and patients.
• The surgeon has no method of estimating the amount of force that is being applied by the tip of the catheter on the walls of the blood vessel. Excessive pressure could rupture the blood vessel with dire consequences. Plaque could also be dislodged which may block blood vessels in the brain or heart and cause a stroke or a myocardial infarction.
• The healthcare provider could have prolonged exposure to radiation or be subjected to a high-level of noise caused by the machinery generating magnetic fields for MRI. These pose danger or discomfort to the healthcare providers who perform the procedure over a prolonged period of time.
• Another problem with the present procedure of Angioplasty is restenosis. The deployment of a stent at the site of a blockage only provides a temporary solution for resuming blood flow. In 40% of the procedures already done, the plaque begins to build up after a couple of years, a process called restenosis. In addition, the stent, made of superelastic Shape Memory Alloy, has a life of only 10 years. For these reasons, there may arise a need to perform repeated angioplasties for a single patient.
• Accordingly it would be advantageous to provide a catheter wherein the movement of the tip can be controlled. Further it would be advantageous to provide a three degree of freedom force sensor that can be attached to a catheter.
SUMMARY OF THE INVENTION
• The present invention describes a system for controlling a thin flexible thermoplastic catheter. The system includes a plurality of shape memory alloy filaments attached to the distal end of the catheter, each filament having a phase and a temperature; a means for receiving a strain value for at least one of the filaments; a means for determining the phase change that will result in the strain value, whereby the phase change is dependent on the temperature and the temperature is dependent on a voltage; and a means for setting the voltage in each filament thereby resulting in the bending of the catheter.
• In another aspect of the invention there is provided a three dimensional force sensor adapted to be attached to the distal end of a catheter. The sensor includes a two-dimensional optical position sensing detector, a spring, an optical fibre embedded within the catheter and a means for determining the position of the beam spot on the detector. The two-dimensional optical position sensing detector is spaced from the distal end of the catheter. The spring attaches the detector to the catheter whereby a force acting on the detector moves the detector relative to the distal end of the catheter. A light beam is emitted from the end of an optical fibre that is embedded within the lumen of the catheter.
• Further features of the invention will be described or will become apparent in the course of the following detailed description.
BRIEF DESCRIPTION OF THE DRAWINGS
• The invention will now be described by way of example only, with reference to the accompanying drawings, in which:
• FIG. 1 is a schematic diagram of the autonomous control of the instrumented catheter constructed in accordance with the present invention;
• FIG. 2 is a schematic diagram of the instrumented catheter constructed in accordance with the present invention and controlled by a haptic device;
• FIG. 3A is a cross sectional view of the three-dimensional force sensor of the present invention;
• FIG. 3B is a perspective view of the three-dimensional force sensor of the present invention;
• FIG. 4 is a three dimensional rendering of the three-dimensional force sensor of the present invention;
• FIG. 5 is a cross sectional view of the active catheter showing multiple shape memory alloy sections;
• FIG. 6 is a cross sectional view of the active catheter showing patterned shape memory alloy sections;
• FIG. 7 is hysteresis characteristics of the shape memory alloy for varying sinusoidal input;
• FIG. 8 is a simulated result showing a closed loop shape memory alloy response to a step reference input using the Linear Quadratic Regulator (LQR) controller;
• FIG. 9 is a simulated result showing a closed loop shape memory alloy response to a sinusoidal reference input using the Linear Quadratic Regulator (LQR) controller;
• FIG. 10 is an experimental result showing a shape memory alloy response to a step reference input using the Linear Quadratic Regulator (LQR) controller;
• FIG. 11 is an experimental result showing a shape memory alloy response to a sinusoidal reference input using the Linear Quadratic Regulator (LQR) controller;
• FIG. 12 is an experimental result showing a closed loop shape memory alloy responses to two sinusoidal reference inputs with different DC offsets using the Linear Quadratic Regulator (LQR) controller;
• FIG. 13 is an experimental result showing a closed loop shape memory alloy responses to two sinusoidal reference inputs with different frequencies using the Linear Quadratic Regulator (LQR) controller;
• FIG. 14 is a simulation result showing a closed loop shape memory alloy responses to a step reference input using the H_∞ loop-shaping controller
• FIG. 15 is a simulation result showing a closed loop shape memory alloy response to a sinusoidal reference input using the H_∞ loop-shaping controller;
• FIG. 16 is an experimental result showing a shape memory alloy response to a step reference input using the H_∞ loop-shaping controller;
• FIG. 17 is an experimental result showing a shape memory alloy response to a sinusoidal reference input using the H_∞ loop-shaping controller;
• FIG. 18 is a block diagram of a force control for autonomous guidance for the catheter of the present invention;
• FIG. 19 is a block diagram of a master-slave configuration for the catheter of the present invention;
• FIG. 20 is a picture of the haptic device for controlling the catheter of the present invention; and
• FIG. 21 is a cross-sectional diagram of the three-dimensional force sensor of the present invention and showing a surgical tool attached thereto.
DETAILED DESCRIPTION OF THE INVENTION
• Researchers in the past have tried to instrument the catheter with Shape Memory Alloy actuators to actively control the tip of the catheter. This way the catheter can be smoothly guided into a branch in the blood vessel. The problem, however, is that the active catheters developed so far consist of actuators without any position feedback. Since most actuators used for the active catheter have a non-linear behavior, the control of the catheter by means of a joy-stick is extremely difficult and inaccurate and could cause damage to the blood vessels.
• In one embodiment, this invention describes a completely autonomous control of a catheter instrumented with Shape Memory Alloy (SMA)actuators12 and a novel three Degree Of Freedom (DOF)force sensor14 on the tip of thecatheter16, as shown inFIG. 1. Theforce sensor14 at the tip measures the magnitude and direction of the force acting on the tip of thecatheter16 due to contact with the blood vessel when used in angioplasty. The force measured by thesensor14 is fed to a control system to provide an input to the SMA actuators to minimize the force acting at the tip of the catheter. This way the force at the tip of the catheter is minimized to significantly reduce the likelihood of damage to epithelial cells of the blood vessel.
• In addition, thecatheter16 instrumented withSMA actuators12 and the 3-DOF force sensor14 can be used as a tool to perform minimally invasive surgery and therapy at remote locations. The instrumentedcatheter12 can be interfaced to ahaptic device20 in a master-slave configuration, as shown inFIG. 2. Thehaptic device20 should have a minimum of 3DOF position output24 and forceinput26. Thesurgeon22, using the haptic device, provides the required position command to actuate theactive catheter16. The force felt at the tip of thecatheter16 is reflected in 3DOFs to thehaptic device20. In this way, thesurgeon22 has access to remote sites in the blood vessel using the instrumentedcatheter16. The active catheter can also be fitted with microtools at the tip to perform on-site diagnosis, biopsies or plaque removal. In addition to application in angioplasty, the catheter can also be used for fetal and gastrointestinal surgeries.
• The sections below describes the design for the 3-DOF force sensor, he modeling of the Shape Memory Alloy Actuators, the LQR and H_∞ loop shaping controller design, the autonomous control of the catheter, and the force reflection to a haptic device.
Sensor
• Referring toFIGS. 3 and 4, thesensor14 designed in accordance with the present invention is a miniaturized 3-DOF force sensor which is attached to the tip of thecatheter16. Thesensor14 is capable of measuring the magnitude and direction of the force acting on the tip. Thesensor14 consists of aspring30 preferably a conical spring (or any suitable elastic flexible material)—one end of which is attached to anoptical fibre32 which remains constrained to thecatheter16, the other end is attached to a 2-d Optical Position Sensing Detector (2d PSD)34. Theoptical fibre32 is coupled to a source ofLASER36 or a high intensity LED. The 2d PSD, which could be a quadrant detector or a lateral effect detector, consists of four independent photodetectors, which measure the intensity of light falling on it. Theposition sensor34 provides four outputs from which the position of the beam spot on the surface of the PSD can be accurately calculated (FIG. 3 andFIG. 4). Theposition sensor34 has ahemispherical tip35 attached thereto.
• A force acting on thesensor14 due to contact with the walls of the blood vessel causes thespring30 to bend. Since theoptical fibre32 and one end of thespring30 are fixed to thecatheter16, the spring bends30 relative to thecatheter16 and as a result the position of the beam spot from theoptical fibre32 shifts along the surface of the2d PSD34. The deviation of the beam spot along the surface of the2d PSD34 can be used to calculate the amount of bending in the spring which is proportional to the force acting on the tip of the catheter, by Hooke's Law. Therefore, the magnitude and direction of the force acting on the sensor can be measured accurately. The total voltage from all four outputs is an indication of the total intensity of light falling laterally on the PSD.
• Theforce sensor14 is mounted on the tip of thecatheter16. Thesensor14 measures the force exerted by the artery walls on the tip of the catheter. This force reading can be used to autonomously guide the catheter to prevent damage to the epithelial cells of the artery or it can be reflected to a haptic device.
Actuators
• Referring toFIGS. 5 and 6, theactuators12 guide thecatheter14 through the lumen of the blood vessel without causing damage to the inner linings of the blood vessel. Theactuators12 also facilitate easy movement of thecatheter14 into one of the blood vessels at a branch. Theactuators12 could be made of Shape Memory Alloy (SMA) or Electro active polymers. SMA filaments (wire40 and tube44) are used to build theactuators12 on thecatheter14. The property of SMAs of demonstrating a change in strain by a change in temperature is due to a solid-state phase change from the Martensite form to the Austenite form. Martensite is relatively softer and easily deformable and exists at a lower temperature. Heating by means of current or hot water results in a phase transition from the Martensite phase to the Austenite phase. In the Austenite phase, the SMA is more structured and is harder and stronger than in the Martensite phase.
• The SMA actuators12 on the instrumentedcatheter16 could be either in the form of a wire, tube or sheet. Since the machinability of SMA is low, the wires orfilaments40 should be microwelded to thestainless steel pads42 which can be glued to the catheter, as illustrated inFIG. 5. A minimum of threewires40 placed at 120° is required to bend thecatheter16 in all directions. The bending can be created by a number of sections of SMA to create variable bending radius and angles. The example shown inFIG. 1 shows one section of the SMA with one set of threewires40 and the example shown inFIG. 5 shows two sections each having threewires40. SMA tubes and sheets can also be used to create the bending portion of the catheter.Patterns44 can be laser-machined on the SMA tubes or sheets and can be directly glued to the catheter, as shown inFIG. 6.
• In order to control the bending angle of the catheter, the amount of strain generated in theSMA actuators12 should be controlled. For this, the physical behavior of the SMA should be described in dynamic equations. The model of the SMA set out below is based on the laws of physics. The modeling equations are as follows:
Modelling of Phase Transformations
• Since an SMA exists in two states, it can be modelled as a two-state system. There is a similarity between this system and an electron, which can exist only in two states—the positive spin or the negative spin. The Fermi-Dirac statistics are used to describe the number of electrons in the two states depending on the energy of the electron. The state of an SMA in the Martensite and Austenite forms using the same statistics has been modeled.
• Two modelling equations were used based on whether the alloy is being heated or cooled due to hysteresis with two different transition temperatures. Since the SMA is in the Martensite form at lower temperatures, the phase transformation equation during heating is described by analogy with the Fermi-Dirac statistics in the form:
• $\begin{array}{cc}\xi =\frac{{\xi }_m}{1+\exp \left(\frac{T_{\mathrm{fa}}-T}{{\sigma }_a}+K_a\sigma \right)}& \left(1\right)\end{array}$
• where ξ is the fraction of the Austenite phase, ξ_mis the fraction of the Martensite phase prior to the present transformation from Martensite to Austenite, T is the temperature, T_fais the transition temperature from Martensite to Austenite, σ_ais an indication of the range of temperature around the transition temperature T_faduring which the phase change occurs, σ is the stress and K_ais the stress curve-fitting parameter which is obtained from the loading plateau of the stress-strain characteristic with no change in temperature.
• On cooling, the Austenite phase gets converted to the Martensite phase and the modelling equation during cooling is described by analogy with the Fermi-Dirac statistics in the form:
• $\begin{array}{cc}\xi =\frac{{\xi }_a}{1+\exp \left(\frac{T_{\mathrm{fm}}-T}{{\sigma }_m}+K_m\sigma \right)}& \left(2\right)\end{array}$
• where ξ_ais the fraction of the Austenite phase prior to the present transformation from Austenite to Martensite, T is the temperature, T_fmis the transition temperature from Austenite to Martensite, σ_mis an indication of the range of temperature around the transition temperature T_fmduring which the phase change occurs, σ is the stress and K_mis the stress curve-fitting parameter which is obtained from the unloading part of the
• Since the SMA is modelled as a two-component system, at any given time, the sum of the mole fractions of the Austenite and Martensite phase is 1, i.e.,
• 
  ξ_a+ξ_m=1   (3)
• The time derivatives of Eqns.(1) and (2) are as follows:
  For heating:
• $\begin{array}{cc}\stackrel{.}{\xi }=\frac{{\xi }^2}{{\xi }_m}\left[\exp \left(\frac{T_{\mathrm{fa}}-T}{{\sigma }_a}+K_a\sigma \right)\right]\left[\frac{\stackrel{.}{T}}{{\sigma }_a}-K_a\stackrel{.}{\sigma }\right]& \left(4\right)\end{array}$
• For cooling:
• $\begin{array}{cc}\stackrel{.}{\xi }=\frac{{\xi }^2}{{\xi }_a}\left[\exp \left(\frac{T_{\mathrm{fm}}-T}{{\sigma }_m}+K_m\sigma \right)\right]\left[\frac{\stackrel{.}{T}}{{\sigma }_m}-K_m\stackrel{.}{\sigma }\right]& \left(5\right)\end{array}$
Modelling of Temperature Dynamics
• An SMA actuator is heated by the process of Joules heating by applying a voltage across the SMA. The loss of heat from the SMA is through natural convection. Mathematically the dynamics of the temperature are given by equation (6).
• $\begin{array}{cc}\stackrel{.}{T}=\frac{1}{mc_p}\left[\frac{V^2}{R}-hA\left(T-T_a\right)\right]& \left(6\right)\end{array}$
• where m is the mass per unit length, c_pis the specific heat capacity, V is the voltage applied across the SMA, R is the resistance per unit length, h is the coefficient of convection, A=πd is the circumferential area of cooling, d is the diameter of the wire, T is the temperature and T_ais the ambient temperature. The coefficient h is assumed to have the characteristics of a second-order polynomial to enhance the rate of convection at higher temperatures as observed in open-loop results:
• 
  h=h₀+h₂T²  (7)
Constitutive Equation
• The constitutive equation relating changes in stress, strain, temperature and mole fraction is given by the following equation:
• 
  {dot over (σ)}=D{dot over (ε)}+θ_t{dot over (T)}+Ω{dot over (ξ)}  (8)
• where σ is the stress in the SMA, D is the average of the Young's Moduli for the Martensite and Austenite phases, ε is the strain, θ_tis the thermal expansion factor and Ω is the phase transformation contribution factor. The model explains the hysteresis as well as the minor loops in the hysteresis, as illustrated inFIG. 7.
• The dynamic characteristics of the SMA are completely defined by equation (4) or (5) (heating or cooling), together with equations (6) and (8). The σ_eis defined as the integral of the error, i.e.,
• 
  {dot over (σ)}_e=ε−ε_ref  (9)
• where ε is the strain of the SMA actuator and ε_refis the reference trajectory. The dynamic equations of the SMA along with equation (9) can be represented in the state-space form:
• $\begin{array}{cc}\stackrel{.}{\stackrel{\_}{z}}=f\left(\stackrel{\_}{z},u,t\right)\mathrm{where}\stackrel{\_}{z}=\left[\begin{array}{c}\epsilon \\ T\\ \xi \\ {\sigma }_e\end{array}\right]& \left(10\right)\end{array}$
• and u is the input voltage to the SMA wire. The nonlinear equations are linearized about a set of operating points (ε₀, T₀, ξ₀, u₀) on the reference trajectory.
• To obtain the operating points, ε₀is chosen as the value of the reference strain at that instant of time. T₀and ξ₀are obtained by integrating equations (4) or (5) and (8), depending on whether the SMA is being heated or cooled, for a given value of ε₀. The value of u₀is obtained from equation (6) for a given value of T₀, assuming steady-state conditions.
• Equation (10) is linearized about the calculated operating points to obtain linear models in the form:
• $\begin{array}{cc}\stackrel{.}{\stackrel{\_}{z}}=A\stackrel{\_}{z}+Bu& \left(11\right)\\ y=C\stackrel{\_}{z}\mathrm{where}A={\left[\frac{\partial f}{\partial \stackrel{\_}{z}}\right]}_{{\mathrm{<figure-callout\; label="\epsilon "\; filenames="US20080051704A1-20080228-D00008.png"\; state="\{\{state\}\}">\epsilon </figure-callout>}}_0,{\mathrm{<figure-callout\; label="T"\; filenames="US20080051704A1-20080228-D00008.png"\; state="\{\{state\}\}">T</figure-callout>}}_0,{\mathrm{<figure-callout\; label="\xi "\; filenames="US20080051704A1-20080228-D00008.png"\; state="\{\{state\}\}">\xi </figure-callout>}}_0,{\mathrm{<figure-callout\; label="u"\; filenames="US20080051704A1-20080228-D00008.png"\; state="\{\{state\}\}">u</figure-callout>}}_0}B={\left[\frac{\partial f}{\partial u}\right]}_{{\mathrm{<figure-callout\; label="\epsilon "\; filenames="US20080051704A1-20080228-D00008.png"\; state="\{\{state\}\}">\epsilon </figure-callout>}}_0,{\mathrm{<figure-callout\; label="T"\; filenames="US20080051704A1-20080228-D00008.png"\; state="\{\{state\}\}">T</figure-callout>}}_0,{\mathrm{<figure-callout\; label="\xi "\; filenames="US20080051704A1-20080228-D00008.png"\; state="\{\{state\}\}">\xi </figure-callout>}}_0,u_0}& \left(12\right)\end{array}$
• For the no-load case, which is the case considered here, σ and {dot over (σ)} are equal to zero. In this case, the model given by (11) is not controllable since the number of independently controllable states is only 2. On removing the uncontrollable modes, the following state space model is obtained:
• $\begin{array}{cc}\stackrel{.}{\stackrel{\_}{x}}=A^{\prime }\stackrel{\_}{x}+B^{\prime }u\mathrm{where}\stackrel{\_}{x}=\left[\begin{array}{c}\epsilon \\ {\sigma }_{\epsilon }\end{array}\right]& \left(13\right)\end{array}$
Controller Design
• The force measured by the sensor at the tip of the catheter is provided as input to the control strategy developed to provide an output to the SMA actuators to minimize the force at the tip of the catheter. The actuators could also be activated by the surgeon to guide the catheter into the branches of the blood vessel. The control scheme should provide fast and accurate control of the strain in the SMA. Two control strategies have been developed—a Gain Scheduled Pi controller and a robust H_∞ loop-shaping controller using normalized coprime stabilization.
Gain Scheduled controller
• For a PI controller, the feedback is of the form,
• 
  u=−K_p(ε−ε_ref)−K_Iσ_e+u₀  (14)
• where K_pis the proportional gain and K, is the integral gain. Writing K=[K_pK_I], the gains are computed such that the resulting controller minimizes the quadratic cost function,
• $\begin{array}{cc}J\left(u\right)={\int }_0^{\infty }\left(x^TQx+u^TRu\right)t& \left(15\right)\end{array}$
• The choice of gains for minimizing J(u) are obtained by first solving for S in the algebraic Riccati equation:
• 
  A′^TS+SA′−SB′R⁻¹B′^TS+Q=0   (16)
• The matrix K is given by:
• 
  K=R⁻¹B′^TS   (17)
• The simulation result for a step reference input is shown inFIG. 8 and that for a sinusoidal input is shown inFIG. 9.
• The experimental verification of the controller was done on a 700 MHz Windows based PC at a sampling rate of about 65 Hz. Since the form of matrices A′ and B′ are fairly simple, the solution for the Ricatti equation was obtained in closed form using MAPLE. This also greatly reduced the computation time by avoiding matrix decompositions like the Schur decomposition. The response of the SMA to a step reference input is shown inFIG. 10.
• From the experimental results, it can be observed that the response time for heating is approximately 1.0 second and for cooling is 2.1 seconds for a 0.012″ diameter SMA wire. The rate of heating and cooling is much higher for a thinner wire since the ratio of surface area to volume increases as the wire diameter is reduced, thereby increasing the rate of convectional cooling.
• The experimental result for a sinusoidal reference input is shown inFIG. 11 which shows an excellent tracking of the reference trajectory by the SMA in a closed loop. The DC offset of the reference was also varied to check the performance of the controller in the entire operating range. The results are shown inFIG. 12. The graphs show an excellent response over the entire operating region. The frequency for the reference input was also varied to check the performance of the controller and the SMA. The results are shown inFIG. 13. The response of the controller is satisfactory for frequencies lower than 0.1 Hz.
Robust Controller
• The LQR based PI controller shows excellent response to a step and a sinusoidal reference input. The controller also has a good stability margin, thereby making the control law fairly insensitive to uncertainties in the parameters of the model. However, the controller is sensitive to unmodelled dynamics at high frequencies and thereby could lead to instability. For this reason, an H_∞ loop-shaping controller using normalized coprime stabilization is designed such that the gains are high when the model describes the SMA accurately and low at higher frequencies when the model is inaccurate. The loop-shaping controller, however, obtains a performance/robust stability tradeoff.
• The gain matrix K=[K_pK_I]of the PI controller, as given by equation (14), is computed by the loop-shaping procedure using normalized coprime stabilization.
• The nominal plant (SMA) is first shaped by using a precompensator W₁so that the singular values of the nominal plant G are shaped to a desired open-loop shape. The compensator also ensures that the shaped plant G_s=GW₁is square and G_shas no hidden modes.
• For ensuring robust stabilization, ε_maxis calculated as,
• $\begin{array}{cc}{\epsilon }_{\max }={\left(\underset{K\mathrm{stabilizing}}{\inf }{\left[\begin{array}{c}I\\ K\end{array}\right]{\left(I+G_sK\right)}^{-1}{\tilde{M}}_s^{-1}}_{\infty }\right)}^{-1}& \left(18\right)\\ =\sqrt{1-{{\tilde{N}}_s{\tilde{M}}_s}_H^2}<1& \left(19\right)\end{array}$
• where {tilde over (M)}_sand Ñ_sdefine the normalized coprime factors of G_s. An ε=ε_maxis selected to form a stabilizing controller K_∞, which satisfies:
• $\begin{array}{cc}{\left[\begin{array}{c}I\\ K_{\infty }\end{array}\right]{\left(I+G_sK_{\infty }\right)}^{-1}{\tilde{M}}_s^{-1}}_{\infty }\le {\epsilon }^{-1}\stackrel{.}{=}\gamma & \left(20\right)\end{array}$
The final feedback controller K is obtained as the product of the precompensator W₁and the H_∞ controller K_∞ i.e.
• 
  K=W₁K_∞  (21)
• The simulation of the robust controller and the SMA was done on MATLAB™. The desired loop shape was chosen so that the gains are high at low frequencies where the model describes the SMA accurately and the gains roll off at −20 dB beyond the corner frequency, thereby ensuring a low gain at high frequencies.
• The value of the corner frequency was chosen to be 1000.424 rad.sec⁻¹. The simulation result for input consisting of step changes is shown inFIG. 14 and for a sinusoidal reference input is shown inFIG. 15.
• The experimental verification of the controller was done on a 700 MHz Windows based PC at a sampling rate of about 100 Hz. The response of the SMA to step changes is shown inFIG. 16 and for sinusoidal reference input is shown inFIG. 17.
Autonomous Control of the Active Catheter
• In the previous section, modeling and robust control of SMA actuators have been described. The simulation and experimental results show excellent tracking response for the SMA, thereby validating both the model and the control scheme. In autonomous control of anactive catheter16 fitted withSMA actuators12 and a 3-DOF force sensor14, the force acting on the tip of thecatheter16 is regulated so that there is no damage to the epithelial cells of the artery. A force control algorithm is implemented to ensure that the tip interaction forces are minimized.
• Theminiature force sensor14 at the tip of the catheter provides the magnitude and direction of the force acting on the tip. This force is provided as a negative feedback signal to the robust control system, as shown inFIG. 18. The sampling rate of the control algorithm is 100 Hz, providing fast response to minimize the error between the force signal F_eand the desired force F_d. The autonomous control of anSMA actuator12 based on the force felt at the tip prevents any damage to the epithelial cells of the arteries. In addition, the catheter autonomously guides itself in the proximity of the blockage to prevent accidental dislodging of the plaque due to application of excessive force. This ensures a higher success rate of performing angioplasty by reducing the risk of rupturing blood vessels or dislodging plaque, which can cause a stroke or myocardial infarction if the plaque lodges itself in minute arteries of the heart or the brain.
Force reflection to a Haptic device
• One of the major problems with the conventional way of performing angioplasty is restenosis. In nearly 40% of procedures, the plaque begins to rebuild at the site after a couple of years, causing an obstruction to blood flow. In addition, the superelastic stent used for reopening the blood vessel has a lifetime of only 10 years. Therefore, the patient would have to undergo multiple angioplasties during his/her life time. A solution to this problem is to perform microsurgery at the site of the blockage.
• Theactive catheter16 fitted withSMA actuators12 and a 3-DOF force sensor14 can be controlled by a masterhaptic device20, as shown inFIG. 2 andFIG. 20. Thehaptic device20 provides position commands to theSMA actuators12. A minimum of 3 DOFs in position and force are required for thehaptic device20. An inverse dynamics control scheme under image guidance is used to control the bending of the active catheter. The force acting on the tip of thecatheter16 is reflected to the stylus of thehaptic device20. The surgeon can, therefore, feel these forces caused by interactions between thecatheter16 and the artery walls. In none of the technologies developed thus far has the surgeon had the capability of feeling the forces exerted on the tip of the catheter at a remote location.Micro tools48 can be attached to the tip of thecatheter16 with 3 DOF position control and the ability to feel forces, also in 3 DOF. Using these tools, the surgeon can gently remove the plaque from the site under image guidance and force feedback, thereby leading to a permanent treatment of plaque buildup. This tool can also be used to perform various procedures like biopsies and on-site diagnosis, as shown inFIG. 21.
• In addition to the application in Angioplasty, the instrumentedcatheter16 can also be used to perform minimally invasive fetal surgery. An abnormality can arise in some fetuses where the diaphragm separating the chest cavity from the abdomen does not completely develop. As a result the intestines begin to grow in the chest cavity, preventing the development of the lungs. This abnormality is called Congenital Diaphragmatic Hernia (CDH). In such cases, the surgeon can perform “fetal endoscopic surgery” for CDH. In this procedure, a catheter along with surgical instruments is inserted through tiny ports in the mother's womb. The surgeon operates on the fetus entirely within the womb. The trachea of the fetus is clipped to block the windpipe. This causes the trapped liquid in the lungs to aid in the growth and expansion of the lungs. This is a high risk operation and requires high dexterity and maneuverability on the part of the surgeon to position the fetus and the mother accurately. In addition, it is extremely difficult to operate on a fetus floating in the amniotic fluid without any force or visual feedback. However, the procedure is safer for the mother and the fetus in comparison to open fetal surgery, where a large incision must be made in the uterus and in the amniotic sac, since it reduces the chances of pre-term labor. A premature birth would reduce the amount of time for the fetus to develop in the womb after prenatal surgery, and reduce its chances of survival.
• An instrumentedcatheter16 as described in the invention would be useful in guiding the catheter into the windpipe of the fetus under image and haptic guidance. This would reduce the time taken and the chances of injury to the fetus in inserting the instrumentedcatheter16 into the fetus, thereby making the procedure more reliable and efficient and reducing the amount of stress to the surgeon.
• The primary application of this invention is to insert a catheter instrumented16 with a 3-DOF force sensor14 and ShapeMemory Alloy actuators12 smoothly into an artery and guide it through the lumen of the blood vessel till thecatheter16 reaches a blockage. The stent is then deployed to resume normal blood flow. The 3-Degrees of Freedom (DOF)force sensor14 measures the magnitude and direction of the force being applied. It could be used for a number of other applications. Specifically it could be used at the tip of a catheter to measure the forces exerted by the walls of the blood vessel on the tip of the catheter. It could also be used as a minimally invasive tool to measure the force exerted by the tool during surgery. As well it could be used to study soft tissue properties by measuring accurately the forces exerted by the soft tissue during deformation of the tissue.
• The ShapeMemory alloy actuators12 have a number of advantages. Specifically, SMA actuators have large plastic deformations which can be recovered; can generate large load forces; have high power to weight ratio; have low driving voltages; have no moving parts; are noiseless; can be used in compact workspace; and are biocompatible and corrosion resistant.
• The modeling and accurate control of theSMA actuators12 enables its use for a number of applications. Specifically it could be used as an active catheter, biomedical forceps, tissue stapler, stent, and minimally invasive tools.
• As used herein, the terms “comprises” and “comprising” are to construed as being inclusive and opened rather than exclusive. Specifically, when used in this specification including the claims, the terms “comprises” and “comprising” and variations thereof mean that the specified features, steps or components are included. The terms are not to be interpreted to exclude the presence of other features, steps or components.
• It will be appreciated that the above description related to the invention by way of example only. Many variations on the invention will be obvious to those skilled in the art and such obvious variations are within the scope of the invention as described herein whether or not expressly described.

Claims (43)

1. A system for controlling a thin flexible thermoplastic catheter comprising:

a plurality of shape memory alloy filaments attached to the distal end of the catheter, each filament having a phase and a temperature;

a means for receiving a strain value for at least one of the filaments;

a means for determining a phase change that will result in the strain value, whereby the phase change is dependent on the temperature and the temperature is dependent on a voltage; and

a means for setting the voltage in each filament thereby resulting in bending of the catheter in three dimensional space.

2. A system as claimed inclaim 1 further including a force sensor attached to the distal end of the catheter and a means for receiving a signal from the force sensor.

3. A system as claimed inclaim 2 wherein the force sensor measures forces in three degrees of freedom.

4. A system as claimed inclaim 3 wherein the signal from the force sensor is used to determine the strain value such that the distal end of the catheter is deflected autonomously to minimize the force on the tip of the catheter when the signal from the force sensor is over a predetermined value.

5. A system as claimed inclaim 1 wherein the strain value is determined by a user.

6. A system as claimed inclaim 3 further including a haptic device operably connected to the plurality of filaments and the force sensor; wherein the haptic device transmits strain values to the plurality of filaments and receives the signal from the force sensor.

7. A system as claimed inclaim 3 wherein the plurality of filaments include three filaments equally spaced around the catheter.

8. A system as claimed inclaim 7 wherein each filament is an elongate filament attached to the catheter at each end thereof.

9. A system as claimed inclaim 8 wherein the plurality of filaments is a first set of filaments and further including a second set of filaments wherein the second set of filaments are adjacent to the first set of filaments along the catheter.

10. A system as claimed inclaim 8 wherein the plurality of filaments is a first set of filaments and further including a second set of filaments wherein the second set of filaments is spaced a predetermined space from the first set of filaments along the catheter.

11. A system as claimed inclaim 7 wherein each filament is micro-machined into a predetermined pattern and each filament is attached to the catheter along its length.

12. A system as claimed inclaim 3 wherein the sensor includes:

a two-dimensional optical position sensing detector located at a predetermined distance from the distal end of the catheter;

a spring attaching the detector to the catheter whereby force on the detector moves the detector relative to the distal end of the catheter;

a light source located at the end of the catheter is coupled to an optical fibre whereby the light source emits a beam spot; and

a means for determining the position of the beam spot on the detector.

13. A system as claimed inclaim 12 further including a hemispherical tip attached to the detector extending outwardly from the detector away from the spring.

14. A system as claimed inclaim 13 further including a means for determining the intensity of the light on the detector.

15. A system as claimed inclaim 12 wherein the two-dimensional optical position detector is a quadrant detector.

16. A system as claimed inclaim 12 wherein the two-dimensional optical position detector is a lateral-effect detector.

17. A system as claimed inclaim 13 wherein the spring is a conical spring.

18. A system as claimed inclaim 12 wherein the light source is a laser.

19. A system as claimed inclaim 12 wherein the light source is a high intensity light emitting diode.

20. A system as claimed inclaim 3 wherein the sensor further includes a surgical tool extending outwardly therefrom.

21. A system as claimed inclaim 3 wherein the catheter is used in one of heart surgery, gastrointestinal surgery and fetal surgery.

22. A system as claimed inclaim 1 wherein the temperature is determined using a formula given by

$\stackrel{.}{T}=\frac{1}{mc_p}\left[\frac{V^2}{R}-hA\left(T-T_a\right)\right]$

where m is the mass per unit length, c_pis the specific heat capacity, V is the voltage applied across the shape memory alloy filament, R is the resistance per unit length, h is the coefficient of convection, A=πd is the circumferential area of cooling, d is the diameter of the filament, T is the temperature and T_ais the ambient temperature.

23. A system as claimed inclaim 22 where the strain is determined using a formula given by

{dot over (σ)}=D{dot over (ε)}+θ_t{dot over (T)}+Ω{dot over (ξ)}

where σ is the stress in the filament, D is the average of the Young's Moduli for the Martensite and Austenite phases, ε is the strain, θ_tis the thermal expansion factor, and Ω is the phase transformation contribution factor.

24. A system as claimed inclaim 23 wherein the mole fraction of shape memory alloy filament in the Austenite phase while heating is determined using a formula given by:

$\stackrel{.}{\xi }=\frac{{\xi }^2}{{\xi }_m}\left[\exp \left(\frac{T_{f_a}-T}{{\sigma }_a}+K_a\sigma \right)\right]\left[\frac{\stackrel{.}{T}}{{\sigma }_a}-K_a\stackrel{.}{\sigma }\right]$

where ξ is the fraction of the Austenite phase, ξ_mis the fraction of the Martensite phase prior to the present transformation from Martensite to Austenite, T is the temperature, T_fais the transition temperature from Martensite to Austenite, σ_ais an indication of the range of temperature around the transition temperature T_faduring which the phase change occurs, σ is the stress and K_ais the stress curve-fitting parameter which is obtained from the loading plateau of the stress-strain characteristic with no change in temperature; and

the mole fraction of shape memory alloy filament in the Austenite phase while cooling is determined using the formula given by:

$\stackrel{.}{\xi }=\frac{{\xi }^2}{{\xi }_a}\left[\exp \left(\frac{T_{f_m}-T}{{\sigma }_m}+K_m\sigma \right)\right]\left[\frac{\stackrel{.}{T}}{{\sigma }_m}-K_m\stackrel{.}{\sigma }\right]$

where ξ_ais the fraction of the Austenite phase prior to the present transformation from Austenite to Martensite, T is the temperature, T_fmis the transition temperature from Austenite to Martensite, σ_mis an indication of the range of temperature around the transition temperature T_fmduring which the phase change occurs, σ is the stress and K_mis the stress curve-fitting parameter which is obtained from the unloading part of the stress-strain characteristic.

25. A system as claimed inclaim 24 wherein the strain in the filament is controlled by means of a gain scheduled Linear Quadratic Regulator control, which is done by optimizing the cost function,

$J\left(u\right)={\int }_0^{\infty }\left(x^TQx+u^TRu\right)t\mathrm{where}\stackrel{\_}{x}=\left[\begin{array}{c}\epsilon \\ {\sigma }_{\epsilon }\end{array}\right]$

where ε is the strain and σ_eis the integral of error, u is the voltage input, Q and R are positive definite matrices.

26. A system as claimed inclaim 24 wherein the strain in the filament is controlled by means of a robust H_∞ loop shaping controller to ensure robust stability and performance.

27. A three-dimensional force sensor adapted to be attached to the distal end of a catheter comprising:

a two-dimensional optical position sensing detector spaced from the distal end of the catheter;

a spring attaching the detector to the catheter whereby force on the detector moves the detector relative to the distal end of the catheter;

a light source located at the end of the catheter is coupled to an optical fibre whereby the light source emits a beam spot; and

a means for determining the position of the beam spot on the detector.

28. A three-dimensional force sensor as claimed inclaim 27 further including a hemispherical tip attached to the detector extending outwardly from the detector away from the spring.

29. A three-dimensional force sensor as claimed inclaim 28 further including a means for determining the intensity of the light on the detector.

30. A three-dimensional force sensor as claimed inclaim 29 wherein the spring is a conical spring.

31. A three-dimensional force sensor as claimed inclaim 27 wherein the detector is a four quadrant photo detector.

32. A three-dimensional force sensor as claimed inclaim 27 wherein the detector is a lateral-effect detector.

33. A three-dimensional force sensor as claimed inclaim 27 wherein the light source is a laser.

34. A three-dimensional force sensor as claimed inclaim 27 wherein the light source is a high intensity light emitting diode.

35. A method of controlling a catheter having a plurality of shape memory alloy filaments at the distal end thereof comprising the steps of:

receiving a strain value for at least one of the filaments;

determining a phase change that will result in the strain value, whereby the phase change is dependent on the temperature and the temperature is dependent on a voltage; and

setting the voltage in each filament thereby resulting in bending of the catheter in three dimensional space

36. A method of controlling a catheter as claimed inclaim 35 wherein the catheter has a force sensor at the distal end thereof and further including the step of receiving a signal from the force sensor.

37. A method of controlling a catheter as claimed inclaim 36 further including the step of determining the strain value from the signal received from the force sensor.

38. A method of controlling a catheter as claimed inclaim 36 wherein the strain value is received from a haptic device and further including the step of transmitting the signal from the force sensor to the haptic device.

39. A method of controlling a catheter as claimed inclaim 37 wherein the temperature is determined using a formula given by

$\stackrel{.}{T}=\frac{1}{mc_p}\left[\frac{V^2}{R}-hA\left(T-T_a\right)\right]$

where m is the mass per unit length, c_pis the specific heat capacity, V is the voltage applied across the shape memory alloy filament, R is the resistance per unit length, h is the coefficient of convection, A=πd is the circumferential area of cooling, d is the diameter of the filament, T is the temperature and T_ais the ambient temperature.

40. A method of controlling a catheter as claimed inclaim 39 wherein the mole fraction of shape memory alloy filament in the Austenite phase while heating is determined using a formula given by:

$\stackrel{.}{\xi }=\frac{{\xi }^2}{{\xi }_m}\left[\exp \left(\frac{T_{f_a}-T}{{\sigma }_a}+K_a\sigma \right)\right]\left[\frac{\stackrel{.}{T}}{{\sigma }_a}-K_a\stackrel{.}{\sigma }\right]$

where ξ is the fraction of the Austenite phase, ξ_mis the fraction of the Martensite phase prior to the present transformation from Martensite to Austenite, T is the temperature, T_fais the transition temperature from Martensite to Austenite, σ_ais an indication of the range of temperature around the transition temperature T_faduring which the phase change occurs, σ is the stress and K_ais the stress curve-fitting parameter which is obtained from the loading plateau of the stress-strain characteristic with no change in temperature; and

the mole fraction of shape memory alloy filament in the Austenite phase while cooling is determined using the formula given by:

$\stackrel{.}{\xi }=\frac{{\xi }^2}{{\xi }_a}\left[\exp \left(\frac{T_{f_m}-T}{{\sigma }_m}+K_m\sigma \right)\right]\left[\frac{\stackrel{.}{T}}{{\sigma }_m}-K_m\stackrel{.}{\sigma }\right]$

where ξ_ais the fraction of the Austenite phase prior to the present transformation from Austenite to Martensite, T is the temperature, T_fmis the transition temperature from Austenite to Martensite, σ_mis an indication of the range of temperature around the transition temperature T_fmduring which the phase change occurs, σ is the stress and K_mis the stress curve-fitting parameter which is obtained from the unloading part of the stress-strain characteristic.

41. A method of controlling a catheter as claimed inclaim 40 where the strain is determined using a formula given by

{dot over (σ)}=D{dot over (ε)}+θ_t{dot over (T)}+Ω{dot over (ξ)}

where σ is the stress in the filament, D is the average of the Young's Moduli for the Martensite and Austenite phases, ε is the strain, θ_tis the thermal expansion factor, and Ω is the phase transformation contribution factor.

42. A method of controlling a catheter as claimed inclaim 41 wherein the strain in the filament is controlled by means of a gain scheduled Linear Quadratic Regulator control, which is done by optimizing the cost function,

$J\left(u\right)={\int }_0^{\infty }\left(x^TQx+u^TRu\right)t\mathrm{where}\stackrel{\_}{x}=\left[\begin{array}{c}\epsilon \\ {\sigma }_{\epsilon }\end{array}\right]$

ε is the strain and σ_eis the integral of error, u is the voltage input, Q and R are positive definite matrices.

43. A method of controlling a catheter as claimed inclaim 42 wherein the strain in the filament is controlled by means of a robust H_∞ loop shaping controller to ensure robust stability and performance.

Images