A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Hill Instability Analysis Of Tlp Tether Subjected To Combined Platform Surge And Heave Motions

This paper presents the Hill instability analysis of Tension Leg Platform (TLP) tether it, deep sea. The 2-D nonlinear beam model which is Undergoing Coupled axial and transverse vibrations, is applied. The governing equations are reduced to nonlinear Hill equation by use of the Galerkin's method and the modes superposition principle. The Hill instability charted Lip to large parameters is obtained. An important parameter M is defined and can he expressed as the functions of tether length, the platform surge and heave motion amplitudes. Some example studies are performed for various environmental conditions. The results demonstrate that the nonlinear coupling between the axial and transverse vibrations has a significant effect on the response of structure.. It needs to be considered for the accurate dynamic analysis of long TLP tether subjected to the combined platform surge and heave motions.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

Micromechanical Modelling of Non-Homogenous Materials by Meshless Methods

A detailed study of bi-material composites, using meshless methods (MMs), is presented in this paper. Firstly, representative volume elements (RVEs) for different bi-material combinations are analysed by the element-free Galerkin (EFG) method in order to confirm the effective properties of heterogeneous material through homogenization. The results are shown to be in good agreement with experimental results and those obtained using the finite element method (FEM) which required a higher node density. Secondly, a functionally graded material (FGM), with a crack, is analysed using the EFG method. This investigation was motivated by the possibility of replacing the distinct fibrematrix interface with a FGM interface. Finally, an illustrative example showing crack propagation, in a two-dimension micro-scale model of a SiC/Al composite is presented. 

Pulse Wave Propagating in Coplanar Waveguides (cPWS) on LTCC Substrates

The propagation of pulse waves in coplanar waveguides (CPWs) is investigated, and these CPWs are assumed to be fabricated on a single -layer low- temperature co-fired ceramic (LTCC) substrate. The input pulse wave can be a Gaussian pulse or a sinusoldally modulated Gaussian pulse. Based on the standard Galerkin 's method in the spectral domain, combined with fast Fourier transform (FFT), the pulse waveform and delay in CPWs are demonstrated and compared for a second plate, oriented orthogonally to the primary planar element, thus producing a crossed planar monopole (CPM), which is simpler to produce and has lower cost than a conical monopole. In this paper, further measurements have been made on this element

On integral equation and least squares methods for scattering by diffraction gratings

In this paper we consider the scattering of a plane acoustic or electromagnetic wave by a one-dimensional, periodic rough surface. We restrict the discussion to the case when the boundary is sound soft in the acoustic case, perfectly reflecting with TE polarization in the EM case, so that the total field vanishes on the boundary. We propose a uniquely solvable first kind integral equation formulation of the problem, which amounts to a requirement that the normal derivative of the Green's representation formula for the total field vanish on a horizontal line below the scattering surface. We then discuss the numerical solution by Galerkin's method of this (ill-posed) integral equation. We point out that, with two particular choices of the trial and test spaces, we recover the so-called SC (spectral-coordinate) and SS (spectral-spectral) numerical schemes of DeSanto et al., Waves Random Media, 8, 315-414 1998. We next propose a new Galerkin scheme, a modification of the SS method that we term the SS* method, which is an instance of the well-known dual least squares Galerkin method. We show that the SS* method is always well-defined and is optimally convergent as the size of the approximation space increases. Moreover, we make a connection with the classical least squares method, in which the coefficients in the Rayleigh expansion of the solution are determined by enforcing the boundary condition in a least squares sense, pointing out that the linear system to be solved in the SS* method is identical to that in the least squares method. Using this connection we show that (reflecting the ill-posed nature of the integral equation solved) the condition number of the linear system in the SS* and least squares methods approaches infinity as the approximation space increases in size. We also provide theoretical error bounds on the condition number and on the errors induced in the numerical solution computed as a result of ill-conditioning. Numerical results confirm the convergence of the SS* method and illustrate the ill-conditioning that arises.

Desenvolvimento de antenas de microfita com Patch em anel utilizando materiais ferrimagnéticos e metamateriais

In general, the materials used as substrates in the project of microstrip antennas are: isotropic, anisotropic dielectrics and ferrimagnetic materials (magnetic anisotropy). The use of ferrimagnetic materials as substrates in microstrip patch antennas has been concentrated on the analysis of antennas with circular and rectangular patches. However, a new class of materials, called metamaterials, has been currently the focus of a great deal of interest. These materials exhibit bianisotropic characteristics, with permittivity and permeability tensors. The main objective of this work is to develop a theoretical and numerical analysis for the radiation characteristics of annular ring microstrip antennas, using ferrites and metamaterials as substrates. The full wave analysis is performed in the Hankel transform domain through the application of the Hertz vector potentials. Considering the definition of the Hertz potentials and imposing the boundary conditions, the dyadic Green s function components are obtained relating the surface current density components at the plane of the patch to the electric field tangential components. Then, Galerkin s method is used to obtain a system of matrix equations, whose solution gives the antenna resonant frequency. From this modeling, it is possible to obtain numerical results for the resonant frequency, radiation pattern, return loss, and antenna bandwidth as a function of the annular ring physical parameters, for different configurations and substrates. The theoretical analysis was developed for annular ring microstrip antennas on a double ferrimagnetic/isotropic dielectric substrate or metamaterial/isotropic dielectric substrate. Also, the analysis for annular ring microstrip antennas on a single ferrimagnetic or metamaterial layer and for suspended antennas can be performed as particular cases

Caracterização de antenas planares com substrato metamaterial

This work presents a theoretical and numerical analysis for the radiation characteristics of rectangular microstrip antenna using metamaterial substrate. The full wave analysis is performed in the Fourier transform domain through the application of the Transverse Transmission Line - TTL method. A study on metamaterial theory was conducted to obtain the constructive parameters, which were characterized through permittivity and permeability tensors to arrive at a set of electromagnetic equations. The general equations for the electromagnetic fields of the antenna are developed using the Transverse Transmission Line - TTL method. Imposing the boundary conditions, the dyadic Green s function components are obtained relating the surface current density components at the plane of the patch to the electric field tangential components. Then, Galerkin s method is used to obtain a system of matrix equations, whose solution gives the antenna resonant frequency. From this modeling, it is possible to obtain numerical results for the resonant frequency and return loss for different configurations and substrates

Antenas de microfita com patch em anel e múltiplas camadas dielétricas anisotrópicas uniaxiais

This work presents an analysis of the annular ring microstrip antennas printed on uniaxial anisotropic substrates and with superstrate.The analysis uses the full-wave formulation by means of the Hertz vector potentials method, in the Hankel transform domain. The definition of the Hertz vector potentials and the application of the appropriate boundary conditions to the structure allow determining the dyadic Green functions, relating the current densities in the conducting patch to the transforms of the tangential electric field components. Galerkin s method is then used to obtain the matrix equation whose nontrivial solution gives the complex resonant frequency of the antenna. From the modeling, it is possible to obtain results for the resonant frequency, bandwidth and quality factor, as a function of several parameters of the antenna, for different configurations. We have considered annular ring microstrip antennas on a single dielectric layer, antennas with two anisotropic dielectric layers, and annular ring microstrip antennas on suspended substrates. Numerical results for the resonant frequency of the these structures printed on isotropic substrates are also presented and compared with those published by other authors, showing a good agreement

Uma contribuição ao estudo do ajuste ótimo de parâmetros para cálculo de campos eletromagnéticos descontínuos com o método EFG

Pós-graduação em Ciência da Computação - IBILCE

Método numérico-analítico generalizado para estimação do campo eletromagnético de linhas de transmissão de energia elétrica utilizando a teoria dos elementos finitos

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

Mixing Snapshots and Fast Time Integration of PDEs

A local proper orthogonal decomposition (POD) plus Galerkin projection method was recently developed to accelerate time dependent numerical solvers of PDEs. This method is based on the combined use of a numerical code (NC) and a Galerkin sys- tem (GS) in a sequence of interspersed time intervals, INC and IGS, respectively. POD is performed on some sets of snapshots calculated by the numerical solver in the INC inter- vals. The governing equations are Galerkin projected onto the most energetic POD modes and the resulting GS is time integrated in the next IGS interval. The major computa- tional e®ort is associated with the snapshots calculation in the ¯rst INC interval, where the POD manifold needs to be completely constructed (it is only updated in subsequent INC intervals, which can thus be quite small). As the POD manifold depends only weakly on the particular values of the parameters of the problem, a suitable library can be con- structed adapting the snapshots calculated in other runs to drastically reduce the size of the ¯rst INC interval and thus the involved computational cost. The strategy is success- fully tested in (i) the one-dimensional complex Ginzburg-Landau equation, including the case in which it exhibits transient chaos, and (ii) the two-dimensional unsteady lid-driven cavity problem

Analysis of Meshfree Methods for Lagrangian Fluid-Structure Interaction

In this dissertation a new numerical method for solving Fluid-Structure Interaction (FSI) problems in a Lagrangian framework is developed, where solids of different constitutive laws can suffer very large deformations and fluids are considered to be newtonian and incompressible. For that, we first introduce a meshless discretization based on local maximum-entropy interpolants. This allows to discretize a spatial domain with no need of tessellation, avoiding the mesh limitations. Later, the Stokes flow problem is studied. The Galerkin meshless method based on a max-ent scheme for this problem suffers from instabilities, and therefore stabilization techniques are discussed and analyzed. An unconditionally stable method is finally formulated based on a Douglas-Wang stabilization. Then, a Langrangian expression for fluid mechanics is derived. This allows us to establish a common framework for fluid and solid domains, such that interaction can be naturally accounted. The resulting equations are also in the need of stabilization, what is corrected with an analogous technique as for the Stokes problem. The fully Lagrangian framework for fluid/solid interaction is completed with simple point-to-point and point-to-surface contact algorithms. The method is finally validated, and some numerical examples show the potential scope of applications.

A Petrov-Galerkin method for a singularly perturbed ordinary differential equation with non-smooth data

In this paper, a singularly perturbed ordinary differential equation with non-smooth data is considered. The numerical method is generated by means of a Petrov-Galerkin finite element method with the piecewise-exponential test function and the piecewise-linear trial function. At the discontinuous point of the coefficient, a special technique is used. The method is shown to be first-order accurate and singular perturbation parameter uniform convergence. Finally, numerical results are presented, which are in agreement with theoretical results.

Meshless Local Petrov-Galerkin Method for Rotating Euler-Bernoulli Beam

Free vibration problem of a rotating Euler-Bernoulli beam is solved with a truly meshless local Petrov-Galerkin method. Radial basis function and summation of two radial basis functions are used for interpolation. Radial basis function satisfies the Kronecker delta property and makes it simpler to apply the essential boundary conditions. Interpolation with summation of two radial basis functions increases the node carrying capacity within the sub-domain of the trial function and higher natural frequencies can be computed by selecting the complete domain as a sub-domain of the trial function. The mass and stiffness matrices are derived and numerical results for frequencies are obtained for a fixed-free beam and hinged-free beam simulating hingeless and articulated helicopter blades. Stiffness and mass distribution suitable for wind turbine blades are also considered. Results show an accurate match with existing literature.