The main topic of this thesis concerns the behavior of concrete pile cap supported by four piles with two varying positions of longitudinal reinforcements. The positions include top of piles and bottom of the pile cap. For this purpose, non-linear finite element models of a pile cap are created using software ATENA 3D. The goal was to observe which position of reinforcement yields the higher bearing capacity and to observe the failure modes in the models. To achieve the above goals, a short review of theoretical background concerning shear phenomena is performed. This, in order to enhance the knowledge regarding shear stresses, shear transfer mechanism, factors affecting shear capacity, modes of shear failure and relate them to the behavior of pile cap. Furthermore, the calculation of shear resistance capacity based on Eurocode 2 using strut and tie method and sectional approach is presented. The numerical analysis started by creating four pile cap models in ATENA 3D. The difference between the models being the position and ratio of longitudinal reinforcement. The purpose behind two reinforcement ratios were to observe the behavior of pile cap model in two cases: a) when failure occurs prior to yielding of reinforcement; b) when failure occurs while reinforcement is yielding. The models are then analyzed using software ATENA Studio. The results revealed that placing the reinforcement on top of piles in case (a) increased the capacity of the model by 23.5 % and in case (b) increased the capacity by 18.5 %. This because the tensile stresses were found to be concentrated on top of piles rather than the bottom of the pile cap. The final failure mode in the model with top reinforcement position was crushing of the inclined compressive strut at the node beneath the column and in the model with bottom reinforcement position, the splitting of the compressive strut due to tensile stresses developed perpendicular to the inclined strut. The potential advantage of placing the reinforcement at the bottom were a better crack control in serviceability limit state and a slightly less fragile failure mode compared to the top position of reinforcement. A parametric study was performed in the model as well to observe the effects of various parameters on the results obtained. It was found that fracture energy had the most significant effect on the results obtained. Finally, a comparison between the results of numerical analysis and analytical design approaches based on strut and tie method and sectional approach was performed. The comparison reveals that the design values obtained based on strut and tie method for the model were very conservative. In particular, the equation for the strength of inclined compressive strut based on Eurocode 2 was very general. / Det huvudsakliga ämnet för den här avhandlingen handlar om beteendet hos pålfundament som stöds av fyra pålar med två olika positioner av längsgående armering. Positionerna inkluderar toppen av pålarna och botten av slagdynan. För detta ändamål skapas icke-linjära finita elementmodeller av en slagdyna med mjukvaran ATENA 3D. Målet var att observera vilket armeringsläge som ger den högre bärkapaciteten och att identifiera brottmekanismen i modellerna. För att uppnå ovanstående mål utförs en kort genomgång av teoretisk bakgrund rörande skjuvningsfenomen. Detta för att förbättra kunskapen om skjuvspänningar, skjuvöverföringsmekanism, faktorer som påverkar skjuvkapacitet, skjuvbrott och relaterar dem till beteendet hos slagdynan. Beräkningen av skjuvmotståndet baserad på Eurocode2 med hjälp av Srut and tie-metod och sektionsmetod. Den numeriska analysen började med att skapa fyra pålfundament i ATENA 3D. Skillnaden mellan modellerna är positionen och förhållandet mellan längsgående armering. Syftet bakom två armeringsförhållanden var att observera beteendet hos slagdynan i två fall: a) när brott inträffar innan armering plasticeras; b) när brott inträffar medan armeringen plasticeras. Modellerna analyseras sedan med hjälp av programvaran ATENA Studio. Resultaten visade att placering av armeringen ovanpå pålarna i fall a) ökade modellens kapacitet med 23,5% och i fall (b) ökade kapaciteten med 18,5%. Detta på grund av att dragspänningarna visade sig vara koncentrerade på toppen av pålarna snarare än på botten av slagdynan. Det slutliga brottet i modellen med topparmeringsposition var krossning av det lutande tryckstaget vid noden under pelaren. I modellen med bottenarmeringsposition delades kompressionsstaget på grund av dragspänningar vinkelrätt mot det lutande staget. The potential advantage of placing the reinforcement at the bottom were a better crack control and slightly less fragile failure mode compared to the top position of reinforcement. En parametrisk studie genomfördes också i modellen för att observera effekterna av olika parametrar på de erhållna resultaten. Det visade sig att brottenergi hade den mest signifikanta effekten på de erhållna resultaten. Slutligen genomfördes en jämförelse mellan resultaten från numerisk analys och analytiska designmetoder baserade på strut and tie-metoden och sektionsmetoden. Jämförelsen avslöjar att de designvärden som erhölls baserat på strut and tie-metoden för modellen var mycket konservativa. I synnerhet var ekvationen för kapaciteten hos det lutande tryckstag baserad på Eurocode 2 mycket generell.

Implementing pre-stressing cables is a viable option aiming at controlling deformation and cracking of concrete flat slabs in serviceability limit state. The pre-stressing cables also contribute to punching shear capacity of the slab when they are located in vicinity of the column. The positive influence of pre-stressing cables on punching capacity of the concrete slabs is mainly due to the vertical component of inclined cables, compressive in-plane stresses and counter acting bending moments near the support region. The method presented in Eurocode 2 to determine the punching capacity of the pre-stressed concrete flat slabs considers the in-plane compressive stresses but totally neglects the effect of counter acting moments. The effect of vertical forces introduced by inclined cables is only considered when they are within the distance 2d from the face of the column. This area is called basic control area in the Eurocode 2. In this master thesis nonlinear finite element analysis is carried out to study the effect of pre-stressing cables on punching shear capacity of concrete slabs respecting the distance of cables from the face of the column. To attain this objective, the concrete damage plasticity model is implemented to model the concrete. The results indicate that until the distance of 6d from the face of the column the contribution of pre-stressing cables in punching shear capacity of slabs is significant. Furthermore, comparing the numerical results with the punching shear capacity of slabs predicted by Eurocode 2 reveals that Eurocode tremendously underestimates the punching shear capacity when the cables are located outside the basic control area.

In Nonlinear Finite Element Analysis (FEA) applied to structures, displacements at which the tangent stiffness matrix KT becomes singular are called critical points, and correspond to instabilities such as buckling or elastoplastic softening (e.g., necking). Prior to the introduction of Arc Length Methods (ALMs), critical points posed severe computational challenges, which was unfortunate since behavior at instabilities is of great interest as a precursor to structural failure. The original ALM was shown to be capable in some circumstances of continued computation at critical points, but limited success and unattractive features of the formulation were noted and addressed in extensive subsequent research. The widely used Crisfield Cylindrical and Spherical ALMs may be viewed as representing the 'state-of-the-art'. The more recent Stiff Arc Length method, which is attractive on fundamental grounds, was introduced in 2004, but without implementation, benchmarking or performance assessment. The present thesis addresses (a) implementation and (b) performance comparisons for the Crisfield and Stiff methods, using simple benchmarks formulated to incorporate elastoplastic softening. It is seen that, in contrast to the Crisfield methods, the Stiff ALM consistently continues accurate computation at, near and beyond critical points.

Global-local finite element analyses were used to study the damage tolerance of curvilinearly stiffened panels; fabricated using the modern additive manufacturing process, the so-called unitized structures, and that of adhesive joints. A damage tolerance study of the unitized structures requires cracks to be defined in the vicinity of the critical stress zone. With the damage tolerance study of unitized structures as the focus, responses of curvilinearly stiffened panels to the combined shear and compression loadings were studied for different stiffeners' height. It was observed that the magnitude of the minimum principal stress in the panel was larger than the magnitudes of the maximum principal and von Mises stresses. It was also observed that the critical buckling load factor increased significantly with the increase of stiffeners' height. To study the damage tolerance of curvilinearly stiffened panels, in the first step, buckling analysis of panels was performed to determine whether panels satisfied the buckling constraint. In the second step, stress distributions of the panel were analyzed to determine the location of the critical stress under the combined shear and compression loadings. Then, the fracture analysis of the curvilinearly stiffened panel with a crack of size 1.45 mm defined at the location of the critical stress, which was the common location with the maximum magnitude of the principal stresses and von Mises stress, was performed under combined shear and tensile loadings. This crack size was used because of the requirement of a sufficiently small crack, if the crack is in the vicinity of any stress raiser. A mesh sensitivity analysis was performed to validate the choice of the mesh density near the crack tip. All analyses were performed using global-local finite element method using MSC. Marc, and global finite element methods using MSC. Marc and ABAQUS. Negligible difference in results and 94% saving in the CPU time was achieved using the global-local finite element method over the global finite element method by using a mesh density of 8.4 element/mm ahead of the crack tip. To study the influence of different loads on basic modes of fracture, the shear and normal (tensile) loads were varied differently. It was observed that the case with the fixed shear load but variable normal loads and the case with the fixed normal load but variable shear loads were Mode-I. Under the maximum combined loading condition, the largest effective stress intensity factor was very smaller than the critical stress intensity factor. Therefore, considering the critical stress intensity factor of the panel with the crack of size 1.45 mm, the design of the stiffened panel was an optimum design satisfying damage tolerance constraints. To acquire the trends in stress intensity factors for different crack lengths under different loadings, fracture analyses of curvilinearly stiffened panels with different crack lengths were performed by using a global-local finite element method under three different load cases: a) a shear load, b) a normal load, and c) a combined shear and normal loads. It was observed that 85% data storage space and the same amount in CPU time requirement could be saved using global-local finite element method compared to the standard global finite element analysis. It was also observed that the fracture mode in panels with different crack lengths was essentially Mode-I under the normal load case; Mode-II under the shear load case; and again Mode-I under the combined load case. Under the combined loading condition, the largest effective stress intensity factor of the panel with a crack of recommended size, if the crack is not in the vicinity of any stress raiser, was very smaller than the critical stress intensity factor. This work also includes the performance evaluation of adhesive joints of two different materials. This research was motivated by our experience of an adhesive joint failure on a test-fixture that we used to experimentally validate the design of stiffened panels under a compression-shear load. In the test-fixture, steel tabs were adhesively bonded to an aluminum panel and this adhesive joint debonded before design loads on the test panel were fully applied. Therefore, the requirement of studying behavior of adhesive joints for assembling dissimilar materials was found to be necessary. To determine the failure load responsible for debonding of adhesive joints of two dissimilar materials, stress distributions in adhesive joints of the nonlinear finite element model of the test-fixture were studied under a gradually increasing compression-shear load. Since the design of the combined load test fixture was for transferring the in-plane shear and compression loads to the panel, in-plane loads might have been responsible for the debonding of the steel tabs, which was similar to the results obtained from the nonlinear finite element analysis of the combined load test fixture. Then, fundamental studies were performed on the three-dimensional finite element models of adhesive lap joints and the Asymmetric Double Cantilever Beam (ADCB) joints for shear and peel deformations subjected to a loading similar to the in-plane loading conditions in the test-fixtures. The analysis was performed using ABAQUS, and the cohesive zone modeling was used to study the debonding growth. It was observed that the stronger adhesive joints could be obtained using the tougher adhesive and thicker adherends. The effect of end constraints on the fracture resistance of the ADCB specimen under compression was also investigated. The numerical observations showed that the delamination for the fixed end ADCB joints was more gradual than for the free end ADCB joints. Finally, both the crack propagation and the characteristics of adhesive joints were studied using a global-local finite element method. Three cases were studied using the proposed global-local finite element method: a) adhesively bonded Double Cantilever Beam (DCB), b) an adhesive lap joint, and c) a three-point bending test specimen. Using global-local methods, in a crack propagation problem of an adhesively bonded DCB, more than 80% data storage space and more than 65% CPU time requirement could be saved. In the adhesive lap joints, around 70% data storage space and 70% CPU time requirement could be saved using the global-local method. For the three-point bending test specimen case, more than 90% for both data storage space and CPU time requirement could be saved using the global-local method. / Ph. D.