Dynamic interface model for masonry walls subjected to high strain rate out-of-plane loads

The present study proposes a dynamic constitutive material interface model that includes non-associated flow rule and high strain rate effects, implemented in the finite element code ABAQUS as a user subroutine. First, the model capability is validated with numerical simulations of unreinforced block work masonry walls subjected to low velocity impact. The results obtained are compared with field test data and good agreement is found. Subsequently, a comprehensive parametric analysis is accomplished with different joint tensile strengths and cohesion, and wall thickness to evaluate the effect of the parameter variations on the impact response of masonry walls.

Evaluation of temperature influence on embedded kinematic system

Dissertação de mestrado integrado em Engenharia Mecânica

A multi–physics approach applied to masonry structures with non–hydraulic lime mortars

Tese de Doutoramento em Engenharia Civil (área de especialização em Engenharia de Estruturas).

Una arqueometria del canvi tecnològic: producció i consum d'àmfores durant el canvi d'era en la zona nord de la costa catalana

L’objectiu de la recerca és definir un marc teòric i metodològic per a l’estudi del canvi tecnològic en Arqueologia. Aquest model posa èmfasi en caracteritzar els compromisos que configuren una tecnologia i avaluar-los en funció dels factors de situació —tècnics, econòmics, polítics, socials i ideològics. S’ha aplicat aquest model a un cas d’estudi concret: la producció d’àmfores romanes durant el canvi d’Era en la província Tarraconensis. L’estudi tecnològic dels envasos s’ha realitzat mitjançant diverses tècniques analítiques: Fluorescència de raigs X (FRX), Difracció de raigs X (DRX), Microscòpia òptica (MO) i Microscòpia electrònica de rastreig (MER). Les dades obtingudes permeten, a més, establir els grups de referència per a cada centre productor d’àmfores i, així, identificar la provinença dels individus recuperats en els centres consumidors. Donat que les àmfores en estudi són artefactes dissenyats específicament per a ser estibats en una nau i servir com a envàs de transport, l’estudi inclou la caracterització de les propietats mecàniques de resistència a la fractura i de tenacitat. En aquest sentit, i per primera vegada, s’ha aplicat l’Anàlisi d’Elements Finits (AEF) per a conèixer el comportament dels diferents dissenys d’àmfora en ésser sotmesos a diverses forces d’ús. L’AEF permet simular per ordinador les activitats en què les àmfores haurien participat durant el seu ús i avaluar-ne el seu comportament tècnic. Els resultats mostren una gran adequació entre les formulacions teòriques i el programa analític implementat per a aquest estudi. Respecte el cas d’estudi, els resultats mostren una gran variabilitat en les eleccions tecnològiques preses pels ceramistes de diferents tallers, però també al llarg del període de funcionament d’un mateix taller. L’aplicació del model ha permès proposar una explicació al canvi de disseny de les àmfores romanes.

Translation of biomechanical concepts in bone tissue engineering: from animal study to revision knee arthroplasty.

Bone defects in revision knee arthroplasty are often located in load-bearing regions. The goal of this study was to determine whether a physiologic load could be used as an in situ osteogenic signal to the scaffolds filling the bone defects. In order to answer this question, we proposed a novel translation procedure having four steps: (1) determining the mechanical stimulus using finite element method, (2) designing an animal study to measure bone formation spatially and temporally using micro-CT imaging in the scaffold subjected to the estimated mechanical stimulus, (3) identifying bone formation parameters for the loaded and non-loaded cases appearing in a recently developed mathematical model for bone formation in the scaffold and (4) estimating the stiffness and the bone formation in the bone-scaffold construct. With this procedure, we estimated that after 3 years mechanical stimulation increases the bone volume fraction and the stiffness of scaffold by 1.5- and 2.7-fold, respectively, compared to a non-loaded situation.

Analysis of a finite volume element method for the Stokes problem

In this paper we propose a stabilized conforming finite volume element method for the Stokes equations. On stating the convergence of the method, optimal a priori error estimates in different norms are obtained by establishing the adequate connection between the finite volume and stabilized finite element formulations. A superconvergence result is also derived by using a postprocessing projection method. In particular, the stabilization of the continuous lowest equal order pair finite volume element discretization is achieved by enriching the velocity space with local functions that do not necessarily vanish on the element boundaries. Finally, some numerical experiments that confirm the predicted behavior of the method are provided.

In vivo loading increases mechanical properties of scaffold by affecting bone formation and bone resorption rates.

A successful bone tissue engineering strategy entails producing bone-scaffold constructs with adequate mechanical properties. Apart from the mechanical properties of the scaffold itself, the forming bone inside the scaffold also adds to the strength of the construct. In this study, we investigated the role of in vivo cyclic loading on mechanical properties of a bone scaffold. We implanted PLA/β-TCP scaffolds in the distal femur of six rats, applied external cyclic loading on the right leg, and kept the left leg as a control. We monitored bone formation at 7 time points over 35 weeks using time-lapsed micro-computed tomography (CT) imaging. The images were then used to construct micro-finite element models of bone-scaffold constructs, with which we estimated the stiffness for each sample at all time points. We found that loading increased the stiffness by 60% at 35 weeks. The increase of stiffness was correlated to an increase in bone volume fraction of 18% in the loaded scaffold compared to control scaffold. These changes in volume fraction and related stiffness in the bone scaffold are regulated by two independent processes, bone formation and bone resorption. Using time-lapsed micro-CT imaging and a newly-developed longitudinal image registration technique, we observed that mechanical stimulation increases the bone formation rate during 4-10 weeks, and decreases the bone resorption rate during 9-18 weeks post-operatively. For the first time, we report that in vivo cyclic loading increases mechanical properties of the scaffold by increasing the bone formation rate and decreasing the bone resorption rate.

Total shoulder arthroplasty: downward inclination of the glenoid component to balance supraspinatus deficiency.

HYPOTHESIS: Supraspinatus deficiency associated with total shoulder arthroplasty (TSA) provokes eccentric loading and may induce loosening of the glenoid component. A downward inclination of the glenoid component has been proposed to balance supraspinatus deficiency. METHODS: This hypothesis was assessed by a numeric musculoskeletal model of the glenohumeral joint during active abduction. Three cases were compared: TSA with normal muscular function, TSA with supraspinatus deficiency, and TSA with supraspinatus deficiency and downward inclination of the glenoid. RESULTS: Supraspinatus deficiency increased humeral migration and eccentric loading. A downward inclination of the glenoid partly balanced the loss of stability, but this potential advantage was counterbalanced by an important stress increase within the glenoid cement. The additional subchondral bone reaming required to incline the glenoid component indeed reduced the bone support, increasing cement deformation and stress. CONCLUSION: Glenoid inclination should not be obtained at the expense of subchondral bone support.

A finite-population revenue management model and a risk-ratio procedure for the joint estimation of population size and parameters

Many dynamic revenue management models divide the sale period into a finite number of periods T and assume, invoking a fine-enough grid of time, that each period sees at most one booking request. These Poisson-type assumptions restrict the variability of the demand in the model, but researchers and practitioners were willing to overlook this for the benefit of tractability of the models. In this paper, we criticize this model from another angle. Estimating the discrete finite-period model poses problems of indeterminacy and non-robustness: Arbitrarily fixing T leads to arbitrary control values and on the other hand estimating T from data adds an additional layer of indeterminacy. To counter this, we first propose an alternate finite-population model that avoids this problem of fixing T and allows a wider range of demand distributions, while retaining the useful marginal-value properties of the finite-period model. The finite-population model still requires jointly estimating market size and the parameters of the customer purchase model without observing no-purchases. Estimation of market-size when no-purchases are unobservable has rarely been attempted in the marketing or revenue management literature. Indeed, we point out that it is akin to the classical statistical problem of estimating the parameters of a binomial distribution with unknown population size and success probability, and hence likely to be challenging. However, when the purchase probabilities are given by a functional form such as a multinomial-logit model, we propose an estimation heuristic that exploits the specification of the functional form, the variety of the offer sets in a typical RM setting, and qualitative knowledge of arrival rates. Finally we perform simulations to show that the estimator is very promising in obtaining unbiased estimates of population size and the model parameters.

Electromagnetic Design of a Solid Steel Rotor Motor for Demanding Operation Environments

The solid-rotor induction motor provides a mechanically and thermally reliable solution for demanding environments where other rotor solutions are prohibited or questionable. Solid rotors, which are manufactured of single pieces of ferromagnetic material, are commonly used in motors in which the rotationspeeds exceed substantially the conventional speeds of laminated rotors with squirrel-cage. During the operation of a solid-rotor electrical machine, the rotor core forms a conductor for both the magnetic flux and the electrical current. This causes an increase in the rotor resistance and rotor leakage inductance, which essentially decreases the power factor and the efficiency of the machine. The electromagnetic problems related to the solid-rotor induction motor are mostly associated with the low performance of the rotor. Therefore, the main emphasis in this thesis is put on the solid steel rotor designs. The rotor designs studied in thisthesis are based on the fact that the rotor construction should be extremely robust and reliable to withstand the high mechanical stresses caused by the rotational velocity of the rotor. In addition, the demanding operation environment sets requirements for the applied materials because of the high temperatures and oxidizing acids, which may be present in the cooling fluid. Therefore, the solid rotors analyzed in this thesis are made of a single piece of ferromagnetic material without any additional parts, such as copper end-rings or a squirrel-cage. A pure solid rotor construction is rigid and able to keep its balance over a large speed range. It also may tolerate other environmental stresses such as corroding substances or abrasive particles. In this thesis, the main target is to improve the performance of an induction motor equipped with a solid steel rotor by traditional methods: by axial slitting of the rotor, by selecting a proper rotor core material and by coating the rotor with a high-resistive stainless ferromagnetic material. In the solid steel rotor calculation, the rotor end-effects have a significant effect on the rotor characteristics. Thus, the emphasis is also put on the comparison of different rotor endfactors. In addition, a corrective slip-dependent end-factor is proposed. The rotor designs covered in this thesis are the smooth solid rotor, the axially slitted solid rotor and the slitted rotor having a uniform ferromagnetic coating cylinder. The thesis aims at design rules for multi-megawatt machines. Typically, mega-watt-size solidrotor machines find their applications mainly in the field of electric-motor-gas-compression systems, in steam-turbine applications, and in various types of largepower pump applications, where high operational speeds are required. In this thesis, a 120 kW, 10 000 rpm solid-rotor induction motor is usedas a small-scale model for such megawatt-range solid-rotor machines. The performance of the 120 kW solid-rotor induction motors is determined by experimental measurements and finite element calculations.

Effect of Mechanical and Geometric Mismatching on Fatigue and Damage of Welded Joints

The future of high technology welded constructions will be characterised by higher strength materials and improved weld quality with respect to fatigue resistance. The expected implementation of high quality high strength steel welds will require that more attention be given to the issues of crack initiation and mechanical mismatching. Experiments and finite element analyses were performed within the framework of continuum damage mechanics to investigate the effect of mismatching of welded joints on void nucleation and coalescence during monotonic loading. It was found that the damage of undermatched joints mainly occurred in the sandwich layer and the damageresistance of the joints decreases with the decrease of the sandwich layer width. The damage of over-matched joints mainly occurred in the base metal adjacent to the sandwich layer and the damage resistance of the joints increases with thedecrease of the sandwich layer width. The mechanisms of the initiation of the micro voids/cracks were found to be cracking of the inclusions or the embrittled second phase, and the debonding of the inclusions from the matrix. Experimental fatigue crack growth rate testing showed that the fatigue life of under-matched central crack panel specimens is longer than that of over-matched and even-matched specimens. Further investigation by the elastic-plastic finite element analysis indicated that fatigue crack closure, which originated from the inhomogeneousyielding adjacent to the crack tip, played an important role in the fatigue crack propagation. The applicability of the J integral concept to the mismatched specimens with crack extension under cyclic loading was assessed. The concept of fatigue class used by the International Institute of Welding was introduced in the parametric numerical analysis of several welded joints. The effect of weld geometry and load condition on fatigue strength of ferrite-pearlite steel joints was systematically evaluated based on linear elastic fracture mechanics. Joint types included lap joints, angle joints and butt joints. Various combinations of the tensile and bending loads were considered during the evaluation with the emphasis focused on the existence of both root and toe cracks. For a lap joint with asmall lack-of-penetration, a reasonably large weld leg and smaller flank angle were recommended for engineering practice in order to achieve higher fatigue strength. It was found that the fatigue strength of the angle joint depended strongly on the location and orientation of the preexisting crack-like welding defects, even if the joint was welded with full penetration. It is commonly believed that the double sided butt welds can have significantly higher fatigue strength than that of a single sided welds, but fatigue crack initiation and propagation can originate from the weld root if the welding procedure results in a partial penetration. It is clearly shown that the fatigue strength of the butt joint could be improved remarkably by ensuring full penetration. Nevertheless, increasing the fatigue strength of a butt joint by increasing the size of the weld is an uneconomical alternative.

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.

Modeling the transport phenomena within arterial wall: porous media approach

The transport of macromolecules, such as low-density lipoprotein (LDL), and their accumulation in the layers of the arterial wall play a critical role in the creation and development of atherosclerosis. Atherosclerosis is a disease of large arteries e.g., the aorta, coronary, carotid, and other proximal arteries that involves a distinctive accumulation of LDL and other lipid-bearing materials in the arterial wall. Over time, plaque hardens and narrows the arteries. The flow of oxygen-rich blood to organs and other parts of the body is reduced. This can lead to serious problems, including heart attack, stroke, or even death. It has been proven that the accumulation of macromolecules in the arterial wall depends not only on the ease with which materials enter the wall, but also on the hindrance to the passage of materials out of the wall posed by underlying layers. Therefore, attention was drawn to the fact that the wall structure of large arteries is different than other vessels which are disease-resistant. Atherosclerosis tends to be localized in regions of curvature and branching in arteries where fluid shear stress (shear rate) and other fluid mechanical characteristics deviate from their normal spatial and temporal distribution patterns in straight vessels. On the other hand, the smooth muscle cells (SMCs) residing in the media layer of the arterial wall respond to mechanical stimuli, such as shear stress. Shear stress may affect SMC proliferation and migration from the media layer to intima. This occurs in atherosclerosis and intimal hyperplasia. The study of blood flow and other body fluids and of heat transport through the arterial wall is one of the advanced applications of porous media in recent years. The arterial wall may be modeled in both macroscopic (as a continuous porous medium) and microscopic scales (as a heterogeneous porous medium). In the present study, the governing equations of mass, heat and momentum transport have been solved for different species and interstitial fluid within the arterial wall by means of computational fluid dynamics (CFD). Simulation models are based on the finite element (FE) and finite volume (FV) methods. The wall structure has been modeled by assuming the wall layers as porous media with different properties. In order to study the heat transport through human tissues, the simulations have been carried out for a non-homogeneous model of porous media. The tissue is composed of blood vessels, cells, and an interstitium. The interstitium consists of interstitial fluid and extracellular fibers. Numerical simulations are performed in a two-dimensional (2D) model to realize the effect of the shape and configuration of the discrete phase on the convective and conductive features of heat transfer, e.g. the interstitium of biological tissues. On the other hand, the governing equations of momentum and mass transport have been solved in the heterogeneous porous media model of the media layer, which has a major role in the transport and accumulation of solutes across the arterial wall. The transport of Adenosine 5´-triphosphate (ATP) is simulated across the media layer as a benchmark to observe how SMCs affect on the species mass transport. In addition, the transport of interstitial fluid has been simulated while the deformation of the media layer (due to high blood pressure) and its constituents such as SMCs are also involved in the model. In this context, the effect of pressure variation on shear stress is investigated over SMCs induced by the interstitial flow both in 2D and three-dimensional (3D) geometries for the media layer. The influence of hypertension (high pressure) on the transport of lowdensity lipoprotein (LDL) through deformable arterial wall layers is also studied. This is due to the pressure-driven convective flow across the arterial wall. The intima and media layers are assumed as homogeneous porous media. The results of the present study reveal that ATP concentration over the surface of SMCs and within the bulk of the media layer is significantly dependent on the distribution of cells. Moreover, the shear stress magnitude and distribution over the SMC surface are affected by transmural pressure and the deformation of the media layer of the aorta wall. This work reflects the fact that the second or even subsequent layers of SMCs may bear shear stresses of the same order of magnitude as the first layer does if cells are arranged in an arbitrary manner. This study has brought new insights into the simulation of the arterial wall, as the previous simplifications have been ignored. The configurations of SMCs used here with elliptic cross sections of SMCs closely resemble the physiological conditions of cells. Moreover, the deformation of SMCs with high transmural pressure which follows the media layer compaction has been studied for the first time. On the other hand, results demonstrate that LDL concentration through the intima and media layers changes significantly as wall layers compress with transmural pressure. It was also noticed that the fraction of leaky junctions across the endothelial cells and the area fraction of fenestral pores over the internal elastic lamina affect the LDL distribution dramatically through the thoracic aorta wall. The simulation techniques introduced in this work can also trigger new ideas for simulating porous media involved in any biomedical, biomechanical, chemical, and environmental engineering applications.

A Thermodynamically Consistent Damage Model for Advanced Composites

A continuum damage model for the prediction of damage onset and structural collapse of structures manufactured in fiber-reinforced plastic laminates is proposed. The principal damage mechanisms occurring in the longitudinal and transverse directions of a ply are represented by a damage tensor that is fixed in space. Crack closure under load reversal effects are taken into account using damage variables established as a function of the sign of the components of the stress tensor. Damage activation functions based on the LaRC04 failure criteria are used to predict the different damage mechanisms occurring at the ply level. The constitutive damage model is implemented in a finite element code. The objectivity of the numerical model is assured by regularizing the dissipated energy at a material point using Bazant’s Crack Band Model. To verify the accuracy of the approach, analyses ofcoupon specimens were performed, and the numerical predictions were compared with experimental data

Polymeeripinnoitteen aiheuttamat itseherätteiset värähtelyt nipillisessä telasysteemissä

Lukuisissa teollisuussovelluksissa materiaalien, kuten paperin ja teräslevyjen, muokkaamiseen käytettävät pyörivät nippitelat kärsivät aina erilaisten herätteiden synnyttämistä mekaanisista värähtelyistä, jotka voivat aiheuttaa virheitä valmistettaviin tuotteisiin. Tässä työssä tutkittiin viskoelastisia polymeerejä ja polymeeripinnoitteen nipilliseen telasysteemiin synnyttämiä haitallisia itseherätteisiä värähtelyjä. Työn polymeerejä käsittelevässä kirjallisuusosassa luotiin katsaus amorfisten polymeerien fysikaalisiin ominaisuuksiin. Kokeellisessa osuudessa tutkittiin tarkemmin kahden amorfisen telapinnoitepolymeerin termoreologisia ja mekaanisia ominaisuuksia suoritettujen DMTA-mittausten perusteella. Sovittamalla toisen polymeerin master-käyrään yleistetty lineaarisen standardiaineen malli saatiin selville polymeerin mekaaniset parametrit ja approksimaatio sen relaksaatiospektrille. Telapinnoitteen nipilliseen systeemiin synnyttämiä itseherätteisiä värähtelyjä ja niiden seurauksia tarkasteltiin kahdelle telalle ja polymeeripinnoitteelle kehitetyn analyyttisen mallin ja numeeristen laskujen avulla. Pinnoite mallinnettiin lineaarisen standardiaineen mukaisesti. Telasysteemin parametrit määritettiin DMTA-mittaustuloksista ja systeemiä vastaavasta koelaitteesta kokeellisella moodianalyysillä ja elementtimenetelmällä. Numeerisesta stabiilisuusanalyysistä ja liikeyhtälöiden integroinneista saadut tulokset kertovat telapinnoitteen aaltomaisista deformaatiomuodoista ja niiden synnyttämistä taajuusalueittain esiintyvistä epästabiileista värähtelyistä. Telasysteemi on epästabiili pinnoitedeformaatiokuvion systeemiin aiheuttaman herätevoiman taajuuden ollessa lähellä systeemin korkeampaa ominaistaajuutta. Numeerisista tuloksista voitiin ennustaa nopean ja hitaan barringin olemassaolo.