|
|||||||
Orthopedics Application - Total Hip Replacement
A partial human model (lower body) is created in this example and augmented with a total hip replacement at the right hip. The hip replacement consists of the geometry files (shell files) for the hip stem and the acetabular liner. These files are imported into the human model and secured to the bones using bushing force elements. A walk-sit-walk simulation is performed using motion capture data. The objective of this exercise is to determine if the hip stem impinges on the acetabular liner, and to determine the forces at the bone/component interfaces for this common everyday activity. A combination inverse-forward dynamics simulation will be performed with full foot contact. Key skills exercised in this tutorial include:
Sections
Generating of the Body Segments Joints, and Motion DataIn this phase, the SLF file is used to create the human body model from measurements, joints from joint data, posture from posture data and motion from recorded motion data. The body segments are created using the parameters stored in the SLF file. This file contains information on the subject name, gender, age, height and weight. LifeMOD™ uses this information to extract body segment measurements and mass properties from the internal anthropometric database. Passive joints are created for the inverse-dynamics phase of the simulation process. For this model passive joints will be created for the inverse-dynamics simulation. The passive joint consists of a tri-axis hinge joint (3 DOF) which includes angulation stops, stiffness and damping torques. This type of joint is used primarily to stabilize the body during the inverse-dynamics simulation. The parameters are included in the SLF file Finally, the motion data (MOCAP) for the walk-sit-walk activity is imported into the model and used to drive the motion agents created on the model. There are two components to the motion agent. A yellow sphere designates the location of the data point and the red sphere designates the marker location on the human model. The yellow sphere is attached to the red sphere via a bushing element with properties designated below. During the inverse dynamic simulation, the yellow sphere will move according to the MOCAP data, while influencing the motion of the red sphere attached to the body. It is during this analysis that muscle contraction histories will be recorded. The motion agent stiffness properties are entered in the panel in Figure 1. The motion trajectory data is included in the SLF file. Step 1: Bring up import panel
Step 2: Import the body, joints, posture and motion
Running the Equilibrium AnalysisIn order to fit the model to the data positions, an equilibrium analysis must be performed. This is a dynamics analysis which holds the positions of the data-driven motion agents (yellow balls) fixed, while finding the minimum energy configuration in the springs of the motion agents. The motion agents with the higher weights will have more influence on the model and the initial configuration. Step 3: Bring up analyze panel
Step 4: Run the equilibrium simulation
Step 5: Update the model configuration with static results
Step 6: Align the body markers with data
Creating Passive Muscle ForcesThe next step in the process is to create soft tissues (muscles) on the model. LifeMOD™ automatically creates a set of basic muscle groups for the body. Muscles consist of recording elements or trained elements. The recording elements are simple data collectors which record the contraction history of the muscle during an activity when the model is moved using external drivers such as motion agents. Trained elements can be either PID closed-loop force actuators or actuation curve open-loop force actuators acting to drive the skeleton's motion. Muscle parameters such as physiological cross sectional area (pCSA) and maximum tissue stress are used to calculate the maximum force potential of the particular muscle. LifeMOD™ contains a database of pCSA values for each muscle and is scaled accordingly based on the input body parameters (ht, wt, gender and age). Further, the force output of the muscle may be scaled from 0% to 200% to change the contributions of each particular muscle. For information on selecting specific model parameters for this section see Choosing Model Parameters. Step 7: Bring up the create muscle-tendon panel
Step 8: Set the fields for the muscle generation
Step 9: Create the muscles
Adding the Hip Replacement DeviceWith the model in the proper position and posture, the joints of the right hip will be deleted and replaced with the total hip replacement. Geometry files of the hip system are imported as shell files. Geometry may also be imported via IGES, Step, Parasolids or any format currently supported by MSC.Adams.
Step 11: Display right hip joint
Step 12: Delete the right hip joint
Step 13: Bring up the single part creation panel
Step 14: Create the acetabular cup part
Step 15: Create a fixation force between the cup and the hip bone
OR Create the bushing using the following ADAMS/View commands: marker cre marker=.World.Diego_Lower_Torso.cup loc=(loc_relative_to({0,0,0},.World.Diego_cup.cm)) ori=(ori_relative_to({0,0,0},.World.Diego_cup.cm))
Step 16: Display the femur bone
Step 17: Bring up the single part creation panel
Step 18: Import the geometry and create the hip stem part
Step 19: Create a joint between the stem and the femur bone
OR Create the bushing enter the following ADAMS/View commands: marker cre ma=.World.Diego_Right_Upper_Leg.stem loc=(loc_relative_to({0,0,0},.World.Diego_stem.cm)) ori=(ori_relative_to({0,0,0},.World.Diego_stem.cm)) force create element bushing bushing=.World.stem_femur i_mark=.World.Diego_stem.cm j_mark=.World.Diego_Right_Upper_Leg.stem damping = 1e8, 1e8, 1e8 stiffness = 1e9, 1e9, 1e9 force_preload = 0.0, 0.0, 0.0 tdamping = 1e8, 1e8, 1e8 tstiffness = 1e9, 1e9, 1e9 torque_preload = 0.0, 0.0, 0.0 Step 20: Create a bushing force between the stem ball and the cup
OR Create the bushing enter the following ADAMS/View commands: marker create marker= .World.Diego_Cup.Stem location = (loc_relative_to({0,0,-8},.world.Diego_cup.cm)) orientation = 0.0, 0.0, 0.0 marker create marker= .World.Diego_Stem.Cup location = (loc_relative_to({0,0,-8},.world.Diego_cup.cm)) orientation = 0.0, 0.0, 0.0 force create element bushing bushing= .World.stem_cup damping = 1e7, 1e7, 1e7 stiffness = 1e8, 1e8, 1e8 force_preload = 0.0, 0.0, 0.0 tdamping = 100,100,100 tstiffness = 0,0,0 torque_preload = 0.0, 0.0, 0.0 i_marker_name = .World.Diego_Cup.Stem j_marker_name = .World.Diego_Stem.cup Creating the Foot Contact ForcesContact forces must be created at the foot to generate a reaction between the feet floor and the pelvis and the stool. When the model segments were created, contact ellipsoids were scaled and positioned at the metatarsal heads and calcaneous. These ellipsoids will be used for the foot-floor contact. Step 21: Bring up the contact create panel
Step 22: Create the ground contact markers
OR Create the marker using the following ADAMS/View command: marker cre marker= .World.ground.CON location = 0,5,0 orientation = 180.0, 90.0, 180.0 relative_to = .World Step 23: Create the foot/floor contact forces
Step 24: Create the stool (chair)
Fix the chair to the ground by selecting the fixed joint icon from the main toolbox. Select the leg of the chair as the first object and the ground as the second. Select .World.ground.CON as the location. Create another joint to fix the seat of the chair to the leg. Select the seat as the first object and the leg as the second. Select .World.ground.CONb as the location. OR Create a marker and simple graphics to designate the location of the stool using the following ADAMS/View commands: marker create marker=.world.ground.CONb location=0,500.7,0 orientation=0,-90,0
Step 25: Create pelvis-stool contact forces
Step 26: Create the specific data requests
Create a reference marker .World.Diego_Stem.Ref. Select Modify, input (0, 0, 50) as the location relative to .World.Diego_Cup.Stem. Create another request named REQ_Stem_Cup_Orientation. Set parameters to figure 12. Select OK. OR Create the requests using the following ADAMS/View commands: output_control create request request_name = REQ_Stem_Cup_Force f1 = "BUSH(.World.Diego_Stem_Cup, 0, 1, 0)" marker create marker=.World.Diego_Stem.Ref loc=(loc_relative_to({0,0,50},.World.Diego_Cup.stem))
Running the Inverse-Dynamics SimulationWith the model in the proper position, the hip components implanted in the model, the motion capture data read in, and the motion agents installed, an inverse-dynamics simulation may be performed. This simulation is performed to record the muscle contractile history for the stair-stepping activity. In the forward-dynamics simulation, to be done later, the muscle contraction histories will be used in the contractile elements in the muscles to produce forces to allow the model the replicate the motion. Step 27: Bring up the analyze panel
Step 28: Run the dynamics simulation
Step 29: Display animation
Preparing the Model for the Forward-Dynamics SimulationWith the muscle contraction history recorded from the inverse-dynamics simulation, it is now used in linear PID-Servo formulation to produce a force to recreate the motion history. The process entails deactivating the Motion Agents and updating the muscles. For information on selecting specific model parameters for this section see Choosing Model Parameters. The motion agents are removed from the model and a "Tracker Agent" is installed. The tracker agent is a motion agent located at the center of the pelvis which provides force-stabilization for the forward-dynamics simulation. During the inverse-dynamics simulation the location and orientation of the frame of the tracker agent is recorded (it is not generating a force during the inverse-dynamics simulation). The location and orientation information may then be used to drive the tracker agent in the forward-dynamics simulation. Usually various degrees-of freedom are specified as "free" to allow for proper dynamical interaction. For this example the freedom in the direction normal to the floor would be specified as free, to allow for proper ground reaction force generation between the feet and the steps. In this example the tracker agent accounts for the fact that the upper body is missing from the model. It compensates for the forces of the arms and upper body transmitted through the trunk. Step 30: Bring up the tissue training panel
Step 31: Install ACTIVE contractile element
Step 32: Set fields and update joints
Step 33: Bring up the motion agent tracker panel
Step 34: Create the tracking agent
Running the Forward-Dynamics SimulationWith the contact forces created, the tracker agent place, and the contractile elements in the muscles to include the motion splines from the inverse-dynamics simulation, the forward-dynamics simulation is ready to be performed. Step 35: Bring up the analysis panel
Step 36: Run the forward dynamics simulation
Interrogating the ResultsWhen the simulation is complete the model may be animated and the results reviewed. Various data may be presented from the forward-dynamics simulation including:
Step 37: Display animation
Step 38: Display simulation with skin/skel model
Step 39: Display simulation close up view
Step 40: Display simulation close up side view
Step 41: Bring up results panel
Step 42: Plot the right foot ground reaction force
Step 43: Plot the left foot ground reaction force
Step 44: Plot the pelvis stool contact forces
Step 45: Animate side view
Step 46: Animate side view close up
Step 47: Plot the stem/cup forces
Step 48: Animate front view
mdi graphic_force object=.World.Diego_stem_cup type=1 Select front view, divide window. Check zoom and enter the center as -270,496,0 and the scale as 1.7. Select PLAY. Step 49: Plot the stem orientation with respect to the cup
Step 50: Animate side view close up
mdi graphic_force object=.World.Diego_stem_cup type=0 Select front view, divide window. Check zoom and enter the center as -797,876,0 and the scale as 6.7. Check camera and enter the lower_torso as the part to fix the camera. Select PLAY Step 52: Plot the right iliacus tension
Step 53: Plot the left iliacus tension
Step 54: Turn muscle graphics scaling on and animate side view
Step 55: Animate side view
Step 56: Animate side view close up
Step 57: DEMO COMPLETE FurtherThis model could be used to explore may involved with the design of total hip replacements including:
AcknowledgementA special thanks for furnishing the data for this model to: Diego Crovato
|
|||||||
|
|||||||