Biomechanical modelling of human knee during living activities

The knee joint is a key structure of the human locomotor system. The knowledge of how each single anatomical structure of the knee contributes to determine the physiological function of the knee, is of fundamental importance for the development of new prostheses and novel clinical, surgical, and rehabilitative procedures. In this context, a modelling approach is necessary to estimate the biomechanic function of each anatomical structure during daily living activities. The main aim of this study was to obtain a subject-specific model of the knee joint of a selected healthy subject. In particular, 3D models of the cruciate ligaments and of the tibio-femoral articular contact were proposed and developed using accurate bony geometries and kinematics reliably recorded by means of nuclear magnetic resonance and 3D video-fluoroscopy from the selected subject. Regarding the model of the cruciate ligaments, each ligament was modelled with 25 linear-elastic elements paying particular attention to the anatomical twisting of the fibres. The devised model was as subject-specific as possible. The geometrical parameters were directly estimated from the experimental measurements, whereas the only mechanical parameter of the model, the elastic modulus, had to be considered from the literature because of the invasiveness of the needed measurements. Thus, the developed model was employed for simulations of stability tests and during living activities. Physiologically meaningful results were always obtained. Nevertheless, the lack of subject-specific mechanical characterization induced to design and partially develop a novel experimental method to characterize the mechanics of the human cruciate ligaments in living healthy subjects. Moreover, using the same subject-specific data, the tibio-femoral articular interaction was modelled investigating the location of the contact point during the execution of daily motor tasks and the contact area at the full extension with and without the whole body weight of the subject. Two different approaches were implemented and their efficiency was evaluated. Thus, pros and cons of each approach were discussed in order to suggest future improvements of this methodologies. The final results of this study will contribute to produce useful methodologies for the investigation of the in-vivo function and pathology of the knee joint during the execution of daily living activities. Thus, the developed methodologies will be useful tools for the development of new prostheses, tools and procedures both in research field and in diagnostic, surgical and rehabilitative fields.

Development and characterization of non-oxide ceramic composites for mechanical and tribological applications

The main reasons for the attention focused on ceramics as possible structural materials are their wear resistance and the ability to operate with limited oxidation and ablation at temperatures above 2000°C. Hence, this work is devoted to the study of two classes of materials which can satisfy these requirements: silicon carbide -based ceramics (SiC) for wear applications and borides and carbides of transition metals for ultra-high temperatures applications (UHTCs). SiC-based materials: Silicon carbide is a hard ceramic, which finds applications in many industrial sectors, from heat production, to automotive engineering and metals processing. In view of new fields of uses, SiC-based ceramics were produced with addition of 10-30 vol% of MoSi2, in order to obtain electro conductive ceramics. MoSi2, indeed, is an intermetallic compound which possesses high temperature oxidation resistance, high electrical conductivity (21·10-6 Ω·cm), relatively low density (6.31 g/cm3), high melting point (2030°C) and high stiffness (440 GPa). The SiC-based ceramics were hot pressed at 1900°C with addition of Al2O3-Y2O3 or Y2O3-AlN as sintering additives. The microstructure of the composites and of the reference materials, SiC and MoSi2, were studied by means of conventional analytical techniques, such as X-ray diffraction (XRD), scanning electron microscopy (SEM) and energy dispersive spectroscopy (SEM-EDS). The composites showed a homogeneous microstructure, with good dispersion of the secondary phases and low residual porosity. The following thermo-mechanical properties of the SiC-based materials were measured: Vickers hardness (HV), Young’s modulus (E), fracture toughness (KIc) and room to high temperature flexural strength (σ). The mechanical properties of the composites were compared to those of two monolithic SiC and MoSi2 materials and resulted in a higher stiffness, fracture toughness and slightly higher flexural resistance. Tribological tests were also performed in two configurations disco-on-pin and slideron cylinder, aiming at studying the wear behaviour of SiC-MoSi2 composites with Al2O3 as counterfacing materials. The tests pointed out that the addition of MoSi2 was detrimental owing to a lower hardness in comparison with the pure SiC matrix. On the contrary, electrical measurements revealed that the addition of 30 vol% of MoSi2, rendered the composite electroconductive, lowering the electrical resistance of three orders of magnitude. Ultra High Temperature Ceramics: Carbides, borides and nitrides of transition metals (Ti, Zr, Hf, Ta, Nb, Mo) possess very high melting points and interesting engineering properties, such as high hardness (20-25 GPa), high stiffness (400-500 GPa), flexural strengths which remain unaltered from room temperature to 1500°C and excellent corrosion resistance in aggressive environment. All these properties place the UHTCs as potential candidates for the development of manoeuvrable hypersonic flight vehicles with sharp leading edges. To this scope Zr- and Hf- carbide and boride materials were produced with addition of 5-20 vol% of MoSi2. This secondary phase enabled the achievement of full dense composites at temperature lower than 2000°C and without the application of pressure. Besides the conventional microstructure analyses XRD and SEM-EDS, transmission electron microscopy (TEM) was employed to explore the microstructure on a small length scale to disclose the effective densification mechanisms. A thorough literature analysis revealed that neither detailed TEM work nor reports on densification mechanisms are available for this class of materials, which however are essential to optimize the sintering aids utilized and the processing parameters applied. Microstructural analyses, along with thermodynamics and crystallographic considerations, led to disclose of the effective role of MoSi2 during sintering of Zrand Hf- carbides and borides. Among the investigated mechanical properties (HV, E, KIc, σ from room temperature to 1500°C), the high temperature flexural strength was improved due to the protective and sealing effect of a silica-based glassy phase, especially for the borides. Nanoindentation tests were also performed on HfC-MoSi2 composites in order to extract hardness and elastic modulus of the single phases. Finally, arc jet tests on HfC- and HfB2-based composites confirmed the excellent oxidation behaviour of these materials under temperature exceeding 2000°C; no cracking or spallation occurred and the modified layer was only 80-90 μm thick.

Caratterizzazione meccanica multiscala di materiali rocciosi strutturalmente complessi

The research work was aimed at studying, with a deterministic approach, the relationships between the rock’s texture and its mechanical properties determined at the laboratory scale. The experimentation was performed on a monomineralic crystalline rock, varying in texture, i.e. grains shape. Multi-scale analysis has been adopted to determine the elasto-mechanical properties of the crystals composing the rock and its strength and deformability at the macro-scale. This let us to understand how the structural variability of the investigated rock affects its macromechanical behaviour. Investigations have been performed on three different scales: nano-scale (order of nm), micro-scale (tens of m) and macro-scale (cm). Innovative techniques for rock mechanics, i.e. Depth Sensing Indentation (DSI), have been applied, in order to determine the elasto-mechanical properties of the calcite grains. These techniques have also allowed to study the influence of grain boundaries on the mechanical response of calcite grains by varying the indents’ sizes and to quantify the effect of the applied load on the hardness and elastic modulus of the grain (indentation size effect, ISE). The secondary effects of static indentation Berkovich, Vickers and Knoop were analyzed by SEM, and some considerations on the rock’s brittle behaviour and the effect of microcracks can be made.

Additive Manufacturing by LPBF and WAAM of metals: correlation between production process, microstructure and mechanical properties

In the most recent years, Additive Manufacturing (AM) has drawn the attention of both academic research and industry, as it might deeply change and improve several industrial sectors. From the material point of view, AM results in a peculiar microstructure that strictly depends on the conditions of the additive process and directly affects mechanical properties. The present PhD research project aimed at investigating the process-microstructure-properties relationship of additively manufactured metal components. Two technologies belonging to the AM family were considered: Laser-based Powder Bed Fusion (LPBF) and Wire-and-Arc Additive Manufacturing (WAAM). The experimental activity was carried out on different metals of industrial interest: a CoCrMo biomedical alloy and an AlSi7Mg0.6 alloy processed by LPBF, an AlMg4.5Mn alloy and an AISI 304L austenitic stainless steel processed by WAAM. In case of LPBF, great attention was paid to the influence that feedstock material and process parameters exert on hardness, morphological and microstructural features of the produced samples. The analyses, targeted at minimizing microstructural defects, lead to process optimization. For heat-treatable LPBF alloys, innovative post-process heat treatments, tailored on the peculiar hierarchical microstructure induced by LPBF, were developed and deeply investigated. Main mechanical properties of as-built and heat-treated alloys were assessed and they were well-correlated to the specific LPBF microstructure. Results showed that, if properly optimized, samples exhibit a good trade-off between strength and ductility yet in the as-built condition. However, tailored heat treatments succeeded in improving the overall performance of the LPBF alloys. Characterization of WAAM alloys, instead, evidenced the microstructural and mechanical anisotropy typical of AM metals. Experiments revealed also an outstanding anisotropy in the elastic modulus of the austenitic stainless-steel that, along with other mechanical properties, was explained on the basis of microstructural analyses.

Numerical and experimental analysis of 3D printed smart structures and systems

Three dimensional (3D) printers of continuous fiber reinforced composites, such as MarkTwo (MT) by Markforged, can be used to manufacture such structures. To date, research works devoted to the study and application of flexible elements and CMs realized with MT printer are only a few and very recent. A good numerical and/or analytical tool for the mechanical behavior analysis of the new composites is still missing. In addition, there is still a gap in obtaining the material properties used (e.g. elastic modulus) as it is usually unknown and sensitive to printing parameters used (e.g. infill density), making the numerical simulation inaccurate. Consequently, the aim of this thesis is to present several work developed. The first is a preliminary investigation on the tensile and flexural response of Straight Beam Flexures (SBF) realized with MT printer and featuring different interlayer fiber volume-fraction and orientation, as well as different laminate position within the sample. The second is to develop a numerical analysis within the Carrera' s Unified Formulation (CUF) framework, based on component-wise (CW) approach, including a novel preprocessing tool that has been developed to account all regions printed in an easy and time efficient way. Among its benefits, the CUF-CW approach enables building an accurate database for collecting first natural frequencies modes results, then predicting Young' s modulus based on an inverse problem formulation. To validate the tool, the numerical results are compared to the experimental natural frequencies evaluated using a digital image correlation method. Further, we take the CUF-CW model and use static condensation to analyze smart structures which can be decomposed into a large number of similar components. Third, the potentiality of MT in combination with topology optimization and compliant joints design (CJD) is investigated for the realization of automated machinery mechanisms subjected to inertial loads.

Elastic Propagation in random media: applications to the imaging of volcano structures

High-frequency seismograms contain features that reflect the random inhomogeneities of the earth. In this work I use an imaging method to locate the high contrast small- scale heterogeneity respect to the background earth medium. This method was first introduced by Nishigami (1991) and than applied to different volcanic and tectonically active areas (Nishigami, 1997, Nishigami, 2000, Nishigami, 2006). The scattering imaging method is applied to two volcanic areas: Campi Flegrei and Mt. Vesuvius. Volcanic and seismological active areas are often characterized by complex velocity structures, due to the presence of rocks with different elastic properties. I introduce some modifications to the original method in order to make it suitable for small and highly complex media. In particular, for very complex media the single scattering approximation assumed by Nishigami (1991) is not applicable as the mean free path becomes short. The multiple scattering or diffusive approximation become closer to the reality. In this thesis, differently from the ordinary Nishigami’s method (Nishigami, 1991), I use the mean of the recorded coda envelope as reference curve and calculate the variations from this average envelope. In this way I implicitly do not assume any particular scattering regime for the "average" scattered radiation, whereas I consider the variations as due to waves that are singularly scattered from the strongest heterogeneities. The imaging method is applied to a relatively small area (20 x 20 km), this choice being justified by the small length of the analyzed codas of the low magnitude earthquakes. I apply the unmodified Nishigami’s method to the volcanic area of Campi Flegrei and compare the results with the other tomographies done in the same area. The scattering images, obtained with frequency waves around 18 Hz, show the presence of high scatterers in correspondence with the submerged caldera rim in the southern part of the Pozzuoli bay. Strong scattering is also found below the Solfatara crater, characterized by the presence of densely fractured, fluid-filled rocks and by a strong thermal anomaly. The modified Nishigami’s technique is applied to the Mt. Vesuvius area. Results show a low scattering area just below the central cone and a high scattering area around it. The high scattering zone seems to be due to the contrast between the high rigidity body located beneath the crater and the low rigidity materials located around it. The central low scattering area overlaps the hydrothermal reservoirs located below the central cone. An interpretation of the results in terms of geological properties of the medium is also supplied, aiming to find a correspondence of the scattering properties and the geological nature of the material. A complementary result reported in this thesis is that the strong heterogeneity of the volcanic medium create a phenomenon called "coda localization". It has been verified that the shape of the seismograms recorded from the stations located at the top of the volcanic edifice of Mt. Vesuvius is different from the shape of the seismograms recorded at the bottom. This behavior is justified by the consideration that the coda energy is not uniformly distributed within a region surrounding the source for great lapse time.

Kinematics of the Sicily and Calabria Subduction System from Elastic Block Modeling of GPS Data

We use data from about 700 GPS stations in the EuroMediterranen region to investigate the present-day behavior of the the Calabrian subduction zone within the Mediterranean-scale plates kinematics and to perform local scale studies about the strain accumulation on active structures. We focus attenction on the Messina Straits and Crati Valley faults where GPS data show extentional velocity gradients of ∼3 mm/yr and ∼2 mm/yr, respectively. We use dislocation model and a non-linear constrained optimization algorithm to invert for fault geometric parameters and slip-rates and evaluate the associated uncertainties adopting a bootstrap approach. Our analysis suggest the presence of two partially locked normal faults. To investigate the impact of elastic strain contributes from other nearby active faults onto the observed velocity gradient we use a block modeling approach. Our models show that the inferred slip-rates on the two analyzed structures are strongly impacted by the assumed locking width of the Calabrian subduction thrust. In order to frame the observed local deformation features within the present- day central Mediterranean kinematics we realyze a statistical analysis testing the indipendent motion (w.r.t. the African and Eurasias plates) of the Adriatic, Cal- abrian and Sicilian blocks. Our preferred model confirms a microplate like behaviour for all the investigated blocks. Within these kinematic boundary conditions we fur- ther investigate the Calabrian Slab interface geometry using a combined approach of block modeling and χ2ν statistic. Almost no information is obtained using only the horizontal GPS velocities that prove to be a not sufficient dataset for a multi-parametric inversion approach. Trying to stronger constrain the slab geometry we estimate the predicted vertical velocities performing suites of forward models of elastic dislocations varying the fault locking depth. Comparison with the observed field suggest a maximum resolved locking depth of 25 km.

Behaviour and applications of elastic waves in structures and metamaterials

The present thesis focuses on elastic waves behaviour in ordinary structures as well as in acousto-elastic metamaterials via numerical and experimental applications. After a brief introduction on the behaviour of elastic guided waves in the framework of non-destructive evaluation (NDE) and structural health monitoring (SHM) and on the study of elastic waves propagation in acousto-elastic metamaterials, dispersion curves for thin-walled beams and arbitrary cross-section waveguides are extracted via Semi-Analytical Finite Element (SAFE) methods. Thus, a novel strategy tackling signal dispersion to locate defects in irregular waveguides is proposed and numerically validated. Finally, a time-reversal and laser-vibrometry based procedure for impact location is numerically and experimentally tested. In the second part, an introduction and a brief review of the basic definitions necessary to describe acousto-elastic metamaterials is provided. A numerical approach to extract dispersion properties in such structures is highlighted. Afterwards, solid-solid and solid-fluid phononic systems are discussed via numerical applications. In particular, band structures and transmission power spectra are predicted for 1P-2D, 2P-2D and 2P-3D phononic systems. In addition, attenuation bands in the ultrasonic as well as in the sonic frequency regimes are experimentally investigated. In the experimental validation, PZTs in a pitch-catch configuration and laser vibrometric measurements are performed on a PVC phononic plate in the ultrasonic frequency range and sound insulation index is computed for a 2P-3D phononic barrier in the sonic frequency range. In both cases the numerical-experimental results comparison confirms the existence of the numerical predicted band-gaps. Finally, the feasibility of an innovative passive isolation strategy based on giant elastic metamaterials is numerically proved to be practical for civil structures. In particular, attenuation of seismic waves is demonstrated via finite elements analyses. Further, a parametric study shows that depending on the soil properties, such an earthquake-proof barrier could lead to significant reduction of the superstructure displacement.

Analytical methods for modeling elastic metasurfaces

This dissertation aims at developing advanced analytical tools able to model surface waves propagating in elastic metasurfaces. In particular, four different objectives are defined and pursued throughout this work to enrich the description of the metasurface dynamics. First, a theoretical framework is developed to describe the dispersion properties of a seismic metasurface composed of discrete resonators placed on a porous medium considering part of it fully saturated. Such a model combines classical elasticity theory, Biot’s poroelasticity and an effective medium approach to describe the metasurface dynamics and its coupling with the poroelastic substrate. Second, an exact formulation based on the multiple scattering theory is developed to extend the two-dimensional classical Lamb’s problem to the case of an elastic half-space coupled to an arbitrary number of discrete surface resonators. To this purpose, the incident wavefield generated by a harmonic source and the scattered field generated by each resonator are calculated. The substrate wavefield is then obtained as solutions of the coupled problem due to the interference of the incident field and the multiple scattered fields of the oscillators. Third, the above discussed formulation is extended to three-dimensional contexts. The purpose here is to investigate the dynamic behavior and the topological properties of quasiperiodic elastic metasurfaces. Finally, the multiple scattering formulation is extended to model flexural metasurfaces, i.e., an array of thin plates. To this end, the resonant plates are modeled by means of their equivalent impedance, derived by exploiting the Kirchhoff plate theory. The proposed formulation permits the treatment of a general flexural metasurface, with no limitation on the number of plates and the configuration taken into account. Overall, the proposed analytical tools could pave the way for a better understanding of metasurface dynamics and their implementation in engineered devices.