Rendering Deformed Speckle Images with a Boolean Model

International audienceRendering speckle images affected by a given deformation field is of primary importance to assess the metrological performance of displacement measurement methods used in experimental mechanics and based on digital image correlation (DIC). This article describes how to render deformed speckle images with a classic model of stochastic geometry, the Boolean model. The advantage of the proposed approach is that it does not depend on any interpolation scheme likely to bias the assessment process, and that it allows the user to render speckle images deformed with any deformation field given by an analytic formula. The proposed algorithm mimics the imaging chain of a digital camera, and its parameters are carefully discussed. A MATLAB software implementation and synthetic ground-truth datasets for assessing DIC software programs are publicly available

Improvements in digital image correlation and application in material mechanical test

Digital image correlation (DIC) is a non-contact full-field optical measurement method. With the advantages of high accuracy, low cost, and simple implementation, it has been widely applied in the area of experimental mechanics. In this study, DIC algorithm has been improved in the aspects of the pixel-level searching method and reference frame update strategy. The feature matching based method is proposed to provide an initial guess for all points of interest with semi-subpixel level accuracy in cases with small or large translation, deformation, or rotation. The bisection searching strategy is presented to automatically adjust the frame step for varying practical circumstances. The improved DIC algorithm is implemented and applied to the miniature tensile test. A convenient experimental method to determine the true stress-strain curve is proposed. The instantaneous cross-section area is estimated by only one camera in aid of DIC method. The derived true stress-strain curves and mechanical parameters of metal material Al6061 and CP-Ti from miniature specimens match well with the results of standard specimens, and no dimension dependence has been observed in the results --Abstract, page iv

A Markov Random Field Based Approach to 3D Mosaicing and Registration Applied to Ultrasound Simulation

A novel Markov Random Field (MRF) based method for the mosaicing of 3D ultrasound volumes is presented in this dissertation. The motivation for this work is the production of training volumes for an affordable ultrasound simulator, which offers a low-cost/portable training solution for new users of diagnostic ultrasound, by providing the scanning experience essential for developing the necessary psycho-motor skills. It also has the potential for introducing ultrasound instruction into medical education curriculums. The interest in ultrasound training stems in part from the widespread adoption of point-of-care scanners, i.e. low cost portable ultrasound scanning systems in the medical community. This work develops a novel approach for producing 3D composite image volumes and validates the approach using clinically acquired fetal images from the obstetrics department at the University of Massachusetts Medical School (UMMS). Results using the Visible Human Female dataset as well as an abdominal trauma phantom are also presented. The process is broken down into five distinct steps, which include individual 3D volume acquisition, rigid registration, calculation of a mosaicing function, group-wise non-rigid registration, and finally blending. Each of these steps, common in medical image processing, has been investigated in the context of ultrasound mosaicing and has resulted in improved algorithms. Rigid and non-rigid registration methods are analyzed in a probabilistic framework and their sensitivity to ultrasound shadowing artifacts is studied. The group-wise non-rigid registration problem is initially formulated as a maximum likelihood estimation, where the joint probability density function is comprised of the partially overlapping ultrasound image volumes. This expression is simplified using a block-matching methodology and the resulting discrete registration energy is shown to be equivalent to a Markov Random Field. Graph based methods common in computer vision are then used for optimization, resulting in a set of transformations that bring the overlapping volumes into alignment. This optimization is parallelized using a fusion approach, where the registration problem is divided into 8 independent sub-problems whose solutions are fused together at the end of each iteration. This method provided a speedup factor of 3.91 over the single threaded approach with no noticeable reduction in accuracy during our simulations. Furthermore, the registration problem is simplified by introducing a mosaicing function, which partitions the composite volume into regions filled with data from unique partially overlapping source volumes. This mosaicing functions attempts to minimize intensity and gradient differences between adjacent sources in the composite volume. Experimental results to demonstrate the performance of the group-wise registration algorithm are also presented. This algorithm is initially tested on deformed abdominal image volumes generated using a finite element model of the Visible Human Female to show the accuracy of its calculated displacement fields. In addition, the algorithm is evaluated using real ultrasound data from an abdominal phantom. Finally, composite obstetrics image volumes are constructed using clinical scans of pregnant subjects, where fetal movement makes registration/mosaicing especially difficult. Our solution to blending, which is the final step of the mosaicing process, is also discussed. The trainee will have a better experience if the volume boundaries are visually seamless, and this usually requires some blending prior to stitching. Also, regions of the volume where no data was collected during scanning should have an ultrasound-like appearance before being displayed in the simulator. This ensures the trainee\u27s visual experience isn\u27t degraded by unrealistic images. A discrete Poisson approach has been adapted to accomplish these tasks. Following this, we will describe how a 4D fetal heart image volume can be constructed from swept 2D ultrasound. A 4D probe, such as the Philips X6-1 xMATRIX Array, would make this task simpler as it can acquire 3D ultrasound volumes of the fetal heart in real-time; However, probes such as these aren\u27t widespread yet. Once the theory has been introduced, we will describe the clinical component of this dissertation. For the purpose of acquiring actual clinical ultrasound data, from which training datasets were produced, 11 pregnant subjects were scanned by experienced sonographers at the UMMS following an approved IRB protocol. First, we will discuss the software/hardware configuration that was used to conduct these scans, which included some custom mechanical design. With the data collected using this arrangement we generated seamless 3D fetal mosaics, that is, the training datasets, loaded them into our ultrasound training simulator, and then subsequently had them evaluated by the sonographers at the UMMS for accuracy. These mosaics were constructed from the raw scan data using the techniques previously introduced. Specific training objectives were established based on the input from our collaborators in the obstetrics sonography group. Important fetal measurements are reviewed, which form the basis for training in obstetrics ultrasound. Finally clinical images demonstrating the sonographer making fetal measurements in practice, which were acquired directly by the Philips iU22 ultrasound machine from one of our 11 subjects, are compared with screenshots of corresponding images produced by our simulator

A Markov Random Field Based Approach to 3D Mosaicing and Registration Applied to Ultrasound Simulation

A novel Markov Random Field (MRF) based method for the mosaicing of 3D ultrasound volumes is presented in this dissertation. The motivation for this work is the production of training volumes for an affordable ultrasound simulator, which offers a low-cost/portable training solution for new users of diagnostic ultrasound, by providing the scanning experience essential for developing the necessary psycho-motor skills. It also has the potential for introducing ultrasound instruction into medical education curriculums. The interest in ultrasound training stems in part from the widespread adoption of point-of-care scanners, i.e. low cost portable ultrasound scanning systems in the medical community. This work develops a novel approach for producing 3D composite image volumes and validates the approach using clinically acquired fetal images from the obstetrics department at the University of Massachusetts Medical School (UMMS). Results using the Visible Human Female dataset as well as an abdominal trauma phantom are also presented. The process is broken down into five distinct steps, which include individual 3D volume acquisition, rigid registration, calculation of a mosaicing function, group-wise non-rigid registration, and finally blending. Each of these steps, common in medical image processing, has been investigated in the context of ultrasound mosaicing and has resulted in improved algorithms. Rigid and non-rigid registration methods are analyzed in a probabilistic framework and their sensitivity to ultrasound shadowing artifacts is studied. The group-wise non-rigid registration problem is initially formulated as a maximum likelihood estimation, where the joint probability density function is comprised of the partially overlapping ultrasound image volumes. This expression is simplified using a block-matching methodology and the resulting discrete registration energy is shown to be equivalent to a Markov Random Field. Graph based methods common in computer vision are then used for optimization, resulting in a set of transformations that bring the overlapping volumes into alignment. This optimization is parallelized using a fusion approach, where the registration problem is divided into 8 independent sub-problems whose solutions are fused together at the end of each iteration. This method provided a speedup factor of 3.91 over the single threaded approach with no noticeable reduction in accuracy during our simulations. Furthermore, the registration problem is simplified by introducing a mosaicing function, which partitions the composite volume into regions filled with data from unique partially overlapping source volumes. This mosaicing functions attempts to minimize intensity and gradient differences between adjacent sources in the composite volume. Experimental results to demonstrate the performance of the group-wise registration algorithm are also presented. This algorithm is initially tested on deformed abdominal image volumes generated using a finite element model of the Visible Human Female to show the accuracy of its calculated displacement fields. In addition, the algorithm is evaluated using real ultrasound data from an abdominal phantom. Finally, composite obstetrics image volumes are constructed using clinical scans of pregnant subjects, where fetal movement makes registration/mosaicing especially difficult. Our solution to blending, which is the final step of the mosaicing process, is also discussed. The trainee will have a better experience if the volume boundaries are visually seamless, and this usually requires some blending prior to stitching. Also, regions of the volume where no data was collected during scanning should have an ultrasound-like appearance before being displayed in the simulator. This ensures the trainee\u27s visual experience isn\u27t degraded by unrealistic images. A discrete Poisson approach has been adapted to accomplish these tasks. Following this, we will describe how a 4D fetal heart image volume can be constructed from swept 2D ultrasound. A 4D probe, such as the Philips X6-1 xMATRIX Array, would make this task simpler as it can acquire 3D ultrasound volumes of the fetal heart in real-time; However, probes such as these aren\u27t widespread yet. Once the theory has been introduced, we will describe the clinical component of this dissertation. For the purpose of acquiring actual clinical ultrasound data, from which training datasets were produced, 11 pregnant subjects were scanned by experienced sonographers at the UMMS following an approved IRB protocol. First, we will discuss the software/hardware configuration that was used to conduct these scans, which included some custom mechanical design. With the data collected using this arrangement we generated seamless 3D fetal mosaics, that is, the training datasets, loaded them into our ultrasound training simulator, and then subsequently had them evaluated by the sonographers at the UMMS for accuracy. These mosaics were constructed from the raw scan data using the techniques previously introduced. Specific training objectives were established based on the input from our collaborators in the obstetrics sonography group. Important fetal measurements are reviewed, which form the basis for training in obstetrics ultrasound. Finally clinical images demonstrating the sonographer making fetal measurements in practice, which were acquired directly by the Philips iU22 ultrasound machine from one of our 11 subjects, are compared with screenshots of corresponding images produced by our simulator

Assisting digital volume correlation with mechanical image-based modeling: application to the measurement of kinematic fields at the architecture scale in cellular materials

La mesure de champs de déplacement et de déformation aux petites échelles dans des microstructures complexes représente encore un défi majeur dans le monde de la mécanique expérimentale. Ceci est en partie dû aux acquisitions d'images et à la pauvreté de la texture à ces échelles. C'est notamment le cas pour les matériaux cellulaires lorsqu'ils sont imagés avec des micro-tomographes conventionnels et qu'ils peuvent être sujets à des mécanismes de déformation complexes. Comme la validation de modèles numériques et l'identification des propriétés mécaniques de matériaux se base sur des mesures précises de déplacements et de déformations, la conception et l'implémentation d'algorithmes robustes et fiables de corrélation d'images semble nécessaire. Lorsque l'on s'intéresse à l'utilisation de la corrélation d'images volumiques (DVC) pour les matériaux cellulaires, on est confronté à un paradoxe: l'absence de texture à l'échelle du constituant conduit à considérer l'architecture comme marqueur pour la corrélation. Ceci conduit à l'échec des techniques ordinaires de DVC à mesurer des cinématiques aux échelles subcellulaires en lien avec des comportements mécaniques locaux complexes tels que la flexion ou le flambement de travées. L'objectif de cette thèse est la conception d'une technique de DVC pour la mesure de champs de déplacement dans des matériaux cellulaires à l'échelle de leurs architectures. Cette technique assiste la corrélation d'images par une régularisation élastique faible en utilisant un modèle mécanique généré automatiquement et basé sur les images. La méthode suggérée introduit une séparation d'échelles au dessus desquelles la DVC est dominante et en dessous desquelles elle est assistée par le modèle mécanique basé sur l'image. Une première étude numérique consistant à comparer différentes techniques de construction de modèles mécaniques basés sur les images est conduite. L'accent est mis sur deux méthodes de calcul particulières: la méthode des éléments finis (FEM) et la méthode des cellules finies (FCM) qui consiste à immerger la géométrie complexe dans une grille régulière de haut ordre sans utiliser de mailleurs. Si la FCM évite une première phase délicate de discrétisation, plusieurs paramètres restent néanmoins délicats à fixer. Dans ce travail, ces paramètres sont ajustés afin d'obtenir (a) la meilleure précision (bornée par les erreurs de pixellisation) tout en (b) assurant une complexité minimale. Pour l'aspect mesure par corrélation d'images régularisée, plusieurs expérimentations virtuelles à partir de différentes simulations numériques (en élasticité, en plasticité et en non-linéarité géométrique) sont d'abord réalisées afin d'analyser l'influence des paramètres de régularisation introduits. Les erreurs de mesures peuvent dans ce cas être quantifiées à l'aide des solutions de référence éléments finis. La capacité de la méthode à mesurer des cinématiques complexes en absence de texture est démontrée pour des régimes non-linéaires tels que le flambement. Finalement, le travail proposé est généralisé à la corrélation volumique des différents états de déformation du matériau et à la construction automatique de la micro-architecture cellulaire en utilisant soit une grille B-spline d'ordre arbitraire (FCM) soit un maillage éléments finis (FEM). Une mise en évidence expérimentale de l'efficacité et de la justesse de l'approche proposée est effectuée à travers de la mesure de cinématiques complexes dans une mousse polyuréthane sollicitée en compression lors d'un essai in situ.Measuring displacement and strain fields at low observable scales in complex microstructures still remains a challenge in experimental mechanics often because of the combination of low definition images with poor texture at this scale. The problem is particularly acute in the case of cellular materials, when imaged by conventional micro-tomographs, for which complex highly non-linear local phenomena can occur. As the validation of numerical models and the identification of mechanical properties of materials must rely on accurate measurements of displacement and strain fields, the design and implementation of robust and faithful image correlation algorithms must be conducted. With cellular materials, the use of digital volume correlation (DVC) faces a paradox: in the absence of markings of exploitable texture on/or in the struts or cell walls, the available speckle will be formed by the material architecture itself. This leads to the inability of classical DVC codes to measure kinematics at the cellular and a fortiori sub-cellular scales, precisely because the interpolation basis of the displacement field cannot account for the complexity of the underlying kinematics, especially when bending or buckling of beams or walls occurs. The objective of the thesis is to develop a DVC technique for the measurement of displacement fields in cellular materials at the scale of their architecture. The proposed solution consists in assisting DVC by a weak elastic regularization using an automatic image-based mechanical model. The proposed method introduces a separation of scales above which DVC is dominant and below which it is assisted by image-based modeling. First, a numerical investigation and comparison of different techniques for building automatically a geometric and mechanical model from tomographic images is conducted. Two particular methods are considered: the finite element method (FEM) and the finite-cell method (FCM). The FCM is a fictitious domain method that consists in immersing the complex geometry in a high order structured grid and does not require meshing. In this context, various discretization parameters appear delicate to choose. In this work, these parameters are adjusted to obtain (a) the best possible accuracy (bounded by pixelation errors) while (b) ensuring minimal complexity. Concerning the ability of the mechanical image-based models to regularize DIC, several virtual experimentations are performed in two-dimensions in order to finely analyze the influence of the introduced regularization lengths for different input mechanical behaviors (elastic, elasto-plastic and geometrically non-linear) and in comparison with ground truth. We show that the method can estimate complex local displacement and strain fields with speckle-free low definition images, even in non-linear regimes such as local buckling. Finally a three-dimensional generalization is performed through the development of a DVC framework. It takes as an input the reconstructed volumes at the different deformation states of the material and constructs automatically the cellular micro-architeture geometry. It considers either an immersed structured B-spline grid of arbitrary order or a finite-element mesh. An experimental evidence is performed by measuring the complex kinematics of a polyurethane foam under compression during an in situ test

Three-dimensional geometry characterization using structured light fields

Tese de doutoramento. Engenharia Mecânica. Faculdade de Engenharia. Universidade do Porto. 200

3D Full Field Displacement and Strain Measurements at the Microscale in Fiber Reinforced Composites Under Transverse Load Using Digital Image Correlation

RÉSUMÉ La rentabilité de l’industrie aérospatiale est fortement liée à une réduction de la consommation de carburant et de la masse des structures utilisées tout en maintenant un niveau de performance et de sécurité similaire. Pour atteindre ces objectifs, les matériaux composites à renfort fibreux (CRF) sont de plus en plus utilisés dans cette industrie. Ces matériaux offrent une résistance spécifique élevée, résistent à la corrosion, sont légers par rapport à leurs propriétés mécaniques et ont de bonnes propriétés de résistance à la fatigue. Ces matériaux sont de plus en plus utilisés mais sont aussi considérés comme une des causes de plusieurs catastrophes aériennes. Les capacités de prédictions de l’état d’endommagement et de rupture de pièces faites en CRF sont encore limitées, tel que démontré par les deux éditions terminées du World Wide Failure Exercise. L’endommagement dans les CRF est caractérisé par une multitude d’événements microscopiques qui se développent puis se regroupent graduellement jusqu’à former un large réseau de micro fissures à travers le matériau. À l’échelle de la fibre, le type d’endommagement le plus critique pour des CRF unidirectionnels est la décohésion inter-faciale entre les fibres et la matrice d’après la littérature. Ce mécanisme commence avec une décohésion inter-faciale en Mode I entre la fibre et la matrice, la fissure inter-faciale se propage ensuite autour de la fibre dans et hors plan. L’interface entre la fibre et la matrice ne peut alors plus transférer de contraintes correctement, ce qui entraîne une augmentation locale de la contrainte autour de la fibre en décohésion. Une des fibres environnantes va ensuite à son tour avoir une décohésion inter-faciale qui se produit dû à l’augmentation de contrainte, et ainsi de suite pour toutes les fibres environnantes. La rupture de l’échantillon se produit éventuellement lorsque les fissures de décohésion commencent à croître dans la matrice, se regroupent et forment un large réseau de fissures qui grandit à travers tout le spécimen. Des observations de décohésion inter-faciale sont disponibles dans la littérature, mais il n’y a pas encore de modèle pour ce mécanisme qui est généralement accepté. La décohésion inter-faciale implique une croissance de la fissure en Mode I, en Mode II et en mode mixte. D’autres mécanismes s’ajoutent aussi, tels que la friction entre la fibre et la matrice, les contraintes résiduelles dues à la cuisson de la matrice et le retrait chimique de la matrice durant la cuisson. Cette combinaison de mécanismes qui participent à la décohésion inter-faciale en font un mécanisme d’endommagement complexe. Des données expérimentales additionnelles, telles que les champs de déplacement ou de déformation in-situ, permettraient de fournir une compréhension plus complète de la décohésion inter-faciale pour différents types de fibres.----------ABSTRACT This thesis aimed at experimentally investigate damage initiation and growth of FRCs under transverse loading at the fiber level, provide in plane and out of plane full field measurements and crack area measurements for different single fiber composites and for a bundle of carbon fibers. Firstly, a single fiber composite specimen is designed and manufactured in such a way that a large fiber, approximately 1 mm in diameter, is under transverse loading during a tensile test. Four specimens were manufactured out of fibers having strong adhesive bonding with epoxies and no adhesive bonding with epoxies combined with an epoxy and a modified epoxy. A stereoscopic Digital Image Correlation (DIC) setup is then used to track 3D displacements and compute in plane strains for a fiber’s free surface and its vicinity. The experimental results showed that inter-facial debonding happened in three steps, an inter-facial crack opened under Mode I at the fiber / matrix interface at the location where an out of plane displacement difference between the fiber and matrix was the highest and where "y was maximum for all specimens. The inter-facial debonding crack then grew under mixed mode around the fiber while it kept protruding out of the matrix. Finally, specimen failure occurred differently for the specimens with no adhesive bonding compared to the strongly bonded ones. For the ones without adhesive bonding, specimen failure was caused by a crack growing under Mode I in the matrix where fibers were horizontally compressed and large out of plane deformation was experienced. Strongly bonded specimens’ failure was also caused by a Mode I crack growing in the matrix but located where the tension is maximal within the inter-facial crack’s free surface. The complete experimental results, containing the stereoscopic full field displacement and strain results for each test and timestep, were provided in a data package for further analysis, or benchmarking of simulation results. DIC provided quantitative information about displacement and strain fields, however, the method has limitations in the vicinity of cracks. In addition, DIC did not provide any quantitative information about cracks themselves. A method using the raw images from the experiment and the DIC results was developed to combine both results and accurately determine the crack area, crack path in the reference coordinates and the exact applied stress on the specimen for the crack to grow through a certain area. This method was applied to the previously mentioned experiments. Results showed that inter-facial crack initiation happens in fact at the same strain value applied on the specimen, whether the fiber has strong or weak adhesive bonding

Tracking and Hands Motion Detection Approach for Monitoring Hand-Hygiene Compliance for Food Handling and Processing Industry

Hand-hygiene is a very critical issue for both food handling and processing industry and health care service providers. Poor hand-hygiene practice can easily lead to foodborne illness or large scale decease transmission. In this research, an automatic tracking and monitoring system was developed that used a 3D camera for hand washing and hands motion detection and a sensor-based monitoring system for hand-hygiene activities evaluation. An active Wi-Fi portable Radio Frequency Identification (RFID) tag was used for personal ID tracking. The effective hand washing time, soaping time were measured based on the hands motion detection and hand movement tracking. Water temperature, water flow, paper towel, soap and hand sanitizer usage were also measured for each hand washing event. All the data were forwarded to a system server for data recording, storage and management. Preliminary test data were collected to evaluate the system performance. The results showed that the system could effectively collect most of the hand-hygiene related factors including hand-hygiene product usage, hand washing time and soap lathering time for hand-hygiene evaluation.Biosystems & Agricultural Engineerin

BIAXIAL STRESS TESTING OF SS-304L MICROTUBES BY AXIAL LOAD AND INTERNAL PRESSURE

The mechanical behavior and material properties of a Stainless Steel SS-304L microtube, with an OD of 2.40 mm and wall thickness of 160 µm, was investigated through uniaxial, isothermal, biaxial, and metallographic testing. The grain structure, microhardness, and tube eccentricity were investigated using optical microscopy. The rate- and temperature-dependence of the material was characterized by isothermal uniaxial tension experiments. A biaxial experimental setup, consisting of a 2 kN electromechanical tensile stage and a 1.4 kbar hydraulic pump, was created to internally pressurize and axially load the microtube in biaxial stress states. Fourteen radial nominal stress path tests were conducted to determine the formability, failure mode, and anisotropy during biaxial stress states. The Yld2000-2D and Yld2004-3D yield functions were fit to the data at the initial yield surface and higher levels of plastic work. The path-dependence of failure stresses and strains was investigated by comparing radial path results to corner paths

Multi-angle valve seat machining: experimental analysis and numerical modelling

Modern automotive manufacturers operate in highly competitive markets, heavily influenced by Government regulation and ever more environmentally conscious consumers. Modern high-temperature, high-pressure engines that use high hardness multi-angle valve seats are an attractive environmental option, but one that manufacturers find requires more advanced materials and tighter geometric tolerances to maintain engine performance.Tool manufacturers meet these increasingly tougher demands by using, higher hardness cutting materials such as polycrystalline cubic boron nitride (pcBN), that on paper, promise to wear at a lower rate, require less coolant and deliver tighter tolerances than their carbide counterparts.The low brittle fracture toughness of pcBN makes tools that use it vulnerable to minute chipping. A review of literature for this work pointed to no clear answer to this problem, although suggestions range from manufacturing defects, dynamic and flexibility problems with the production line machinery and fixtures, and radial imbalances in the cutting loads.This work set about experimentally investigating those potential explanations, coming to the conclusion that the high radial imbalance of the cutting loads is responsible for pcBN cutting insert failure during multi-angle valve seat machining, and that by simply relocating the cutting inserts around the multi angle cutting tool, the imbalance can be reduced, thus extending the life of the cutting inserts.It is not always easy to predict the imbalance due to the multiple flexibilities in the system, and simulating such a system in 3D with all its associated cutting phenomena such as friction, thermal expansion, chip flow and shearing, would call upon extraordinary computational power and extremely precise experimental inputs to reduce cumulative error.This thesis proves that such a 3D simulation can be made, that runs in exceptionally short durations compared to traditional methods, by making a number of simplifications.MSC Marc was used to host the simulation, with a parametric script written in Python responsible for generating the model geometry and cutter layout. A Fortran program was developed that is called upon by Marc to calculate the required cutting load outputs and generate new workpiece meshes as material is removed.</div