Computation of worst geometric imperfection profiles of composite cylindrical shell panels by minimizing the non-linear buckling load

In this article the “most unfavorable” shape of initial geometric imperfection profile for laminated cylindrical shell panel is obtained analytically by minimizing the limit point load. The partial differential equations governing the shell stability problem are reduced to a set of non-linear algebraic equations using Galerkin's technique. The non-linear equilibrium path is traced by employing Newton–Raphson method in conjunction with the Riks approach. A double Fourier series is used to represent the initial geometric imperfection profile for the cylindrical shell panel. The optimum values of these Fourier coefficients are determined by minimizing the limit point load using genetic algorithm. The results are determined for simply supported composite cylindrical shell panel. Numerical results show that more number of terms is needed in Fourier series representation to obtain the “worst” geometric imperfection profile which gives lower limit load compared to single term representation of imperfection. We have incorporated constraints on the shape of imperfection to avoid unrealistic limit point loads (due to imperfection shape) as we have assumed that the imperfection is due to machining/manufactuting.     

Analysis of the effects of a pulsed electromagnetic field on the dynamic response of electrically conductive composites

The effects of pulsed electromagnetic fields on the dynamic mechanical response of electrically conductive anisotropic plates are studied. The analysis is based on the simultaneous solving of the system of nonlinear partial differential equations that include equations of motion and Maxwell’s equations. Physics-based hypotheses for electro-magneto-mechanical coupling in anisotropic composite plates and dimension reduction solution procedures for the nonlinear system of the governing equations are presented. A numerical solution procedure for the resulting two-dimensional nonlinear system of the governing equations has been developed and consists of the sequential application of time and spatial integration and quasilinearization. The developed methodology is applied to the problem of the dynamic response of a long current-carrying unidirectional carbon fiber polymer matrix composite plate subjected to transverse impact load and in-plane pulsed electromagnetic load. The interacting effects of the pulsed electric current, external magnetic field, and mechanical load are studied.     

A new approach to the stabilization of a fluid–structure interaction model

This article is a continuation of our work on a linear fluid–structure interaction model [Grobbelaar-Van Dalsen, On a fluid–structure model in which the dynamics of the structure involves the shear stress due to the fluid, J. Math. Fluid Mech. 10(3) (2008), pp. 388–401; Grobbelaar-Van Dalsen, Strong stability for a fluid––structure model, Math. Methods Appl. Sci., 32(2009) pp. 1452–1466]. The model describes the interaction between a 3-D incompressible fluid and a 2-D plate, the interface, which coincides with a flat flexible part of the surface of the vessel containing the fluid. The mathematical model comprises the Stokes equations and the equations for the longitudinal deflections of the plate with the inclusion of the shear stress that the fluid exerts on the plate. A dissipative damping mechanism of Kelvin–Voigt type is applied to the interior of the plate. While our earlier work shows that weak solutions in a space of finite energy are strongly asymptotically stable under no-slip transmission conditions at the interface with uniform exponential stability only attainable under an additional domination condition, the present research is directed at achieving uniform exponential stability of weak solutions without imposing the domination condition. Using energy methods we establish uniform exponential decay under a modified transmission condition at the interface. This condition entails that the fluid velocity at the interface is coupled to a linear combination of the plate velocity and displacement.     

Mathematical problems of liquid crystal flows

This article deals with flow phenomena in liquid crystal materials. The model that I analyse is that due to Ericksen for liquid crystals with variable degree of orientation. One of the goals is understanding the role of the order parameter in the model. I examine the multiplicity of solutions, phases, that occur in shear flow regimes and analyse their stability under physically realistic conditions. In particular, I show that ellipticity of the linearized system of governing equations is a consequence of the Clausius–Duhem inequality.     

Stability of hydrostatic-pressured FGM thick rings with material nonlinearity

Buckling behaviors of elastoplastic ceramic/metallic functionally graded material (FGM) rings are investigated by using the first order shear deformation theory. The hydrostatic-pressured rings are assumed to be in both the plane-stress case and the plane-strain case, which lead respectively to a uniaxial and a biaxial elastoplastic stress states in prebuckling stage. A uniform strain hypothesis helps to deal with the elastoplastic stress states. By introducing in the graded material properties, the constitutive model of FGMs is formulated under the framework of J₂ deformation theory. By considering the kinetic relations of von-Kárman type and employing the principle of virtual displacement, the equilibrium equations and the buckling governing equations of FGM circular rings are formulated, and the analytical solution of the anisotropic rings is obtained. Finally, the elastoplastic buckling problem is numerically solved through a semi-analytical method, which is proposed to seek the real circumferential strain of FGM rings at the buckling point and determinate the elastoplastic buckling critical hydrostatic pressure. The effects of the inhomogeneous and geometrical parameters on the buckling critical load and the position of the elastoplastic interface are discussed. Results show that, in both the plane-stress and the plane-strain cases, the elastoplastic critical loads are generally lower than their elastic counterparts due to material flow, and the plane-strain critical load is generally larger than the plane-stress one. The elastoplastic critical load does not always decrease monotonously with the increase of the inhomogeneous parameters, which is quite different from their elastic counterparts.     

1D gas dynamic modelling of mass conservation in engine duct systems with thermal contact discontinuities

A detailed analysis of mass non-conservation in the proximity of thermal contact discontinuities, when solving 1-D gas dynamic flow equations with finite difference numerical methods, is carried out in this paper. A wide spectrum of finite difference numerical methods has been applied to solve such conditions. Thermal contact discontinuities are very common in current diesel engines due to back-flow in the intake valves during the valve overlap period. Every method has been shown to be incapable of correctly solving the problem raised, displaying (or revealing) a different behavior. Taking as base line these analyses a study regarding mesh size reduction in ducts has been also performed. This solution becomes suitable since it leads to making mass conservation problems disappear. Nevertheless, most extended calculation structure in 1D gas dynamic models is not advised due to the increase of computational effort required. Thus, a new calculation structure for solving governing equations in ducts is suggested. This proposed calculation structure is based on independent time discretisation of every duct according to its CFL stability criterion. Its application to thermal contact discontinuities points out its advantages with regard to computational demand as the calculation time of every duct is adapted to its mesh size.     

Elasticity solutions for functionally graded rectangular plates with two opposite edges simply supported

England (2006) [13] proposed a novel method to study the bending of isotropic functionally graded plates subject to transverse biharmonic loads. His method is extended here to functionally graded plates with materials characterizing transverse isotropy. Using the complex variable method, the governing equations of three plate displacements appearing in the expansions of displacement field are formulated based on the three-dimensional theory of elasticity for a transverse load satisfying the biharmonic equation. The solution may be expressed in terms of four analytic functions of the complex variable, in which the unknown constants can be determined from the boundary conditions similar to that in the classical plate theory. The elasticity solutions of an FGM rectangular plate with opposite edges simply supported under 12 types of biharmonic polynomial loads are derived as appropriate sums of the general and particular solutions of the governing equations. A comparison of the present results for a uniform load with existing solutions is made and good agreement is observed. The influence of boundary conditions, material inhomogeneity, and thickness to length ratio on the plate deflection and stresses for the load x²yq are studied numerically.     

Nonlinear response of a beam under distributed moving contact load

In this paper, the nonlinear behavior of a one-dimensional model of the disc brake pad is examined. The contact normal force between the disc brake pad lining and rotor is represented by a second order polynomial of the relative displacement between the two elastic bodies. The frictional force due to the sliding motion of the rotor against the stationary pad is modeled as a distributed follower-type axial load with time-dependent terms. By Galerkin discretization, the equation governing the transverse motion of the beam model is reduced to a set of extended Duffing system with quasi-periodically modulated excitations. Retaining the first two vibration modes in the governing equations, frequency response curves are obtained by applying a two-dimensional spectral balance method. For the first time, it is predicted that nonlinearity resulting from the contact mechanics between the disc brake pad lining and rotor can lead to a possible irregular motion (chaotic vibration) of the pad in the neighborhood of simple and parametric resonance. This chaotic behavior is identified and quantitatively measured by examining the Poincaré maps, Fourier spectra, and Lyapunov exponents. It is also found that these chaotic motions emerge as a result of successive Hopf bifurcations characterized by the torus breakdown and torus doubling routes as the excitation frequency varies. Various aspects of the numerical difficulties in the solution of the nonlinear equations are also discussed.     

On boundary layer flow of a Sisko fluid over a stretching sheet

Abstract

In this paper, the steady boundary layer flow of a non-Newtonian fluid over a nonlinear stretching sheet is investigated. The Sisko fluid model, which is combination of power-law and Newtonian fluids in which the fluid may exhibit shear thinning/thickening behaviors, is considered. The boundary layer equations are derived for the two-dimensional flow of an incompressible Sisko fluid. Similarity transformations are used to reduce the governing nonlinear equations and then solved analytically using the homotopy analysis method. In addition, closed form exact analytical solutions are provided for n = 0 and n = 1. Effects of the pertinent parameters on the boundary layer flow are shown and solutions are contrasted with the power-law fluid solutions.     

State of stress in a plat circular ring with a crack : PMM vol. 39, n 6, 1975, pp. 1118–1128

The stress distribution in a circular isotropie ring with a crack on part of the concentric circle is investigated. A system of functional equations governing the coefficients of the complex Fourier series expansion of the stresses acting on the circle on which the crack is located is obtained. The solution of the mentioned system of equations is obtained by using a factorization method, which permitted reduction of the initial system of equations to two coupled infinite systems of algebraic equations. The possibility of using the method of truncation to solve these systems is proved. The singularity originating in the neighborhood of ends of the crack in the formulas governing the stresses is isolated. Stress intensity coefficients for the effect of a uniform load on the external contour are presented.     

A direct method for solving an anisotropic mean curvature flow of plane curves with an external force

A new method for solution of the evolution of plane curves satisfying the geometric equation v=β(x,k,ν), where v is the normal velocity, k and ν are the curvature and tangential angle of a plane curve Γ ? ?² at the point x∈Γ, is proposed. We derive a governing system of partial differential equations for the curvature, tangential angle, local length and position vector of an evolving family of plane curves and prove local in time existence of a classical solution. These equations include a non‐trivial tangential velocity functional governing a uniform redistribution of grid points and thus preventing numerically computed solutions from forming various instabilities. We discretize the governing system of equations in order to find a numerical solution for 2D anisotropic interface motions and image segmentation problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Exact boundary controllability of nodal profile for 1-D quasilinear wave equations

Based on the theory of semi-global C ² solution for 1-D quasilinear wave equations, the local exact boundary controllability of nodal profile for 1-D quasilinear wave equations is obtained by a constructive method, and the corresponding global exact boundary controllability of nodal profile is also obtained under certain additional hypotheses.     

Unified analysis of <Emphasis Type="Italic">BMAP</Emphasis>/<Emphasis Type="Italic">G</Emphasis>/1 cyclic polling models

In this paper we present a unified analysis of the BMAP/G/1 cyclic polling model and its application to the gated and exhaustive service disciplines as examples. The applied methodology is based on the separation of the analysis into service discipline independent and dependent parts. New expressions are derived for the vector-generating function of the stationary number of customers and for its mean in terms of vector quantities depending on the service discipline. They are valid for a broad class of service disciplines and both for zero- and nonzero-switchover-times polling models. We present the service discipline specific solution for the nonzero-switchover-times model with gated and exhaustive service disciplines. We set up the governing equations of the system by using Kronecker product notation. They can be numerically solved by means of a system of linear equations. The resulting vectors are used to compute the service discipline specific vector quantities.     

Influence of variable heat flux on natural convection along a corrugated wall in porous media

In this work the coupled non-linear partial differential equations, governing the free convection from a wavy vertical wall under a power law heat flux condition, are solved numerically. For both Darcy and Forchheimer extended non-Darcy models, a wavy to flat surface transformation is applied and the governing equations are reduced to boundary layer equations. A finite difference scheme based on the Keller Box approach has been used in conjunction with a block tri-diagonal solver for obtaining the solution. Detailed simulations are carried out to investigate the effect of varying parameters such as power law heat flux exponent m, wavelength–amplitude ratio a and the transformed Grashof number Gr′. Both surface undulations and inertial forces increase the temperature of the vertical surface while increasing m reduces it. The wavy pattern observed in surface temperature plots, become more prominent with increasing m or a but reduces as Gr′ increases.     

Local radial basis function collocation method for bending analyses of quasicrystal plates

The local radial basis function collocation method (LRBFCM) is proposed for plate bending analysis in orthorhombic quasicrystals (QCs) under static and transient dynamic loads. Three common types of the plate bending problems are considered: (1) QC plates resting on Winkler foundation (2) QC plates with variable thickness and (3) QC plates under a transient dynamic load. According to the Reissner–Mindlin plate bending theory, there is allowed to simulate the behavior of the two excitations in QC plates, phonon and phason, by 2D strong formulations for the system of governing equations. The governing equations, which describe the phason displacements, are based on Agiasofitou and Lazar elastodynamic model. Numerical results demonstrate the effect of the elastic foundation, as well as plate thickness on the phonon and phason characteristics in this paper. For the transient dynamic analysis, the influence of the phason friction coefficients on the responses of QC plate to transient dynamic loads is also studied.     

Performance modeling of state dependent system with mixed standbys and two modes of failure

The present investigation deals with a multicomponent repairable system with state dependent rates. For smooth functioning of the system, mixed standbys (warm and cold) are provided so that the failed units are immediately replaced by standbys if available. To prevent congestion in the system due to failure of units, permanent along with additional repairmen are provided to restore the failed units. It is assumed that the units may fail in two modes. The units have exponential life time and repair time distributions. The failed unit may balk in case of heavy load of failed units. The failed units may also wait in the queue and renege on finding the repairmen busy according to a pre-specified rule. The Chapman–Kolmogorov equations, governing the model in the form of matrix are constructed using transition flow rates of different states. The steady state solution of queue size distribution is derived using product formula. A cost function is suggested to determine the optimal number of warm and cold standbys units required for the desired level of quality of service. The numerical illustrations are carried out to explore the effect of different parameters on performance measures.     

THE VANISHING PRESSURE LIMIT OF SOLUTIONS TO THE SIMPLIFIED EULER EQUATIONS FOR ISENTROPIC FLUIDS

In this paper,the Riemann problem of the 1-D reduced model for the 2-D Euler equations is considered and the Riemann solutions are obtained.It is proved that,as the pressure vanishes,they converge to two kinds of Riemann solutions to the 1D reduced model for the 2-D transport equations:one contains δ-shocks,the other contains vacuum.     

Convergence of the point vortex method for the 3-D euler equations

We prove consistency, stability, and convergence of a point vortex approximation to the 3-D incompressible Euler equations with smooth solutions. The 3-D algorithm we consider here is similar to the corresponding 3-D vortex blob algorithm introduced by Beale and Majda; see [3]. We first show that the discretization error is second-order accurate. Then we show that the method is stable in l^p norm for the particle trajectories and in w^?1.p norm for discrete vorticity. Consequently, the method converges up to any time for which the Euler equations have a smooth solution. One immediate application of our convergence result is that the vortex filament method without smoothing also converges.     

Effect of pressure work on free convection flow along a vertical circular cone with suction and variable surface temperature

Free convection laminar flow from a vertical circular cone maintained a variable surface temperature with suction and pressure work effects has been investigated. The governing boundary layer equations are transformed into a non-dimensional form and the resulting nonlinear system of partial differential equations are reduced to local non-similarity equations. The governing non-similarity equations are then solved numerically by implicit finite difference method together with Keller box scheme. Numerical results are presented in terms of velocity and temperature profiles of the fluid as well as the local skin-friction coefficients and the local heat transfer rate for different values of Prandtl number Pr, suction parameter ξ, temperature gradient parameter n and the pressure work parameter ∈. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

动力学方程的积分型直接积分法

提出了一个求解动力学问题的新方法(DIM-IM).将动力学方程化成积分方程的形式,借助于该方程构造出了具有显式预测-校正的单步、自起动和四阶精度的积分型直接积分算法.理论分析和算例指出,这一方法较中心差分法、Houbolt法、Newmark法和Wilson-θ法都有较高的精度.本方法适用于强非线性,非保守系统.