Estudamos a Eletrodinâmica de Wheeler-Feynman usando um princípio variacional para um funcional de ação finito acoplado a um problema de valor na fronteira. Para trajetórias C2 por trechos, a condição de ponto crítico desse funcional fornece as equações de movimento de Wheeler-Feynman mais uma condição de continuidade dos momentos parciais e energias parciais, conhecida como condição de quina de Weierstrass-Erdmann. Estudamos em detalhe um sub-caso mais simples, onde os dados de fronteira têm um comprimento mínimo. Nesse caso, mostramos que a condição de extremo se reduz a um problema de valor na chegada para uma equação diferencial com retardo misto dependente do estado e do tipo neutro. Resolvemos numericamente esse problema usando um método de shooting e um método de Runge-Kutta de quarta ordem. Para os casos em que as fronteiras mínimas têm velocidades descontínuas, elaboramos uma técnica para resolver as condições de quina de Weierstrass-Erdmann junto com o problema de valor na chegada. As trajetórias com velocidades descontínuas previstas pelo método variacional foram verificadas por experimentos numéricos. Em um segundo desenvolvimento, para o caso mais difícil de fronteiras de comprimento arbitrário, implementamos um método de minimização com gradiente fraco para o princípio variacional e problema de fronteira acima citado. Elaboramos dois métodos numéricos, ambos implementados em MATLAB, para encontrar soluções do problema eletromagnético de dois corpos. O primeiro combina o método de elementos finitos com o método de Newton para encontrar as soluções que anulam o gradiente fraco do funcional para fronteiras genéricas. O segundo usa o método do declive máximo para encontrar as soluções que minimizam a ação. Nesses dois métodos as trajetórias são aproximadas dentro de um espaço de dimensão finita gerado por uma Galerkiana que suporta velocidades descontínuas. Foram realizados diversos testes e experimentos numéricos para verificar a convergência das trajetórias calculada numericamente; também comparamos os valores do funcional calculados numericamente com alguns resultados analíticos sobre órbitas circulares. / We study the Wheeler-Feynman electrodynamics using a variational principle for an action functional coupled to a finite boundary value problem. For piecewise C2 trajectories, the critical point condition for this functional gives the Wheeler-Feynman equations of motion in addition to a continuity condition of partial moments and partial energies, known as the Weierstrass-Erdmann corner conditions. In the simplest case, for the boundary value problem of shortest length, we show that the critical point condition reduces to a two-point boundary value problem for a state-dependent mixed-type neutral differential-delay equation. We solve this special problem numerically using a shooting method and a fourth order Runge-Kutta. For the cases where the boundary segment has discontinuous velocities we developed a technique to solve the Weierstrass-Erdmann corner conditions and the two-point boundary value problem together. The trajectories with discontinuous velocities presupposed by the variational method were verified by numerical experiments. In a second development, for the harder case with boundaries of arbitrary length, we implemented a method of minimization with weak gradient for the variational principle quoted above. Two numerical methods were implemented in MATLAB to find solutions of the two-body electromagnetic problem. The first combines the finite element method with Newtons method to find the solutions that vanish the weak gradient. The second uses the method of steepest descent to find the solutions that minimize the action. In both methods the trajectories are approximated within a finite-dimensional space generated by a Galerkian that supports discontinuous velocities. Many tests and numerical experiments were performed to verify the convergence of the numerically calculated trajectories; also were compared the values of the functional computed numerically with some known analytical results on circular orbits.

Wave propagation is a phenomenon resulting from high transient loadings where the duration of the load is in µ seconds range. In aerospace and space craft industries it is important to gain knowledge about the high frequency characteristics as it aids in structural health monitoring, wave transmission/attenuation for vibration and noise level reduction. The wave propagation problem can be approached by the conventional Finite Element Method(FEM); but at higher frequencies, the wavelengths being small, the size of the finite element is reduced to capture the response behavior accurately and thus increasing the number of equations to be solved, leading to high computational costs. On the other hand such problems are handled in the frequency domain using Fourier transforms and one such method is the Spectral Finite Element Method(SFEM). This method is introduced first by Doyle ,for isotropic case and later popularized in developing specific purpose elements for structural diagnostics for inhomogeneous materials, by Gopalakrishnan. The general approach in this method is that the partial differential wave equations are reduced to a set of ordinary differential equations(ODEs) by transforming these equations to another space(transformed domain, say Fourier domain). The reduced ODEs are usually solved exactly, the solution of which gives the dynamic shape functions. The interpolating functions used here are exact solution of the governing differential equations and hence, the exact elemental dynamic stiffness matrix is derived. Thus, in the absence of any discontinuities, one element is sufficient to model 1-D waveguide of any length. This elemental stiffness matrix can be assembled to obtain the global matrix as in FEM, but in the transformed space. Thus after obtaining the solution, the original domain responses are obtained using the inverse transform. Both the above mentioned manuscripts present the Fourier transform based spectral finite element (FSFE), which has the inherent aliasing problem that is persistent in the application of the Fourier series/Fourier transforms. This is alleviated by using an additional throw-off element and/or introducing slight damping in to the system. More recently wave let transform based spectral finite element(WSFE) has been formulated which alleviated the aliasing problem; but has a limitation in obtaining the frequency characteristics, like the group speeds are accurate only up-to certain fraction of the Nyquist(central frequency). Currently in this thesis Laplace transform based spectral finite elements(LSFE) are developed for sandwich members. The advantages and limitations of the use of different transforms in the spectral finite element framework is presented in detail in Chapter-1. Sandwich structures are used in the space craft industry due to higher stiffness to weight ratio. Many issues considered in the design and analysis of sandwich structures are discussed in the well known books(by Zenkert, Beitzer). Typically the main load bearing structures are modeled as beam sand plates. Plate structures with kh<1 is analysed based on the Kirch off plate theory/Classical Plate Theory(CPT) and when the bending wavelength is small compared to the plate thickness, the effect of shear deformation and rotary inertia needs to be included where, k is the wave number and h is the thickness of the plate. Many works regarding the wave propagation in sandwich structures has been published in the past literature for wave propagation in infinite sandwich structure and giving the complete description of dispersion relation with no restriction on frequency and wavelength. More recently exact analytical solution or simply supported sandwich plate has been derived. Also it is seen by comparison of dispersion curves obtained with exact (3D formulation of theory of elasticity) and simplified theories (2D formulation as generalization of Timoshenko theory) made on infinite domain and concluded that the simplified theory can be reliably used to assess the waveguide properties of sandwich plate in the frequency range of interest. In order to approach the problems with finite domain and their implementation in the use of general purpose code; finite degrees of freedom is enforced. The concept of displacement based theories provides the flexibility in assuming different kinematic deformations to approach these problems. Many of the displacement based theories incorporate the Equivalent Single Layer(ESL) approach and these can capture the global behavior with relative ease. Chapter-2 presents the Laplace spectral finite element for thick beams based on the First order Shear Deformation Theory (FSDT). Here the effect of different choices of the real part of the Laplace variable is demonstrated. It is shown that the real part of the Laplace variable acts as a numerical damping factor. The spectrum and dispersion relations are obtained and the use of these relations are demonstrated by an example. Here, for sandwich members based on FSDT, an appropriate choice of the correction factor ,which arises due to the inconsistency between the kinematic hypothesis and the desired accuracy is presented. Finally the response obtained by the use of the element is validated with experimental results. For high shock loading cases, the core flexibility induces local effects which are very predominant and this can lead to debonding of face sheets. The ESL theories mentioned above cannot capture these effects due to the computation of equivalent through the thickness section properties. Thus, higher order theories such as the layer-wise theories are required to capture the local behaviour. One such theory for sandwich panels is the Higher order Sandwich Plate theory (HSaPT). Here, the in-plane stress in the core has been neglected; but gives a good approximation for sandwich construction with soft cores. Including the axial inertial terms of the core will not yield constant shear stress distribution through the height of the core and hence more recently the Extended Higher order Sandwich Plate theory (EHSaPT) is proposed. The LSFE based on this theory has been formulated and is presented in Chapter-4. Detailed 3D orthotropic properties of typical sandwich construction is considered and the core compressibility effect of local behavior due to high shock loading is clearly brought out. As detailed local behavior is sought the degrees of freedom per element is high and the specific need for such theory as compared with the ESL theories is discussed. Chapter-4 presents the spectral finite element for plates based on FSDT. Here, multi-transform method is used to solve the partial differential equations of the plate. The effect of shear deformation is brought out in the spectrum and dispersion relations plots. Response results obtained by the formulated element is compared and validated with many different experimental results. Generally structures are built-up by connecting many different sub-structures. These connecting members, called joints play a very important role in the wave transmission/attenuation. Usually these joints are modeled as rigid joints; but in reality these are flexible and exhibits non-linear characteristics and offer high damping to the energy flow in the connected structures. Chapter-5 presents the attenuation and transmission of wave energy using the power flow approach for rigid joints for different configurations. Later, flexible spectral joint model is developed and the transmission/attenuation across the flexible joints is studied. The thesis ends with conclusion and highlighting futures cope based on the developments reported in this thesis.

Problems governed by partial differential equations (PDEs) in deformable domains, t Rd; d = 2; 3; are of fundamental importance in science and engineering. They are of particular relevance in the design of many engineering systems e.g., aircrafts and bridges as well as to the analysis of several biological phenomena e.g., blood ow in arteries. However, developing numerical scheme for such problems is still very challenging even when the deformation of the boundary of domain is prescribed a priori. Possibility of excessive mesh distortion is one of the major challenge when solving such problems with numerical methods using boundary tted meshes. The arbitrary Lagrangian- Eulerian (ALE) approach is a way to overcome this difficulty. Numerical simulations of convection-dominated problems have for long been the subject to many researchers. Galerkin formulations, which yield the best approximations for differential equations with high diffusivity, tend to induce spurious oscillations in the numerical solution of convection dominated equations. Though such spurious oscillations can be avoided by adaptive meshing, which is computationally very expensive on ne grids. Alternatively, stabilization methods can be used to suppress the spurious oscillations. In this work, the considered equation is designed within the framework of ALE formulation. In the first part, Streamline Upwind Petrov-Galerkin (SUPG) finite element method with conservative ALE formulation is proposed. Further, the first order backward Euler and the second order Crank-Nicolson methods are used for the temporal discretization. It is shown that the stability of the semi-discrete (continuous in time) ALE-SUPG equation is independent of the mesh velocity, whereas the stability of the fully discrete problem is unconditionally stable for implicit Euler method and is only conditionally stable for Crank-Nicolson time discretization. Numerical results are presented to support the stability estimates and to show the influence of the SUPG stabilization parameter in a time-dependent domain. In the second part of this work, SUPG stabilization method with non-conservative ALE formulation is proposed. The implicit Euler, Crank-Nicolson and backward difference methods are used for the temporal discretization. At the discrete level in time, the ALE map influences the stability of the corresponding discrete scheme with different time discretizations, and it leads to schemes where conservative and non-conservative formulations are no longer equivalent. The stability of the fully discrete scheme, irrespective of the temporal discretization, is only conditionally stable. It is observed from numerical results that the Crank-Nicolson scheme induces high oscillations in the numerical solution compare to the implicit Euler and the backward difference time discretiza-tions. Moreover, the backward difference scheme is more sensitive to the stabilization parameter k than the other time discretizations. Further, the difference between the solutions obtained with the conservative and non-conservative ALE forms is significant when the deformation of domain is large, whereas it is negligible in domains with small deformation. Finally, the local projection stabilization (LPS) and the higher order dG time stepping scheme are studied for convection dominated problems. The analysis is based on the quadrature formula for approximating the integrals in time. We considered the exact integration in time, which is impractical to implement and the Radau quadrature in time, which can be used in practice. The stability and error estimates are shown for the mathematical basis of considered numerical scheme with both time integration methods. The numerical analysis reveals that the proposed stabilized scheme with exact integration in time is unconditionally stable, whereas Radau quadrature in time is conditionally stable with time-step restriction depending on the ALE map. The theoretical estimates are illustrated with appropriate numerical examples with distinct features. The second order dG(1) time discretization is unconditionally stable while Crank-Nicolson gives the conditional stable estimates only. The convergence order for dG(1) is two which supports the error estimate.

A mecânica celular jaz nas propriedades materiais da membrana plasmática, fundamentalmente uma bicamada fosfolipídica com espessura de dimensões moleculares. Além de forças elásticas, tal material bidimensional também experimenta tensões viscosas devido ao seu comportamento fluido (presumivelmente newtoniano) na direção tangencial. A despeito da notável concordância entre teoria e experimentos biofísicos sobre a geometria de membranas celulares, ainda não se faz presente um método computacional para simulação de sua (real) dinâmica viscosa governada pela lei de Boussinesq-Scriven. Assim sendo, introduzimos uma formulação variacional mista de três campos para escoamentos viscosos de superfícies fechadas curvas. Nela, o fluido circundante é levado em conta considerando-se uma restrição de volume interior, ao passo que uma restrição de área corresponde à inextensibilidade. As incógnitas são a velocidade, o vetor curvatura e a pressão superficial, todas estas interpoladas com elementos finitos lineares contínuos via estabilização baseada na projeção do gradiente de pressão. O método é semi-implícito e requer a solução de apenas um único sistema linear por passo de tempo. Outro ingrediente numérico proposto é uma força que mimetiza a ação de uma pinça óptica, permitindo interação virtual com a membrana, onde a qualidade e o refinamento de malha são mantidos por remalhagem adaptativa automática. Extensivos experimentos numéricos de dinâmica de relaxação são apresentados e comparados com soluções quasi-analíticas. É observada estabilidade temporal condicional com uma restrição de passo de tempo que escala como o quadrado do tamanho de malha. Reportamos a convergência e os limites de estabilidade de nosso método e sua habilidade em predizer corretamente o equilíbrio dinâmico de compridas e finas elongações cilíndricas (tethers) que surgem a partir de pinçamentos membranais. A dependência de forma membranal com respeito a uma velocidade imposta de pinçamento também é discutida, sendo que há um valor limiar de velocidade abaixo do qual um tether não se forma de início. Testes adicionais ilustram a robustez do método e a relevância dos efeitos viscosos membranais quando sob a ação de forças externas. Sem dúvida, ainda há um longo caminho a ser trilhado para o entendimento completo da mecânica celular (há de serem consideradas outras estruturas tais como citoesqueleto, canais iônicos, proteínas transmembranares, etc). O primeiro passo, porém, foi dado: a construção de um esquema numérico variacional capaz de simular a intrigante dinâmica das membranas celulares. / Cell mechanics lies on the material properties of the plasmatic membrane, fundamentally a two-molecule-thick phospholipid bilayer. Other than bending elastic forces, such a two-dimensional interfacial material also experiences viscous stresses due to its (presumably Newtonian) surface fluid tangential behaviour. Despite the remarkable agreement on membrane shapes between theory and biophysical experiments, there is no computational method for simulating its (actual) viscous dynamics governed by the Boussinesq- Scriven law. Accordingly, we introduce a mixed three-field variational formulation for viscous flows of closed curved surfaces. In it, the bulk fluid is taken into account by means of an enclosed-volume constraint, whereas an area constraint stands for the membranes inextensible character. The unknowns are the velocity, vector curvature and surface pressure fields, all of which are interpolated with linear continuous finite elements by means of a pressure-gradient-projection scheme. The method is semi-implicit and it requires the solution of a single linear system per time step. Another proposed ingredient is a numerical force that emulates the action of an optical tweezer, allowing for virtual interaction with the membrane, where mesh quality and refinement are maintained by adaptively remeshing. Extensive relaxation experiments are reported and compared with quasi-analytical solutions. Conditional time stability is observed, with a time step restriction that scales as the square of the mesh size. We discuss both convergence and the stability limits of our method, its ability to correctly predict the dynamical equilibrium of the tether due to tweezing. The dependence of the membrane shape on imposed tweezing velocities is also addressed. Basically, there is a threshold velocity value below which the tethers shape does not arise at first. Further tests illustrate the robustness of the method and show the significance of viscous effects on membranes deformation under external forces. Undoubtedly, there is still a long way to track toward the understanding of celullar mechanics (one still has to account other structures such as cytoskeleton, ion channels, transmembrane proteins, etc). The first step has given, though: the design of a numerically robust variational scheme capable of simulating the engrossing dynamics of fluid cell membranes.

This work deals with issue of human wrist and appropriate total wrist implant allowing a restoration of hand mobility approaching physiological condition after traumatic and degenerative diseases. Treating these diseases are very complex. These issues including a biological and medical issues. To determine the appropriate treatment method and select right total wrist implant is important to know the behavior the human wrist at all stages in terms of medical and biomechanical. For this reason, it was developed a biomechanical study including computation model of human wrist allowing solution of strain and stress of hand in physiological and pathological conditions and condition after total wrist implant. The frost remodeling of bone tissue was used for analysis of human wrist bone tissues and bone tissues after application of total wrist implant RE-MOTION™ Total Wrist.

Medical devices have made a big step forward in the past decades. One of the most noticeable medical events of the twenties century was the development of long-lasting, wireless electronic implants such as identification tags, pacemakers and neuronal stimulators. These devices were only made possible after the development of small scale radio frequency electronics. Small radio electronic circuits provided a way to operate in both transmission and reception mode allowing an implant to communicate with an external world from inside a living organism. Bidirectional communication is a vital feature that has been increasingly implemented in similar systems to continuously record biological parameters, to remotely configure the implant, or to wirelessly stimulate internal organs. Further miniaturisation of implantable devices to make the operation of the device more comfortable for the patient requires rethinking of the whole radio system concept making it both power efficient and of high performance. Nowadays, high data throughput, large bandwidth, and long term operation requires new radio systems to operate at UHF (ultra-high frequency) bands as this is the most suitable for implantable applications. For instance, the MICS (Medical Implant Communication System) band was introduced for the communication with implantable devices. However, this band could only enable communication at low data rates. This was acceptable for the transmission of telemetry data such as heart beat rate, respiratory and temperature with sub Mbps rates. Novel developments such as neuronal and prosthetic implants require significantly higher data rates more than 10 Mbps that can be achieved with large bandwidth communicating systems operating at higher frequencies in a GHz range. Higher operating frequency would also resolve a strong issue of MICS devices, namely the scale of implants defined by dimensions of antennas used at this band. Operation at 2.4 GHz ISM band was recognized to be the most adequate as it has a moderate absorption in the human body providing a compromise between an antenna/implant scale and a total power efficiency of the communicating system. This thesis addresses a key challenge of implantable radio communicating systems namely an efficient and small scale antenna design which allows a high yield fabrication in a microelectronic fashion. It was demonstrated that a helical antenna design allows the designer to precisely tune the operating frequency, input impedance, and bandwidth by changing the geometry of a self-assembled 3D structure defined by an initial 2D planar layout. Novel stimuli responsive materials were synthesized, and the rolled-up technology was explored for fabrication of 5.5-mm-long helical antenna arrays operating in ISM bands at 5.8 and 2.4 GHz. Characterization and various applications of the fabricated antennas are successfully demonstrated in the thesis.

The objective of this thesis was to identify the mechanisms associated with the recurrence of rock-slope instability along the R518 and R523 roads in Limpopo. Advanced geotechnical methods and ASTER imagery were used for the purpose while a predictive rockfall hazard rating matrix chart and rock slope stability charts for unsaturated sensitive clay soil and rock slopes were to be developed. The influence of extreme rainfall on the slope stability of the sensitive clay soil was also evaluated. To achieve the above, field observations, geological mapping, kinematic analysis, and limit equilibrium were performed. The latter involved toppling, transitional and rotational analyses. Numerical simulation was finally resorted to. The following software packages were employed: SWEDGE, SLIDE, RocData, RocFall, DIPS, RocPlane, and Phase 2. The simulation outputs were analyzed in conjunction with ASTER images. The advanced remote sensing data paved the way for landslide susceptibility analysis. From all the above, rockfall hazard prediction charts and slope stability prediction charts were developed. Several factors were also shown by numerical simulation to influence slope instability in the area of study, i.e. sites along the R518 and R523 roads in the Thulamela Municipality. The most important factors are extreme rainfall, steep slopes, geological features and water streams in the region, and improper road construction. Owing to the complexity of the failure mechanisms in the study area, it was concluded that both slope stability prediction charts and rock hazard matrix charts are very useful. They indeed enable one to characterize slope instability in sensitive clay soils as well as rockfall hazards in the study area. It is however recommended that future work is undertaken to explore the use of sophisticated and scientific methods. This is instrumental in the development of predictive tools for rock deformation and displacement in landslide events. / Electrical and Mining Engineering / D. Phil. (Mining Engineering)

Infimalfolgen nichtkonvexer Variationsprobleme haben aufgrund feiner Oszillationen häufig keinen starken Grenzwert in Sobolevräumen. Diese Oszillationen haben eine physikalische Bedeutung; Finite-Element-Approximationen können sie jedoch im Allgemeinen nicht auflösen. Relaxationsmethoden ersetzen die nichtkonvexe Energie durch ihre (semi)konvexe Hülle. Das entstehende makroskopische Modell ist degeneriert: es ist nicht strikt konvex und hat eventuell mehrere Minimalstellen. Die fehlende Kontrolle der primalen Variablen führt zu Schwierigkeiten bei der a priori und a posteriori Fehlerschätzung, wie der Zuverlässigkeits- Effizienz-Lücke und fehlender starker Konvergenz. Zur Überwindung dieser Schwierigkeiten erweitern Stabilisierungstechniken die relaxierte Energie um einen diskreten, positiv definiten Term. Bartels et al. (IFB, 2004) wenden Stabilisierung auf zweidimensionale Probleme an und beweisen dabei starke Konvergenz der Gradienten. Dieses Ergebnis ist auf glatte Lösungen und quasi-uniforme Netze beschränkt, was adaptive Netzverfeinerungen ausschließt. Die vorliegende Arbeit behandelt einen modifizierten Stabilisierungsterm und beweist auf unstrukturierten Netzen sowohl Konvergenz der Spannungstensoren, als auch starke Konvergenz der Gradienten für glatte Lösungen. Ferner wird der sogenannte Fluss-Fehlerschätzer hergeleitet und dessen Zuverlässigkeit und Effizienz gezeigt. Für Interface-Probleme mit stückweise glatter Lösung wird eine Verfeinerung des Fehlerschätzers entwickelt, die den Fehler der primalen Variablen und ihres Gradienten beschränkt und so starke Konvergenz der Gradienten sichert. Der verfeinerte Fehlerschätzer konvergiert schneller als der Fluss- Fehlerschätzer, und verringert so die Zuverlässigkeits-Effizienz-Lücke. Numerische Experimente mit fünf Benchmark-Tests der Mikrostruktursimulation und Topologieoptimierung ergänzen und bestätigen die theoretischen Ergebnisse. / Infimising sequences of nonconvex variational problems often do not converge strongly in Sobolev spaces due to fine oscillations. These oscillations are physically meaningful; finite element approximations, however, fail to resolve them in general. Relaxation methods replace the nonconvex energy with its (semi)convex hull. This leads to a macroscopic model which is degenerate in the sense that it is not strictly convex and possibly admits multiple minimisers. The lack of control on the primal variable leads to difficulties in the a priori and a posteriori finite element error analysis, such as the reliability-efficiency gap and no strong convergence. To overcome these difficulties, stabilisation techniques add a discrete positive definite term to the relaxed energy. Bartels et al. (IFB, 2004) apply stabilisation to two-dimensional problems and thereby prove strong convergence of gradients. This result is restricted to smooth solutions and quasi-uniform meshes, which prohibit adaptive mesh refinements. This thesis concerns a modified stabilisation term and proves convergence of the stress and, for smooth solutions, strong convergence of gradients, even on unstructured meshes. Furthermore, the thesis derives the so-called flux error estimator and proves its reliability and efficiency. For interface problems with piecewise smooth solutions, a refined version of this error estimator is developed, which provides control of the error of the primal variable and its gradient and thus yields strong convergence of gradients. The refined error estimator converges faster than the flux error estimator and therefore narrows the reliability-efficiency gap. Numerical experiments with five benchmark examples from computational microstructure and topology optimisation complement and confirm the theoretical results.