Mathematical musculoskeletal modeling and simulation provide a means for proactive injury prevention. To be effective, these models must physiologically replicate shoulder function. Although several muscle force prediction (MFP) shoulder models exist, few have attempted to integrate the force contributions of ligaments, especially the glenohumeral ligaments. The purpose of the current study was to integrate seven shoulder ligaments into an existing computational shoulder model, and analyze both individual ligament characteristics and the influence on the model outputs. Using data from the literature, seven shoulder ligaments were integrated into the model: the costoclavicular, conoid, trapezoid, coracohumeral, superior glenohumeral, middle glenohumeral, and inferior glenohumeral. 10 subjects performed isometric exertions in 56 posture-force combinations. Upper body posture and hand force collected were used as inputs for three different model versions; No-Ligaments (NL) included, Glenohumeral-Ligaments (GH) included, and All-Ligaments (AL) included. Electromyographic (EMG) signals from 11 muscle sites were used for comparison with model MFPs. The primary analysis focused on the differences between the GH and NL versions. Normalized EMG amplitudes were plotted against normalized MFPs from both models. Ligament effects on model outputs were measured by comparing changes in correlation between EMG and MFP, changes in slopes regression lines relating EMG to MFP, and the frequency of zero-force prediction by the model. Paired Student’s t-tests were used to measure significant differences. Results showed significant correlations (Pearson product) between EMG amplitude and MFP in the lower trapezius and infraspinatus muscles (p<0.01). No significant differences were found in r-values for these muscles between the NL and GH model. Slopes of regression lines decreased when GH ligaments were added, while the change in zero-force predictions varied by muscle. This study highlights the sensitivity of musculoskeletal models to the inclusion of ligament forces. Though correlations did not change, decreases in slope indicate increased force prediction by the GH model. Though zero-force predictions for some muscles increased, the results from those that decreased suggest muscles are active in postures where they were originally believed to be inactive. This finding suggests that inclusion of GH ligaments into our model may help predict antagonist muscle activity. However, further research is required.

Mathematical musculoskeletal modeling and simulation provide a means for proactive injury prevention. To be effective, these models must physiologically replicate shoulder function. Although several muscle force prediction (MFP) shoulder models exist, few have attempted to integrate the force contributions of ligaments, especially the glenohumeral ligaments. The purpose of the current study was to integrate seven shoulder ligaments into an existing computational shoulder model, and analyze both individual ligament characteristics and the influence on the model outputs. Using data from the literature, seven shoulder ligaments were integrated into the model: the costoclavicular, conoid, trapezoid, coracohumeral, superior glenohumeral, middle glenohumeral, and inferior glenohumeral. 10 subjects performed isometric exertions in 56 posture-force combinations. Upper body posture and hand force collected were used as inputs for three different model versions; No-Ligaments (NL) included, Glenohumeral-Ligaments (GH) included, and All-Ligaments (AL) included. Electromyographic (EMG) signals from 11 muscle sites were used for comparison with model MFPs. The primary analysis focused on the differences between the GH and NL versions. Normalized EMG amplitudes were plotted against normalized MFPs from both models. Ligament effects on model outputs were measured by comparing changes in correlation between EMG and MFP, changes in slopes regression lines relating EMG to MFP, and the frequency of zero-force prediction by the model. Paired Student’s t-tests were used to measure significant differences. Results showed significant correlations (Pearson product) between EMG amplitude and MFP in the lower trapezius and infraspinatus muscles (p<0.01). No significant differences were found in r-values for these muscles between the NL and GH model. Slopes of regression lines decreased when GH ligaments were added, while the change in zero-force predictions varied by muscle. This study highlights the sensitivity of musculoskeletal models to the inclusion of ligament forces. Though correlations did not change, decreases in slope indicate increased force prediction by the GH model. Though zero-force predictions for some muscles increased, the results from those that decreased suggest muscles are active in postures where they were originally believed to be inactive. This finding suggests that inclusion of GH ligaments into our model may help predict antagonist muscle activity. However, further research is required.

Spinal injections for back-pain management are frequently carried out in hospitals and radiological clinics. Currently, these procedures are performed under fluoroscopy or CT guidance in specialized interventional radiology facilities, and thus incur a major financial burden on the healthcare system. Additionally, the current practice exposes patients and surgeons to X-ray radiation. The use of US for image guided navigation of the spine would greatly reduce the exposure of both the patient and the physician to ionizing radiation and allow the procedure to be performed outside of a specialized facility. However, US as the sole guidance modality has its own challenges. In particular, due to the significant level of occlusion in spinal US images, it can be difficult to accurately identify the appropriate injection site. Here, a groupwise US to CT registration algorithm for guiding percutaneous spinal interventions is presented. In our registration methodology, each vertebra in CT is treated as a sub-volume and transformed individually. A biomechanical model is used to constrain the displacement of the vertebrae relative to one another. The sub-volumes are then reconstructed into a single volume. In each iteration of registration, an US image is simulated from the reconstructed CT volume and an intensity-based similarity metric with the real US image is calculated. Validation studies are performed on datasets from a lamb cadaver, five patient-based phantoms designed to preserve realistic curvatures of the spine and a sixth patient-based phantom where the curvature of the spine is changed between preoperative and intraoperative imaging. For datasets where the spine curve between two imaging modalities was artificially perturbed, the proposed methodology was able to register initial misalignments of up to 20 mm with a success rate of 95%. For the phantom with a physical change in the curvature of the spine introduced between the US and CT datasets, the registration success rate was 98.5%. Finally, the registration success rate for the lamb cadaver with soft tissue information was 87%. The results demonstrate that our algorithm robustly registers US and CT datasets of the spine, regardless of a change in the patients pose between preoperative and intraoperative image acquisitions. / Thesis (Master, Computing) -- Queen's University, 2009-11-27 13:48:33.288

No description available.

Loading vitrifiable concentrations of cryoprotective agents is an important step for cryopreservation of biological tissues by vitrification for research and transplantation purposes. This may be done by immersing the tissue in a cryoprotective agent (CPA) solution, and increasing the concentration, continuously or in multiple steps, and simultaneously decreasing the temperature to decrease the toxicity effects of the cryoprotective agent on the tissue cellular system. During cryoprotective agent loading, osmotic water movement from the tissue to the surrounding solution, and the resultant tissue shrinkage and stress-strain in the tissue matrix as well as on the cellular system can significantly alter the outcome of the cryopreservation protocol. In this thesis, a biomechanical model for articular cartilage is developed to account for the transport of the cryoprotective agent, the nonideal-nondilute properties of the vitrifiable solutions, the osmotic water movement and the resultant tissue shrinkage and stress-strain in the tissue matrix, and the osmotic volume change of the chondrocytes, during cryoprotective agent loading in the cartilage matrix. Four essential transport parameters needed for the model were specified, the values of which were obtained uniquely by fitting the model to experimental data from porcine articular cartilage. Then, it was shown that using real nonuniform initial distributions of water and fixed charges in cartilage, measured separately in this thesis using MRI, in the model can significantly affect the model predictions. The model predictions for dimethyl sulfoxide diffusion in porcine articular cartilage were verified by comparing to spatially and temporally resolved measurements of dimethyl sulfoxide concentration in porcine articular cartilage using a spectral MRI technique, developed for this purpose and novel to the field of cryobiology. It was demonstrated in this thesis that the developed mathematical model provides a novel tool for studying transport phenomena in cartilage during cryopreservation protocols, and can make accurate predictions for the quantities of interest for applications in the cryopreservation of articular cartilage. / Chemical Engineering

Biomedical image processing is a very important research area. Image analysis is one of the most important techniques in studies related to heart functions. The clinical assessment of LV function is very important to evaluate the heart function for patients or suspected heart disease sufferers. 2D echocardiography allows us to study the dynamic analysis of the heart which results in obtaining the quantitative and qualitative analysis of the LV. Cardiac function quantitative analysis depends on the heart’s shape characteristics like the enclosed area and heart wall thickness. The segmentation of medical images and obtaining the traces of the LV boundaries is an essential procedure to get the quantitative and qualitative analysis. Yet, in clinical procedure, this task depends on manual tracing which is slow, tedious and time consuming job. Hence, automating this clinical procedure during the cardiac cycle is of great importance. The aim of this thesis is to automate the manual process of detecting and tracking the LV boundaries of 2D echocardiographic image sequence. Instead of depending only on the imaging based techniques, the designed and implemented framework utilizes the LV mechanics beside the imaging based techniques. When it comes to information extraction from patterns which have been classified, it has been proved that the different contour detection methods complement each other. As a result, efficient combination of different contour detectors is expected to achieve better contour detection than if only one detector is used. This combination of contour detectors produces incremental gains in overall performance. In the first framework, the detection and tracking are accomplished by employing the extended Kalman filter framework to combine the contours estimated by the biomechanical model and the contours extracted using the deformable models. An alternative framework is used by employing averaging fusion followed by level set method. A gold standard is created from three manual outlines and utilized in the experimental results to evaluate the automated results. The tracking and segmentation of LV during the cardiac cycle was accomplished successfully in all cases. The results showed limits of agreement for an average perpendicular distance of 1.277 ±0.252 mm versus the created gold standard. This proved that this framework achieved better performance in tracking and segmenting the LV through the cardiac cycle.

No description available.

Filamentous fungi are important in an enormous variety of ways to our life, with examples ranging from bioremediation, through the food and drinks industry to human health. These organisms can form huge networks stretching metres and even kilometres. However, their mode of growth is by the extension of individual hyphal tips only a few microns in diameter. Tip growth is mediated by the incorporation of new wall building materials at the soft apex. Just how this process is controlled (in fungi and in cell elongation in other organisms) has been the subject of intense study over many years and has attracted considerable attention from mathematical modellers. In this thesis, we consider mathematical models of fungal tip growth that can be classified as either geometrical or biomechanical. In every model we examine, a 2-D axisymmetric semihemisphere-like curve represents half the medial section of fungal tip geometry. A geometrical model for the role of the Spitzenkorper in the tip growth was proposed by Bartnicki-Garcia et al (1989), where a number of problems with the mathematical derivation were pointed out by Koch (2001). A suggestion is given as an attempt to revise the derivation by introducing a relationship between arc length of a growing tip, deposition of wall-building materials and tip curvature. We also consider two types of geometrical models as proposed by Goriely et al (2005). The first type considers a relationship between the longitudinal curvature and the function used to model deposition of wall-building materials. For these types of models, a generalized formulae for the tip shape is introduced, which allows localization of deposition of wall-building materials to be examined. The second type considers a relationship between longitudinal and latitudinal curvatures and the function used to model deposition of wall-building materials. For these types of models, a new formulation of the function used to model deposition of wall-building materials is introduced. Finally, a biomechanical model as proposed by Goriely et al (2010). Varying arc length of the stretchable region on the tip suggests differences in geometry of tip shape and the effective pressure profile. The hypothesis of orthogonal growth is done by focusing only on the apex of a "germ tube". Following that, it suggests that material points on the tip appear to move in a direction perpendicular to the tip either when surface friction is increased or decreased.

Contexte: La mammographie est une modalité d’imagerie médicale à faible dose permettant la détection du cancer mammaire à un stade précoce. Lors de l'examen, le sein est comprimé entre deux plaques afin d'uniformiser son épaisseur et d'étaler les tissus. Cette compression améliore la qualité clinique de l'examen mais elle est également source d'inconfort chez les patientes. Bien que la mammographie soit la méthode de dépistage la plus efficace du cancer du sein, l’inconfort ressenti peut dissuader les femmes de passer cet examen. Par conséquent, une technique alternative de compression du sein prenant en compte le confort de la patiente, en plus de l’amélioration de la qualité d'image, présente un grand intérêt.Méthodes: Dans ce travail, nous avons proposé un nouvel environnement de simulation permettant l'évaluation de différentes techniques de compression du sein. La qualité de la compression a été caractérisée en termes de confort de la patiente, de la qualité d'image et de la dose glandulaire moyenne délivrée. Afin d'évaluer la déformation du sein lors de la compression, un modèle biomécanique par éléments finis du sein a été développé. Ce dernier a été calibré et évalué en utilisant des volumes IRM d'une volontaire dans trois configurations différentes (sur le dos, le ventre et de côté). Par ailleurs, la qualité d'image a été évaluée en utilisant un environnement de simulation d'imagerie auparavant validé pour la simulation de l'acquisition d'images en mammographie.Résultats: La capacité de notre modèle biomécanique à reproduire les déformations réelles des tissus a été évaluée. Tout d'abord, la géométrie du sein dans les trois configurations a été estimée en utilisant des matériaux Néo-Hookeens pour la modélisation des tissus mous. Les propriétés mécaniques des différents constituants du sein ont été estimés afin que les géométries du sein dans les positions couchée sur le ventre et couchée soient le plus proches possibles des mesures. La distance de Hausdorff entre les données estimées et les données mesurées est égale à 2.17 mm en position couché sur le ventre et 1.72 mm en position couché sur le dos. Le modèle a ensuite été évalué dans une troisième configuration sur le côté, avec une distance de Hausdorff étant alors égale à 6.14 mm. Cependant, nous avons été montré que le modèle Néo-Hookeen ne peut pas décrire intégralement le comportement mécanique riche des tissus mous. Par conséquent, nous avons introduit d'autres modèles de matériaux basés sur la fonction d'énergie de Gent. Cette dernière hypothèse a permis de réduire l'erreur maximale dans la configuration couchée sur le ventre et dos incliné d’environ 10 mm.Le couplage entre la simulation de la mécanique du sein et la simulation d'aquisition d'image nous ont permis d'effectuer deux études préliminaires. Dans la première étude, les différences entre les pelotes de compression standard rigide et flex ont été évaluées. Selon les simulations effectuées, l'utilisation de la pelote flex pour la compression du sein a le potentiel d'améliorer le confort de la patiente sans affecter la qualité de l'image ou la dose glandulaire moyenne.Dans la seconde étude, l'impact du positionnement du sein sur la mécanique globale de la compression mammaire a été étudié. Nos simulations confirment que rapprocher la pelote de compression de la cage thoracique peut augmenter l'inconfort de la patiente. Selon les données estimées, pour une même épaisseur du sein sous compression, la force appliqée au sein peut être s'accroitre de 150%.Conclusion: L'estimation réaliste de la géométrie du sein pour différentes configurations sous l'effet de la gravité, ainsi que les résultats conformes aux descriptions cliniques sur la compression du sein, ont confirmé l'interêt d'un environnement de simulation dans le cadre de nos études. / Background: Mammography is a specific type of breast imaging that uses low-dose X-rays to detect breast cancer at early stage. During the exam, the women breast is compressed between two plates in order to even out the breast thickness and to spread out the soft tissues. This compression improves the exam quality but can also be a source of discomfort. Though mammography is the most effective breast cancer screening method, the discomfort perceived during the exam might deter women from getting the test. Therefore, an alternative breast compression technique considering the patient comfort in addition to an improved clinical image quality is of large interest.Methods: In this work, a simulation environment allowing the evaluation of different breast compression techniques was put forward. The compression quality was characterized in terms of patient comfort, image quality and average glandular dose. To estimate the breast deformation under compression, a subject-specific finite element biomechanical model was developed. The model was calibrated and evaluated using MR images in three different breast configurations (supine, prone and supine tilted). On the other hand, image quality was assessed by using an already validated simulation framework. This framework was largely used to mimic image acquisitions in mammography.Findings: The capability of our breast biomechanical model to reproduce the real breast deformations was evaluated. To this end, the geometry estimates of the three breast configurations were computed using Neo-Hookean material models. The subject specific mechanical properties of each breast's structures were assessed, such as the best estimates of the supine and prone configurations were obtained. The Hausdorff distances between the estimated and the measured geometries were equal to 2.17 mm and 1.72 mm respectively. Then, the model was evaluated using a supine tilted configuration; with a Hausdorff distance of 6.14 mm was obtained in that case. However, we have showed that the Neo-Hookean strain energy function cannot totally describe the rich mechanical behavior of breast soft tissues. Therefore, alternative material models based on the Gent strain energy function were proposed. The latter assumption reduced the maximal error in supine tilted breast configuration by about 10 mm.The coupling between the simulations of the breast mechanics and the X-ray simulations allowed us to run two preliminary studies. In the first study, the differences between standard rigid and flex compression paddles were assessed. According to the performed simulations, using the flex paddle for breast compression may improve the patient comfort without affecting the image quality and the delivered average glandular dose.In the second study, the impact of breast positioning on the general compression mechanics was described. Our simulations confirm that positioning the paddle closer to the chest wall is suspected to increase the patient discomfort. Indeed, based on the estimated data, for the same breast thickness under compression, the force applied to the breast may increase by 150%.Conclusion: The good results we get for the estimation of breast deformation under gravity, as well as the conforming results on breast compression quality with the already published clinical statements, have shown the feasibility of such studies by the means of a simulation framework.

The ability to estimate body mass from human skeletal remains with a high degree of accuracy would be significant for the identification of identifying unknown individuals in a forensic anthropology context, documenting secular change in modern populations, and evaluating any prevalence in prehistoric populations. Modern research investigating body mass incorporates one of two models: morphometric and biomechanical. The morphometric model views the body as a cylinder, where weight estimates are gathered from extreme points such as the breadth of the pelvis. In contrast, the biomechanical model incorporates engineering principles and biology to understand the effects of mass on the human skeleton. Only the biomechanical model can accommodate extremes in body mass, such as those exhibited by modern populations. This study examined the accuracy of estimating body mass (obesity in particular) from human skeletal remains using a suite of traits shown to be significant in previous studies, including documented biomechanical analysis of obese individuals involving gait and sit-to-stand (STS) movements. It was hypothesized that using a combination of methods, body mass could be estimated with a high degree of accuracy. Using a large skeletal sample (n = 191), composed of male and females with documented age, weight, and height, the following three variables were examined: (1) the spinal manifestation of diffuse idiopathic skeletal hyperostosis (DISH), (2) osteoarthritis (OA) of the tibiae, and (3) external femoral dimensions. These were then subject to statistical tests. Spearman's rank-order correlation and Mann-Whitney U tests showed significant relationships between DISH and obesity in females (p<.05), but not for males. The presence and severity of OA of the medial condyles were also significantly related to BMI in females (p<.05). In males, the relationship between BMI and OA was only significant on the condyles of the right tibiae (p<.05). Finally, ANOVA and Pearson's product-moment correlation tests were performed to evaluate the cross-sectional dimensions of the femur. The effect of age, stature, and BMI were also examined. ANOVA results showed a significant effect between BMI and M-L cross-sectional dimensions among both sexes (p<.05). Initial Pearson's tests performed separately on males and females showed no significant correlations; however, after the sexes were pooled, small to moderate negative correlations between the M-L/A-P ratio along the diaphysis of the femur and BMI were found. Finally, multiple regression analyses were performed. The models for both sexes with all ten variables was statistically significant for BMI. The final accuracy rate was 78.48% for females and 84.37% for males. The primary goal of this study was to evaluate Moore's (2008) body mass estimation study. In this investigation, however, all dimensions of the femur were performed using an osteometric board and sliding calipers following the guidelines used by Agostini and Ross (2011). The results of this study paralleled many of the observations seen in previous studies, particularly the M-L lateral widening of the femur. Future research should continue to examine the relationship of DISH and OA with body mass, particularly regarding the varying manifestations between the sexes and confounding factors such as age.