Towards the stabilization of the low density elements in topology optimization with large deformation

This work addresses the treatment of lower density regions of structures undergoing large deformations during the design process by the topology optimization method (TOM) based on the finite element method. During the design process the nonlinear elastic behavior of the structure is based on exact kinematics. The material model applied in the TOM is based on the solid isotropic microstructure with penalization approach. No void elements are deleted and all internal forces of the nodes surrounding the void elements are considered during the nonlinear equilibrium solution. The distribution of design variables is solved through the method of moving asymptotes, in which the sensitivity of the objective function is obtained directly. In addition, a continuation function and a nonlinear projection function are invoked to obtain a checkerboard free and mesh independent design. 2D examples with both plane strain and plane stress conditions hypothesis are presented and compared. The problem of instability is overcome by adopting a polyconvex constitutive model in conjunction with a suggested relaxation function to stabilize the excessive distorted elements. The exact tangent stiffness matrix is used. The optimal topology results are compared to the results obtained by using the classical Saint Venant–Kirchhoff constitutive law, and strong differences are found.

The evolution of high-temperature mylonitic microfabrics: evidence from simple shearing of a quartz analogue (norcamphor)

Plane strain simple shearing of norcamphor (C7H10O) in a see-through deformation rig to a shear strain of γ = 10.5 at a homologous temperature of Th = 0.81 yields a microfabric similar to that of quartz in amphibolite facies mylonite. Synkinematic analysis of the norcamphor microfabric reveals that the development of a steady-state texture is linked to changes in the relative activities of several grain-scale mechanisms. Three stages of textural and microstructural evolution are distinguished: (1) rotation and shearing of the intracrystalline glide planes are accommodated by localized deformation along three sets of anastomozing microshears. A symmetrical c-axis girdle reflects localized pure shear extension along the main microshear set (Sa) oblique to the bulk shear zone boundary (abbreviated as SZB); (2) progressive rotation of the microshears into parallelism with the SZB increases the component of simple shear on the Sa microshears. Grain-boundary migration recrystallization favours the survival of grains with slip systems oriented for easy glide. This is associated with a textural transition towards two stable c-axis point maxima whose skeletal outline is oblique with respect to the Sa microshears and the SZB; and (3) at high shear strains (γ > 8), the microstructure, texture and mechanism assemblage are strain invariant, but strain continues to partition into rotating sets of microshears. Steady state is therefore a dynamic, heterogeneous condition involving the cyclic nucleation, growth and consumption of grains.

Two-dimensional mesh optimization in the finite element method

The solution to the problem of finding the optimum mesh design in the finite element method with the restriction of a given number of degrees of freedom, is an interesting problem, particularly in the applications method. At present, the usual procedures introduce new degrees of freedom (remeshing) in a given mesh in order to obtain a more adequate one, from the point of view of the calculation results (errors uniformity). However, from the solution of the optimum mesh problem with a specific number of degrees of freedom some useful recommendations and criteria for the mesh construction may be drawn. For 1-D problems, namely for the simple truss and beam elements, analytical solutions have been found and they are given in this paper. For the more complex 2-D problems (plane stress and plane strain) numerical methods to obtain the optimum mesh, based on optimization procedures have to be used. The objective function, used in the minimization process, has been the total potential energy. Some examples are presented. Finally some conclusions and hints about the possible new developments of these techniques are also given.

Analytical bearing capacity of strip footing in weightless materials with power-law failure criteria

Sokolovskii’s method of characteristics is extended to provide analytical solutions for the ultimate load at the moment of plastic failure under plane-strain conditions of shallow strip foundations on weightless rigid-plastic media with a noncohesive power-law failure envelope. The formulation is made parametrically in terms of the instantaneous friction angle, and the key idea to obtain the bearing capacity is that information can be transmitted from the free surface (where external loads are known) to the contact plane of the foundation. The methodology can consider foundations adjacent to a slope, external surcharges at the free surface, and inclined loads (both on the slope and on the foundation). Sensitivity analyses illustrate the influence on bearing capacity of changes in the different geometrical parameters involved. An application example is presented and design plots are provided, and model predictions are compared with results of bearing capacity tests under low gravity.

Benchmarking analogue models of brittle thrust wedges

We performed a quantitative comparison of brittle thrust wedge experiments to evaluate the variability among analogue models and to appraise the reproducibility and limits of model interpretation. Fifteen analogue modeling laboratories participated in this benchmark initiative. Each laboratory received a shipment of the same type of quartz and corundum sand and all laboratories adhered to a stringent model building protocol and used the same type of foil to cover base and sidewalls of the sandbox. Sieve structure, sifting height, filling rate, and details on off-scraping of excess sand followed prescribed procedures. Our analogue benchmark shows that even for simple plane-strain experiments with prescribed stringent model construction techniques, quantitative model results show variability, most notably for surface slope, thrust spacing and number of forward and backthrusts. One of the sources of the variability in model results is related to slight variations in how sand is deposited in the sandbox. Small changes in sifting height, sifting rate, and scraping will result in slightly heterogeneous material bulk densities, which will affect the mechanical properties of the sand, and will result in lateral and vertical differences in peak and boundary friction angles, as well as cohesion values once the model is constructed. Initial variations in basal friction are inferred to play the most important role in causing model variability. Our comparison shows that the human factor plays a decisive role, and even when one modeler repeats the same experiment, quantitative model results still show variability. Our observations highlight the limits of up-scaling quantitative analogue model results to nature or for making comparisons with numerical models. The frictional behavior of sand is highly sensitive to small variations in material state or experimental set-up, and hence, it will remain difficult to scale quantitative results such as number of thrusts, thrust spacing, and pop-up width from model to nature.

Finite elastic analysis of an infinite plate with an elliptic hole.

1. Plane strain.

The effect of the reduction of carbon content on the toughness of high chromium white irons in the as-cast state

Three high chromium white cast irons were examined in the as-cast state to determine the effect of the carbon content on the fracture toughness. The plane strain fracture toughness K-Ic and the fracture strength were measured for each alloy. X-ray mapping was used to identify the phases on the fracture surfaces. Scanning electron fractography and optical microscopy were used to determine the volume fraction of each phase on the fracture surfaces. It was found that most fracture occurred in the eutectic carbides, but that for the alloys with a reduced volume fraction of eutectic carbides, a small amount of crack propagation occurred in the austenitic dendrites. This change in crack path correlated with an increase in fracture toughness. The Ritchie-Knott-Rice model of brittle fracture was applied. It was found to sensibly predict the critical length for fracture for each alloy. Deep etching was employed to examine the distribution of eutectic carbides. It was found that the eutectic carbides formed a continuous network in each case. (C) 2004 Kluwer Academic Publishers.

On the stress analysis of cracked bodies by means of finite elements

The work described in this thesis deals with the development and application of a finite element program for the analysis of several cracked structures. In order to simplify the organisation of the material presented herein, the thesis has been subdivided into two Sections : In the first Section the development of a finite element program for the analysis of two-dimensional problems of plane stress or plane strain is described. The element used in this program is the six-mode isoparametric triangular element which permits the accurate modelling of curved boundary surfaces. Various cases of material aniftropy are included in the derivation of the element stiffness properties. A digital computer program is described and examples of its application are presented. In the second Section, on fracture problems, several cracked configurations are analysed by embedding into the finite element mesh a sub-region, containing the singularities and over which an analytic solution is used. The modifications necessary to augment a standard finite element program, such as that developed in Section I, are discussed and complete programs for each cracked configuration are presented. Several examples are included to demonstrate the accuracy and flexibility of the technique.

Creep crack growth in cast steam turbine casing steels

The creep rupture properties of cast ½Cr½Mo¼V and 1Cr1Mo¼V alloy steel used in the manufacture of power station steam generating plant. have been investigated. The effects of constraint and geometry on the creep rupture properties are also considered. The validity of various criteria controlling macroscopic creep crack growth in cast CrMoV alloys has been examined. It is found that neither the stress intensity factor nor reference stress correlate satisfactorily the creep crack growth rates at the test temperature of 550°C. Certain minimum displacements must be achieved for crack initiation and propagation. It is found that this displacement as measured by crack opening displacement or crack aspect ratio, is the same in both compact tension and centre-cracked panel geometries, is invariant with crack length and decreases with increasing constraint. The effect of constraint on creep crack growth rate in the two geometries is less conclusive. A new model describing creep crack growth in cast CrMoV alloy steels has been developed. The model is based on the results from a numerical finite element creep analysis of the relaxation and redistribution of stress ahead of an incubating creep crack . It is found that macroscopic creep crack growth in a material undergoing either plane stress or plane strain deformation can be described by a fracture stress which is based on the Von Mises equivalent stress. It has been shown that this model is capable of rationalising all of the experimental crack velocity data from the cast CrMoV alloys. The resultant degree of data correlation is far superior to that obtained when using the stress intensity factor or reference stress. A cumulative damage creep fracture model based upon the results from the numerical analysis has been developed. It is found that the model is capable of predicting the behaviour of propagating creep cracks in cast CrMoV alloys from smooth bar creep rupture data.

Ductile fracture mechanisms in a structural steel

Microstructural fracture processes in a BS4360 Grade 50D structural steel with lower sulphur content were studied in smooth tensile specimen tests and Charpy-size bend bar tests. Based on the experimental analysis, an experimental void growth relation with the plastic strain and stress triaxiality and multiplying factor on void growth were determined. Experimental results show that the void growth relation can be reasonably used to estimate the constraint in the specimens containing the notch or crack, also they can be used to evaluate the variations of the stress triaxiality in front of the notch and crack tip under general yielding condition. Side-grooves obviously increase the constraint of the CVN specimens. Strain hardening leads to increasing the stress triaxiality, and decelerating the net void growth. This is especially true for the values of stress triaxiality more than about one. Additionally, the effect of the stress triaxiality on the critical void growth corresponding to the onset of ductile tearing was preliminarily investigated. In this work, a large number of smaller specimens were tested to investigate the ductile-brittle transition behaviour of the structural steel. A void growth rate explanation was suggested for evaluating the temperature transition behaviour. The elastic-plastic fracture tough-ness values based on small specimen tests, such as pre-cracked side-grooved bending specimen and short bar tensile specimen, may give large overestimates of the plane strain fracture toughness.

The effect of microstructure on the fracture toughness of a metal-matrix composite

The fracture behaviour and plane strain fracture toughness, KIC, of four 8090-based metal-matrix composites containing 20 weight % SiC particles, 3, 6 and 23 μm in diameter, has been evaluated as a function of matrix ageing condition. Toughness values are found to be almost independent of reinforcement size. Ageing at 170°C results in a monotonic decrease in toughness with increasing strength up to the peak condition, with no subsequent recovery in toughness on overageing. However, unlike reinforced 8090, the composites are not found to be susceptible to intergranular embrittlement on overageing. The observed trends are found to be independent of reinforcement size. These findings are explained in terms of the strength, work hardening behaviour and nature and distribution of void-nucleating particles in the matrix. © 1993.

Crack tip stress study for elastic-perfectly plastic materials with some applications

The problems of plasticity and non-linear fracture mechanics have been generally recognized as the most difficult problems of solid mechanics. The present dissertation is devoted to some problems on the intersection of both plasticity and non-linear fracture mechanics. The crack tip is responsible for the crack growth and therefore is the focus of fracture science. The problem of crack has been studied by an army of outstanding scholars and engineers in this century, but has not, as yet, been solved for many important practical situations. The aim of this investigation is to provide an analytical solution to the problem of plasticity at the crack tip for elastic-perfectly plastic materials and to apply the solution to a classical problem of the mechanics of composite materials.^ In this work, the stresses inside the plastic region near the crack tip in a composite material made of two different elastic-perfectly plastic materials are studied. The problems of an interface crack, a crack impinging an interface at the right angle and at arbitrary angles are examined. The constituent materials are assumed to obey the Huber-Mises yielding condition criterion. The theory of slip lines for plane strain is utilized. For the particular homogeneous case these problems have two solutions: the continuous solution found earlier by Prandtl and modified by Hill and Sokolovsky, and the discontinuous solution found later by Cherepanov. The same type of solutions were discovered in the inhomogeneous problems of the present study. Some reasons to prefer the discontinuous solution are provided. The method is also applied to the analysis of a contact problem and a push-in/pull-out problem to determine the critical load for plasticity in these classical problems of the mechanics of composite materials.^ The results of this dissertation published in three journal articles (two of which are under revision) will also be presented in the Invited Lecture at the 7$\rm\sp{th}$ International Conference on Plasticity (Cancun, Mexico, January 1999). ^

Análise numérica de dutos enterrados sujeitos a elevação e recalque localizado

Problems associated to longitudinal interactions in buried pipelines are characterized as three-dimensional and can lead to different soil-pipe issues. Despite the progress achieved in research on buried pipelines, little attention has been given to the three-dimensional nature of the problem throughout the last decades. Most of researches simplify the problem by considering it in plane strain condition. This dissertation aims to present a study on the behavior of buried pipelines under local settlement or elevation, using three-dimensional simulations. Finite element code Plaxis 3D was used for the simulations. Particular aspects of the numerical modeling were evaluated and parametric analyzes were performed, was investigated the effects of soil arching in three-dimensional form. The main variables investigated were as follows: relative density, displacement of the elevation or settlement zone, elevated zone size, height of soil cover and pipe diameter/thickness ratio. The simulations were performed in two stages. The first stage was involved the validation of the numerical analysis using the physical models put forward by Costa (2005). In the second stage, numerical analyzes of a full-scale pipeline subjected to a localized elevation were performed. The obtained results allowed a detailed evaluation of the redistribution of stresses in the soil mass and the deflections along the pipe. It was observed the reduction of stresses in the soil mass and pipe deflections when the height of soil cover was decreased on regions of the pipe subjected to elevation. It was also shown for the analyzed situation that longitudinal thrusts were higher than vi circumferential trusts and exceeded the allowable stresses and deflections. Furthermore, the benefits of minimizing stress with technical as the false trench, compressible cradle and a combination of both applied to the simulated pipeline were verified

Seismic Response Analysis of a Full-Scale Base-Isolated Structure via Measurements and Modeling

The full-scale base-isolated structure studied in this dissertation is the only base-isolated building in South Island of New Zealand. It sustained hundreds of earthquake ground motions from September 2010 and well into 2012. Several large earthquake responses were recorded in December 2011 by NEES@UCLA and by GeoNet recording station nearby Christchurch Women's Hospital. The primary focus of this dissertation is to advance the state-of-the art of the methods to evaluate performance of seismic-isolated structures and the effects of soil-structure interaction by developing new data processing methodologies to overcome current limitations and by implementing advanced numerical modeling in OpenSees for direct analysis of soil-structure interaction.

This dissertation presents a novel method for recovering force-displacement relations within the isolators of building structures with unknown nonlinearities from sparse seismic-response measurements of floor accelerations. The method requires only direct matrix calculations (factorizations and multiplications); no iterative trial-and-error methods are required. The method requires a mass matrix, or at least an estimate of the floor masses. A stiffness matrix may be used, but is not necessary. Essentially, the method operates on a matrix of incomplete measurements of floor accelerations. In the special case of complete floor measurements of systems with linear dynamics, real modes, and equal floor masses, the principal components of this matrix are the modal responses. In the more general case of partial measurements and nonlinear dynamics, the method extracts a number of linearly-dependent components from Hankel matrices of measured horizontal response accelerations, assembles these components row-wise and extracts principal components from the singular value decomposition of this large matrix of linearly-dependent components. These principal components are then interpolated between floors in a way that minimizes the curvature energy of the interpolation. This interpolation step can make use of a reduced-order stiffness matrix, a backward difference matrix or a central difference matrix. The measured and interpolated floor acceleration components at all floors are then assembled and multiplied by a mass matrix. The recovered in-service force-displacement relations are then incorporated into the OpenSees soil structure interaction model.

Numerical simulations of soil-structure interaction involving non-uniform soil behavior are conducted following the development of the complete soil-structure interaction model of Christchurch Women's Hospital in OpenSees. In these 2D OpenSees models, the superstructure is modeled as two-dimensional frames in short span and long span respectively. The lead rubber bearings are modeled as elastomeric bearing (Bouc Wen) elements. The soil underlying the concrete raft foundation is modeled with linear elastic plane strain quadrilateral element. The non-uniformity of the soil profile is incorporated by extraction and interpolation of shear wave velocity profile from the Canterbury Geotechnical Database. The validity of the complete two-dimensional soil-structure interaction OpenSees model for the hospital is checked by comparing the results of peak floor responses and force-displacement relations within the isolation system achieved from OpenSees simulations to the recorded measurements. General explanations and implications, supported by displacement drifts, floor acceleration and displacement responses, force-displacement relations are described to address the effects of soil-structure interaction.

Experimental study of forces and energies during shearing of steel sheet with angled tools

Shearing is a fast and inexpensive method to cut sheet metal that has been used since the beginning of the industrialism. Consequently, published experimental studies of shearing can be found from over a century back in time. Recent studies, however, are due to the availability of low cost digital computation power, mostly based on finite element simulations that guarantees quick results. Still, for validation of models and simulations, accurate experimental data is a requisite. When applicable, 2D models are in general desirable over 3D models because of advantages like low computation time and easy model formulation. Shearing of sheet metal with parallel tools is successfully modelled in 2D with a plane strain approximation, but with angled tools the approximation is less obvious. Therefore, plane strain approximations for shearing with angled tools were evaluated by shear experiments of high accuracy. Tool angle, tool clearance, and clamping of the sheet were varied in the experiments. The results showed that the measured forces in shearing with angled tools can be approximately calculated using force measurements from shearing with parallel tools. Shearing energy was introduced as a quantifiable measure of suitable tool clearance range. The effects of the shearing parameters on forces were in agreement with previous studies. Based on the agreement between calculations and experiments, analysis based on a plane strain assumption is considered applicable for angled tools with a small (up to 2 degrees) rake angle.