A Practical Micro Mechanical Model Based Local Approach Methodology for the Analysis of Ductile Fracture of Welded T-Joints

This thesis concentrates on developing a practical local approach methodology based on micro mechanical models for the analysis of ductile fracture of welded joints. Two major problems involved in the local approach, namely the dilational constitutive relation reflecting the softening behaviour of material, and the failure criterion associated with the constitutive equation, have been studied in detail. Firstly, considerable efforts were made on the numerical integration and computer implementation for the non trivial dilational Gurson Tvergaard model. Considering the weaknesses of the widely used Euler forward integration algorithms, a family of generalized mid point algorithms is proposed for the Gurson Tvergaard model. Correspondingly, based on the decomposition of stresses into hydrostatic and deviatoric parts, an explicit seven parameter expression for the consistent tangent moduli of the algorithms is presented. This explicit formula avoids any matrix inversion during numerical iteration and thus greatly facilitates the computer implementation of the algorithms and increase the efficiency of the code. The accuracy of the proposed algorithms and other conventional algorithms has been assessed in a systematic manner in order to highlight the best algorithm for this study. The accurate and efficient performance of present finite element implementation of the proposed algorithms has been demonstrated by various numerical examples. It has been found that the true mid point algorithm (a = 0.5) is the most accurate one when the deviatoric strain increment is radial to the yield surface and it is very important to use the consistent tangent moduli in the Newton iteration procedure. Secondly, an assessment of the consistency of current local failure criteria for ductile fracture, the critical void growth criterion, the constant critical void volume fraction criterion and Thomason's plastic limit load failure criterion, has been made. Significant differences in the predictions of ductility by the three criteria were found. By assuming the void grows spherically and using the void volume fraction from the Gurson Tvergaard model to calculate the current void matrix geometry, Thomason's failure criterion has been modified and a new failure criterion for the Gurson Tvergaard model is presented. Comparison with Koplik and Needleman's finite element results shows that the new failure criterion is fairly accurate indeed. A novel feature of the new failure criterion is that a mechanism for void coalescence is incorporated into the constitutive model. Hence the material failure is a natural result of the development of macroscopic plastic flow and the microscopic internal necking mechanism. By the new failure criterion, the critical void volume fraction is not a material constant and the initial void volume fraction and/or void nucleation parameters essentially control the material failure. This feature is very desirable and makes the numerical calibration of void nucleation parameters(s) possible and physically sound. Thirdly, a local approach methodology based on the above two major contributions has been built up in ABAQUS via the user material subroutine UMAT and applied to welded T joints. By using the void nucleation parameters calibrated from simple smooth and notched specimens, it was found that the fracture behaviour of the welded T joints can be well predicted using present methodology. This application has shown how the damage parameters of both base material and heat affected zone (HAZ) material can be obtained in a step by step manner and how useful and capable the local approach methodology is in the analysis of fracture behaviour and crack development as well as structural integrity assessment of practical problems where non homogeneous materials are involved. Finally, a procedure for the possible engineering application of the present methodology is suggested and discussed.

Surgically-treated Hip Fracture in older People. With Special Emphasis on Mortality Analysis

Hip fractures are associated with significant morbidity and mortality. Cervical and trochanteric fractures have a different morphometry, surgical treatment, and outcome. Polypharmacy, common in older people, is associated with increased mortality. The risk factors for mortality can be identified based on cause-of-death analysis. In this population-based study, 461 older, surgically in 1999-2000 treated hip fracture patients were enrolled. Incidence, morphometry, medication, mortality, and cause-of-death were analysed. Hip fractures were most commonly sustained by women, occurred mostly indoors, and often in institutions. One in four patients had sustained a previous fracture. Routine clinical radiographs revealed no differences in the hip geometry between hip fracture types. Age-adjusted mortality was higher in men than in women during the follow-up. Chronic lung disease and male sex were predictors of mortality after cervical fracture. In men, potent anticholinergics were associated with excess age-adjusted mortality. Men were more likely to die from circulatory disease and dementia after hip fracture than women. Mortality after hip fracture was 3-fold higher than that of the general population, including every cause-of-death class. Fracture prevention in institutions and homes, indoor safety measures, and treatment of chronic lung diseases should be encouraged. Hip morphometry analyses require more accurate measures than that provided by routine radiographs. Careful use of potent anticholinergics may reduce mortality. Compared to the general population, excess mortality after hip fracture was evident up to 9 years after hip fracture. Cause-of-death analysis indicates that all major comorbidities require optimal treatment after hip fracture surgery.

Fatigue Analysis of Welded Cruciform Joints with Cold Laps

Työssä on tutkittu vetojännityskuormituksen alaisena olevien hitsattujen kuormaa kantamattomien X-liitosten hitsin paikallisen geometrian variaation vaikutusta väsymislujuuteen. Muuttujina olivat reunan pyöristyssäde, kylmäjuoksun suuruus ja kylkikulma. Geometristen muuttujien parametrinen riippuvuussuhde on analysoitu usealla elementtimallilla. Väsymistarkastelu on suoritettu käyttämällä lineaaris-elastista murtumismekaniikkaa (LEFM) tasovenymätilassa ja materiaalina terästä. Särönkasvun suunnan ennustamisessaon käytetty maksimipääjännityskriteeriä sekä jännitysintensiteettikertoimet on määritetty J-integraalilla. Särön ydintymisvaihetta ei ole otettu huomioon. Rakenteen on oletettu olevan hitsatussa tilassa ja jännitysheilahdus on kokonaan tehollinen. Särön kasvunopeuden ennustamiseen on käytetty Paris'n lakia. Väsymislujuustulokset on esitetty karakteristisina väsymisluokkina (FAT) ja sovitettu parametriseksi yhtälöksi. Lopuksi väsymisanalyysin ennustamia tuloksia on verrattu saatavilla oleviin väsytystestituloksiin.

Multivariable regression analysis and hot-spot stress approach into fatigue life estimations of welded joints

Kolmen eri hitsausliitoksen väsymisikä arvio on analysoitu monimuuttuja regressio analyysin avulla. Regression perustana on laaja S-N tietokanta joka on kerätty kirjallisuudesta. Tarkastellut liitokset ovat tasalevy liitos, krusiformi liitos ja pitkittäisripa levyssä. Muuttujina ovat jännitysvaihtelu, kuormitetun levyn paksuus ja kuormitus tapa. Paksuus effekti on käsitelty uudelleen kaikkia kolmea liitosta ajatellen. Uudelleen käsittelyn avulla on varmistettu paksuus effektin olemassa olo ennen monimuuttuja regressioon siirtymistä. Lineaariset väsymisikä yhtalöt on ajettu kolmelle hitsausliitokselle ottaen huomioon kuormitetun levyn paksuus sekä kuormitus tapa. Väsymisikä yhtalöitä on verrattu ja keskusteltu testitulosten valossa, jotka on kerätty kirjallisuudesta. Neljä tutkimustaon tehty kerättyjen väsymistestien joukosta ja erilaisia väsymisikä arvio metodeja on käytetty väsymisiän arviointiin. Tuloksia on tarkasteltu ja niistä keskusteltu oikeiden testien valossa. Tutkimuksissa on katsottu 2mm ja 6mm symmetristäpitkittäisripaa levyssä, 12.7mm epäsymmetristä pitkittäisripaa, 38mm symmetristä pitkittäisripaa vääntökuormituksessa ja 25mm/38mm kuorman kantavaa krusiformi liitosta vääntökuormituksessa. Mallinnus on tehty niin lähelle testi liitosta kuin mahdollista. Väsymisikä arviointi metodit sisältävät hot-spot metodin jossa hot-spot jännitys on laskettu kahta lineaarista ja epälineaarista ekstrapolointiakäyttäen sekä paksuuden läpi integrointia käyttäen. Lovijännitys ja murtumismekaniikka metodeja on käytetty krusiformi liitosta laskiessa.

Use of field measurements in fatigue analysis of forestry cranes

Tässä työssä on esitetty väsyttävän kuormituksen mittaamiseen ja mittausdatan jälkikäsittelyyn sekä väsymismitoitukseen liittyviä menetelmiä. Menetelmien sovelluskohteena oli metsäkoneen kuormain, joka on väsyttävästi kuormitettu hitsattu rakenne. Teoriaosassa on kuvattu väsymisilmiötä ja väsymismitoitusmenetelmiä sekä kuormitusten tunnistamiseen ja mittausten jälkikäsittelyyn liittyviä menetelmiä. Yleisimmin käytettyjen väsymismitoitusmenetelmien rinnalle on esitetty luotettavuuteen perustuvaa väsymismitoitusmenetelmää. Kuormainten suunnittelussa on keveys- j a kestoikävaatimusten takia erityisen suuri merkitys väsymisen huomioimisella. Rakenteille on ominaista tietyt toiminnan kannalta välttämättömät hitsatut yksityiskohdat, jotka usein määräävät koko rakenteen kestoiän. Koska nämä ongelmakohdat pystytään useimmiten tunnistamaan jo suunnitteluvaiheessa, voidaan yksityiskohtien muotoilulla usein parantaa huomattavasti koko rakenteen kestoikää. Näiden yksityiskohtien optimointi on osittain mahdollista toteuttaa ilman kuormituskertymätietoa, mutta useimmiten kuormitusten tunnistaminen on edellytys parhaan ratkaisun löytymiselle. Tällöin toistaiseksi paras keino todellisen väsyttävän kuormituksen tunnistamiseksi on pitkäaikaiset kenttämittaukset. Kenttämittauksilla selvitetään rakenteeseen kohdistuvat kuormitukset venymäliuskojen avulla. Kuormitusten tunnistamisella on erityisen suuri merkitys kun halutaan määrittää rakenteen kestoikä. Väsyminen ja väsyttävä kuormitus ovat kuitenkin tilastollisia muuttujia j a yksittäiselle rakenteelle ei ole mahdollista määrittää tarkkaa k estoikää. Tilastollisia menetelmiä käyttäen on kuitenkin mahdollista määrittää rakenteen vaurioitumisriski. Laskettaessa vaurioitumisriskiä suurelle määrälle yksittäisiä rakenteita voidaan muodostaa tarkkojakin ennusteita mahdollisten vaurioiden lukumäärästä. Tällöin kuormituskertymätiedosta voi olla tavanomaisen suunnittelun lisäksi laajempaa hyötyä esimerkiksi takuukäsittelyssä. Tässä työssä on sovellettu esitettyjä teorioita käytännössä metsäkoneen harvesterin puomiston väsymistarkasteluun. Kyseisen rakenteen kuormituksia mitattiin kahden viikon aikana yhteensä 35 tuntia, jonka perusteella laskettiin väsyttävän kuormituksen tilastollinen jakauma esimerkkitapaukselle. Mittauksen perusteella ei voitu tehdä kuitenkaan johtopäätöksiä tuotteen koko elinkaaren kuormituksista eikä muiden samanlaisten tuotteiden kuormituksista, koska mitattu otos oli suhteellisen lyhyt ja rajoittui vain yhteen käyttäjään ja muutamaan käyttökohteeseen. Menetelmien testaamiseksi kyseinen otos oli kuitenkin riittävä. Kuormituskertymätietoa käytettiin hyväksi myös laatumääritysten muodostamisessaesimerkkitapaukselle. Murtumismekaniikkaan perustuvalla menetelmällä arvioitiinharvesteripilarin valun mahdollisten valuvirheiden suurin sallittu koko. Luotettavuuteen pohjautuvan mitoitusmenettelyn tarve näyttää olevanlisääntymässä, joten pitkäaikaisten kenttämittausten tehokas hyödyntäminen tulee olemaan keskeinen osa väsymismitoitusta lähitulevaisuudessa. Menetelmiä olisi mahdollista tehostaa yhdistämällä kuormituskertymään erilaisia kuormitusten suhteen riippuvia tunnettuja suureita kuten käsiteltävän puun halkaisija. Todellisettuotekohtaiset tilastolliset jakaumat kuormituksista voitaisiin muodostaa mahdollisesti tehokkaammin, jos esimerkiksi kuormitusten riippuvuus metsätyypistä pystyttäisiin ensin määrittämään.

Flexible Multibody Simulation Approach in the Dynamic Analysis of Bone Strains Activity

The objective of this study is to show that bone strains due to dynamic mechanical loading during physical activity can be analysed using the flexible multibody simulation approach. Strains within the bone tissue play a major role in bone (re)modeling. Based on previous studies, it has been shown that dynamic loading seems to be more important for bone (re)modeling than static loading. The finite element method has been used previously to assess bone strains. However, the finite element method may be limited to static analysis of bone strains due to the expensive computation required for dynamic analysis, especially for a biomechanical system consisting of several bodies. Further, in vivo implementation of strain gauges on the surfaces of bone has been used previously in order to quantify the mechanical loading environment of the skeleton. However, in vivo strain measurement requires invasive methodology, which is challenging and limited to certain regions of superficial bones only, such as the anterior surface of the tibia. In this study, an alternative numerical approach to analyzing in vivo strains, based on the flexible multibody simulation approach, is proposed. In order to investigate the reliability of the proposed approach, three 3-dimensional musculoskeletal models where the right tibia is assumed to be flexible, are used as demonstration examples. The models are employed in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model is developed using the actual geometry of the subject’s tibia, which is obtained from 3 dimensional reconstruction of Magnetic Resonance Images. Inverse dynamics simulation based on motion capture data obtained from walking at a constant velocity is used to calculate the desired contraction trajectory for each muscle. In the forward dynamics simulation, a proportional derivative servo controller is used to calculate each muscle force required to reproduce the motion, based on the desired muscle contraction trajectory obtained from the inverse dynamics simulation. Experimental measurements are used to verify the models and check the accuracy of the models in replicating the realistic mechanical loading environment measured from the walking test. The predicted strain results by the models show consistency with literature-based in vivo strain measurements. In conclusion, the non-invasive flexible multibody simulation approach may be used as a surrogate for experimental bone strain measurement, and thus be of use in detailed strain estimation of bones in different applications. Consequently, the information obtained from the present approach might be useful in clinical applications, including optimizing implant design and devising exercises to prevent bone fragility, accelerate fracture healing and reduce osteoporotic bone loss.

Fatigue analysis of brackets with and without scallop for bridge and other deck structure using effective notch stress.

Due to functional requirement of a structural detail brackets with and without scallop are frequently used in bridges, decks, ships and offshore structure. Scallops are designed to serve as passage way for fluids, to reduce weld length and plate distortions. Moreover, scallops are used to avoid intersection of two or more welds for the fact that there is the presence of inventible inherent initial crack except for full penetrated weld and the formation of multi-axial stress state at the weld intersection. Welding all around the scallop corner increase the possibility of brittle fracture even for the case the bracket is not loaded by primary load. Avoiding of scallop will establish an initial crack in the corner if bracket is welded by fillet welds. If the two weld run pass had crossed, this would have given a 3D residual stress situation. Therefore the presences and absence of scallop necessitates the 3D FEA fatigue resistance of both types of brackets using effective notch stress approach ( ). FEMAP 10.1 with NX NASTRAN was used for the 3D FEA. The first and main objective of this research was to investigate and compare the fatigue resistance of brackets with and without scallop. The secondary goal was the fatigue design of scallops in case they cannot be avoided for some reason. The fatigue resistance for both types of brackets was determined based on approach using 1 mm fictitiously rounded radius based on IIW recommendation. Identical geometrical, boundary and loading conditions were used for the determination and comparison of fatigue resistance of both types of brackets using linear 3D FEA. Moreover the size effect of bracket length was also studied using 2D SHELL element FEA. In the case of brackets with scallop the flange plate weld toe at the corner of the scallop was found to exhibit the highest and made the flange plate weld toe critical for fatigue failure. Whereas weld root and weld toe at the weld intersections were the highly stressed location for brackets without scallop. Thus weld toe for brackets with scallop, and weld root and weld toe for brackets without scallop were found to be the critical area for fatigue failure. Employing identical parameters on both types of brackets, brackets without scallop had the highest except for full penetrated weld. Furthermore the fatigue resistance of brackets without scallop was highly affected by the lack of weld penetration length and it was found out that decreased as the weld penetration was increased. Despite the fact that the very presence of scallop reduces the stiffness and also same time induce stress concentration, based on the 3D FEA it is worth concluding that using scallop provided better fatigue resistance when both types of brackets were fillet welded. However brackets without scallop had the highest fatigue resistance when full penetration weld was used. This thesis also showed that weld toe for brackets with scallop was the only highly stressed area unlike brackets without scallop in which both weld toe and weld root were the critical locations for fatigue failure when different types of boundary conditions were used. Weld throat thickness, plate thickness, scallop radius, lack of weld penetration length, boundary condition and weld quality affected the fatigue resistance of both types of brackets. And as a result, bracket design procedure, especially welding quality and post weld treatment techniques significantly affect the fatigue resistance of both type of brackets.

The Use of Radiostereometric Analysis in Fractures of the Distal Radius: From Phantom Models to Clinical Application

Radiostereometric analysis (RSA) is a highly accurate method for the measurement of in vivo micromotion of orthopaedic implants. Validation of the RSA method is a prerequisite for performing clinical RSA studies. Only a limited number of studies have utilised the RSA method in the evaluation of migration and inducible micromotion during fracture healing. Volar plate fixation of distal radial fractures has increased in popularity. There is still very little prospective randomised evidence supporting the use of these implants over other treatments. The aim of this study was to investigate the precision, accuracy, and feasibility of using RSA in the evaluation of healing in distal radius fractures treated with a volar fixed-angle plate. A physical phantom model was used to validate the RSA method for simple distal radius fractures. A computer simulation model was then used to validate the RSA method for more complex interfragmentary motion in intra-articular fractures. A separate pre-clinical investigation was performed in order to evaluate the possibility of using novel resorbable markers for RSA. Based on the validation studies, a prospective RSA cohort study of fifteen patients with plated AO type-C distal radius fractures with a 1-year follow-up was performed. RSA was shown to be highly accurate and precise in the measurement of fracture micromotion using both physical and computer simulated models of distal radius fractures. Resorbable RSA markers demonstrated potential for use in RSA. The RSA method was found to have a high clinical precision. The fractures underwent significant translational and rotational migration during the first two weeks after surgery, but not thereafter. Maximal grip caused significant translational and rotational interfragmentary micromotion. This inducible micromotion was detectable up to eighteen weeks, even after the achievement of radiographic union. The application of RSA in the measurement of fracture fragment migration and inducible interfragmentary micromotion in AO type-C distal radius fractures is feasible but technically demanding. RSA may be a unique tool in defining the progress of fracture union.