Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.

Integral Transformations of Volterra Gaussian Processes

We study integral representations of Gaussian processes with a pre-specified law in terms of other Gaussian processes. The dissertation consists of an introduction and of four research articles. In the introduction, we provide an overview about Volterra Gaussian processes in general, and fractional Brownian motion in particular. In the first article, we derive a finite interval integral transformation, which changes fractional Brownian motion with a given Hurst index into fractional Brownian motion with an other Hurst index. Based on this transformation, we construct a prelimit which formally converges to an analogous, infinite interval integral transformation. In the second article, we prove this convergence rigorously and show that the infinite interval transformation is a direct consequence of the finite interval transformation. In the third article, we consider general Volterra Gaussian processes. We derive measure-preserving transformations of these processes and their inherently related bridges. Also, as a related result, we obtain a Fourier-Laguerre series expansion for the first Wiener chaos of a Gaussian martingale. In the fourth article, we derive a class of ergodic transformations of self-similar Volterra Gaussian processes.

Investigation of the cold crucible melting process: experimental and numerical study

The dynamic process of melting different materials in a cold crucible is being studied experimentally with parallel numerical modelling work. The numerical simulation uses a variety of complementing models: finite volume, integral equation and pseudo-spectral methods combined to achieve the accurate description of the dynamic melting process. Results show the temperature history of the melting process with a comparison of the experimental and computed heat losses in the various parts of the equipment. The free surface visual observations are compared to the numerically predicted surface shapes.

PERIODIC PERTURBATION of QUADRATIC SYSTEMS WITH TWO INFINITE HETEROCLINIC CYCLES

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

Time-periodic perturbation of a Lienard equation with an unbounded homoclinic loop

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

Periodic perturbations of quadratic planar polynomial vector fields

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Effective potential in curved space and cut-off regularizations

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Periodic orbits for a class of reversible quadratic vector field on R-3

For a class of reversible quadratic vector fields on R-3 we study the periodic orbits that bifurcate from a heteroclinic loop having two singular points at infinity connected by an invariant straight line in the finite part and another straight line at infinity in the local chart U-2. More specifically, we prove that for all n is an element of N, there exists epsilon(n) > 0 such that the reversible quadratic polynomial differential systemx = a(0) + a(1y) + a(3y)(2) + a(4Y)(2) + epsilon(a(2x)(2) + a(3xz)),y = b(1z) + b(3yz) + epsilon b(2xy),z = c(1y) +c(4az)(2) + epsilon c(2xz)in R-3, with a(0) < 0, b(1)c(1) < 0, a(2) < 0, b(2) < a(2), a(4) > 0, c(2) < a(2) and b(3) is not an element of (c(4), 4c(4)), for epsilon is an element of (0, epsilon(n)) has at least n periodic orbits near the heteroclinic loop. (c) 2007 Elsevier B.V. All rights reserved.

On the existence of infinite heteroclinic cycles in polynomial systems and its dynamic consequences

In this work we consider the dynamic consequences of the existence of infinite heteroclinic cycle in planar polynomial vector fields, which is a trajectory connecting two saddle points at infinity. It is stated that, although the saddles which form the cycle belong to infinity, for certain types of nonautonomous perturbations the perturbed system may present a complex dynamic behavior of the solutions in a finite part of the phase plane, due to the existence of tangencies and transversal intersections of their stable and unstable manifolds. This phenomenon might be called the chaos arising from infinity. The global study at infinity is made via the Poincare Compactification and the argument used to prove the statement is the Birkhoff-Smale Theorem. (c) 2004 WILEY-NCH Verlag GmbH & Co. KGaA, Weinheim.

Noncoaxial multivortices in the complex sine-Gordon theory on the plane

We construct explicit multivortex solutions for the complex sine-Gordon equation (the Lund-Regge model) in two Euclidean dimensions. Unlike the previously found (coaxial) multivortices, the new solutions comprise n single vortices placed at arbitrary positions (but confined within a finite part of the plane.) All multivortices, including the single vortex, have an infinite number of parameters. We also show that, in contrast to the coaxial complex sine-Gordon multivortices, the axially-symmetric solutions of the Ginzburg-Landau model (the stationary Gross-Pitaevskii equation) do not belong to a broader family of noncoaxial multivortex configurations.

Zweischleifenkorrekturen zum leptonischen Zerfall des Z-Bosons

Im Rahmen der vorliegenden Dissertation wurde, basierend auf der Parallel-/Orthogonalraum-Methode, eine neue Methode zur Berechnung von allgemeinen massiven Zweischleifen-Dreipunkt-Tensorintegralen mit planarer und gedrehter reduzierter planarer Topologie entwickelt. Die Ausarbeitung und Implementation einer Tensorreduktion fuer Integrale, welche eine allgemeine Tensorstruktur im Minkowski-Raum besitzen koennen, wurde durchgefuehrt. Die Entwicklung und Implementation eines Algorithmus zur semi-analytischen Berechnung der schwierigsten Integrale, die nach der Tensorreduktion verbleiben, konnte vollendet werden. (Fuer die anderen Basisintegrale koennen wohlbekannte Methoden verwendet werden.) Die Implementation ist bezueglich der UV-endlichen Anteile der Masterintegrale, die auch nach Tensorreduktion noch die zuvor erwaehnten Topologien besitzen, abgeschlossen. Die numerischen Integrationen haben sich als stabil erwiesen. Fuer die verbleibenden Teile des Projektes koennen wohlbekannte Methoden verwendet werden. In weiten Teilen muessen lediglich noch Links zu existierenden Programmen geschrieben werden. Fuer diejenigen wenigen verbleibenden speziellen Topologien, welche noch zu implementieren sind, sind (wohlbekannte) Methoden zu implementieren. Die Computerprogramme, die im Rahmen dieses Projektes entstanden, werden auch fuer allgemeinere Prozesse in das xloops-Projekt einfliessen. Deswegen wurde sie soweit moeglich fuer allgemeine Prozesse entwickelt und implementiert. Der oben erwaehnte Algorithmus wurde insbesondere fuer die Evaluation der fermionischen NNLO-Korrekturen zum leptonischen schwachen Mischungswinkel sowie zu aehnlichen Prozessen entwickelt. Im Rahmen der vorliegenden Dissertation wurde ein Grossteil der fuer die fermionischen NNLO-Korrekturen zu den effektiven Kopplungskonstanten des Z-Zerfalls (und damit fuer den schachen Mischungswinkel) notwendigen Arbeit durchgefuehrt.

Sigmoidal cosine series on the interval

We construct a set of functions, say, psi([r])(n) composed of a cosine function and a sigmoidal transformation gamma(r) of order r > 0. The present functions are orthonormal with respect to a proper weight function on the interval [-1, 1]. It is proven that if a function f is continuous and piecewise smooth on [-1, 1] then its series expansion based on psi([r])(n) converges uniformly to f so long as the order of the sigmoidal transformation employed is 0 < r

Fermat's Equation in Matrices

The Fermat equation is solved in integral two by two matrices of determinant one as well as in finite order integral three by three matrices.

The Behavior of compressively loaded stiffener runout specimens - Part II: Finite element analysis

The termination of stiffeners in composite aircraft structures give rise to regions of high interlaminar shear and peel stresses as the load in the stiffener is diffused into the skin. This is of particular concern in co-cured composite stiffened structures where there is a relatively low resistance to through-thickness stress components at the skin-stiffener interface. In Part I, experimental results of tested specimens highlighted the influence of local design parameters on their structural response. Indeed some of the observed behavior was unexpected. There is a need to be able to analyse a range of changes in geometry rapidly to allow the analysis to form an integral part of the structural design process.

This work presents the development of a finite element methodology for modelling the failure process of these critical regions. An efficient thick shell element formulation is presented and this element is used in conjuction with the Virtual Crack Closure Technique (VCCT) to predict the crack growth characteristics of the modelled specimens. Three specimens were modelled and the qualitative aspects of crack growth were captured successfully. The shortcomings in the quantitative correlation between the predicted and observed failure loads are discussed. There was evidence to suggest that high through-thickness compressive stresses enhanced the fracture toughness in these critical regions.

Efficient solution of multiple cracks in great number using eigen COD boundary integral equations with iteration procedure

A newly developed computational approach is proposed in the paper for the analysis of multiple crack problems based on the eigen crack opening displacement (COD) boundary integral equations. The eigen COD particularly refers to a crack in an infinite domain under fictitious traction acting on the crack surface. With the concept of eigen COD, the multiple cracks in great number can be solved by using the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix to determine all the unknown CODs step by step. To deal with the interactions among cracks for multiple crack problems, all cracks in the problem are divided into two groups, namely the adjacent group and the far-field group, according to the distance to the current crack in consideration. The adjacent group contains cracks with relatively small distances but strong effects to the current crack, while the others, the cracks of far-field group are composed of those with relatively large distances. Correspondingly, the eigen COD of the current crack is computed in two parts. The first part is computed by using the fictitious tractions of adjacent cracks via the local Eshelby matrix derived from the traction boundary integral equations in discretized form, while the second part is computed by using those of far-field cracks so that the high computational efficiency can be achieved in the proposed approach. The numerical results of the proposed approach are compared not only with those using the dual boundary integral equations (D-BIE) and the BIE with numerical Green's functions (NGF) but also with those of the analytical solutions in literature. The effectiveness and the efficiency of the proposed approach is verified. Numerical examples are provided for the stress intensity factors of cracks, up to several thousands in number, in both the finite and infinite plates.