Bifurcation theory of nonlinear boundary value problems

The theory of bifurcation of solutions to two-point boundary value problems is developed for a system of nonlinear first order ordinary differential equations in which the bifurcation parameter is allowed to appear nonlinearly. An iteration method is used to establish necessary and sufficient conditions for bifurcation and to construct a unique bifurcated branch in a neighborhood of a bifurcation point which is a simple eigenvalue of the linearized problem. The problem of bifurcation at a degenerate eigenvalue of the linearized problem is reduced to that of solving a system of algebraic equations. Cases with no bifurcation and with multiple bifurcation at a degenerate eigenvalue are considered.

The iteration method employed is shown to generate approximate solutions which contain those obtained by formal perturbation theory. Thus the formal perturbation solutions are rigorously justified. A theory of continuation of a solution branch out of the neighborhood of its bifurcation point is presented. Several generalizations and extensions of the theory to other types of problems, such as systems of partial differential equations, are described.

The theory is applied to the problem of the axisymmetric buckling of thin spherical shells. Results are obtained which confirm recent numerical computations.

Part I. A study of the velocity structure of the earth by the use of core phases. Part II. The 1971 San Fernando earthquake series focal mechanisms and tectonics

The initial objective of Part I was to determine the nature of upper mantle discontinuities, the average velocities through the mantle, and differences between mantle structure under continents and oceans by the use of P'dP', the seismic core phase P'P' (PKPPKP) that reflects at depth d in the mantle. In order to accomplish this, it was found necessary to also investigate core phases themselves and their inferences on core structure. P'dP' at both single stations and at the LASA array in Montana indicates that the following zones are candidates for discontinuities with varying degrees of confidence: 800-950 km, weak; 630-670 km, strongest; 500-600 km, strong but interpretation in doubt; 350-415 km, fair; 280-300 km, strong, varying in depth; 100-200 km, strong, varying in depth, may be the bottom of the low-velocity zone. It is estimated that a single station cannot easily discriminate between asymmetric P'P' and P'dP' for lead times of about 30 sec from the main P'P' phase, but the LASA array reduces this uncertainty range to less than 10 sec. The problems of scatter of P'P' main-phase times, mainly due to asymmetric P'P', incorrect identification of the branch, and lack of the proper velocity structure at the velocity point, are avoided and the analysis shows that one-way travel of P waves through oceanic mantle is delayed by 0.65 to 0.95 sec relative to United States mid-continental mantle.

A new P-wave velocity core model is constructed from observed times, dt/dΔ's, and relative amplitudes of P'; the observed times of SKS, SKKS, and PKiKP; and a new mantle-velocity determination by Jordan and Anderson. The new core model is smooth except for a discontinuity at the inner-core boundary determined to be at a radius of 1215 km. Short-period amplitude data do not require the inner core Q to be significantly lower than that of the outer core. Several lines of evidence show that most, if not all, of the arrivals preceding the DF branch of P' at distances shorter than 143° are due to scattering as proposed by Haddon and not due to spherically symmetric discontinuities just above the inner core as previously believed. Calculation of the travel-time distribution of scattered phases and comparison with published data show that the strongest scattering takes place at or near the core-mantle boundary close to the seismic station.

In Part II, the largest events in the San Fernando earthquake series, initiated by the main shock at 14 00 41.8 GMT on February 9, 1971, were chosen for analysis from the first three months of activity, 87 events in all. The initial rupture location coincides with the lower, northernmost edge of the main north-dipping thrust fault and the aftershock distribution. The best focal mechanism fit to the main shock P-wave first motions constrains the fault plane parameters to: strike, N 67° (± 6°) W; dip, 52° (± 3°) NE; rake, 72° (67°-95°) left lateral. Focal mechanisms of the aftershocks clearly outline a downstep of the western edge of the main thrust fault surface along a northeast-trending flexure. Faulting on this downstep is left-lateral strike-slip and dominates the strain release of the aftershock series, which indicates that the downstep limited the main event rupture on the west. The main thrust fault surface dips at about 35° to the northeast at shallow depths and probably steepens to 50° below a depth of 8 km. This steep dip at depth is a characteristic of other thrust faults in the Transverse Ranges and indicates the presence at depth of laterally-varying vertical forces that are probably due to buckling or overriding that causes some upward redirection of a dominant north-south horizontal compression. Two sets of events exhibit normal dip-slip motion with shallow hypocenters and correlate with areas of ground subsidence deduced from gravity data. Several lines of evidence indicate that a horizontal compressional stress in a north or north-northwest direction was added to the stresses in the aftershock area 12 days after the main shock. After this change, events were contained in bursts along the downstep and sequencing within the bursts provides evidence for an earthquake-triggering phenomenon that propagates with speeds of 5 to 15 km/day. Seismicity before the San Fernando series and the mapped structure of the area suggest that the downstep of the main fault surface is not a localized discontinuity but is part of a zone of weakness extending from Point Dume, near Malibu, to Palmdale on the San Andreas fault. This zone is interpreted as a decoupling boundary between crustal blocks that permits them to deform separately in the prevalent crustal-shortening mode of the Transverse Ranges region.

Earthquake response of steel braces and braced steel frames

This thesis consists of three parts. Chapter 2 deals with the dynamic buckling behavior of steel braces under cyclic axial end displacement. Braces under such a loading condition belong to a class of "acceleration magnifying" structural components, in which a small motion at the loading points can cause large internal acceleration and inertia. This member-level inertia is frequently ignored in current studies of braces and braced structures. This chapter shows that, under certain conditions, the inclusion of the member-level inertia can lead to brace behavior fundamentally different from that predicted by the quasi-static method. This result is to have significance in the correct use of the quasi-static, pseudo-dynamic and static condensation methods in the simulation of braces or braced structures under dynamic loading. The strain magnitude and distribution in the braces are also studied in this chapter.

Chapter 3 examines the effect of column uplift on the earthquake response of braced steel frames and explores the feasibility of flexible column-base anchoring. It is found that fully anchored braced-bay columns can induce extremely large internal forces in the braced-bay members and their connections, thus increasing the risk of failures observed in recent earthquakes. Flexible braced-bay column anchoring can significantly reduce the braced bay member force, but at the same time also introduces large story drift and column uplift. The pounding of an uplifting column with its support can result in very high compressive axial force.

Chapter 4 conducts a comparative study on the effectiveness of a proposed non-buckling bracing system and several conventional bracing systems. The non-buckling bracing system eliminates buckling and thus can be composed of small individual braces distributed widely in a structure to reduce bracing force concentration and increase redundancy. The elimination of buckling results in a significantly more effective bracing system compared with the conventional bracing systems. Among the conventional bracing systems, bracing configurations and end conditions for the bracing members affect the effectiveness.

The studies in Chapter 3 and Chapter 4 also indicate that code-designed conventionally braced steel frames can experience unacceptably severe response under the strong ground motions recorded during the recent Northridge and Kobe earthquakes.

A retrofitting framework for pre-Northridge steel moment-frame buildings

In the 1994 M_w 6.7 Northridge and 1995 M_w 6.9 Kobe earthquakes, steel moment-frame buildings were exposed to an unexpected flaw. The commonly utilized welded unreinforced flange, bolted web connections were observed to experience brittle fractures in a number of buildings, even at low levels of seismic demand. A majority of these buildings have not been retrofitted and may be susceptible to structural collapse in a major earthquake.

This dissertation presents a case study of retrofitting a 20-story pre-Northridge steel moment-frame building. Twelve retrofit schemes are developed that present some range in degree of intervention. Three retrofitting techniques are considered: upgrading the brittle beam-to-column moment resisting connections, and implementing either conventional or buckling-restrained brace elements within the existing moment-frame bays. The retrofit schemes include some that are designed to the basic safety objective of ASCE-41 Seismic Rehabilitation of Existing Buildings.

Detailed finite element models of the base line building and the retrofit schemes are constructed. The models include considerations of brittle beam-to-column moment resisting connection fractures, column splice fractures, column baseplate fractures, accidental contributions from ``simple'' non-moment resisting beam-to-column connections to the lateral force-resisting system, and composite actions of beams with the overlying floor system. In addition, foundation interaction is included through nonlinear translational springs underneath basement columns.

To investigate the effectiveness of the retrofit schemes, the building models are analyzed under ground motions from three large magnitude simulated earthquakes that cause intense shaking in the greater Los Angeles metropolitan area, and under recorded ground motions from actual earthquakes. It is found that retrofit schemes that convert the existing moment-frames into braced-frames by implementing either conventional or buckling-restrained braces are effective in limiting structural damage and mitigating structural collapse. In the three simulated earthquakes, a 20% chance of simulated collapse is realized at PGV of around 0.6 m/s for the base line model, but at PGV of around 1.8 m/s for some of the retrofit schemes. However, conventional braces are observed to deteriorate rapidly. Hence, if a braced-frame that employs conventional braces survives a large earthquake, it is questionable how much service the braces provide in potential aftershocks.

Imperfection Insensitive Thin Shells

The buckling of axially compressed cylindrical shells and externally pressurized spherical shells is extremely sensitive to even very small geometric imperfections. In practice this issue is addressed by either using overly conservative knockdown factors, while keeping perfect axial or spherical symmetry, or adding closely and equally spaced stiffeners on shell surface. The influence of imperfection-sensitivity is mitigated, but the shells designed from these approaches are either too heavy or very expensive and are still sensitive to imperfections. Despite their drawbacks, these approaches have been used for more than half a century.

This thesis proposes a novel method to design imperfection-insensitive cylindrical shells subject to axial compression. Instead of following the classical paths, focused on axially symmetric or high-order rotationally symmetric cross-sections, the method in this thesis adopts optimal symmetry-breaking wavy cross-sections (wavy shells). The avoidance of imperfection sensitivity is achieved by searching with an evolutionary algorithm for smooth cross-sectional shapes that maximize the minimum among the buckling loads of geometrically perfect and imperfect wavy shells. It is found that the shells designed through this approach can achieve higher critical stresses and knockdown factors than any previously known monocoque cylindrical shells. It is also found that these shells have superior mass efficiency to almost all previously reported stiffened shells.

Experimental studies on a design of composite wavy shell obtained through the proposed method are presented in this thesis. A method of making composite wavy shells and a photogrametry technique of measuring full-field geometric imperfections have been developed. Numerical predictions based on the measured geometric imperfections match remarkably well with the experiments. Experimental results confirm that the wavy shells are not sensitive to imperfections and can carry axial compression with superior mass efficiency.

An efficient computational method for the buckling analysis of corrugated and stiffened cylindrical shells subject to axial compression has been developed in this thesis. This method modifies the traditional Bloch wave method based on the stiffness matrix method of rotationally periodic structures. A highly efficient algorithm has been developed to implement the modified Bloch wave method. This method is applied in buckling analyses of a series of corrugated composite cylindrical shells and a large-scale orthogonally stiffened aluminum cylindrical shell. Numerical examples show that the modified Bloch wave method can achieve very high accuracy and require much less computational time than linear and nonlinear analyses of detailed full finite element models.

This thesis presents parametric studies on a series of externally pressurized pseudo-spherical shells, i.e., polyhedral shells, including icosahedron, geodesic shells, and triambic icosahedra. Several optimization methods have been developed to further improve the performance of pseudo-spherical shells under external pressure. It has been shown that the buckling pressures of the shell designs obtained from the optimizations are much higher than the spherical shells and not sensitive to imperfections.

Pulsed neutron measurements in two adjacent finite media

The pulsed neutron technique has been used to investigate the decay of thermal neutrons in two adjacent water-borated water finite media. Experiments were performed with a 6x6x6 inches cubic assembly divided in two halves by a thin membrane and filled with pure distilled water on one side and borated water on the other side.

The fundamental decay constant was measured versus the boric acid concentration in the poisoned medium. The experimental results showed good agreement with the predictions of the time dependent diffusion model. It was assumed that the addition of boric acid increases the absorption cross section of the poisoned medium without affecting its diffusion properties: In these conditions, space-energy separability and the concept of an “effective” buckling as derived from diffusion theory were introduced. Their validity was supported by the experimental results.

Measurements were performed with the absorption cross section of the poisoned medium increasing gradually up to 16 times its initial value. Extensive use of the IBM 7090-7094 Computing facility was made to analyze properly the decay data (Frantic Code). Attention was given to the count loss correction scheme and the handling of the statistics involved. Fitting of the experimental results into the analytical form predicted by the diffusion model led to

Ʃ_av = 4721 sec^-1 (±150)

D_o = 35972 cm²sec^-1 (±800) for water at 21˚C

C (given) = 3420 cm⁴sec^-1

These values, when compared with published data, show that the diffusion model is adequate in describing the experiment.

Strength of thin-walled elliptical cylinders supported at the minor axis

In this investigation it was found that the instability failure of curved sheet is nearly independent of the type of loading and is primarily a function of the maximum stress, radius-thickness ration and modulus of elasticity. A method of correlating the critical stress of thin sheet under several different types of loading is given. An explanation for the experimental critical stress of thin walled cylinders under bending being greater than that for pure compression is given. The strength of unstiffened thin walled circular nose sections under pure bending was found to be controlled by local instability of the section, rather than a large scale instability. The equation of local instability of curved sheet gives values which are in fair agreement with those found experimentally.

The strength of elliptical cylinders supported at the minor axis under bending plus shear loads is governed primarily by the bending strength, and is little effected by the sheer force unless the amount of shear is quite large with respect to the moment. The effect of increasing the amount of elliptically greatly reduces the bending and shear strength of nose sections. Under torsional loads the stress at buckling falls off as the ration of the major to minor axis increases but the failure stress decreases at a slower rate than the buckling stress. The length effect of semi-circular sections under torsion is similar to that of a circular tube, and can be obtained by Donnell's theoretical equation.

Avaliação numérica de estabilidade lateral de vigas casteladas

Restrições de espaço e altura são frequentemente impostas às edificações residenciais, comerciais, industriais, depósitos e galpões com um ou diversos pavimentos em função de aspectos de regulamentos regionais, técnicos, econômicos ou ainda de natureza estética. A fim de proporcionar a passagem de tubulações e dutos de grande diâmetro sob vigas de aço, grandes alturas são normalmente requeridas, demandando por vezes, magnitudes de altura inviáveis entre pavimentos de edificações. Diversas soluções estruturais podem ser utilizadas para equacionar tais obstáculos, onde dentre outras, pode-se citar as vigas com inércia variável, stub-girders, treliças mistas, vigas misuladas e vigas com uma ou múltiplas aberturas na alma com geometrias variadas. No que tange às vigas casteladas, solução estrutural pautada neste estudo, a estabilidade é sempre um motivo de preocupação tipicamente durante a construção quando os contraventamentos laterais ainda não estão instalados. De qualquer forma, o comprimento destravado em geral alcançado pelos vãos destas vigas, são longos o suficiente para que a instabilidade ocorra. Todavia, o acréscimo substancial da resistência à flexão de tais membros devido ao aumento da altura oriundo de seu processo fabril em relação ao perfil matriz, aliada a economia de material e utilidade fim de serviço, garante a atratividade no aproveitamento destas, para grandes vãos junto aos projetistas. Não obstante, este aumento proporcional no comprimento dos vãos faz com que a instabilidade lateral ganhe importância especial. Neste contexto, o presente trabalho tem por objetivo desenvolver um modelo numérico que permita a realização de uma avaliação paramétrica a partir da calibração do modelo com resultados experimentais, efetuar a análise do comportamento de vigas casteladas e verificar seus mecanismos de falha, considerando comportamento elasto-plástico, além das não-linearidades geométricas. Também é objetivo deste trabalho, avaliar, quantificar e determinar a influência das diferenças geométricas características das vigas casteladas em relação às vigas maciças com as mesmas dimensões, analisando e descrevendo o comportamento estrutural destas vigas de aço para diversos comprimentos de vãos. A metodologia empregada para tal estudo baseou-se em uma análise paramétrica com o auxílio do método numérico dos elementos finitos.

Diseño y fabricación de un instrumento coaxial para cirugía laparoscópica.

[ES]En este trabajo se ha diseñado un producto secundario con una fijación extra en formato trocar para operaciones percutáneas dirigidas por sistemas de imagen. Este instrumento que hasta ahora era inexistente en el mercado actuará como túnel, permitiendo la introducción de agujas de radiofrecuencia y biopsia por su interior. La particularidad de este producto es el diseño del sistema de sujeción, que garantizará al cirujano que el punto a intervenir sea el mismo durante toda la operación. Este diseño se ha realizado después de hacer un estudio en profundidad de las alternativas existentes de los sistemas de fijación en el campo de la medicina oncológica y el sector industrial. Además se detalla el material de cada elemento del producto en base a la norma sanitaria correspondiente, también se explica el proceso más adecuado para la fabricación de dichos elementos y por último se realiza un cálculo sencillo del pandeo de la aguja, que garantice un uso óptimo del instrumento.

Comportamento estrutural de vigas casteladas em aço inoxidável.

As vigas casteladas já são utilizadas em diversos tipos de estruturas para se vencer grandes vãos. Há uma nova realidade para o uso de perfis castelados, agora em aço inoxidável. Este trabalho apresenta uma análise comparativa entre vigas casteladas e vigas de alma cheia em perfis IPE, baseadas na norma europeia, no método da resistência contínua; e, também, através de um modelo em elementos finitos desenvolvido no programa Ansys. Este trabalho tem como objetivo verificar o comportamento e a resistência à flexão das vigas casteladas em aço inoxidável. As não-linearidades físicas e geométricas foram incorporadas aos modelos, a fim de se mobilizar totalmente a capacidade resistente desta estrutura. A não-linearidade do material foi considerada através do critério de plastificação de von Mises e da lei constitutiva tensão versus deformação multi-linear, de forma a exibir um comportamento elasto-plástico com encruamento. A não-linearidade geométrica foi considerada a partir de uma formulação de Lagrange atualizado. Os resultados numéricos das vigas em estudo foram avaliados quanto aos modos de falha e da distribuição de tensões. Os momentos resistentes analíticos foram comparados com os resultados do modelo numérico apresentando valores satisfatórios e a favor da segurança.

Observation of large-scale features on graphite by scanning tunnelling microscopy

Superlattice structures and rippling fringes were imaged on two separate pieces of graphite (HOPG) by scanning tunnelling microscopy (STM). We observed the corrugation conservation phenomenon on one of the superlattice structures where an overlayer does not attenuate the corrugation amplitude of the superlattice. Such a phenomenon may illustrate an implication that nanoscale defects a few layers underneath the surface may propagate through many layers without decay and form the superlattice structure on the topmost surface. Some rippling fringes with periodicities of 20 nm and 30 nm and corrugations of 0.1 nm and 0.15nm were observed in the superlattice area and in nearby regions. Such fringes are believed to be due to physical buckling of the surface. The stress required to generate such structures is estimated, and a possible cause is discussed. An equation relating the attenuation factor to the number of overlayers is proposed. © 2005 The Japan Society of Applied Physics.

The use of object oriented modelling and simulation tools in constraint satisfaction for embodiment design

The software package Dymola, which implements the new, vendor-independent standard modelling language Modelica, exemplifies the emerging generation of object-oriented modelling and simulation tools. This paper shows how, in addition to its simulation capabilities, it may be used as an embodiment design tool, to size automatically a design assembled from a library of generic parametric components. The example used is a miniature model aircraft diesel engine. To this end, the component classes contain extra algebraic equations calculating the overload factor (or its reciprocal, the safety factor) for all the different modes of failure, such as buckling or tensile yield. Thus the simulation results contain the maximum overload or minimum safety factor for each failure mode along with the critical instant and the device state at which it occurs. The Dymola "Initial Conditions Calculation" function, controlled by a simple software script, may then be used to perform automatic component sizing. Each component is minimised in mass, subject to a chosen safety factor against failure, over a given operating cycle. Whilst the example is in the realm of mechanical design, it must be emphasised that the approach is equally applicable to the electrical or mechatronic domains, indeed to any design problem requiring numerical constraint satisfaction.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Pipeline uplift response in thawed sandy backfills: Preliminary findings

Offshore and onshore buried pipelines under high operating temperature and pressures may lead to upheaval buckling (UHB) if sufficient soil cover is not present to prevent the upward movement of the pipeline. In regions where seasonal changes involve ground soil undergoing freezing-thawing cycles, the uplift resistance from soil cover may be minimum when the soil is undergoing thawing. This paper presents the results from 2 directly-comparable minidrum centrifuge tests conducted at the Schofield Centre, University of Cambridge, to investigate the difference in uplift resistance responses between fully-saturated and thawed sandy backfill conditions. Both tests were conducted drained at 30g using an 8.6 mm diameter aluminium model pipe, corresponding to a prototype pipe diameter of 258 mm. The soil cover/pipe diameter ratio, H/D, was kept at 1. Fraction E fine silica sand was used as the backfill. Preliminary experimental results indicated that the ultimate uplift resistance of a thawing sand backfill to be lower than that of a fully saturated sand backfill. This suggests that in regions where backfill soil undergoes freeze-thaw cycles, the thawing backfill may be more critical than fully saturated backfill for uplift resistance. The 2-dimensional displacement field during the experiment was accurately measured and analysed using the Particle Image Velocimetry technique. Copyright © 2011 by the International Society of Offshore and Polar Engineers (ISOPE).