 That’s me in the photo on the left during open heart surgery to replace my native aortic valve.
Years ago I was diagnosed with an “organic murmur” and after following it closely with my cardiologist we decided it was a good time to have a replacement of my aortic valve. I was able to plan for and have exactly the surgery I wanted and I couldn’t be happier with the result!
I look forward to sharing more about my surgery, my new valve and especially what I learned in preparing for my surgery and recovery.
Thanks to everyone who wished me well and helped my through this!
 An important input to FEA is material characterization and limit information. In our laboratory we use a Bose 3230 test system for achieving test speeds up to 200Hz. Depending on the project, we may test complete devices, sub-components, coupon samples and/or simple material characterization samples. Many times, we will test more than one configuration to make sure our data are consistent and ultimately demonstrates the validation of an implantable medical device’s safety and efficacy.
Many factors are involved in deciding what testing would be most efficient to support a product’s validation effort. Characterization of material properties is typically the first step. For medical implants that utilize advanced materials, fine geometric features and deliberate metallurgical processing, it is important to capture these in the specimen design. Fabrication of appropriate samples takes the most effort compared to the relatively short test times and the simplicity of material characterization tests.
After material characterization, testing to validate the performance features of an implant are typically necessary. These tests may also require the fabrication of appropriate samples but can additionally involve the reproduction of complex in vivo load states. Besides verifying that the test setup reproduces the desired load state, performance testing is short term and individual samples can be tested in series to study the effect of loading magnitudes, variations in processing, load combinations, etc. Many times, these studies to quantify the impact of such variables play an important role in a validation strategy and can be efficiently performed since only one or a few cycles of loading are required.
Beyond material characterization and performance testing, material limit properties are the next category of data necessary for input to FEA models. When the required number of cycles is low, say 10 million or fewer, accelerated testing can be performed on individual samples tested in series or multiple samples tested simultaneously. The advantage of testing individual samples was mentioned above while testing multiple samples simultaneously produces stronger statistics more quickly.
When the required number of cycles reaches 100 million or more, a balance between testing enough samples for the desired confidence level and testing the samples to run out at the high cycle limit must be reached. A typical strategy might involve a combination of testing individual samples to 10 million cycles and then selecting one or moreloading conditions/levels for testing multiple samples simultaneously. Such a scenario takes advantage of the relatively short testing time frame for a 10 million cycle test (on the order of 1-2 days, compared to 4-6 weeks for testing to 400 million cycles) and the statistical advantage of testing multiple samples.
In the photo above, a multiple sample fixture capable of testing 12 fatigue coupon samples in a controlled temperature chamber is shown. Such a fixture must be fabricated and aligned very carefully to insure identical loading conditions are achieved for all of the samples. Stroboscopic verification, randomized experimental design, fracture surface investigation and other experimental mechanics techniques also help to build confidence in a validation test methodology.
 I just returned from a workshop sponsored by the FDA (and the NHLBI and NSF) on Computer Methods for Cardiovascular Devices. It was an excellent workshop providing an audience of regulatory, academic and industrial interests a chance to get caught up on the state-of-the-art, trends and, in general, to exchange ideas on the issues of using computational methods to support regulatory filings for medical devices. A general theme emerged for me during the workshop that I’d like to discuss in this article:
We are not providing the FDA with adequate validation of our computational models!
For years now I’ve been helping companies demonstrate the safety and effectiveness of their products. I’ve written many FEA reports that have been reviewed and accepted by the agency, including cases where we’ve argued to forgo expensive and time consuming durability testing in lieu of providing computational results to support safety claims. It has been my experience that the FDA has been very open to such an approach, provided there was an adequate demonstration of the validity of the FEA models.
From what I heard from reviewers at the workshop, however, the typical submission of FEA results does not include adequate validation. I don’t know if it is because companies don’t know how or what to provide for validation of their FEA models or if they are reluctant to share testing or data that the FDA has not specifically asked for or if they have unreasonable expectations about what computational models can replace in terms of physical testing, but it is clear to me that if we want to leverage FEA to streamline the development and approval process then we need to take a proactive role at demonstrating how well our models describe our products.
It is far less expensive and time consuming to perform carefully designed bench tests to validate computational results than it is to run long term durability tests on our devices and hope they pass. Not to mention far less risky from a product development perspective. It was clear in the workshop that the FDA understands this and that they too are motivated to see a better balance between physical testing and computational modeling in a submission.
In my over ten years experience in the field, I have yet to run across a device or a specified test or loading scenario that I could not analyze using Abaqus and achieve excellent agreement between experiment and computer simulation. Many times the endeavor to match experiment and analysis reveals critical insight into the mechanics of the product involved or nuances associated with the loading conditions that lead to important improvements. With advanced contact, strong nonlinear capabilities and the extensibility of user subroutines, Abaqus provides a platform to model almost any physical scenario giving the engineer and product designer a more than ample toolkit for validating any device.
Still, as open and receptive as the FDA may be, they are not in a position to advise on how best to perform the appropriate validation. As engineers we need to establish the validity of our computational models and we need to do so BEFORE we submit results to the FDA. In fact, we need to begin this effort early in the development process before we start making decisions based on our computational data. Otherwise, how can we expect the FDA to accept that our results have emerged from a rigorous engineering methodology?
How much and what type of validation is necessary in any given case depends on how a model is going to be used. Conversely, the confidence we have in a computational model depends on how extensively it has been applied and shown to agree with reality. There are numerous opportunities we have during the development process to establish the validity and range of our computational models. Radial force testing of different stent designs for example provides an excellent opportunity to confirm our models ability to predict reality.
In summary, the time is right for advancing the use of computational models for demonstrating the safety of our products. But we need to be proactive and utilize models that are well grounded in experimental data. How far we are able to leverage these results with the FDA will depend on how good of a job we do at convincing them that they represent actual experience.
 For medical devices, rarely does “one size fit all”. The human body is extremely variable and patient populations, especially those with disease present a wide range of differences that the medical device design engineer must consider.
Identifying the worst case size for a product family is an important part of the validation process for implantable medical devices. Understanding your product and especially the in vivo loading conditions are essential for engineering the structure and material specifications for each size over the intended product range.
Typically, there are numerous sources of nonlinearity associated with implantable medical device design. The materials we use, the geometries and especially the physiology we treat all respond in ways that are difficult to describe in simple terms. It is tempting to design a product by considering an idealized patient population and then simply “scale” that design to smaller and larger sizes. But this approach can can result in a poorly optimized product family. Furthermore, when one considers device/lumen ineraction and teh resulting compliance under physiological loads, identifying the worst case loading condition is not a straightforward activity.
The figure above illustrates how the alternating fatigue strain for a stent-like product can vary for deployment to different diameters. It is based on an analytic model of lumen compliance and finite element analysis models of the two device sizes. Clearly, the results indicate a highly nonlinear system that precludes the selection of a single worst case device size and implant condition based on the “four corners” approach. Assuming that the largest device put into the smallest lumen will result in the most challenged loading condition is niave.
When it is possible, it is preferred to model ALL device sizes to determine the worst case size. It is also advantageous to develop and validate a simulated model of the intended physiology for the implant and use this model to verify the performance of each size in a design family. Parameters such as radial force, anchoringg, dynamic compliance, vessel tortuosity and fatigue loading conditions can then all be evaluated for each device size and safety established for the complete instructions for use (IFU) for the product.
 My first job out of graduate school and post-doctorate work was with the footwear company, NIKE. I was fortunate to hire into the heart of the most advanced footwear technology department in the world and I could not have had a more exciting job description. I began by learning how to design NIKE AIR cushioning bladders. After becoming familiar with the design and development process I was able to make my first contribution, which was to pioneer the use of Finite Element Analysis to predict the inflated geometries of NIKE AIR bladders.
The engineers in the NIKE AIR Technology Lab were very proficient at design and had established an efficient process for prototyping new designs-at least for components based on flat film and lay-flat tubing. However, blow molded AIR bladders which were becoming quite popular required building mold tooling prior to prototyping. A design tool that allowed the refinement of a design prior to cutting tools was needed. These molds were expensive and involved several weeks of lead time. The biggest benefit of an FEA method, however, was the ability to predict the geometry of the bladders so that midsole tooling could be made for production. The current practice was to inflate ten AIR bladders, set them aside, allow them to equilibrate for 6-8 weeks and then measure them with calipers. Several mechanisms affected their change in shape but viscoelastic creep was primarily responsible for the time dependent shape change.
There had been several previous attempts to predict the inflated geometries of NIKE AIR bladders prior to my working at the company. It was generally believed that FEA would not be able to predict inflated geometries with sufficient accuracy. Furthermore, the viscoelastic behavior of the urethane elastomers was considered too complicated to characterize. Having been an ardent experimentalist challenging FEA codes for advanced elastic-plastic fracture and rock mechanics problems, I thought for sure that we would be able to utilize FEA for NIKE AIR bladders.
Force deflection data was available for the NIKE AIR bladder material but before I embarked on fitting that data to a hyperelastic model, I was curious to know what strain levels were associated with an inflated bladder. One of the first things I did was to scrounge up some polarizing filters and rig up a crude polariscope to estimate the strain in the film of an inflated bag. My estimate surprised me—less than 1% strain. I verified my measurement using a grid method and sure enough, the strains were really that low. (I also explored using mercury strain gages for directly measuring strains with some success.)
Armed with this simplification, some basic creep data for urethane, and the help of some software vendors, I was able to demonstrate excellent predictive ability of FEA in determining inflated AIR bladder shape and growth over time. The application of FEA helped reduce product development lead times and significantly decreased the cost and time to develop complex blow molded designs. This last figure illustrates how FEA could help identify “hot spots” that would not inflate uniformly. I have since left NIKE and the AIR Technology Lab but not before seeing a whole department grow up to apply FEA to cushioning technology and many other aspects of footwear and athletic equipment.
 There are many circumstances when it is necessary to model contact and contact interactions when analyzing medical implants. Such circumstances arise when devices interact with tooling, catheters, vessels and when self-contact occurs. This latter situation arises for example when a stent is crimped down tightly to catheter dimensions.
ABAQUS/Standard provides a range of methods for defining and managing contact in FEA simulations. One first defines a master surface and a slave surface which can be based on node or elements sets and/or analytically defined surfaces. Then one must select a contact interaction model and ascribe specific properties to that model. The image on the left shows both node and element sets used in a typical contact definition for self-contact of a stent.
Considerable care must be used to define appropriate surfaces. Contact interactions are computationally expensive and although ABAQUS efficiently manages their computational cost, considerations such as smoothness, consistency and other geometric factors need to be considered to maintain solution stability.
ABAQUS provides a variety of contact interaction models for mechanical contact. These include “hard” contact, “softened contact”, “no separation” to name just a few. Hard contact is the most simple as it simulates no reaction pressure until contact occurs then a quick ramp up as “penetration” increases. Softened contact provides a small, exponentially increasing contact pressure just prior to penetration and increasing pressure thereafter. The no separation model provides increasing pressure as penetration is increased and maintains contact with negative pressure if the surfaces attempt to move away from one another. The behavior of softened contact and no separation contact models are illustrated in the plots below.
Each of these various models have utility in simulating implantable medical devices. For example, the hard contact model and occasionally the no separation model are useful in simulating the interaction of devices with tooling during processing and manufacture and also modeling the interaction of the device with a catheter. The softened contact model is useful for simulating a device when interacting with tissue. Regardless of the model chosen, it is important to experiment with and evaluate the effects of the chosen parameters for the models. The selection of these parameters will affect the solution and therefore it is important to study such effects and verify that the results make sense for the problem being considered.
The Table below summarizes one such study using different models and sets of parameters for typical stent analysis scenarios. For this Table, a two-strut model is crimped to catheter dimensions using different contact interaction models. As much as 0.2% strain difference occurs between the hardest and the softest of the various models considered.
 Finite Element Analysis (FEA) can produce an enormous amount of data as output. Solution variables such as stress and strain are computed throughout an analysis for each increment and at each location within the model. These solution variables are computed at what are called “integration points”. These locations ARE NOT the same as the nodes of an element and it is important when post-processing FEA results to understand how the actual solution data is used to create contour plots, and how to extract accurate data representing the solution of the problem you are modeling. As is generally the case, it is the responsibility of the analyst to make sure the engineering is consistent with the problem at hand.
During a typical non-linear FEA solution process, numerous increments are taken to establish an equilibrated solution for the given applied loading. Stress equilibrium, strain compatibility and other mechanics equations are simultaneously “solved” by adjusting local solution variables such as stress and strain throughout the entire model being analyzed. These equations are written in such a way as to ultimately satisfy equilibrium conditions at an elemental level—the details of which depend on the shape function of the element being used. Regardless of element type and formulation, the shape functions determine how the discrete solutions for field variables like stress and strain are represented throughout an element, and compared with forces and displacements to evaluate the various equations governing the solution of the problem.
What is important to understand is that while forces and displacements are computed at nodal locations, stresses and strains are computed at integration points. The figure above for a typical 8-noded, linear, three-dimensional element illustrates the difference between nodal locations and integration points. The nodal locations are at the corners of the cube while the integration points are located within the element.
In order to generate contour plots from FEA results, it is necessary to extrapolate the stress or strain values from the integration points to the nodal locations. Since each node in a model will generally be shared by more than one element these extrapolated values will also have to be averaged in order to produce a smooth contour plot. While there are parameters that can be set to control how this process is handled by your post-processing software, there are many factors which can affect the accuracy of the extrapolated values.
The degree to which the contour plot data differs from the actual data depends on the element type/formulation and particularly on the steepness of the underlying field variable gradient and to a large extent the quality of the underlying finite element mesh. Large gradients and a coarse mesh will obviously produce a greater degree of difference between the integration point data and extrapolated nodal averages.
Relying upon discrete data taken from contour plots for stress or strain can be misleading and inconsistent between various analysis runs. One reliable method for extracting specific stress and strain values is to use “integration point” data. In this way, you are assured of getting an accurate stress or strain value that relates directly to the solution of the underlying mechanics equations. However, in the case of performing a fatigue analysis, for example, stresses and strains at the surface of a component may be more relevant, as it is well known that cracks and defects generally initiate on the surface a part. In such a case, it is the analysts responsibility to understand the impact of using either integration point data or averaged nodal data.
The contour plot below on the left was produced with the nodal averaging turned off. You can see that when the nodal data is extrapolated for each element that the result produces discontinuities between adjacent elements. This is inevitable because the variation of stress and strain over an element is much more complicated than a simple linear relationship. The contour plot on the right was produced by turning nodal averaging back on. It is much smoother, but that smoothness comes with a loss of information.
 
 Finite Element Analysis is routinely used to evaluate the performance and durability of medical devices. When using an implicit method such as ABAQUS/Standard, FEA can also be used to evaluate the geometric stability of a proposed design. Such issues could arise for example when designing stents for large vessels, such as the aorta, when the large diameters exceed the length of the stent. When the length/diameter aspect ratio of a stent exceeds one, stability of a design could be an issue, with unstable deformation occurring causing the stent to dislodge and turn sideways in the vessel.
We routinely evaluate stent components for stability and occasionally are required to make suggestions for improving a component’s stability. In the plot above showing lateral force as a function of diameter, the original (in blue) and improved design (in red) demonstrate the tendency of a device to remain stable. The further into negative territory of lateral force for any given design directly correlates to a device becoming unstable. The design represented by the blue curve quickly turns into negative territory and remains there throughout most of the loading. The curve in red shows a dramatic improvement with only a portion of the loading curve becoming negative and then with only a small magnitude.
The images below illustrate buckling behavior of a stent with the image on the right becoming unstable.
A typical process for making stents is electropolishing to remove material and processing defects and to impart a smooth finish with rounded corners. It is well known that a polished part will have superior fatigue resistance, but how much does the rounding of edges affect the stress and strain in the stent when loaded?
We studied the effect of material removal using a generic Nitinol two-strut stent model. In the figure to the left are four models: From left to right, a typical stent geometry meshed with no consideration for edge detail; a model with a 50 micron chamfer; a model with a 100 micron chamfer; and finally, a model with a 100 micron radius. These four geometries were analyzed by simulating catheter crimp loading. When crimped to minimum dimensions, the strains were as follows:
As can be seen, there is very little impact of the chamfer or the radius on the predicted strain values. These results give credence to the standard practice of neglecting the effect of edge features. It is common, however, to account for the effects of electropolishing by removing an appropriate amount from the as-cut laser geometry. As will be discussed in subsequent posts, verification of actual component geometry is a very important part of the FEA process.<-->
 Finite Element Analysis is an excellent engineering tool for optimizing medical implants. It provides a fast and reliable method for evaluating both the performance and safety issues associated with any given product. A variety of design concepts can be evaluated and relative safety factors determined for the full size range of an intended product.
Parametric studies are conveniently carried out and relative improvements in performance are clearly distinguished in objective engineering terms. Interactions of the device with the human body can be performed and can incorporate bench, animal and human data to improve predictive capabilities. Self-consistent methodology provides a rational foundation to make improvements on existing technology.
Advances in mechanical testing, materials selection and processing continue to provide more choices for the design and implementation of existing and new concepts for minimally invasive medical therapies. Finite Element Analysis provides the capability to keep up with these changes and offers an ideal methodology to minimize risk, shorten development time and increase confidence in the development of life saving technology.
|
|
That’s me in the photo on the left during open heart surgery to replace my native aortic valve.
