Computation of stress intensity factors in functionally graded materials using partition-of-unity meshfee method

This paper describes the computation of stress intensity factors (SIFs) for cracks in functionally graded materials (FGMs) using an extended element-free Galerkin (XEFG) method. The SIFs are extracted through the crack closure integral (CCI) with a local smoothing technique, non-equilibrium and incompatibility formulations of the interaction integral and the displacement method. The results for mode I and mixed mode case studies are presented and compared with those available in the literature. They are found to be in good agreement where the average absolute error for the CCI with local smoothing, despite its simplicity, yielded a high level of accuracy.

Element-Free Galerkin modelling of composite damage

Delamination and matrix cracking are routine damage mechanisms, observed by post-mortem analysis of laminated structures containing geometrical features such as notches or bolts. Current finite element tools cannot explicitly model an intralaminar matrix microcrack, except if the location of the damage is specified a priori. In this work, a meshless technique, the Element-Free Galerkin (EFG) method, is utilized for the first time to simulate delamination (interlaminar) and intralaminar matrix microcracking in composite laminates.

Element-Free Galerkin Modelling of Crack Migration in Composite Laminates

The development of a virtual testing environment, as a cost-effective industrial design tool in the design and analysis of composite structures, requires the need to create models efficiently, as well as accelerate the analysis by reducing the number of degrees of freedom, while still satisfying the need for accurately tracking the evolution of a debond, delamination or crack front. The eventual aim is to simulate both damage initiation and propagation in components with realistic geometrical features, where crack propagation paths are not trivial. Meshless approaches, and the Element-Free Galerkin (EFG) method, are particularly suitable for problems involving changes in topology and have been successfully applied to simulate damage in homogeneous materials and concrete. In this work, the method is utilized to model initiation and mixed-mode propagation of cracks in composite laminates, and to simulate experimentally-observed crack migration which is difficult to model using standard finite element analysis. N

Nichtüberlappende Gebietszerlegungsmethoden für lineare und quasilineare (monotone und nichtmonotone) Probleme

In dieser Arbeit werden nichtüberlappende Gebietszerlegungsmethoden einerseits hinsichtlich der zu lösenden Problemklassen verallgemeinert und andererseits in bisher nicht untersuchten Kontexten betrachtet. Dabei stehen funktionalanalytische Untersuchungen zur Wohldefiniertheit, eindeutigen Lösbarkeit und Konvergenz im Vordergrund. Im ersten Teil werden lineare elliptische Dirichlet-Randwertprobleme behandelt, wobei neben Problemen mit dominantem Hauptteil auch solche mit singulärer Störung desselben, wie konvektions- oder reaktionsdominante Probleme zugelassen sind. Der zweite Teil befasst sich mit (gleichmäßig) monotonen koerziven quasilinearen elliptischen Dirichlet-Randwertproblemen. In beiden Fällen wird das Lipschitz-Gebiet in endlich viele Lipschitz-Teilgebiete zerlegt, wobei insbesondere Kreuzungspunkte und Teilgebiete ohne Außenrand zugelassen sind. Anschließend werden Transmissionsprobleme mit frei wählbaren $L^{\infty}$-Parameterfunktionen hergeleitet, wobei die Konormalenableitungen als Funktionale auf geeigneten Funktionenräumen über den Teilrändern ($H_{00}^{1/2}(\Gamma)$) interpretiert werden. Die iterative Lösung dieser Transmissionsprobleme mit einem Ansatz von Deng führt auf eine Substrukturierungsmethode mit Robin-artigen Transmissionsbedingungen, bei der eine Auswertung der Konormalenableitungen aufgrund einer geschickten Aufdatierung der Robin-Daten nicht notwendig ist (insbesondere ist die bekannte Robin-Robin-Methode von Lions als Spezialfall enthalten). Die Konvergenz bezüglich einer partitionierten $H^1$-Norm wird für beide Problemklassen gezeigt. Dabei werden keine über $H^1$ hinausgehende Regularitätsforderungen an die Lösungen gestellt und die Gebiete müssen keine zusätzlichen Glattheitsvoraussetzungen erfüllen. Im letzten Kapitel werden nichtmonotone koerzive quasilineare Probleme untersucht, wobei das Zugrunde liegende Gebiet nur in zwei Lipschitz-Teilgebiete zerlegt sein soll. Das zugehörige nichtlineare Transmissionsproblem wird durch Kirchhoff-Transformation in lineare Teilprobleme mit nichtlinearen Kopplungsbedingungen überführt. Ein optimierungsbasierter Lösungsansatz, welcher einen geeigneten Abstand der rücktransformierten Dirichlet-Daten der linearen Teilprobleme auf den Teilrändern minimiert, führt auf ein optimales Kontrollproblem. Die dabei entstehenden regularisierten freien Minimierungsprobleme werden mit Hilfe eines Gradientenverfahrens unter minimalen Glattheitsforderungen an die Nichtlinearitäten gelöst. Unter zusätzlichen Glattheitsvoraussetzungen an die Nichtlinearitäten und weiteren technischen Voraussetzungen an die Lösung des quasilinearen Ausgangsproblems, kann zudem die quadratische Konvergenz des Newton-Verfahrens gesichert werden.

Method of plant tissue culture and regeneration

Plants may be regenerated from stomatal cells or protoplasts of such cells. Prior to regeneration the cells or protoplasts may be genetically transformed by the introduction of hereditary material most preferably by a DNA construct which is free of genes which specify resistance to antibiotics. The regeneration step may include callus formation on a hormone-free medium. The method is particularly suitable for sugar beet.

Newton`s iterates can converge to non-stationary points

In this note we discuss the convergence of Newton`s method for minimization. We present examples in which the Newton iterates satisfy the Wolfe conditions and the Hessian is positive definite at each step and yet the iterates converge to a non-stationary point. These examples answer a question posed by Fletcher in his 1987 book Practical methods of optimization.

Improving ultimate convergence of an augmented Lagrangian method

Optimization methods that employ the classical Powell-Hestenes-Rockafellar augmented Lagrangian are useful tools for solving nonlinear programming problems. Their reputation decreased in the last 10 years due to the comparative success of interior-point Newtonian algorithms, which are asymptotically faster. In this research, a combination of both approaches is evaluated. The idea is to produce a competitive method, being more robust and efficient than its `pure` counterparts for critical problems. Moreover, an additional hybrid algorithm is defined, in which the interior-point method is replaced by the Newtonian resolution of a Karush-Kuhn-Tucker (KKT) system identified by the augmented Lagrangian algorithm. The software used in this work is freely available through the Tango Project web page:http://www.ime.usp.br/similar to egbirgin/tango/.

Computational study of the step size parameter of the subgradient optimization method

The subgradient optimization method is a simple and flexible linear programming iterative algorithm. It is much simpler than Newton's method and can be applied to a wider variety of problems. It also converges when the objective function is non-differentiable. Since an efficient algorithm will not only produce a good solution but also take less computing time, we always prefer a simpler algorithm with high quality. In this study a series of step size parameters in the subgradient equation is studied. The performance is compared for a general piecewise function and a specific p-median problem. We examine how the quality of solution changes by setting five forms of step size parameter.

Notes on Newton-Euler formulation of robotic manipulators

The equations corresponding to Newton-Euler iterative method for the determination of forces and moments acting on the rigid links of a robotic manipulator are given a new treatment using composed vectors for the representation of both kinematical and dynamical quantities. It is shown that Lagrange equations for the motion of a holonomic system are easily found from the composed vectors defined in this note. Application to a simple model of an industrial robot shows that the method developed in these notes is efficient in solving the dynamics of a robotic manipulator. An example is developed, where it is seen that with the application of appropriate control moments applied to each arm of the robot, starting from a given initial position, it is possible to reach equilibrium in a final pre-assigned position.

Determinação do ponto de máximo carregamento de sistemas de potência utilizando o fluxo de carga desacoplado rápido parametrizado

Os métodos de fluxo de carga por Newton-Raphson e fluxo de carga desacoplado rápido convencionais são considerados inadequados para a obtenção do ponto de máximo carregamento de sistemas de potência, devido à problemas de mal-condicionamento neste ponto crítico e na sua vizinhança. Neste ponto a matriz Jacobiana do método de Newton-Raphson torna-se singular e considera-se que não são mais válidas as hipóteses de desacoplamento P-V e Q-teta utilizadas para a formulação do método fluxo de carga desacoplado rápido. No entanto, mostra-se neste trabalho, que com pequenas modificações, as versões XB e BX do fluxo de carga desacoplado rápido tornam-se adequadas para a obtenção do ponto de máximo carregamento. Estas novas versões modificadas são comparadas entre si com o intuito de explicitar suas características, assim como da influência da atuação dos limites de geração de potência reativa e de tap's de transformadores. Os resultados obtidos para os sistemas testes do IEEE (14, 30, 57 e 118 barras) mostram que as características de convergência das versões originais são preservadas. Além disso, durante o traçado das curvas PV, os diversos métodos podem ser comutados entre si possibilitando o cálculo de todos os pontos da curva com um número reduzido de iterações.

Fluxo de carga continuado utilizando estratégia de parametrização geométrica baseada em parâmetros físicos

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Continuation load flow using automatically determined branch megawatt losses as parameters

The conventional Newton's method is considered to be inadequate for the computation of the maximum loading point (MLP) of power systems since: (i) it encounters difficulties in the vicinity of the MLP: and (ii) the load flow Jacobian matrix becomes singular at the MLP. It is well known that continuation methods are powerful and useful tools that are able to trace the solution PV curve without experiencing such diffculties. However, continuation methods require a parameterisation so that a modified, well conditioned set of load flow equations is obtained. In particular, the Jacobian matrix associated with this modified set of equations should not be singular at the MLP. The authors propose that the actual power losses in transmission branches (lines and transformers) are used to parameterise the approach. Specific procedures for the automatic determination of the most appropriate parameter (branch) are proposed. Such procedures include the utilisation of fast voltage-stability indices. Simulation results are presented to show that the proposed method is able to trace the whole solution PV curve very efficiently.

An efficient geometric parameterization technique for the continuation power flow

Continuation methods have been shown as efficient tools for solving ill-conditioned cases, with close to singular Jacobian matrices, such as the maximum loading point of power systems. Some parameterization techniques have been proposed to avoid matrix singularity and successfully solve those cases. This paper presents a new geometric parameterization scheme that allows the complete tracing of the P-V curves without ill-conditioning problems. The proposed technique associates robustness to simplicity and, it is of easy understanding. The Jacobian matrix singularity is avoided by the addition of a line equation, which passes through a point in the plane determined by the total real power losses and loading factor. These two parameters have clear physical meaning. The application of this new technique to the IEEE systems (14, 30, 57, 118 and 300 buses) shows that the best characteristics of the conventional Newton's method are not only preserved but also improved. (C) 2006 Elsevier B.V. All rights reserved.

Label-free DNA detection based on modified conducting polypyrrole films at microelectrodes

A label-free electrochemical detection method for DNA hybridization based on electrostatic modulation of the ion-exchange kinetics of a polypyrrole film deposited at microelectrodes is reported. Synthetic single-stranded 27-mer oligonucleotides (probe) have been immobilized at 2,5-bis(2-thienyl)-N-(3-phosphorylpropyl)pyrrole film formed by electropolymerization on the previously formed polypyrrole layer. The 27- or 18-mer target oligonucleotides were monitored via the electrochemically driven anion exchange of the inner polypyrrole film. The performance of the miniaturized DNA biosensor system was studied in respect to selectivity, sensitivity, reproducibility, and regeneration of the sensor. Control experiments were performed with a noncomplementary target of 27-mer DNA and 12 base-pair mismatched 18-mer sequences, respectively, and did not show any unspecific binding. Under optimized experimental conditions, the label-free electrochemical biosensor enabled the detection limits of 0.16 and 3.5 fmol for the 18- and 2 7-mer DNA strand, respectively. Furthermore, we demonstrate reusability of the electrochemical DNA biosensor after successful recovery of up to 100% of the original signal by regenerating the DNA label-free electrode with 50 mM HCl at room temperature.

Variational iterative method for scattering problems

A parameter-free variational iterative method is proposed for scattering problems. The present method yields results that are far better, in convergence, stability and precision, than any other momentum space method. Accurate result is obtained for the atomic exponential (Yukawa) potential with an estimated error of less than 1 in 1015 (1010) after some 13 (10) iterations.