The failure mechanism in ductile materials involves void nucleation, their growth and subsequent coalescence to form the fracture surface. The voids are generated due to fracture or debonding of second phase particles or at slip band intersections. The triaxial stress field prevailing around a crack tip and in the necking region strongly influences the growth of these voids. In the initial stages of deformation, these microscale voids are often sufficiently small so that they exist entirely within a single grain of a polycrystalline material. Further, single crystals are used in high technology applications like turbine blades. This motivates the need to study void growth in a single crystal while investigating ductile fracture. Thus, the objectives of this work are to analyze the interaction between a notch tip and void as well as the growth and coalescence of a periodic array of voids under different states of stress in ductile FCC single crystals. First, the growth of a cylindrical void ahead of a notch tip in ductile FCC single crystals is studied. To this end, 2D plane strain finite element simulations are carried out under mode I, small scale yielding conditions, neglecting elastic anisotropy. In most of these computations, the orientation of the FCC single crystal is chosen so that notch lies in the (010) plane, with notch front along the [101] direction and potential crack growth along [101]. This orientation has been frequently observed in experimental studies on fracture of FCC single crystals. Three equivalent slip systems are considered which are deduced by combining three pairs of 3D conjugate slip systems producing only in-plane deformation. Attention is focused on the effects of crystal hardening, ratio of void diameter to spacing from the notch on plastic flow localization in the ligament connecting the notch and the void as well as their growth. The results show strong interaction between slip shear bands emanating from the notch and angular sectors of single slip forming around the void leading to intense plastic strain development in the ligament. However, the ductile fracture processes are retarded by increase in hardening of the single crystal and decrease in ratio of void diameter to spacing from the notch. In order to examine the effect of crystal orientation, computations are performed with an orientation wherein the three effective slip systems are rotated about the normal to the plane of deformation. A strong influence of crystal orientation on near-tip void growth and plastic slip band development is observed. Further, in order to study the synergistic, cooperative growth of multiple voids ahead of the notchtip, an analysis is performed by considering a series of voids located ahead of the tip. It is found that enhanced void growth occurs at higher load levels as compared to the single void model. Next, the growth and coalescence of a periodic array of cylindrical voids in a FCC single crystal is analyzed under different stress states by employing a 2D plane strain, unit cell approach. The orientation of the crystal studied here considers [101] and [010] crystal directions along the minor and major principal stress directions, respectively. Three equivalent slip systems, similar to those in the notch and void simulations are taken into account. Fringe contours of plastic slip and evolution of macroscopic hydrostatic stress and void volume fraction are examined. A criterion for unstable void growth which leads to onset of void coalescence is established. The effects of various stress triaxialities, initial void volume fraction and hardening on void growth and coalescence is assessed. It is observed that plastic slip activity around the void intensifies with increase in stress triaxiality. The macroscopic hydrostatic stress increases with deformation, reaches a peak value and subsequently decreases rapidly. An increase in stress triaxiality enhances the macroscopic hydrostatic stress sustained by the unit cell and promotes void coalescence. The stress triaxiality also has a profound effect on the shape of the void profile. The values of critical void volume fraction and critical strain, which mark onset of void coalescence, decrease within crease in stress triaxiality. However, the onset of void coalescence is delayed by increase in hardening and decrease initial void volume fraction.

Magnesium alloys, owing to their low density and high specific strength, are potential candidates for structural applications in automotive and aerospace industry. While considerable research effort has been devoted in recent years to understand deformation twinning in these alloys and Mg single crystals, only few studies have been conducted on their fracture behavior. This issue assumes importance since some investigations have shown that Mg alloys may possess low fracture toughness (less than Al alloys). Therefore, a combined experimental and numerical study of fracture in Mg single crystals under mode-I loading is performed in this work. The fracture experiments are conducted using three point bend(TPB) specimens inside a scanning electron microscope(SEM) stage equipped with specially designed fixtures. Three crystallographic orientations are considered where c-axis [0001] is along the normal to the flat surface of the notch in the first two orientations, while in the third it is aligned with the notch front. In-situ electron back scattered diffraction (EBSD) observations are made in the region around the notch root to monitor the evolution of tensile twinning on the specimen free surface. Along with EBSD, optical metallography, fractography and surface profilometry are also performed on the specimens to obtain a comprehensive understanding on the micromechanics of fracture in Mg single crystals. From the EBSD data, it is noticed that all the orientations show profuse tensile twinning of {1012}-type. Further, in the first two orientations, basal and prismatic slip traces are identified along with secondary basal slip inside the twins. The growth of the most prominent twin is monitored as a function of load and it is found that its width saturates at around 120 -150 μm, while twins continue to nucleate farther away to accommodate plastic deformation. The 3D nature of twinning is examined by comparing distribution of twin traces and the average twin volume fraction at the free surface and the mid-plane. It is noted that in all the orientations crack initiation occurs before the attainment of peak load and the crack grows stably along twin-matrix interface. Further, zigzaging of the crack path occurs due to deflection of the crack at the twin-twin intersections. It is found that profuse tensile twinning is an important energy dissipating mechanism that enhances the toughness of the material. Indeed, the experimental results show that the energy release rate J versus load histories corroborate with evolution of average twin volume fraction around the notch root. In order to gain further insights on the mechanics of fracture in Mg single crystals, 3D finite element simulations are carried out using a crystal plasticity framework, which includes crystallographic slip and twinning. The predicted load-displacement curves, slip traces and tensile twinning activity from finite element analysis are in good corroboration with the experimental observations. The numerical results are used to understand the 3D nature of the crack tip stress, plastic slip and twin volume fraction distributions near the notch root. The occurrence of tensile twinning in all three orientations is rationalized from the distribution of nor-mal stress ahead of the notch tip. In particular, compressive normal stress beyond the plastic hinge point causes out-of-plane bulging that is accompanied by tensile twinning for the third orientation in which the c-axis is aligned along the specimen thickness. The above behavior emphasizes the importance of tensile twinning since this orientation has relevance to polycrystalline Mg alloys that have a basal texture.

Οι κατασκευές εδαφικής αντιστήριξης εξακολουθούν να βρίσκονται σε ευρύτατη χρήση, με διαρκώς αυξανόμενο ενδιαφέρον λόγω των απαιτήσεων των σύγχρονων έργων υποδομής αλλά και των αναγκών δόμησης σε πυκνό αστικό περιβάλλον. Το ενδιαφέρον εστιάζεται σε κατασκευαστικές λύσεις και μεθόδους σχεδιασμού που συνδυάζουν ασφάλεια και οικονομία. Η ανάλυση των συγκεκριμένων κατασκευών αντιμετωπίζει πλήθος δυσεπίλυτων προβλημάτων στο αντικείμενο της αλληλεπίδρασης εδάφους-κατασκευής που συχνά καθορίζουν τη συμπεριφορά του έργου. Η κατανόηση αυτών των μηχανισμών επιτρέπει το σχεδιασμό με μικρότερα περιθώρια αβεβαιότητας που οδηγούν σε οικονομικότερες και ορθολογικότερες λύσεις. Στην κατεύθυνση αυτή συμβάλει η παρούσα Διατριβή, με την ανάπτυξη αναλυτικών εργαλείων και θεωρητικών ευρημάτων που βοηθούν στην κατανόηση των μηχανισμών της αλληλεπίδρασης και στην εκτίμηση της συμπεριφοράς των τοίχων αντιστήριξης υπό συνδυασμένη βαρυτική και σεισμική φόρτιση. Έμφαση δίνεται στην παραγωγή απλών κλειστών λύσεων και μεθοδολογιών για τον υπολογισμό των εδαφικών ωθήσεων και τη στατική ανάλυση του συστήματος τοίχου εδάφους. Συγκεκριμένα, παράγονται λύσεις άνω και κάτω ορίου για ενδόσιμους τοίχους, οι οποίες, παρότι προσεγγιστικές, πλεονεκτούν έναντι των κλασικών εξισώσεων Coulomb και Mononobe-Okabe τις οποίες μπορούν να αντικαταστήσουν. Σε ειδικές περιπτώσεις, όπως η περίπτωση τοίχων προβόλων με πεπλατυσμένο πέλμα, οι προτεινόμενες λύσεις οδηγούν σε ακριβή αποτελέσματα που βασίζονται σε ένα γενικευμένο πεδίο τάσεων Rankine. Επίσης παρουσιάζονται επεκτάσεις τους οι οποίες επιτρέπουν τον υπολογισμό μη-υδροστατικών κατανομών ωθήσεων γαιών λαμβάνοντας υπόψη την κυματική διάδοση της σεισμικής διέγερσης στο επίχωμα, σύμφωνα με μια ορθότερη παραλλαγή της ιδέας των Steedman & Zeng και τις διαφορετικές κινηματικές συνθήκες που προέρχονται από την απόκριση του τοίχου με περιστροφή περί την κορυφή ή τη βάση σύμφωνα με την τεχνική της Dubrova. Για την περίπτωση ανένδοτων τοίχων παρουσιάζεται μεθοδολογία για τη δραστική απλοποίηση των διαθέσιμων ελαστοδυναμικών, κυματικών λύσεων, όπως αυτή των Veletsos & Younan, η οποία καταλήγει σε κλειστές μαθηματικές εκφράσεις για τον υπολογισμό των ωθήσεων. Τέλος, παρουσιάζονται νέα ευρήματα στην κατεύθυνση της μαθηματικής αντιμετώπισης του δυσεπίλυτου προβλήματος της οριακής ισορροπίας ριπιδίου τάσεων σε εδαφικό μέσο στο οποίο ενεργούν βαρυτικές και αδρανειακές δυνάμεις πεδίου. Η παρούσα εργασία συμβάλλει στην περαιτέρω διερεύνηση του προβλήματος το οποίο θεμελίωσαν θεωρητικά οι Levy, Boussinesq, von Karman και Caquot, μέσω της δραστικής (αλλά ακριβούς) απλοποίησης του σε μία μη-γραμμική συνήθη διαφορική εξίσωση, η οποία επιτρέπει την επίλυση με απλές αριθμητικές και ημιαναλυτικές τεχνικές. Πέρα από τα ακριβή αριθμητικά αποτελέσματα, η προτεινόμενη ανάλυση προσφέρει μια βαθύτερη εποπτεία στο πρόβλημα και ανοίγει το δρόμο για περαιτέρω διερεύνηση ή και επέκταση της μεθόδου πέρα από τα όρια της κλασικής οριακής ανάλυσης. Η αξιοπιστία των προτεινόμενων λύσεων ελέγχεται μέσω συγκρίσεων με καθιερωμένες λύσεις και πειραματικά δεδομένα από τη βιβλιογραφία, αλλά και πρόσφατα πειραματικά αποτελέσματα που παρήχθησαν από τον συγγραφέα και ερευνητές στη σεισμική τράπεζα του Πανεπιστημίου του Bristol του Ηνωμένου Βασιλείου. / Earth retaining structures are still in widespread use, with growing interest due to the demands of modern infrastructure and building needs in a dense urban environment. Building solutions and design methodologies that combine safety and economy are the objectives of modern research. Significant difficulties in the analysis of retaining structures arise from the soil-structure interaction nature of the problem that often prescribes its behavior. Understanding these mechanisms allows design under smaller uncertainties, leading to economical and rational solutions. The contribution of the present thesis consists of the development of analytical tools and theoretical findings, helpful in understanding the mechanisms of interaction and the behavior of walls under combined gravity and seismic loading. Emphasis is given to the derivation of simple closed-form solutions and methodologies for the calculation of earth pressures and the static analysis of wall-soil system. Specifically, approximate Lower and Upper Bound solutions are produced for the case of yielding walls, which are advantageous compared to the classical equations Coulomb and Mononobe-Okabe. In special cases, such as the L-shaped cantilever walls, these solutions lead to exact results, pertaining to a generalized Rankine stress field. Extensions of the above solutions are presented allowing the calculation of non-hydrostatic earth pressure distributions, due to the wave propagation of the seismic excitation in the backfill, according to a better variant of the Steedman & Zeng approach and different kinematic conditions of the wall rotating around the top or bottom, according to the technique of Dubrova. For the case of non-yielding walls, a new methodology for the drastic simplification of available wave solutions, such as the Veletsos & Younan, is presented which leads to closed-form expressions for the dynamic pressure calculation. Finally, new theoretical findings are presented for the mathematical treatment of the intractable problem of plastic limit equilibrium in soil medium subjected to gravitational and inertial forces field. This work contributes to the further investigation of the problem which is founded theoretically by Levy, Boussinesq, von Karman and Caquot, through the significant (but accurate) simplification to a single, non-linear ordinary differential equation, easier to handle by simple numerical and semi-analytical techniques. Apart from the exact numerical results, the proposed analysis provides a deeper physical insight, leading the way to further investigation or extension of the method beyond the classical limit analysis assumptions. The reliability of the proposed solutions is checked through comparisons with established solutions and experimental data from the literature and recent experimental results obtained by the author and researchers in the shake table laboratory of the University of Bristol, UK.

The introductions of large digital computers in the field of engineering have rendered possible the solution of a wide variety of problems without the need to violate the equilibrium and compatibility. The major requirement for such analysis is a good constitutive model that represents the stress strain behaviour of the materials in an accurate way. Nowadays for most of the geotechnical engineering applications the elastoplastic models like Mohr Coulomb model are widely used. All the existing constitutive models which represent the plastic behaviour of soil are developed from the fundamentals of classical theory of plasticity. The classical theory of plasticity is always associated with the concept of yield surface and potential surface to represent the plastic behaviour. The definition of yield surface depends on the location of the yield point. But in practical sense it is very difficult to find out the exact yield point for a material. The expression for yield and potential surfaces are simply mathematical expressions formulated for computational efficiency. Experimentally it is very difficult to find out the yield surface in the case of three dimensional stress spaces. The classical theory of plasticity is developed based on the mechanical process. It is believed that a theory which violates the thermodynamic principle is not able to represent the material behaviour accurately. the initial stage and combined to give the final state of stress. It was proved that the equation proposed by Wu and Wang (1983) can be used to represent the triaxial behaviour of sand very well. The dilation and densification behaviour can be predicted very well with the endochronic constitutive equations. The principal aim of this work is to implement the endochronic constitutive equation in the FLAC3D model library like any other constitutive model and validate it with the triaxial test data. After implementation and validation, the application of the particular constitutive model is extended to some practical geotechnical engineering problems like the stresses and displacements around an underground opening such as tunnels, surface settlement due to shallow tunneling, stress distribution below the footing, settlement analysis of footing on various foundation beds such as sand, clay and sand overlying clay bed, lateral displacement of the secant pile wall due to excavation and the force developed in the horizontal support etc. All the three problems validate the model with the analytical, experimental and field data respectively. The equation proposed by Wu and Wang (1983) is used for the present study. In order to validate the equation proposed by Wu and Wang (1983), MatLab programming is used. The hydrostatic, deviatoric and volumetric behaviour is obtained separately using the concerned equations. The equation is coded in the MatLab and analysis is done for a triaxial element test. Both drained and undrained analyses were done in order to study the behaviour. The pore pressure developed is captured separately with the equation proposed by Geoffrey et al (1975). The results obtained from the analysis of the MatLab were compared with that of the experimental results. The analysis shows that the equation captures the least plastic behaviour well for the triaxial element test. The dilation and densification behaviour obtained using the respective equation shows that it matches well with the experimental results. A parametric study is also conducted in MatLab to see how the parameters affect the stress strain and volumetric behaviour of the sand. The parametric study conducted with the MatLab shows that most of the parameters involved in the equation affect the plastic part of the stress strain curve rather than the initial elastic part. User defined constitutive model was written in visual C++ and compiled as DLL (Dynamic Link Library) files that will be loaded whenever it is needed in FLAC3D. In visual C++, header and source files were written by incorporating the constitutive equation proposed by Wu and Wang (1983), defining the variables and other functions etc, and a dynamic link library is created, which is then integrated into the 3D finite difference code FLAC3D using the CPPUDM module to simulate the stress strain behaviour of the materials. CPPUDM module is an additional option in FLAC3D to implement the user defined constitutive models. The visual C++ code was written in the form of incremental stress strain relationship. The model acts like any other constitutive model in the FLAC3D model library and can be loaded whenever it is required. For the validation of the model in FLAC3D, the data for the MatLab simulation were used. Both drained and undrained tests were simulated with the model. The results obtained from the analysis shows that by suitably selecting the parameters the model can simulate the stress strain behaviour of sand very well. The volumetric and deviatoric behaviour were observed and is matching well with the experimental data. In the case of the undrained test the pore pressure generation is well captured by the equation proposed by Geoffrey et al (1975). In urban areas the construction of shallow tunnels results in excessive settlements of the ground surface and thereby causes damage to the existing above ground structures. In order to minimize the settlements and to reduce the impact due to that, a prior analysis of the displacements and stresses around the opening is very important. Nowadays numerical analysis is widely used for the analysis of such structures. The most important requirement of such analysis is a constitutive model that can represent the unloading behaviour around the tunnel opening of sand very well. Here the endochronic constitutive model implemented in the FLAC3D model library is used to evaluate the stresses and displacements around the tunnel. In the analysis the tunnel is simulated as a cylindrical hole in an infinite medium with the in situ stress. The stresses at the springing line was observed and compared with the analytical solution. The results show that the results are matching well with the analytical results. The comparison of the results with that obtained using the Mohr Coulomb model is also done to see how the model differs from a widely used plastic model. By slightly adjusting the parameters the results obtained from both the models are in well agreement. The strain softening effect which is predominant around an underground opening due to the loosening of soil mass is well captured by the endochronic model compared to the Mohr Coulomb model. The settlement analysis shows that the model is almost in close agreement with the closed form solution proposed by Oteo & Sagaseta (1982) and the results obtained with the Mohr Coulomb model. The settlement trough formed for various shapes is wider and deeper than the Mohr Coulomb model. The vertical stress distribution around the opening of the tunnel is studied with varying the shape of the openings using the proposed constitutive model. The results obtained were compared with that of the Mohr Coulomb model. The slightly higher values in the case of endochronic model are basically due to its plastic nature. The displacement and stresses in the axial direction (along the excavation) is observed with the model. In the case of shallow tunnel the surface get influenced by the loosening of the soil mass which necessitates the use of the support system. The study shows that the model can be used for the simulation of underground opening like tunnel and will capture the behaviour well. Footings are structures used to support the buildings constructed above the ground. The settlement analysis of footings is very important when we consider the stability of the structures supported by it. The vertical stress distribution below the footing is studied using the endochronic constitutive model and compared with the analytical solution proposed by Boussinesq (1885). In the elastic range the model shows matching results with the Boussinesq’s solution. The settlement analysis of footing on various foundation beds such as sand, clay and sand overlying the clay bed were studied using the endochronic constitutive model implemented in the FLAC3D model library. The experimental data conducted in our lab (Sireesh (2006)) was used for the study. The results show that with the chosen parameters the results obtained with the endochronic model are in good agreement with the experimental data. The Mohr Coulomb model over predicts the results. This shows higher modulus value for the Mohr Coulomb model. By conducting the parametric study it was seen that by reducing the value of modulus for the Mohr Coulomb model, the results are in good agreement with the experimental value. The displacement and stress contours obtained for the two models were compared. By analyzing the displacement contours it is seen that the Mohr Coulomb model shows uniform settlement. In the case of endochronic model uniform settlement is observed for about 5% settlement that is in the elastic range. After a certain strain level the displacement contours are tilted more towards one side showing the rotational failure. Here the endochronic model captures the anisotropic behaviour associated with the materials like sand at higher strain level. This result is a concrete evidence that the model can capture the realistic behaviour very well compared to any other model. Even though the model developed is for sand its application can be extended to clay also. The size and shape of the footing is varied to study its effect on the pressure settlement curve. The analysis is done with square shape of 150mm side and circular shape of 150mm diameter. As there is not much variation in the area of influence, the shape has little influence on the pressure settlement curve. As the size of the footing increases the settlement increases for a given pressure. A parametric study is conducted by varying the modulus value used. The study shows that as the modulus value increases, the settlement reduces for a given bearing pressure. The endochronic model analysed with the lower modulus value shows that the model predicts the perfectly plastic behaviour, here the settlement increases for low value of bearing pressure. The application of endochronic model for the simulation of complex geotechnical engineering problems like footings is highly explored in the study. Nowadays most of the infrastructure facilities are concentrated towards the underground space. The excavation and construction of such structures in the urban areas results in damage to the existing above ground structures if the construction is done in close proximity to the structures. In the present study a staged excavation of an underground construction for the Bangalore metro project is simulated with the endochronic constitutive model. In the Bangalore metro project the excavation for the underground station is done at the cricket stadium site. At the site there are two major buildings such as the six storied Hindustan Aeronautical Limited building and 100 years old BSNL masonry building. To minimize the impact on these structures were a major concern during the construction of the work. A robust support system consists of secant pile walls, soldier piles and horizontal struts are installed at the site. The OSV method known as the Onsite Visualization and monitoring is conducted to minimize the damage to the existing structures and the accidents at the construction site. Sensors are connected to LEDs which show change in color when the displacements and forces cross the triggered value. The field instrumentation is done with inclinometers, tilt meters and strain gauges connected to the sensors to observe the lateral deformation of the secant pile wall, tilt of the HAL building and the forces developed in the horizontal struts. The monitoring of field data is done for a period of five months from July to November. From the analysis of the field observed data it is clear that the support system provided were strong enough to resist the forces due to unloading. The lateral deformation of the secant pile wall and the forces developed in the strut were numerically analysed using the endochronic constitutive model and the results were compared with the field monitored data. The results show that the model captures the behaviour very close to the field data for a slightly higher modulus than that reported in the geotechnical report (BMRC report). This may be due to the fact that the value of modulus calculated experimentally might have some error. The analysis with the Mohr Coulomb model shows that the model over predicts the results very close to the surface of the excavation. This indicates that the influence of load is more on that particular depth for the Mohr Coulomb model. But the stiffness of the lateral support system is uniform throughout the depth; the endochronic model predicts the result more accurately than the Mohr Coulomb model. The strut forces developed in the horizontal support system is observed using the two models. The strut forces in the field is affected by so many factors such as the temperature variation, stages of excavation and other live loads acting on the site, so an exact comparison with the field data is quite difficult. The analysis shows that even though it is difficult to simulate the exact three dimensional nature of the problem in the present study the endochronic constitutive model captures the behaviour very well. The results obtained shows that the endochronic constitutive model implemented in the FLAC3D model library provides a very promising solution like any other constitutive model. As the theory is based on the irreversible law of thermodynamics and the formulation of the constitutive equation are based on the internal energy concept it can represent the material behaviour in accordance with the laws of continuum mechanics. The anisotropic behaviour of soil at higher strain level is well represented through the footing problem. The endochronic constitutive model is a very simple one to simulate the stress strain behaviour of the materials without the concept of yield surface; the parameters used in the equation can be obtained directly from a single triaxial stress strain plot. This study highlights the importance of a model without the concept of yield surface to capture the stress strain behaviour of any materials. Since the model is of completely plastic nature it has its own uniqueness in capturing the material behaviour due to loading and unloading.