This study presents an approach based on the floating frame of reference method to model complex three-dimensional bodies in a multibody system. Unlike most of the formulations based on the floating frame of reference method, which assume small or moderate deformations, the present formulation allows large elastic deformations within each frame by using the co-rotational form of the updated Lagrangian description of motion. The implicit integration scheme is based on the Generalized-alpha method, and kinematic joints are invoked in the formulation through the coordinate partitioning method. The resulting numerical scheme permits the usage of relatively large time steps even though the flexible bodies may experience large elastic deformations. A triangular element, based on the first order shear deformable theory, has been developed specifically for folded plate and shell structures. The plate element does not suffer from either shear or aspect-ratio locking under transverse and membrane bending, respectively. A stiffened plate element has been developed that combines a shear deformable plate with a Timoshenko beam. A solid element, that utilized the isoparametric formulation along with incompatible modes, and one-dimensional elements are also included in the element library. The tools developed in the present work are then utilized for detailed rotorcraft applications. As opposed to the conventional approach of using beam elements to represent the rotor blade, the current approach focuses on detailed modeling of the blade using plate and solid elements. A quasi-steady model based on lifting line theory is utilized to compute the aerodynamic loads on the rotor blade in order to demonstrate the capabilities of the proposed tool to model rotorcraft aeroelasticity.

Shear failure in concrete beams are often brittle in nature and potentially dangerous without adequatereinforcing measures. In design of concrete, it is commonly recommended to install transversalreinforcement along the shear span to induce a more ductile structural response, improving the shearcapacity all together and providing sufficient warning prior to collapse. However, it is more frequentlybeing assessed whether analogous performance can be achieved in fiber reinforced concrete beamswithout stirrups, and multiple attempts in literature confirm that it is possible. This alternative technologyintroduces need for better understanding of the modeling aspects of FRC in numerical simulations, as it isbecoming more common for engineers to resort to the finite element method in quality assurance ofstructures.In this thesis, the possibility of predicting shear failure numerically in simply supported fiber reinforcedconcrete beams with flexural bar reinforcement but without stirrups was investigated by means ofnonlinear finite element analysis, using the software package ATENA 2D Engineering. The ultimate aimwas to, as accurately as possible by means of numerical analyses on representative FE-models, replicatethe results from physical three-point-bending tests on simply supported FRC beams of various sizesperformed by Minelli et al. (2014). These beams were merely equipped with flexural reinforcement andexhibited shear failure.This thesis revolved around development and comparative assessment of material models for FRC basedon the smeared crack approach, adopting two different strategies: (1) The first strategy was to calibratematerial parameters based on results from 3PBT on notched FRC beams that were carried out prior totesting of the reinforced FRC beams, as reported by Minelli et al. (2014). Nonlinear finite element analysiswas used on representative FE-models for the notched 3PBT specimens, from which material parameterswere obtained iteratively by employing inverse analysis methods proposed by Červenka Consulting s.r.o.(2). The second strategy comprised of utilizing recommended constitutive relations from designrecommendations in SS812310 and RILEM TC 162-TDF. All of the constructed material models werefinally coupled with the FE-models that represented the beams with flexural reinforcement for evaluationof their performance based on their consistency with experiment data.It was found that the material models that were generated from inverse analysis in general would haveyielded successful predictions for the occurrence of shear failure in the reinforced FRC beams, providedthat the governing post-cracking residual tensile parameters were processed with respect to relevantassumptions as to describe uniaxial tensile behavior. However, although it was possible to utilize theproposed calibration method to replicate the load-displacement data for the notched 3PBT specimens withsufficient conformity, it was not possible to arrive at only one unique solution. Instead, multiple outcomescould be obtained based on the initial choice for the input value of the uniaxial tensile strength, leading tothe conclusion that experience and the engineering judgment of the user is of high importance whenadopting this method.Regarding the material models that were derived from constitutive relations in design recommendations,satisfactory estimates for the shear capacity could be obtained from the FE-models that were based onrecommendations by RILEM. The models that were based on SS812310, on the other hand, demonstratedover-stiff behavior and they were unable to provide accurate graphical visualizations of characteristicshear cracking, although the obtained load bearing capacity overall matched the experiment data in caseswhen size effects seemingly had a minor influence. An important observation from the comparison ofthese material models was that the initial drop in tensile strength during crack initiation within an elementis crucial in modeling of FRC, as it accounts for a more realistic behavior through a gradual transitionfrom aggregate bridging mechanisms of PC to the added fiber bridging mechanisms of FRC. Forsituations with high residual tensile strengths in relation to tensile strength at crack initiation, theguidelines in SS812310 become less practical for predicting shear failure by means of NLFEA.

No description available.

Non-ductile Reinforced Concrete Frames (RCF) with and without Unreinforced Masonry (URM) infills can be found in many places around the world including the Western United States, Eastern Europe, Asia and Latin America. These structures can have an unsatisfactory seismic performance which may even lead to collapse due to brittle failure modes. Furthermore, the effect of the infills on the seismic response of the structural system is not always accounted for in analysis and design. At present, there is no consensus on whether masonry infills are beneficial (by increasing the resistance of the system) or detrimental (by leading to brittle failure modes) for RCF construction. This study focuses on the development of a simplified modeling approach for non-ductile RCF with URMI that combines the simplicity of strut-and-tie models with the accuracy of Nonlinear Finite Element Analysis (NLFEA). Despite the fact that NLFEA procedures are the most advanced way to address the structural analysis of RCF with URM infills, their conceptual complexity and computational cost may hinder their widespread adoption as an analysis and design tool. At the same time, simplified methods, such as those based on the equivalent strut concept, may be overly crude and neglect essential aspects of the nonlinear response. To address the need for an adequately accurate, but computationally and conceptually efficient analysis method, this study establishes a novel method for planar RCF with URM infills subjected to lateral loads. The method, which is based on the Nonlinear Truss Analogy (NLTA) is shown to have an accuracy comparable to that of NLFEA. Specifically, the method is shown to adequately capture the strength and stiffness degradation and the damage patterns while entailing a reduced computational cost (compared to that of NLFEA). The proposed method is expected to bridge the gap between overly crude equivalent strut models and computationally expensive NLFEA. / Ph. D.

No description available.

<p>Breast cancer affects 12.5\% of women over their life time and tissue expansion (TE) is the most common technique for breast reconstruction after mastectomy. However, the rate of complications with TE can be as high as 15\%. Even though the first documented case of TE happened in 1957, there has yet to be a standardized procedure established due to the variations among patients and the TE protocols are currently designed based on surgeon's experience. There are several studies of computational and theoretical framework modeling skin growth in TE but these tools are not used in the clinical setting. This dissertation focuses on bridging the gap between the already existing skin growth modeling efforts and it's potential application in the clinical setting.</p> <p><br></p> <p>We started with calibrating a skin growth model based on porcine skin expansions data. We built a predictive finite element model of tissue expansion. Two types of model were tested, isotropic and anisotropic models. Calibration was done in a probabilistic framework, allowing us to capture the inherent biological uncertainty of living tissue. We hypothesized that the skin growth rate was proportional to stretch. Indeed, the Bayesian calibration process confirmed that this conceptual model best explained the data. </p> <p><br></p> <p>Although the initial model described the macroscale response, it did not consider any activity on the cellular level. To account for the underlying cellular mechanisms at the microscopic scale, we have established a new system of differential equations that describe the dynamics of key mechanosensing pathways that we observed to be activated in the porcine model. We calibrated the parameters of the new model based on porcine skin data. The refined model is still able to reproduce the observed macroscale changes in tissue growth, but now based on mechanistic knowledge of the cell mechanobiology.  </p> <p><br></p> <p>Lastly, we demonstrated how our skin growth model can be used in a clinical setting. We created TE simulations matching the protocol used in human patients and compared the results with clinical data with good agreement. Then we established a personalized model built from 3D scans of a patient unique geometry. We verified our model by comparing the skin growth area with the area of the skin harvested in the procedure, again with good agreement.</p> <p><br></p> <p>Our work shows that skin growth modeling can be a powerful tool to aid surgeons design TE procedures before they are actually performed. The simulations can help with optimizing the protocol to guarantee the correct amount of skin is growth in the shortest time possible without subjecting the skin to deformations that can compromise the procedure.</p>

There exists several different methods that can be used to implement a level of safety when performing non-linear finite element analysis of a structure. These methods are called safety formats and they estimate safety by different means and formulas which are partly discussed further in this thesis. The aim of this master thesis is to evaluate a model uncertainty factor for one safety format method called the estimation of coefficient of variation method (ECOV) since it is suggested to be included in the next version of Eurocode. The ECOV method will also be compared with the most common and widely used safety format which is the partial factor method (PF). The first part of this thesis presents the different safety formats more thoroughly followed by a theoretical part. The theory part aims to provide a deeper knowledge for the finite element method and non-linear finite element analysis together with some beam theory that explains shear mechanism in different beam types. The study was conducted on six beams in total, three deep beams and three slender beams. The deep beams were previously tested in the 1970s and the slender beams were previously tested in the 1990s, both test series were performed in a laboratory. All beams failed due to shear in the experimental tests. A detailed description of the beams are presented in the thesis. The simulations of the beams were all performed in the FEM- programme ATENA 2D to obtain high resemblance to the experimental test. In the results from the simulations it could be observed that the ECOV method generally got a higher capacity than the PF method. For the slender beams both methods received rather high design capacities with a mean of about 82% of the experimental capacity. For the deep beams both method reached low design capacities with a mean of around 46% of the experimental capacity. The results regarding the model uncertainty factor showed that the mean value for slender beams should be around 1.06 and for deep beams it should be around 1.25.

Reinforced concrete (RC) moment frame structures make up a notable proportion of buildings in earthquake-prone regions in the United States and throughout the world. The beam-column (BC) joints are the most crucial regions in a RC moment frame structure as any deterioration of strength and/or stiffness in these areas can lead to global collapse of the structure. Thus, accurate simulations of the joint behavior are important for assessment of the local and global performance of both one-way and two-way interior BC joints. Such simulations can be used to study the flexural-shear-bond interaction, the failure modes, and sensitivity of various parameters of structural elements. Most of the existing analytical approaches for interior BC joints have either failed to account for the cyclic bond-slip behavior and the triaxial compressive state of confined concrete in the joint correctly or require so many calibrations on parameters as to render them impractical. The core motivation for this study is the need to develop robust models to test current design recommendations for 3D beam-column-slab subassemblies subjected to large drifts. The present study aims to first evaluate the flexural-shear-bond interactive behavior of two-way beam-column-slab interior connections by both finite element and nonlinear truss methodologies. The local performance such as bond-slip and strain history of reinforcing steel are compared with the experimental results for the first time. The reliability of applied finite element approach is evaluated against a series of one-way interior BC joints and a two-way interior beam-column-slab joint. The accuracy and efficiency of the nonlinear truss methodology is also evaluated by the same series of joints. Results show good agreement for finite element method against both global and local response, including hysteretic curve, local bond-slip development and beam longitudinal bar stress/strain distributions. The nonlinear truss model is also capable in obtaining satisfactory global response, especially in capturing large shear cracks. A parametric study is exhibited for a prototype two-way interior beam-column-slab joint described in an example to ACI 352R-02, to quantify several non-consensus topics in the design of interior BC connections, such as the joint shear force subjected to bidirectional cyclic loading, the development of bond-slip behavior, and the failure modes of two-way interior joints with slab. Results from connections with different levels of joint shear force subjected to unidirectional loading show that meeting the requirements from ACI 352 is essential to maintain the force transfer mechanism and the integrity of the joint. The connections achieved satisfactory performance under unidirectional loading, while the bidirectional monotonic loading decreases the joint shear force calculated by ACI 352 by 10%~26% based on current results. Poorer performance is obtained for wider beams and connections fail by shear in the joint rather than bond-slip behavior when subjected to bidirectional cyclic loading. In general, the study indicates that the ACI352-02 design methodology generally results in satisfactory performance when applied to 2D joints (planar) under monotonic and cyclic loads. Less satisfactory performance was found for cases of 3D joints with slabs. / Doctor of Philosophy / Reinforced concrete (RC) moment frames are one of the most popular structure types because of their economical construction and adaptable spaces. Moment frames consist of grid-like assemblages of vertical columns and horizontal beams joined by cruciform connections commonly labelled as beam-column joints. Because of the regularity of the grid and the ability to have long column spacing, moment frames are easy to form and cast and result in wide open bays that can be adapted and readapted to many uses. In RC structures, steel bars embedded in the concrete are used to take tensile forces, as concrete is relatively weak when loaded in tension. Forces are transferred between the steel and concrete components by so-called “bond” forces at the perimeter of the bars. The proper modeling of the behavior of bond forces inside the beam-column joints of reinforced concrete moment frames is the primary objective of this dissertation. Reinforced concrete moment frames constitute a notable proportion of the existing buildings in earthquake-prone regions in the United States and throughout the world. The beam-column joints are the most crucial elements in a RC moment frame structure as any deterioration of strength and/or stiffness in these areas can lead to global collapse of the structure. Physical experimentation is the most reliable means of studying the performance of beam-column joints. However, experimental tests are expensive and time-consuming. This is why computational simulation must always be used as a supplemental tool. Accurate simulations of the behavior of beam-column joints is important for assessment of the local and global behavior of beam-column joints. However, most of the existing analytical approaches for interior beam-column joints have either failed to account for the bond-slip behavior and the triaxial compressive state of confined concrete in the joint correctly or require so many calibration parameters as to render them impractical. The present study aims to provide reliable numerical methods for evaluating the behavior of two-way beam-column-slab interior joints. Two methods are developed. The v first method is a complex finite element model in which the beam-column joint is subdivided into many small 3D parts with the geometrical and material characteristics of each part carefully defined. Since the number of parts may be in the hundreds of thousands and the geometry and material behavior highly non-linear, setting up the problem and its solution of this problem requires large effort on the part of the structural engineer and long computation times in supercomputers. Finite element models of this type are generally accurate and are used to calibrate simpler models. The second method developed herein is a nonlinear truss analogy model. In this case the structure is modelled as nonlinear truss elements, or elements carrying only axial forces. When properly calibrated, this method can produce excellent results especially in capturing large shear cracks. To evaluate the accuracy and to quantify the current seismic design procedure for beam-column joints, a prototype two-way interior beam-column-slab joint described in an example to ACI 352R-02, the current design guide used for these elements in the USA, is analytically studied by the finite element methodology. The study indicates that the ACI352-02 design methodology generally results in satisfactory performance when applied to one-way (planar) joints under monotonic and cyclic loads. Less satisfactory performance was found for cases of three-dimensional (3D) joints with slabs.

Two types of connection are generally considered in the design of steel structures in practice. These are classified as completely rigid (moment) and simple (shear) connections. In theory, completely rigid connections can not undergo rotation and simple connections can not transfer moment. However, in reality rigid connections have a relative flexibility which makes them to rotate and simple connections have some reserve capacity to transfer moments. In many modern design specifications, this fact is realized and another type which is called partially restrained or semi-rigid connection is introduced. These types of connections have got the transfer of some beam moment to column together with shear. However, there is a lack of information on the amount of moment transferred and rotation of connection during the action of the moment transfer. The only way to quantify the moment and rotation of the partially restrained connections is to draw momentrotation curves. Nevertheless, drawing such curves requires great amount of expenses for experiments. Taking these into account, the use of finite elements with the help of increased computational power is one way to obtain moment-rotation curves of connections. Available test results guides the finite element analysis for justifications. So these analyses can be further implemented into design functions. This thesis is intended to conduct 3-D non-linear finite element analyses to compliment with tests results for different types of semi-rigid connections with angles and compare them with mathematical models developed by different researchers.

With increasing age of the concrete road bridges, the highly topical question is to determine their reliability and load-bearing capacity level required for the residual life of the structure. Doctoral thesis presents a comprehensive methodology for assessing the reliability of reinforced and prestressed concrete bridges based on non-linear finite element method damage and failure virtual simulations at both deterministic and stochastic levels. Load-bearing capacity values are specified by the structure’s design load capacity estimation by global safety factor methods or they are based on a fully probabilistic load capacity analysis using the direct resistance estimation. For the fully probabilistic calculations, the simulation technique Latin Hypercube Sampling is used.