Numerical Electromagnetic Analysis of Current and Voltage Distribution along Down Conductors and Towers Hit by Direct Lightning

In order to answer the practically important question of whether the down conductors of lightning protection systems to tall towers and buildings can be electrically isolated from the structure itself, this work is conducted. As a first step in this regard, it is presumed that the down conductor placed on metallic tower will be a pessimistic representation of the actual problem. This opinion was based on the fact that the proximity of heavy metallic structure will have a large damping effect. The post-stroke current distributions along the down conductors and towers, which can be quite different from that in the lightning channel, govern the post-stroke near field and the resulting gradient in the soil. Also, for a reliable estimation of the actual stroke current from the measured down conductor currents, it is essential to know the current distribution characteristics along the down conductors. In view of these, the present work attempts to deduce the post-stroke current and voltage distribution along typical down conductors and towers. A solution of the governing field equations on an electromagnetic model of the system is sought for the investigation. Simulation results providing the spatio-temporal distribution of the post-stroke current and voltage has provided very interesting results. It is concluded that it is almost impossible to achieve electrical isolation between the structure and the down conductor. Furthermore, there will be significant induction into the steel matrix of the supporting structure.

STOKES' SYSTEM IN A DOMAIN WITH OSCILLATING BOUNDARY: HOMOGENIZATION AND ERROR ANALYSIS OF AN INTERIOR OPTIMAL CONTROL PROBLEM

Homogenization and error analysis of an optimal interior control problem in the framework of Stokes' system, on a domain with rapidly oscillating boundary, are the subject matters of this article. We consider a three dimensional domain constituted of a parallelepiped with a large number of rectangular cylinders at the top of it. An interior control is applied in a proper subdomain of the parallelepiped, away from the oscillating volume. We consider two types of functionals, namely a functional involving the L-2-norm of the state variable and another one involving its H-1-norm. The asymptotic analysis of optimality systems for both cases, when the cross sectional area of the rectangular cylinders tends to zero, is done here. Our major contribution is to derive error estimates for the state, the co-state and the associated pressures, in appropriate functional spaces.

A FRAMEWORK FOR THE ERROR ANALYSIS OF DISCONTINUOUS FINITE ELEMENT METHODS FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS AND APPLICATIONS TO C-0 IP METHODS

In this article, an abstract framework for the error analysis of discontinuous Galerkin methods for control constrained optimal control problems is developed. The analysis establishes the best approximation result from a priori analysis point of view and delivers a reliable and efficient a posteriori error estimator. The results are applicable to a variety of problems just under the minimal regularity possessed by the well-posedness of the problem. Subsequently, the applications of C-0 interior penalty methods for a boundary control problem as well as a distributed control problem governed by the biharmonic equation subject to simply supported boundary conditions are discussed through the abstract analysis. Numerical experiments illustrate the theoretical findings.

Approximate Controllability Of Semilinear Systems Using Integral Contractors

In this article, a non-autonomous (time-varying) semilinear system is considered and its approximate controllability is investigated. The notion of 'bounded integral contractor', introduced by Altman, has been exploited to obtain sufficient conditions for approximate controllability. This condition is weaker than Lipschitz condition. The main theorems of Naito [11, 12] are obtained as corollaries of our main results. An example is also given to show how our results weaken the conditions assumed by Sukavanam[17].

Approximate solution of u'=-A2u using a relation between exp(-tA2) and exp(itA)

Let A be a positive definite operator in a Hilbert space and consider the initial value problem for u(t) = -A(2)u. Using a representation of the semigroup exp(-A(2)t) in terms of the group exp(iAt) we express u in terms of the solution of the standard heat equation w(t) = W-yy, with initial values v solving the initial value problem for v(y) = iAv. This representation is used to construct a method for approximating u in terms of approximations of v. In the case that A is a 2(nd) order elliptic operator the method is combined with finite elements in the spatial variable and then reduces the solution of the 4(th) order equation for u to that of the 2(nd) order equation for v, followed by the solution of the heat equation in one space variable.

Optimality Conditions and Duality for Semi-Infinite Mathematical Programming Problem with Equilibrium Constraints

This article considers a semi-infinite mathematical programming problem with equilibrium constraints (SIMPEC) defined as a semi-infinite mathematical programming problem with complementarity constraints. We establish necessary and sufficient optimality conditions for the (SIMPEC). We also formulate Wolfe- and Mond-Weir-type dual models for (SIMPEC) and establish weak, strong and strict converse duality theorems for (SIMPEC) and the corresponding dual problems under invexity assumptions.

Numerical analysis of geocell reinforced ballast overlying soft clay subgrade

Geotextiles and geogrids have been in use for several decades in variety of geo-structure applications including foundation of embankments, retaining walls, pavements. Geocells is one such variant in geosynthetic reinforcement of recent years, which provides a three dimensional confinement to the infill material. Although extensive research has been carried on geocell reinforced sand, clay and layered soil subgrades, limited research has been reported on the aggregates/ballast reinforced with geocells. This paper presents the behavior of a railway sleeper subjected to monotonic loading on geocell reinforced aggregates, of size ranging from 20 to 75 mm, overlying soft clay subgrades. Series of tests were conducted in a steel test tank of dimensions 700 mm x 300 mm x 700 mm. In addition to the laboratory model tests, numerical simulations were performed using a finite difference code to predict the behavior of geocell reinforced ballast. The results from numerical simulations were compared with the experimental data. The numerical and experimental results manifested the importance that the geocell reinforcement has a significant effect on the ballast behaviour. The results depicted that the stiffness of underlying soft clay subgrade has a significant influence on the behavior of the geocell-aggregate composite material in redistributing the loading system.

Three-dimensional numerical analysis of geocell-reinforced soft clay beds by considering the actual geometry of geocell pockets

Due to its complex honeycomb structure, the numerical modeling of the geocell has always been a big challenge. Generally, the equivalent composite approach is used to model the geocells. In the equivalent composite approach, the geocellsoil composite is treated as the soil layer with improved strength and stiffness values. Though this approach is very simple, it is unrealistic to model the geocells as the soil layer. This paper presents a more realistic approach of modeling the geocells in three-dimensional (3D) framework by considering the actual curvature of the geocell pocket. A square footing resting on geocell reinforced soft clay bed was modeled using the ``fast Lagrangian analysis of continua in 3D'' (FLAC(3D)) finite difference package. Three different material models, namely modified Cam-clay, Mohr-Coulomb, and linear elastic were used to simulate the behaviour of foundation soil, infill soil and the geocell, respectively. It was found that the geocells distribute the load laterally to the wider area below the footing as compared to the unreinforced case. More than 50% reduction in the stress was observed in the clay bed in the presence of geocells. In addition to geocells, two other cases, namely, only geogrid and geocell with additional basal geogrid cases were also simulated. The numerical model was systematically validated with the results of the physical model tests. Using the validated numerical model, parametric studies were conducted to evaluate the influence of various geocell properties on the performance of reinforced clay beds.

Numerical simulation of conductive particle behaviour at the surface of a plate electrode affected by a DC corona field

The electric field in certain electrostatic devices can be modeled by a grounded plate electrode affected by a corona discharge generated by a series of parallel wires connected to a DC high-voltage supply. The system of differential equations that describe the behaviour (i.e., charging and motion) of the conductive particle in such an electric field has been numerically solved, using several simplifying assumptions. Thus, it was possible to investigate the effect of various electrical and mechanical factors on the trajectories of conductive particles. This model has been employed to study the behaviour of coalparticles in fly-ash corona separators.

On an analytical-numerical procedure for the analysis of cylindrical shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.

Macromechanical behaviour of composite laminates

An improved higher order transverse shear deformation theory is employed to arrive at modified constitutive relations which can be used in the flexural, buckling and vibration analysis of laminated plates and shells. The strain energy for such systems is then expressed in terms of the displacements and the rotations for ready reference and use. Numerical values of vibration frequencies are obtained using this formulation employing Ritz's method of analysis. The results are compared with those available in the literature to validate the analysis presented.

On an analytical-numerical procedure for the analysis of cylindrical-shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.

On a finite element model for the analysis of through cracks in laminated anisotropic cylindrical shells

A finite element model for the analysis of laminated composite cylindrical shells with through cracks is presented. The analysis takes into account anisotropic elastic behaviour, bending-extensional coupling and transverse shear deformation effects. The proposed finite element model is based on the approach of dividing a cracked configuration into triangular shaped singular elements around the crack tip with adjoining quadrilateral shaped regular elements. The parabolic isoparametric cylindrical shell elements (both singular and regular) used in this model employ independent displacement and rotation interpolation in the shell middle surface. The numerical comparisons show the evidence to the conclusion that the proposed model will yield accurate stress intensity factors from a relatively coarse mesh. Through the analysis of a pressurised fibre composite cylindrical shell with an axial crack, the effect of material orthotropy on the crack tip stress intensity factors is shown to be quite significant.

Stress and strength analysis of pin joints in laminated anisotropic plates

Mechanical joints in composites can be tailored to achieve improved performance and better life by appropriately selecting the laminate parameters. In order to gain the best advantage of this possibility of tailoring the laminate, it is necessary to understand the influence of laminate parameters on the behaviour of joints in composites. Most of the earlier studies in this direction were based on simplified assumptions regarding load transfer at the pin-plate interface and such studies were only carried out on orthotropic and quasi-isotropic laminates. In the present study, a more rigorous analysis is carried out to study pin joints in laminates with anisotropic properties. Two types of laminates with (0/ + ?4/90)s and (0/ ± ?2/90)s layups made out of graphite epoxy T300/5208 material system are considered. The analysis mainly concentrates on clearance fit in which the pin is of smaller diameter compared to the hole. The main aspect of the analysis of pin joints is the changing contact between the pin and the plate with increasing load levels. The analysis is carried out by an iterative finite element technique and a computationally efficient routine is developed for this purpose. Numerical studies indicate that the location and magnitude of the peak stresses along the hole boundary are functions of fibre angle and the overall anisotropic properties. It is also shown that the conventional assumption of cosine distribution for the contact pressure between pin and the plate in the analysis lead to underestimation of bearing failure load and overestimation of shear and tensile failure loads in typical (0/905)s cross-ply laminates.