Análise da influência da concentração de feixeis de fibras de bambu em ensaios de piso

Bamboo has one of the highest growth rates among plants, however, its lignifications (which confers resistance) takes around a few years and, therefore, certain physical characteristics and mechanical, that depend on this process will only be acquired between the three to six years old. In addition, bamboo also has significant density variations in different parts of the stem, both in the radial direction as the axial. In particular the radial direction, where the density found in the inner and outer (near the bark) of a single stem can range on average from 0.5 to 0.9 g/cm3. Thus, the application of bamboo as a floor, there to examine whether both sides of the bamboo (internal and external), provide resistance properties required for that purpose. In this study sought to characterize and quantify the influence of the concentration of fiber bundles in the inner and outer sides of rules or bamboo strips of bamboo flooring through testing service. Analyses performed were based on ASTM D 2394- 83 for wooden floors and derivatives. This was necessary because of the absence of a specific prohibition of the use and testing of floors made of bamboo and its products. The data were analyzed by ball indentation test shooting, test for resistance to abrasion, indentation test for stress / load treadmill test and by indentation loads applied to small areas - test the jump. The results of the tests were extremely friendly bamboo, even this presents considerable differences between the resistances obtained from assay of the cover of the inner and outer face, being comparable with those of many commonly used to manufacture wood flooring. This comparison was made possible by information from technical trials of several floors made with wood

Caracterização de cerâmicas a base de nitreto de silício para aplicações estruturais

The need for development of new materials is a natural process in the companies’ technological point of view, seeking improvements in materials and processes. Specifically, among the materials, ceramic exhibit valuable properties, especially the covalent ceramics which have excellent properties for applications which requires the abrasion resistance, hardness, high temperatures, resistence, etc. being a material that has applications in several areas. Most studies are related to improvement of properties, specially fracture toughness that allows the expansion of its application. Among the most promising ceramic materials are silicon nitride (Si3N4) which has excellent properties. The goal of this work was the development and caracterization of Si3N4-based ceramics, doped with yttrium oxide (Y2O3), rar earth concentrate (CTR2O3) and cerium oxide (CeO2) in the same proportion for the evaluation of properties. The powders' mixtures were homogenized, dried and compressed under pressure uniaxial and isostatic. Sintering was carried out in 1850 ⁰C under pressure of 0,1MPa N2 for 1 h with a heating rate of 25 ⁰C / min and cooling in the furnace inertia. The characterizations were performed using Archimedes principle to relative density, weight loss by measuring before and after sintering, phase analysis by X-ray diffraction, microstructure by scanning electron microscope (SEM), hardness and fracture toughness by the method Vickers indentation. The results obtained showed relative density of 97-98%, Vickers hardness 17 to 19 GPa, fracture toughness 5.6 to 6.8 MPa.m1/2, with phases varying from α-SiAlON and β-Si3N4 depending the types of additives used. The results are promising for tribological applications and can be defined according to the types of additives to be used

Reproducibility of microhandeness measurements in different surfices

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Effect of partially demineralized dentin beneath the hybrid layer on dentin-adhesive interface micromechanics

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Avaliação da tenacidade à fratura de diferentes sistemas cerâmicos

Although ceramics present high compressive strength, they are brittle materials due to their low tensile strength so they have lower capacity to absorb shocks. This study evaluated the fracture toughness of different ceramic systems, which refers to the ability of a friable material to absorb defformation energy. Three ceramic systems were investigated. Ten cylindrical samples (5,0mm x 3,0mm), were obtained from each ceramic material as follows: G1- 10 samples of Vitadur Alpha (Vita-Zahnfabrik); G2- 10 samples of IPS Empress2 (Ivoclar-Vivadent); G3- 10 samples of In-Ceram Alumina (Vita-Zahnfabrik). Fracture toughness values were collected upon indentation tests that were performed under a heavy load. A microhardness tester (Digital Microhardness Tester FM) utilized a 500gf load cell during 10seconds to perform four impressions on each sample. Statistically significant results were observed (ANOVA and Kruskal-Wallis tests). In-Ceram Alumina presented the highest median toughness values (2,96N/m3/2), followed by Vitadur Alpha (2,08N/m3/2) and IPS Empress2 (1,05N/m3/2). It may be concluded that different ceramic systems present distinct fracture toughness values, thus In-Ceram is capable of absorbing superior stress when compared to Vitadur Alpha and IPS Empress2.

A comparison of microhardness of indirect composite restorative materials

The purpose of this study was to compare the microhardness of four indirect composite resins. Forty cylindrical samples were prepared according to the manufacturer s recommendations using a Teflon mold. Ten specimens were produced from each tested material, constituting four groups (n=10) as follows: G1 - Artglass; G2 - Sinfony; G3 - Solidex; G4 - Targis. Microhardness was determined by the Vickers indentation technique with a load of 300g for 10 seconds. Four indentations were made on each sample, determining the mean microhardness values for each specimen. Descriptive statistics data for the experimental conditions were: G1 - Artglass (mean ±standard deviation: 55.26 ± 1.15HVN; median: 52.6); G2 - Sinfony (31.22 ± 0.65HVN; 31.30); G3 - Solidex (52.25 ± 1.55HVN; 52.60); G4 - Targis (72.14 ± 2.82HVN; 73.30). An exploratory data analysis was performed to determine the most appropriate statistical test through: (I) Levene's for homogeneity of variances; (II) ANOVA on ranks (Kruskal-Wallis); (III) Dunn's multiple comparison test (0.05). Targis presented the highest microhardness values while Sinfony presented the lowest. Artglass and Solidex were found as intermediate materials. These results indicate that distinct mechanical properties may be observed at specific materials. The composition of each material as well as variations on polymerization methods are possibly responsibles for the difference found in microhardness. Therefore, indirect composite resin materials that guarantee both good esthetics and adequate mechanical properties may be considered as substitutes of natural teeth.

Desenvolvimento de um durômetro portátil para madeiras com o uso de transdutor de deslocamento

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

Sinais digitais de um durômetro portátil para a estimativa de propriedades de madeiras para dormentes

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

Relationship between elastic and mechanical properties of dental ceramics and their index of brittleness

The purpose of the study was to verify the effects of a number of materials' parameters (crystalline content; Young's modulus, E; biaxial flexure strength, sigma(i); Vickers hardness, VH; fracture toughness, K-Ic; fracture surface energy, gamma(f); and index of brittleness, B) on the brittleness of dental ceramics. Five commercial dental ceramics with different contents of glass phase and crystalline particles were studied: a vitreous porcelain (VM7/V), a porcelain with 16 vol% leucite particles (d.Sign/D), a glass-ceramic with 29 vol% leucite particles (Empress/E1), a glass-ceramic with 58 vol% lithium-disilicate needle-like particles (Empress 2/E2), and a glass-infiltrated alumina composite with 65 vol% crystals (In-Ceram Alumina/IC). Discs were constructed according to manufacturers' instructions, ground and polished to final dimensions (12 mm x 1.1 mm). Elastic constants were determined by ultrasonic pulse-echo method. sigma(i) was determined by piston-on-3-balls method in inert condition. VH was determined using 19.6 N load and K-Ic was determined by indentation strength method. gamma(f) was calculated from the Griffith-Irwin relation and B by the ratio of HV to K-Ic. IC and E2 showed higher values of sigma(i), E, K-Ic and gamma(f), and lower values of B compared to leucite-based glass-ceramic and porcelains. Positive correlations were observed for sigma(i) versus K-Ic, and K-Ic versus E-1/2, however, E did not show relationship with HV and B. The increase of crystalline phase content is beneficial to decrease the brittleness of dental ceramics by means of both an increase in fracture surface energy and a lowering in index of brittleness. (C) 2012 Elsevier Ltd and Techna Group Sri. All rights reserved.

Analysis, manufacture and characterization of Ni/Cu functionally graded structures

In this work, an experimental and numerical analysis and characterization of functionally graded structures (FGSs) is developed. Nickel (Ni) and copper (Cu) materials are used as basic materials in the numerical modeling and experimental characterization. For modeling, a MATLAB finite element code is developed, which allows simulation of harmonic and modal analysis considering the graded finite element formulation. For experimental characterization, Ni-Cu FGSs are manufactured by using spark plasma sintering technique. Hardness and Young's modulus are found by using microindentation and ultrasonic measurements, respectively. The effective gradation of Ni/Cu FGS is addressed by means of optical microscopy, energy dispersive spectrometry, scanning electron microscopy and hardness testing. For the purpose of comparing modeling and experimental results, the hardness curve, along the gradation direction, is used for identifying the gradation profile; accordingly, the experimental hardness curve is used for approximating the Young's modulus variation and the graded finite element modeling is used for verification. For the first two resonance frequency values, a difference smaller than 1% between simulated and experimental results is obtained. (C) 2012 Elsevier Ltd. All rights reserved.

Evaluation of the micro-hardness and fracture toughness of amorphous and partially crystallized 3CaO center dot P2O5-SiO2-MgO bioglasses

In this work, the effect of the indentation load on the results of hardness and fracture toughness, determined by Vickers micro-hardness measurements, of some glasses and glass-ceramics has been investigated. Furthermore, in order to verify the effect of crystallinity on the results, glasses of composition 52.75 wt.% 3CaO center dot P2O5, 30 wt.% SiO2 and 17.25 wt.% MgO were fused at 1600 degrees C for 4 h and annealed at 700 degrees C for 2h, and further heat-treated at 700, 775, 800 and 900 degrees C for 4h. The obtained materials were analyzed by high resolution X-ray diffraction, HRXRD, to determine the crystallization degree in function of the heat-treatment temperature. The hardness of the different specimens was determined by Vickers' micro-hardness measurements under various loads. It has been observed that with increasing crystallization of the materials their hardness increased. Furthermore, it has been possible to verify the so-called indentation size effect (ISE), i.e. hardness decreases as the indentation depth, under higher loads, increases. This effect has been more pronounced in the glass-ceramic samples. Fracture toughness has been determined by the crack length induced by the Vickers indentations and relating them to the applied loads. Glass materials presented a fracture pattern with characteristics of cleavage, forming cracks of the half-penny shaped type, while the glass-ceramic materials exhibited crack bridging effects and Palmqvist type cracks. (C) 2011 Elsevier B.V. All rights reserved.

Dependence of brittle-to-ductile transition on crystallographic direction in diamond turning of single-crystal silicon

The objective of this paper is to show the dependence relationship between the crystallographic orientations upon brittle-to-ductile transition during diamond turning of monocrystalline silicon. Cutting tests were performed using a -5 degrees rake angle round nose diamond tool at different machining scales. At the micrometre level, the feedrate was kept constant at 2.5 micrometres per revolution (mu m/r), and the depth of cut was varied from 1 to 5 mu m. At the submicrometre level, the depth of cut was kept constant at 500 nm and the feedrate varied from 5 to 10 mu m/r. At the micrometre level, the uncut shoulder generated with an interrupted cutting test procedure provided a quantitative measurement of the ductile-to-brittle transition. Results show that the critical chip thickness in silicon for ductile material removal reaches a maximum of 285 nm in the [100] direction and a minimum of 115 nm in the [110] direction, when the depth of cut was 5 mu m. It was found that when a submicrometre depth of cut was applied, microcracks were revealed in the [110] direction, which is the softer direction in silicon. Micro Raman spectroscopy was used to estimate surface residual stress after machining. Compressive residual stress in the range 142 MPa and smooth damage free surface finish was probed in the [100] direction for a depth of cut of 5 mu m, whereas residual stresses in the range 350 MPa and brittle damage was probed in the [110] direction for a depth of cut of 500 nm.

Preparation and characterisation of refractory whiskers and selected alumina composites

A whisker is a common name of single crystalline inorganic fibre of small dimensions, typically 0.5-1 μm in diameter and 20-50 μm in length. Whiskers are mainly used as reinforcement of ceramics. This work describes the synthesis and characterisation of new whisker types. Ti0.33Ta0.33Nb0.33CxN1-x, TiB2, B4C, and LaxCe1-xB6 have been prepared by carbothermal vapour–liquid–solid (CTR-VLS) growth mechanisms in the temperature range 900-1800°C, in argon or nitrogen. Generally, carbon and different suitable oxides were used as whisker precursors. The oxides reacted via a carbothermal reduction process. A halogenide salt was added to form gaseous metal halogenides or oxohalogenides and small amount of a transition metal was added to catalyse the whisker growth. In this mechanism, the whisker constituents are dissolved into the catalyst, in liquid phase, which becomes supersaturated. Then a whisker could nucleate and grow out under continuous feed of constituents. The syntheses of TiC, TiB2, and B4C were followed at ordinary synthesis conditions by means of mass spectrometry (MS), thermogravimetry (TG), differential thermal analysis (DTA) and quenching. The main reaction starting temperatures and reaction time for the different mixtures was revealed, and it was found that the temperature inside the crucible during the reactions was up to 100°C below the furnace set-point, due to endothermic nature of the reactions. Quench experiments showed that whiskers were formed already when reaching the temperature plateau, but the yield increased fast with the holding time and reached a maximum after about 20-30 minutes. Growth models for whisker formation have been proposed. Alumina based composites reinforced by (2-5 vol.%) TiCnano and TiNnano and 25 vol.% of carbide, and boride phases (whiskers and particulates of TiC, TiN, TaC, NbC, (Ti,Ta)C, (Ti,Ta,Nb)C, SiC, TiB2 and B4C) have been prepared by a developed aqueous colloidal processing route followed by hot pressing for 90 min at 1700°C, 28 MPa or SPS sintering for 5 minutes at 1200-1600°C and 75 MPa. Vickers indentation measurements showed that the lowest possible sintering temperature is to prefer from mechanical properties point of view. In the TiNnano composites the fracture mode was typically intergranular, while it was transgranular in the SiCnano composites. The whisker and particulate composites have been compared in terms of e.g. microstructure and mechanical properties. Generally, additions of whiskers yielded higher fracture toughness compared to particulates. Composites of commercially available SiC whiskers showed best mechanical properties with a low spread but all the other whisker phases, especially TiB2, exhibited a great potential as reinforcement materials.

Some issues concerning the dynamic response and damage of composite laminates subjected to low velocity impact

Abstract. This thesis presents a discussion on a few specific topics regarding the low velocity impact behaviour of laminated composites. These topics were chosen because of their significance as well as the relatively limited attention received so far by the scientific community. The first issue considered is the comparison between the effects induced by a low velocity impact and by a quasi-static indentation experimental test. An analysis of both test conditions is presented, based on the results of experiments carried out on carbon fibre laminates and on numerical computations by a finite element model. It is shown that both quasi-static and dynamic tests led to qualitatively similar failure patterns; three characteristic contact force thresholds, corresponding to the main steps of damage progression, were identified and found to be equal for impact and indentation. On the other hand, an equal energy absorption resulted in a larger delaminated area in quasi-static than in dynamic tests, while the maximum displacement of the impactor (or indentor) was higher in the case of impact, suggesting a probably more severe fibre damage than in indentation. Secondly, the effect of different specimen dimensions and boundary conditions on its impact response was examined. Experimental testing showed that the relationships of delaminated area with two significant impact parameters, the absorbed energy and the maximum contact force, did not depend on the in-plane dimensions and on the support condition of the coupons. The possibility of predicting, by means of a simplified numerical computation, the occurrence of delaminations during a specific impact event is also discussed. A study about the compressive behaviour of impact damaged laminates is also presented. Unlike most of the contributions available about this subject, the results of compression after impact tests on thin laminates are described in which the global specimen buckling was not prevented. Two different quasi-isotropic stacking sequences, as well as two specimen geometries, were considered. It is shown that in the case of rectangular coupons the lay-up can significantly affect the damage induced by impact. Different buckling shapes were observed in laminates with different stacking sequences, in agreement with the results of numerical analysis. In addition, the experiments showed that impact damage can alter the buckling mode of the laminates in certain situations, whereas it did not affect the compressive strength in every case, depending on the buckling shape. Some considerations about the significance of the test method employed are also proposed. Finally, a comprehensive study is presented regarding the influence of pre-existing in-plane loads on the impact response of laminates. Impact events in several conditions, including both tensile and compressive preloads, both uniaxial and biaxial, were analysed by means of numerical finite element simulations; the case of laminates impacted in postbuckling conditions was also considered. The study focused on how the effect of preload varies with the span-to-thickness ratio of the specimen, which was found to be a key parameter. It is shown that a tensile preload has the strongest effect on the peak stresses at low span-to-thickness ratios, leading to a reduction of the minimum impact energy required to initiate damage, whereas this effect tends to disappear as the span-to-thickness ratio increases. On the other hand, a compression preload exhibits the most detrimental effects at medium span-to-thickness ratios, at which the laminate compressive strength and the critical instability load are close to each other, while the influence of preload can be negligible for thin plates or even beneficial for very thick plates. The possibility to obtain a better explanation of the experimental results described in the literature, in view of the present findings, is highlighted. Throughout the thesis the capabilities and limitations of the finite element model, which was implemented in an in-house program, are discussed. The program did not include any damage model of the material. It is shown that, although this kind of analysis can yield accurate results as long as damage has little effect on the overall mechanical properties of a laminate, it can be helpful in explaining some phenomena and also in distinguishing between what can be modelled without taking into account the material degradation and what requires an appropriate simulation of damage. Sommario. Questa tesi presenta una discussione su alcune tematiche specifiche riguardanti il comportamento dei compositi laminati soggetti ad impatto a bassa velocità. Tali tematiche sono state scelte per la loro importanza, oltre che per l’attenzione relativamente limitata ricevuta finora dalla comunità scientifica. La prima delle problematiche considerate è il confronto fra gli effetti prodotti da una prova sperimentale di impatto a bassa velocità e da una prova di indentazione quasi statica. Viene presentata un’analisi di entrambe le condizioni di prova, basata sui risultati di esperimenti condotti su laminati in fibra di carbonio e su calcoli numerici svolti con un modello ad elementi finiti. È mostrato che sia le prove quasi statiche sia quelle dinamiche portano a un danneggiamento con caratteristiche qualitativamente simili; tre valori di soglia caratteristici della forza di contatto, corrispondenti alle fasi principali di progressione del danno, sono stati individuati e stimati uguali per impatto e indentazione. D’altro canto lo stesso assorbimento di energia ha portato ad un’area delaminata maggiore nelle prove statiche rispetto a quelle dinamiche, mentre il massimo spostamento dell’impattatore (o indentatore) è risultato maggiore nel caso dell’impatto, indicando la probabilità di un danneggiamento delle fibre più severo rispetto al caso dell’indentazione. In secondo luogo è stato esaminato l’effetto di diverse dimensioni del provino e diverse condizioni al contorno sulla sua risposta all’impatto. Le prove sperimentali hanno mostrato che le relazioni fra l’area delaminata e due parametri di impatto significativi, l’energia assorbita e la massima forza di contatto, non dipendono dalle dimensioni nel piano dei provini e dalle loro condizioni di supporto. Viene anche discussa la possibilità di prevedere, per mezzo di un calcolo numerico semplificato, il verificarsi di delaminazioni durante un determinato caso di impatto. È presentato anche uno studio sul comportamento a compressione di laminati danneggiati da impatto. Diversamente della maggior parte della letteratura disponibile su questo argomento, vengono qui descritti i risultati di prove di compressione dopo impatto su laminati sottili durante le quali l’instabilità elastica globale dei provini non è stata impedita. Sono state considerate due differenti sequenze di laminazione quasi isotrope, oltre a due geometrie per i provini. Viene mostrato come nel caso di provini rettangolari la sequenza di laminazione possa influenzare sensibilmente il danno prodotto dall’impatto. Due diversi tipi di deformate in condizioni di instabilità sono stati osservati per laminati con diversa laminazione, in accordo con i risultati dell’analisi numerica. Gli esperimenti hanno mostrato inoltre che in certe situazioni il danno da impatto può alterare la deformata che il laminato assume in seguito ad instabilità; d’altra parte tale danno non ha sempre influenzato la resistenza a compressione, a seconda della deformata. Vengono proposte anche alcune considerazioni sulla significatività del metodo di prova utilizzato. Infine viene presentato uno studio esaustivo riguardo all’influenza di carichi membranali preesistenti sulla risposta all’impatto dei laminati. Sono stati analizzati con simulazioni numeriche ad elementi finiti casi di impatto in diverse condizioni di precarico, sia di trazione sia di compressione, sia monoassiali sia biassiali; è stato preso in considerazione anche il caso di laminati impattati in condizioni di postbuckling. Lo studio si è concentrato in particolare sulla dipendenza degli effetti del precarico dal rapporto larghezza-spessore del provino, che si è rivelato un parametro fondamentale. Viene illustrato che un precarico di trazione ha l’effetto più marcato sulle massime tensioni per bassi rapporti larghezza-spessore, portando ad una riduzione della minima energia di impatto necessaria per innescare il danneggiamento, mentre questo effetto tende a scomparire all’aumentare di tale rapporto. Il precarico di compressione evidenzia invece gli effetti più deleteri a rapporti larghezza-spessore intermedi, ai quali la resistenza a compressione del laminato e il suo carico critico di instabilità sono paragonabili, mentre l’influenza del precarico può essere trascurabile per piastre sottili o addirittura benefica per piastre molto spesse. Viene evidenziata la possibilità di trovare una spiegazione più soddisfacente dei risultati sperimentali riportati in letteratura, alla luce del presente contributo. Nel corso della tesi vengono anche discussi le potenzialità ed i limiti del modello ad elementi finiti utilizzato, che è stato implementato in un programma scritto in proprio. Il programma non comprende alcuna modellazione del danneggiamento del materiale. Viene però spiegato come, nonostante questo tipo di analisi possa portare a risultati accurati soltanto finché il danno ha scarsi effetti sulle proprietà meccaniche d’insieme del laminato, esso possa essere utile per spiegare alcuni fenomeni, oltre che per distinguere fra ciò che si può riprodurre senza tenere conto del degrado del materiale e ciò che invece richiede una simulazione adeguata del danneggiamento.

Modellazione multibody di veicoli ferroviari: sviluppo ed implementazione di modelli innovativi per l'analisi del contatto ruota - rotaia

The wheel - rail contact analysis plays a fundamental role in the multibody modeling of railway vehicles. A good contact model must provide an accurate description of the global contact phenomena (contact forces and torques, number and position of the contact points) and of the local contact phenomena (position and shape of the contact patch, stresses and displacements). The model has also to assure high numerical efficiency (in order to be implemented directly online within multibody models) and a good compatibility with commercial multibody software (Simpack Rail, Adams Rail). The wheel - rail contact problem has been discussed by several authors and many models can be found in the literature. The contact models can be subdivided into two different categories: the global models and the local (or differential) models. Currently, as regards the global models, the main approaches to the problem are the so - called rigid contact formulation and the semi – elastic contact description. The rigid approach considers the wheel and the rail as rigid bodies. The contact is imposed by means of constraint equations and the contact points are detected during the dynamic simulation by solving the nonlinear algebraic differential equations associated to the constrained multibody system. Indentation between the bodies is not permitted and the normal contact forces are calculated through the Lagrange multipliers. Finally the Hertz’s and the Kalker’s theories allow to evaluate the shape of the contact patch and the tangential forces respectively. Also the semi - elastic approach considers the wheel and the rail as rigid bodies. However in this case no kinematic constraints are imposed and the indentation between the bodies is permitted. The contact points are detected by means of approximated procedures (based on look - up tables and simplifying hypotheses on the problem geometry). The normal contact forces are calculated as a function of the indentation while, as in the rigid approach, the Hertz’s and the Kalker’s theories allow to evaluate the shape of the contact patch and the tangential forces. Both the described multibody approaches are computationally very efficient but their generality and accuracy turn out to be often insufficient because the physical hypotheses behind these theories are too restrictive and, in many circumstances, unverified. In order to obtain a complete description of the contact phenomena, local (or differential) contact models are needed. In other words wheel and rail have to be considered elastic bodies governed by the Navier’s equations and the contact has to be described by suitable analytical contact conditions. The contact between elastic bodies has been widely studied in literature both in the general case and in the rolling case. Many procedures based on variational inequalities, FEM techniques and convex optimization have been developed. This kind of approach assures high generality and accuracy but still needs very large computational costs and memory consumption. Due to the high computational load and memory consumption, referring to the current state of the art, the integration between multibody and differential modeling is almost absent in literature especially in the railway field. However this integration is very important because only the differential modeling allows an accurate analysis of the contact problem (in terms of contact forces and torques, position and shape of the contact patch, stresses and displacements) while the multibody modeling is the standard in the study of the railway dynamics. In this thesis some innovative wheel – rail contact models developed during the Ph. D. activity will be described. Concerning the global models, two new models belonging to the semi – elastic approach will be presented; the models satisfy the following specifics: 1) the models have to be 3D and to consider all the six relative degrees of freedom between wheel and rail 2) the models have to consider generic railway tracks and generic wheel and rail profiles 3) the models have to assure a general and accurate handling of the multiple contact without simplifying hypotheses on the problem geometry; in particular the models have to evaluate the number and the position of the contact points and, for each point, the contact forces and torques 4) the models have to be implementable directly online within the multibody models without look - up tables 5) the models have to assure computation times comparable with those of commercial multibody software (Simpack Rail, Adams Rail) and compatible with RT and HIL applications 6) the models have to be compatible with commercial multibody software (Simpack Rail, Adams Rail). The most innovative aspect of the new global contact models regards the detection of the contact points. In particular both the models aim to reduce the algebraic problem dimension by means of suitable analytical techniques. This kind of reduction allows to obtain an high numerical efficiency that makes possible the online implementation of the new procedure and the achievement of performance comparable with those of commercial multibody software. At the same time the analytical approach assures high accuracy and generality. Concerning the local (or differential) contact models, one new model satisfying the following specifics will be presented: 1) the model has to be 3D and to consider all the six relative degrees of freedom between wheel and rail 2) the model has to consider generic railway tracks and generic wheel and rail profiles 3) the model has to assure a general and accurate handling of the multiple contact without simplifying hypotheses on the problem geometry; in particular the model has to able to calculate both the global contact variables (contact forces and torques) and the local contact variables (position and shape of the contact patch, stresses and displacements) 4) the model has to be implementable directly online within the multibody models 5) the model has to assure high numerical efficiency and a reduced memory consumption in order to achieve a good integration between multibody and differential modeling (the base for the local contact models) 6) the model has to be compatible with commercial multibody software (Simpack Rail, Adams Rail). In this case the most innovative aspects of the new local contact model regard the contact modeling (by means of suitable analytical conditions) and the implementation of the numerical algorithms needed to solve the discrete problem arising from the discretization of the original continuum problem. Moreover, during the development of the local model, the achievement of a good compromise between accuracy and efficiency turned out to be very important to obtain a good integration between multibody and differential modeling. At this point the contact models has been inserted within a 3D multibody model of a railway vehicle to obtain a complete model of the wagon. The railway vehicle chosen as benchmark is the Manchester Wagon the physical and geometrical characteristics of which are easily available in the literature. The model of the whole railway vehicle (multibody model and contact model) has been implemented in the Matlab/Simulink environment. The multibody model has been implemented in SimMechanics, a Matlab toolbox specifically designed for multibody dynamics, while, as regards the contact models, the CS – functions have been used; this particular Matlab architecture allows to efficiently connect the Matlab/Simulink and the C/C++ environment. The 3D multibody model of the same vehicle (this time equipped with a standard contact model based on the semi - elastic approach) has been then implemented also in Simpack Rail, a commercial multibody software for railway vehicles widely tested and validated. Finally numerical simulations of the vehicle dynamics have been carried out on many different railway tracks with the aim of evaluating the performances of the whole model. The comparison between the results obtained by the Matlab/ Simulink model and those obtained by the Simpack Rail model has allowed an accurate and reliable validation of the new contact models. In conclusion to this brief introduction to my Ph. D. thesis, we would like to thank Trenitalia and the Regione Toscana for the support provided during all the Ph. D. activity. Moreover we would also like to thank the INTEC GmbH, the society the develops the software Simpack Rail, with which we are currently working together to develop innovative toolboxes specifically designed for the wheel rail contact analysis.