1
|
Kadar A, Haddara MM, Fan S, Chinchalkar S, Ferreira LM, Suh N. Use of Thermoplastic Rings Following Venting of Flexor Tendon Pulleys: A Biomechanical Analysis. J Hand Surg Am 2021; 46:485-492. [PMID: 33357985 DOI: 10.1016/j.jhsa.2020.11.003] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/05/2019] [Revised: 09/08/2020] [Accepted: 11/06/2020] [Indexed: 02/02/2023]
Abstract
PURPOSE Normal digital flexion relies on flexor tendon pulleys to convert linear muscular force to angular digital motion. However, there is a growing trend to vent them partially during flexor tendon repair. The objective of this study was to examine the effects of a thermoplastic ring, acting as an external pulley, on flexor tendon biomechanics and finger range of motion (ROM) after pulley venting. METHODS We tested 15 cadaveric digits using an in vitro active finger motion simulator. We measured loads induced by flexor digitorum superficialis (FDS) and flexor digitorum profundus (FDP) as well as joint ROM with sequential sectioning of the A2, A3, and A4 pulleys compared with an intact pulley condition. At each stage, external thermoplastic pulley rings were applied snugly over the proximal and middle phalanges to recreate A2 and A4 function, respectively. RESULTS After complete venting of the A2, A3, and A4 pulleys, proximal interphalangeal joint ROM significantly decreased by 13.4° ± 2.7° and distal interphalangeal joint ROM decreased by 15.8° ± 2.1°. Application of external rings over the proximal and middle phalanx resulted in a residual ROM decrease of 8.3° ± 1.9° at the proximal interphalangeal joint and 7.9° ± 2.1° at the distal interphalangeal joint, nearly restoring ROM. Similarly, complete pulley venting resulted in reduced FDS load by 37% and FDP load by 50% compared with intact pulleys. After application of external rings, loads were restored almost to normal, with a 9% reduction for FDS load and 9% reduction for FDP load compared with intact pulleys. CONCLUSIONS The application of thermoplastic rings acting as external pulleys is an effective, noninvasive, and reproducible approach to restore flexor tendon biomechanics and digit ROM after pulley venting. CLINICAL RELEVANCE Thermoplastic rings may be a useful therapeutic adjunct in restoring joint ROM and flexor tendon loads after surgical venting of the pulleys.
Collapse
Affiliation(s)
- Assaf Kadar
- Division of Orthopedic Surgery, University of Western Ontario, Roth-McFarlane Hand and Upper Limb Center, St Joseph's Health Care, London, Ontario, Canada; Department of Orthopedics, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Mohammad M Haddara
- Division of Plastic and Reconstructive Surgery, University of Western Ontario, Roth-McFarlane Hand and Upper Limb Center, St Joseph's Health Care, London, Ontario, Canada
| | - Stacy Fan
- Department of Biomedical Engineering, University of Western Ontario, Roth-McFarlane Hand and Upper Limb Center, St Joseph's Health Care, London, Ontario, Canada
| | - Shrikant Chinchalkar
- Division of Orthopedic Surgery, University of Western Ontario, Roth-McFarlane Hand and Upper Limb Center, St Joseph's Health Care, London, Ontario, Canada
| | - Louis M Ferreira
- Division of Plastic and Reconstructive Surgery, University of Western Ontario, Roth-McFarlane Hand and Upper Limb Center, St Joseph's Health Care, London, Ontario, Canada
| | - Nina Suh
- Division of Orthopedic Surgery, University of Western Ontario, Roth-McFarlane Hand and Upper Limb Center, St Joseph's Health Care, London, Ontario, Canada.
| |
Collapse
|
2
|
Snoeck O, Coupier J, Beyer B, Salvia P, Lefèvre P, Van Sint Jan S, Rooze M, Feipel V. The biomechanical role of the lacertus fibrosus of the biceps brachii Muscle. Surg Radiol Anat 2021; 43:1587-1594. [PMID: 33751178 DOI: 10.1007/s00276-021-02739-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Accepted: 03/13/2021] [Indexed: 11/28/2022]
Abstract
PURPOSE The lacertus fibrosus (LF) is involved in various surgeries. However, the biomechanical contribution of the LF remains unclear. The aim of this study was to determine the role of the lacertus fibrosus on the elbow and forearm kinematics and on the biceps brachii muscle lever arms. METHODS This biomechanical study was performed on seven fresh-frozen upper limbs of cadavers. Elbow flexion, forearm supination, and biceps brachii muscle lever arms were analyzed in the intact conditions (I) and after superficial (R) and deep part (R2) of the lacertus fibrosus release, respectively. RESULTS Elbow flexion shows a significant difference (p < 0.0001) between I, R, R2. Abduction/adduction shows a significant difference between I-R (p < 0.0001) and I-R2 (p < 0.0001). Supination does not show a significant difference in mean maximum amplitude, but between 40 and 70%, there are significant differences. There is a significant mean decrease of lever arm in flexion (28%) and supination (50%) after superficial and deep part of the lacertus fibrosus release. CONCLUSION The results of this study show that the lacertus fibrosus increases the lever arm during flexion and supination. It limits the flexion and abduction of the elbow and supination of the forearm. Lacertus fibrosus maintains the rhythmicity between the elbow flexion and supination of the forearm. LEVEL OF EVIDENCE Basic science study, biomechanics.
Collapse
Affiliation(s)
- Olivier Snoeck
- Laboratory of Anatomy, Biomechanics and Organogenesis (LABO), Faculty of Medicine, Université Libre de Bruxelles ULB, Brussels, Belgium. .,Laboratory for Functional Anatomy, Faculty of Motor Sciences, Université Libre de Bruxelles ULB, Brussels, Belgium.
| | - Jérôme Coupier
- Laboratory of Anatomy, Biomechanics and Organogenesis (LABO), Faculty of Medicine, Université Libre de Bruxelles ULB, Brussels, Belgium.,Laboratory for Functional Anatomy, Faculty of Motor Sciences, Université Libre de Bruxelles ULB, Brussels, Belgium
| | - Benoît Beyer
- Laboratory of Anatomy, Biomechanics and Organogenesis (LABO), Faculty of Medicine, Université Libre de Bruxelles ULB, Brussels, Belgium.,Laboratory for Functional Anatomy, Faculty of Motor Sciences, Université Libre de Bruxelles ULB, Brussels, Belgium
| | - Patrick Salvia
- Laboratory of Anatomy, Biomechanics and Organogenesis (LABO), Faculty of Medicine, Université Libre de Bruxelles ULB, Brussels, Belgium.,Laboratory for Functional Anatomy, Faculty of Motor Sciences, Université Libre de Bruxelles ULB, Brussels, Belgium
| | - Philippe Lefèvre
- Laboratory of Anatomy, Biomechanics and Organogenesis (LABO), Faculty of Medicine, Université Libre de Bruxelles ULB, Brussels, Belgium.,Laboratory for Functional Anatomy, Faculty of Motor Sciences, Université Libre de Bruxelles ULB, Brussels, Belgium
| | - Serge Van Sint Jan
- Laboratory of Anatomy, Biomechanics and Organogenesis (LABO), Faculty of Medicine, Université Libre de Bruxelles ULB, Brussels, Belgium.,Laboratory for Functional Anatomy, Faculty of Motor Sciences, Université Libre de Bruxelles ULB, Brussels, Belgium
| | - Marcel Rooze
- Laboratory of Anatomy, Biomechanics and Organogenesis (LABO), Faculty of Medicine, Université Libre de Bruxelles ULB, Brussels, Belgium.,Laboratory for Functional Anatomy, Faculty of Motor Sciences, Université Libre de Bruxelles ULB, Brussels, Belgium
| | - Véronique Feipel
- Laboratory of Anatomy, Biomechanics and Organogenesis (LABO), Faculty of Medicine, Université Libre de Bruxelles ULB, Brussels, Belgium.,Laboratory for Functional Anatomy, Faculty of Motor Sciences, Université Libre de Bruxelles ULB, Brussels, Belgium
| |
Collapse
|
3
|
The Effect of Wrist Position on Finger Tendon Loads Following Pulley Sectioning and Operative Reconstruction. JOURNAL OF HAND SURGERY GLOBAL ONLINE 2019. [DOI: 10.1016/j.jhsg.2019.04.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
|
4
|
Rahman M, Cil A, Stylianou AP. Medial Collateral Ligament Deficiency of the Elbow Joint: A Computational Approach. Bioengineering (Basel) 2018; 5:bioengineering5040084. [PMID: 30308994 PMCID: PMC6316890 DOI: 10.3390/bioengineering5040084] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2018] [Revised: 10/02/2018] [Accepted: 10/08/2018] [Indexed: 11/16/2022] Open
Abstract
Computational elbow joint models, capable of simulating medial collateral ligament deficiency, can be extremely valuable tools for surgical planning and refinement of therapeutic strategies. The objective of this study was to investigate the effects of varying levels of medial collateral ligament deficiency on elbow joint stability using subject-specific computational models. Two elbow joint models were placed at the pronated forearm position and passively flexed by applying a vertical downward motion on humeral head. The models included three-dimensional bone geometries, multiple ligament bundles wrapped around the joint, and the discretized cartilage representation. Four different ligament conditions were simulated: All intact ligaments, isolated medial collateral ligament (MCL) anterior bundle deficiency, isolated MCL posterior bundle deficiency, and complete MCL deficiency. Minimal kinematic differences were observed for isolated anterior and posterior bundle deficient elbows. However, sectioning the entire MCL resulted in significant kinematic differences and induced substantial elbow instability. Joint contact areas were nearly similar for the intact and isolated posterior bundle deficiency. Minor differences were observed for the isolated anterior bundle deficiency, and major differences were observed for the entire MCL deficiency. Complete elbow dislocations were not observed for any ligament deficiency level. As expected, during isolated anterior bundle deficiency, the remaining posterior bundle experiences higher load and vice versa. Overall, the results indicate that either MCL anterior or posterior bundle can provide anterior elbow stability, but the anterior bundle has a somewhat bigger influence on joint kinematics and contact characteristics than posterior one. A study with a larger sample size could help to strengthen the conclusion and statistical significant.
Collapse
Affiliation(s)
- Munsur Rahman
- Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, 5110 Rockhill Road, Kansas City, MO 64110, USA.
| | - Akin Cil
- Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, 5110 Rockhill Road, Kansas City, MO 64110, USA.
- Department of Orthopaedic Surgery, University of Missouri-Kansas City, 2411 Holmes Street, Kansas City, MO 64108, USA.
- Department of Orthopaedics, Truman Medical Centers, 2301 Holmes Street, Kansas City, MO 64108, USA.
| | - Antonis P Stylianou
- Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, 5110 Rockhill Road, Kansas City, MO 64110, USA.
| |
Collapse
|
5
|
West EA, Knowles NK, Athwal GS, Ferreira LM. A 3D comparison of humeral head retroversion by sex and measurement technique. Shoulder Elbow 2018; 10:192-200. [PMID: 29796107 PMCID: PMC5960870 DOI: 10.1177/1758573217711897] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/13/2017] [Accepted: 05/01/2017] [Indexed: 11/16/2022]
Abstract
BACKGROUND Accurate humeral head reconstruction during shoulder arthroplasty is partially dependent on correctly estimating and replicating native version. The present study evaluated the effects of sex and measurement technique on three-dimensional (3D) humeral version measurements made using the transepicondylar, forearm and flexion-extension axes. METHODS Fifty-two full-arm computed tomography scans were converted to 3D models and geometry extracted to define landmarks and coordinate systems. An anatomic humeral head osteotomy plane was used to measure version relative to the three measurement techniques and compare between sexes. RESULTS The measurement technique used had a significant affect (p < 0.001) on the resulting version measurement. The forearm axis technique consistently resulted in higher measured version compared to either the flexion-extension [mean (SD) males 9° (4°), females 13° (5°), p < 0.001] or the transepicondylar axes [mean (SD) males 8° (4°), females 11° (4°), p < 0.001]. Version in males was 7° greater than females when referencing either the flexion-extension [p = 0.029; mean (SD) males 37.7° (11°), females 30.4° (13°)] or transepicondylar axes [p = 0.045; mean (SD) males 39° (11°), females 32° (12°)]. CONCLUSIONS The choice of measurement technique can affect the humeral version angle. These results are important because measuring version using the epicondyles pre-operatively, and subsequently the forearm intra-operatively, will result in approximately 10° under-retroverted osteotomy. For example, 0° neutral version cut during reverse arthroplasty measured referencing the forearm results in 10° anteverted osteotomy when referencing the distal humerus.
Collapse
Affiliation(s)
- Emily A. West
- Department of Biomedical Engineering, The University of Western Ontario, London, ON, Canada,Roth
- McFarlane Hand and Upper Limb Centre, St Joseph's Hospital, London, ON, Canada
| | - Nikolas K. Knowles
- Department of Biomedical Engineering, The University of Western Ontario, London, ON, Canada,Roth
- McFarlane Hand and Upper Limb Centre, St Joseph's Hospital, London, ON, Canada
| | - George S. Athwal
- Department of Surgery, The University of Western Ontario, London, ON, Canada,Roth
- McFarlane Hand and Upper Limb Centre, St Joseph's Hospital, London, ON, Canada
| | - Louis M. Ferreira
- Department of Biomedical Engineering, The University of Western Ontario, London, ON, Canada,Roth
- McFarlane Hand and Upper Limb Centre, St Joseph's Hospital, London, ON, Canada,Louis M. Ferreira, 1151 Richmond St, London ON, N6A 3K7, Canada.
| |
Collapse
|
6
|
Rahman M, Cil A, Bogener JW, Stylianou AP. Lateral collateral ligament deficiency of the elbow joint: A modeling approach. J Orthop Res 2016; 34:1645-55. [PMID: 26763801 DOI: 10.1002/jor.23165] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/23/2015] [Accepted: 01/04/2016] [Indexed: 02/04/2023]
Abstract
A computational model capable of predicting the effects of lateral collateral ligament deficiency of the elbow joint would be a valuable tool for surgical planning and prediction of the long-term consequences of ligament deficiency. The purpose of this study was to simulate lateral collateral ligament deficiency during passive flexion using a computational multibody elbow joint model and investigate the effects of ligament insufficiency on the kinematics, ligament loads, and articular contact characteristics (area, pressure). The elbow was placed initially at approximately 20° of flexion and a 345 mm vertical downward motion profile was applied over 40 s to the humerus head. The vertical displacement induced flexion from the initial position to a maximum flexion angle of 135°. The study included simulations for intact, radial collateral ligament deficient, lateral ulnar collateral ligament deficient, and combined radial and lateral ulnar collateral ligament deficient elbow. For each condition, relative bone kinematics, contact pressure, contact area, and intact ligament forces were predicted. Intact and isolated radial collateral ligament deficient elbow simulations were almost identical for all observed outcomes. Minor differences in kinematics, contact area and pressure were observed for the isolated lateral ulnar collateral ligament deficient elbow compared to the intact elbow, but no elbow dislocation was detected. However, sectioning both ligaments together induced substantial differences in kinematics, contact area, and contact pressure, and caused complete dislocation of the elbow joint. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1645-1655, 2016.
Collapse
Affiliation(s)
- Munsur Rahman
- Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, 5110 Rockhill Road, Kansas City 64110, Missouri
| | - Akin Cil
- Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, 5110 Rockhill Road, Kansas City 64110, Missouri.,Department of Orthopaedic Surgery, University of Missouri-Kansas City, 5110 Rockhill Road, Kansas City 64110, Missouri
| | - James W Bogener
- Department of Orthopaedic Surgery, University of Missouri-Kansas City, 5110 Rockhill Road, Kansas City 64110, Missouri
| | - Antonis P Stylianou
- Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, 5110 Rockhill Road, Kansas City 64110, Missouri
| |
Collapse
|
7
|
Prediction of elbow joint contact mechanics in the multibody framework. Med Eng Phys 2016; 38:257-66. [DOI: 10.1016/j.medengphy.2015.12.012] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2015] [Revised: 12/08/2015] [Accepted: 12/31/2015] [Indexed: 11/21/2022]
|
8
|
McLachlin SD, Bailey CS, Dunning CE. Use of the alpha shape to quantify finite helical axis dispersion during simulated spine movements. J Biomech 2015; 49:112-118. [PMID: 26653673 DOI: 10.1016/j.jbiomech.2015.06.033] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2014] [Revised: 06/25/2015] [Accepted: 06/30/2015] [Indexed: 10/23/2022]
Abstract
In biomechanical studies examining joint kinematics the most common measurement is range of motion (ROM), yet other techniques, such as the finite helical axis (FHA), may elucidate the changes in the 3D motion pathology more effectively. One of the deficiencies with the FHA technique is in quantifying the axes generated throughout a motion sequence. This study attempted to solve this issue via a computational geometric technique known as the alpha shape, which bounds a set of point data within a closed boundary similar to a convex hull. The purpose of this study was to use the alpha shape as an additional tool to visualize and quantify FHA dispersion between intact and injured cadaveric spine movements and compare these changes to the gold-standard ROM measurements. Flexion-extension, axial rotation, and lateral bending were simulated with five C5-C6 motion segments using a spinal loading simulator and Optotrak motion tracking system. Specimens were first tested intact followed by a simulated injury model. ROM and the FHAs were calculated post-hoc, with alpha shapes and convex hulls generated from the anatomic planar intercept points of the FHAs. While both ROM and the boundary shape areas increased with injury (p<0.05), no consistent geometric trends in the alpha shape growth were identified. The alpha shape area was sensitive to the alpha value chosen and values examined below 2.5 created more than one closed boundary. Ultimately, the alpha shape presents as a useful technique to quantify sequences of joint kinematics described by scatter plots such as FHA intercept data.
Collapse
Affiliation(s)
- Stewart D McLachlin
- Department of Mechanical and Materials Engineering, Western University, London, Ontario, Canada
| | - Christopher S Bailey
- Department of Mechanical and Materials Engineering, Western University, London, Ontario, Canada
| | - Cynthia E Dunning
- Department of Mechanical and Materials Engineering, Western University, London, Ontario, Canada.
| |
Collapse
|
9
|
The effect of radial head implant shape on radiocapitellar kinematics during in vitro forearm rotation. J Shoulder Elbow Surg 2015; 24:258-64. [PMID: 25457776 DOI: 10.1016/j.jse.2014.09.019] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/09/2014] [Revised: 09/09/2014] [Accepted: 09/12/2014] [Indexed: 02/01/2023]
Abstract
BACKGROUND A number of radial head implants are in clinical use for the management of radial head fractures and their sequelae. However, the optimal shape of a radial head implant to ensure proper tracking relative to the capitellum has not been established. This in vitro biomechanical study compared radiocapitellar joint kinematics for 3 radial head implant designs as well as the native head. METHODS Eight cadaveric upper extremities were tested using a forearm rotation simulator with the elbow at 90° of flexion. Motion of the radius relative to the capitellum was optically tracked. A stem was navigated into a predetermined location and cemented in place. Three unipolar implant shapes were tested: axisymmetric, reverse-engineered patient-specific, and population-based quasi-anatomic. The patient-specific and quasi-anatomic implants were derived from measurements performed on computed tomography models. RESULTS Medial-lateral and anterior-posterior translation of the radial head with respect to the capitellum varied with forearm rotation and radial head condition. A significant difference in medial-lateral (P = .03) and anterior-posterior (P = .03) translation was found between the native radial head and the 3 implants. No differences were observed among the radial head conditions except for a difference in medial-lateral translation between the axisymmetric and patient-specific implants (P = .04). CONCLUSIONS Radiocapitellar kinematics of the tested radial head implants were similar in all but one comparison, and all had different kinematics from the native radial head. Patient-specific radial head implants did not prove advantageous relative to conventional implant designs. The shape of the fixed stem unipolar radial head implants had little influence on radiocapitellar kinematics when optimally positioned in this testing model.
Collapse
|
10
|
McLachlin SD, Ferreira LM, Dunning CE. A Refined Technique to Calculate Finite Helical Axes From Rigid Body Trackers. J Biomech Eng 2014; 136:124506. [DOI: 10.1115/1.4028413] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2014] [Accepted: 08/27/2014] [Indexed: 11/08/2022]
Abstract
Finite helical axes (FHAs) are a potentially effective tool for joint kinematic analysis. Unfortunately, no straightforward guidelines exist for calculating accurate FHAs using prepackaged six degree-of-freedom (6DOF) rigid body trackers. Thus, this study aimed to: (1) describe a protocol for calculating FHA parameters from 6DOF rigid body trackers using the screw matrix and (2) to maximize the number of accurate FHAs generated from a given data set using a moving window analysis. Four Optotrak® Smart Markers were used as the rigid body trackers, two moving and two fixed, at different distances from the hinge joint of a custom-machined jig. 6DOF pose information was generated from 51 static positions of the jig rotated and fixed in 0.5 deg increments up to 25 deg. Output metrics included the FHA direction cosines, the rotation about the FHA, the translation along the axis, and the intercept of the FHA with the plane normal to the jig's hinge joint. FHA metrics were calculated using the relative tracker rotation from the starting position, and using a moving window analysis to define a minimum acceptable rotational displacement between the moving tracker data points. Data analysis found all FHA rotations calculated from the starting position were within 0.15 deg of the prescribed jig rotation. FHA intercepts were most stable when determined using trackers closest to the hinge axis. Increasing the moving window size improved the FHA direction cosines and center of rotation accuracy. Window sizes larger than 2 deg had an intercept deviation of less than 1 mm. Furthermore, compared to the 0 deg window size, the 2 deg window had a 90% improvement in FHA intercept precision while generating almost an equivalent number of FHA axes. This work identified a solution to improve FHA calculations for biomechanical researchers looking to describe changes in 3D joint motion.
Collapse
Affiliation(s)
- Stewart D. McLachlin
- Department of Mechanical and Materials Engineering,Western University, London N6A 5B9 ON, Canada
| | - Louis M. Ferreira
- Department of Mechanical and Materials Engineering, Western University, London N6A 5B9 ON, Canada
| | - Cynthia E. Dunning
- Department of Mechanical and Materials Engineering, Western University, London N6A 5B9 ON, Canada e-mail:
| |
Collapse
|