Transient effect on ideally plastic stationary crack-tip fields

This paper analyses the transient effect on ideally plastic stationary crack tip fields under mode I plane strain conditions, when the inertial forces are not negligible. It is shown that the governing equation for such a problem can be expressed in formal simplicity when referred to a system of moving curvilinear coordinates, which is a generalization of the system defined by the slip-line field in quasi-static plasticity. A perturbation method of solving the equations is described and illustrated by application to problems of ideally plastic stationary crack tip fields when the inertial forces are not negligible.

Where Tori fear to tread : hypermassive neutron star remnants and absolute event horizons Or topics in computational general relativity

Computational general relativity is a field of study which has reached maturity only within the last decade. This thesis details several studies that elucidate phenomena related to the coalescence of compact object binaries. Chapters 2 and 3 recounts work towards developing new analytical tools for visualizing and reasoning about dynamics in strongly curved spacetimes. In both studies, the results employ analogies with the classical theory of electricity and magnitism, first (Ch. 2) in the post-Newtonian approximation to general relativity and then (Ch. 3) in full general relativity though in the absence of matter sources. In Chapter 4, we examine the topological structure of absolute event horizons during binary black hole merger simulations conducted with the SpEC code. Chapter 6 reports on the progress of the SpEC code in simulating the coalescence of neutron star-neutron star binaries, while Chapter 7 tests the effects of various numerical gauge conditions on the robustness of black hole formation from stellar collapse in SpEC. In Chapter 5, we examine the nature of pseudospectral expansions of non-smooth functions motivated by the need to simulate the stellar surface in Chapters 6 and 7. In Chapter 8, we study how thermal effects in the nuclear equation of state effect the equilibria and stability of hypermassive neutron stars. Chapter 9 presents supplements to the work in Chapter 8, including an examination of the stability question raised in Chapter 8 in greater mathematical detail.

Lyapunov theory for zeno stability

Zeno behavior is a dynamic phenomenon unique to hybrid systems in which an infinite number of discrete transitions occurs in a finite amount of time. This behavior commonly arises in mechanical systems undergoing impacts and optimal control problems, but its characterization for general hybrid systems is not completely understood. The goal of this paper is to develop a stability theory for Zeno hybrid systems that parallels classical Lyapunov theory; that is, we present Lyapunov-like sufficient conditions for Zeno behavior obtained by mapping solutions of complex hybrid systems to solutions of simpler Zeno hybrid systems defined on the first quadrant of the plane. These conditions are applied to Lagrangian hybrid systems, which model mechanical systems undergoing impacts, yielding simple sufficient conditions for Zeno behavior. Finally, the results are applied to robotic bipedal walking. © 2012 IEEE.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.

Unitarity of spin-2 theories with linearized Weyl symmetry in D=2+1 dimensions

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Boundary identification criteria and geometry analysis of Al 7050 stretch zone from elevation profiles

This article describes the development of a method for analysis of the shape of the stretch zone surface based on parallax measurement theory and using digital image processing techniques. Accurate criteria for the definition of the boundaries of the stretch zone are established from profiles of fracture surfaces obtained from crack tip opening displacement tests on Al-7050 alloy samples. The elevation profiles behavior analysis is based on stretch zone width and height parameters. It is concluded that the geometry of the stretch zone profiles under plane strain conditions can be described by a semi-parabolic relationship. (C) Elsevier B.V., 1999. All rights reserved.

Trypanosoma vivax infection dynamics in a cattle herd maintained in a transition area between Pantanal lowlands and highlands of Mato Grosso do Sul, Brazil

Trypanosoma vivax outbreaks in beef cattle in the Pantanal region of Mato Grosso do Sul state, Brazil, causes relevant economical impact due to weight loss, abortion and mortality. Cattle moved from the Pantanal to adjacent areas of this ecosystem for breeding and fattening is a common feature. Therefore an epidemiological study on breeding cows in the transition area between Pantanal lowland and adjacent highlands of Mato Grosso do Sul was performed to determine the T. vivax infection dynamics and outbreak risk. Three experimental groups were formed: Group 1 consisted of cows parasitologically negative by the Woo test and in the enzyme-linked immunosorbent assay for T. vivax antibody detection (Tv-ELISA-Ab); Group 2 parasitologically negative and positive in the Tv-ELISA-Ab; and in Group 3 cows were parasitologically positive and with positive reactions in the Tv-ELISA-Ab. During 24 months, the cows' dislodgment between the above established groups was monitored by Woo test and Tv-ELISA-Ab exams. The tabanid population was also monitored and the highest number occurred during the rainy season. Although parasitemias were detected only in the first four samplings of the experimental period, the cows could be considered as trypanotolerant, because no clinical signs were observed. Despite the higher T. vivax incidence during the dry season, no disease symptoms were seen. Even though T. vivax epidemiological situation in the herd was characterized as endemic with seasonal variation, the probability of outbreaks was null within the conditions of the study.

Relativistic quark spin coupling effects in the nucleon electromagnetic form factors

We investigate the effect of different forms of relativistic spin coupling of constituent quarks in the nucleon electromagnetic form factors. The four-dimensional integrations in the two-loop Feynman diagram are reduced to the null-plane, such that the light-front wave function is introduced in the computation of the form factors. The neutron charge form factor is very sensitive to different choices of spin coupling schemes, once its magnetic moment is fitted to the experimental value. The scalar coupling between two quarks is preferred by the neutron data, when a reasonable fit of the proton magnetic momentum is found. (C) 2000 Elsevier Science B.V.

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.

Modelo de flujo en un acuífero volcánico : el Acuífero de la Aldea (Gran Canaria)

RESUMEN. Se ha llevado a cabo un modelo de flujo de aguas subterráneas en la desembocadura del Barranco de La Aldea (Gran Canaria). El área fue discretizada tridimensionalmente en celdas de 50x50 m considerando 3 capas. La capa superior está constituida por materiales sedimentarios y volcánicos (aluvial, derrubios de ladera y basaltos alterados) y las capas intermedia e inferior por basaltos. Se ha realizado un modelo en régimen estacionario simulando el año hidrológico medio 1991/92 y transitorio para el período de 1991/92-1998/99. Los límites norte, sur y este se han definido como bordes impermeables, la línea de costa se ha definido como nivel constante y la cabecera del barranco se ha simulado mediante un tramo de caudal prefijado representando el aporte de la cuenca alta del barranco. Las entradas en la zona son: recarga por lluvia, retornos de riego, pérdidas en la red de abastecimiento, cabecera del barranco principal y desde la zona de intra-caldera. Las salidas son: bombeos y descarga al mar. El borde inferior se define por el flujo nulo en el contacto entre los basaltos alterados y sin alterar. En el cauce de los barrancos se ha impuesto una condición de dren y las extracciones se han localizado según los datos obtenidos de las captaciones de la zona. Los parámetros resultantes de la calibración del modelo, en particular, la transmisividad, son del mismo orden que los obtenidos en estudios previos modelo. Por otro lado, los niveles calculados y observados presentan un buen ajuste y el balance hídrico resulta consistente. ABSTRACT. A groundwater flow model in La Aldea ravine lower part (Gran Canaria) has been developed. The zone has been tridimensionally discretized as cells of 50 x 50 m considering 3 layers. The superficial layer is formed by sedimentary and volcanic materials (Alluvial, screes and altered basalts) and the intermediate and lower layers are basalts. The model has been developed in stationary state for the average hydrologic year 1991/92 and in transitory state for the period of 1991/92-1998/99. The North, South and East limits have been defined as null flow boundary conditions, the coast line has been defined as constant level and the ravine bed at the east has been defined as a constant flow, representing the contribution from the upper ravine basin. Recharge is a result of rainfall, irrigation returns, supply network leaks and inflow from the intra-caldera zone. Discharge takes place by pumping wells and flows towards the sea. The bottom surface is defined as a null flow condition in the limit between altered and unaltered basalts. A drain condition has been imposed in the ravine and the pumping wells extraction has been located. The simulation results indicate that the transmisivities obtained in previous works present the same order of magnitude than the obtained in the model and the calculated levels are in good agreement with the observed levels measurements.

A new insight in stress recovery tecniques and hessian reconstruction methods

Stress recovery techniques have been an active research topic in the last few years since, in 1987, Zienkiewicz and Zhu proposed a procedure called Superconvergent Patch Recovery (SPR). This procedure is a last-squares fit of stresses at super-convergent points over patches of elements and it leads to enhanced stress fields that can be used for evaluating finite element discretization errors. In subsequent years, numerous improved forms of this procedure have been proposed attempting to add equilibrium constraints to improve its performances. Later, another superconvergent technique, called Recovery by Equilibrium in Patches (REP), has been proposed. In this case the idea is to impose equilibrium in a weak form over patches and solve the resultant equations by a last-square scheme. In recent years another procedure, based on minimization of complementary energy, called Recovery by Compatibility in Patches (RCP) has been proposed in. This procedure, in many ways, can be seen as the dual form of REP as it substantially imposes compatibility in a weak form among a set of self-equilibrated stress fields. In this thesis a new insight in RCP is presented and the procedure is improved aiming at obtaining convergent second order derivatives of the stress resultants. In order to achieve this result, two different strategies and their combination have been tested. The first one is to consider larger patches in the spirit of what proposed in [4] and the second one is to perform a second recovery on the recovered stresses. Some numerical tests in plane stress conditions are presented, showing the effectiveness of these procedures. Afterwards, a new recovery technique called Last Square Displacements (LSD) is introduced. This new procedure is based on last square interpolation of nodal displacements resulting from the finite element solution. In fact, it has been observed that the major part of the error affecting stress resultants is introduced when shape functions are derived in order to obtain strains components from displacements. This procedure shows to be ultraconvergent and is extremely cost effective, as it needs in input only nodal displacements directly coming from finite element solution, avoiding any other post-processing in order to obtain stress resultants using the traditional method. Numerical tests in plane stress conditions are than presented showing that the procedure is ultraconvergent and leads to convergent first and second order derivatives of stress resultants. In the end, transverse stress profiles reconstruction using First-order Shear Deformation Theory for laminated plates and three dimensional equilibrium equations is presented. It can be seen that accuracy of this reconstruction depends on accuracy of first and second derivatives of stress resultants, which is not guaranteed by most of available low order plate finite elements. RCP and LSD procedures are than used to compute convergent first and second order derivatives of stress resultants ensuring convergence of reconstructed transverse shear and normal stress profiles respectively. Numerical tests are presented and discussed showing the effectiveness of both procedures.

Konstruktion und Analyse einer funktionalen Renormierungsgruppengleichung für Gravitation im Einstein-Cartan-Zugang

Während das Standardmodell der Elementarteilchenphysik eine konsistente, renormierbare Quantenfeldtheorie dreier der vier bekannten Wechselwirkungen darstellt, bleibt die Quantisierung der Gravitation ein bislang ungelöstes Problem. In den letzten Jahren haben sich jedoch Hinweise ergeben, nach denen metrische Gravitation asymptotisch sicher ist. Das bedeutet, daß sich auch für diese Wechselwirkung eine Quantenfeldtheorie konstruieren läßt. Diese ist dann in einem verallgemeinerten Sinne renormierbar, der nicht mehr explizit Bezug auf die Störungstheorie nimmt. Zudem sagt dieser Zugang, der auf der Wilsonschen Renormierungsgruppe beruht, die korrekte mikroskopische Wirkung der Theorie voraus. Klassisch ist metrische Gravitation auf dem Niveau der Vakuumfeldgleichungen äquivalent zur Einstein-Cartan-Theorie, die das Vielbein und den Spinzusammenhang als fundamentale Variablen verwendet. Diese Theorie besitzt allerdings mehr Freiheitsgrade, eine größere Eichgruppe, und die zugrundeliegende Wirkung ist von erster Ordnung. Alle diese Eigenschaften erschweren eine zur metrischen Gravitation analoge Behandlung.rnrnIm Rahmen dieser Arbeit wird eine dreidimensionale Trunkierung von der Art einer verallgemeinerten Hilbert-Palatini-Wirkung untersucht, die neben dem Laufen der Newton-Konstante und der kosmologischen Konstante auch die Renormierung des Immirzi-Parameters erfaßt. Trotz der angedeuteten Schwierigkeiten war es möglich, das Spektrum des freien Hilbert-Palatini-Propagators analytisch zu berechnen. Auf dessen Grundlage wird eine Flußgleichung vom Propertime-Typ konstruiert. Zudem werden geeignete Eichbedingungen gewählt und detailliert analysiert. Dabei macht die Struktur der Eichgruppe eine Kovariantisierung der Eichtransformationen erforderlich. Der resultierende Fluß wird für verschiedene Regularisierungsschemata und Eichparameter untersucht. Dies liefert auch im Einstein-Cartan-Zugang berzeugende Hinweise auf asymptotische Sicherheit und damit auf die mögliche Existenz einer mathematisch konsistenten und prädiktiven fundamentalen Quantentheorie der Gravitation. Insbesondere findet man ein Paar nicht-Gaußscher Fixpunkte, das Anti-Screening aufweist. An diesen sind die Newton-Konstante und die kosmologische Konstante jeweils relevante Kopplungen, wohingegen der Immirzi-Parameter an einem Fixpunkt irrelevant und an dem anderen relevant ist. Zudem ist die Beta-Funktion des Immirzi-Parameters von bemerkenswert einfacher Form. Die Resultate sind robust gegenüber Variationen des Regularisierungsschemas. Allerdings sollten zukünftige Untersuchungen die bestehenden Eichabhängigkeiten reduzieren.

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.

Advanced thermal analysis of microelectronics using spreading resistance models

Thermal analysis of electronic devices is one of the most important steps for designing of modern devices. Precise thermal analysis is essential for designing an effective thermal management system of modern electronic devices such as batteries, LEDs, microelectronics, ICs, circuit boards, semiconductors and heat spreaders. For having a precise thermal analysis, the temperature profile and thermal spreading resistance of the device should be calculated by considering the geometry, property and boundary conditions. Thermal spreading resistance occurs when heat enters through a portion of a surface and flows by conduction. It is the primary source of thermal resistance when heat flows from a tiny heat source to a thin and wide heat spreader. In this thesis, analytical models for modeling the temperature behavior and thermal resistance in some common geometries of microelectronic devices such as heat channels and heat tubes are investigated. Different boundary conditions for the system are considered. Along the source plane, a combination of discretely specified heat flux, specified temperatures and adiabatic condition are studied. Along the walls of the system, adiabatic or convective cooling boundary conditions are assumed. Along the sink plane, convective cooling with constant or variable heat transfer coefficient are considered. Also, the effect of orthotropic properties is discussed. This thesis contains nine chapters. Chapter one is the introduction and shows the concepts of thermal spreading resistance besides the originality and importance of the work. Chapter two reviews the literatures on the thermal spreading resistance in the past fifty years with a focus on the recent advances. In chapters three and four, thermal resistance of a twodimensional flux channel with non-uniform convection coefficient in the heat sink plane is studied. The non-uniform convection is modeled by using two functions than can simulate a wide variety of different heat sink configurations. In chapter five, a non-symmetrical flux channel with different heat transfer coefficient along the right and left edges and sink plane is analytically modeled. Due to the edge cooling and non-symmetry, the eigenvalues of the system are defined using the heat transfer coefficient on both edges and for satisfying the orthogonality condition, a normalized function is calculated. In chapter six, thermal behavior of two-dimensional rectangular flux channel with arbitrary boundary conditions on the source plane is presented. The boundary condition along the source plane can be a combination of the first kind boundary condition (Dirichlet or prescribed temperature) and the second kind boundary condition (Neumann or prescribed heat flux). The proposed solution can be used for modeling the flux channels with numerous different source plane boundary conditions without any limitations in the number and position of heat sources. In chapter seven, temperature profile of a circular flux tube with discretely specified boundary conditions along the source plane is presented. Also, the effect of orthotropic properties are discussed. In chapter 8, a three-dimensional rectangular flux channel with a non-uniform heat convection along the heat sink plane is analytically modeled. In chapter nine, a summary of the achievements is presented and some systems are proposed for the future studies. It is worth mentioning that all the models and case studies in the thesis are compared with the Finite Element Method (FEM).

Numerical modelling of non-load bearing light gauge cold-formed steel frame walls under fire conditions

Recently an innovative composite panel system was developed, where a thin insulation layer was used externally between two plasterboards to improve the fire performance of light gauge cold-formed steel frame walls. In this research, finite-element thermal models of both the traditional light gauge cold-formed steel frame wall panels with cavity insulation and the new light gauge cold-formed steel frame composite wall panels were developed to simulate their thermal behaviour under standard and realistic fire conditions. Suitable apparent thermal properties of gypsum plasterboard, insulation materials and steel were proposed and used. The developed models were then validated by comparing their results with available fire test results. This article presents the details of the developed finite-element models of small-scale non-load-bearing light gauge cold-formed steel frame wall panels and the results of the thermal analysis. It has been shown that accurate finite-element models can be used to simulate the thermal behaviour of small-scale light gauge cold-formed steel frame walls with varying configurations of insulations and plasterboards. The numerical results show that the use of cavity insulation was detrimental to the fire rating of light gauge cold-formed steel frame walls, while the use of external insulation offered superior thermal protection to them. The effects of real fire conditions are also presented.