No description available.

No description available.

Cette thèse traite du développement d’un élément continu de plaques orthotropes, sandwichs et multicouches. La démarche consiste dans un premier temps à établir la matrice de raideur dynamique de plaques orthotropes pour des conditions aux limites naturelles à partir d’une reformulation des éléments de plaques isotropes développés au laboratoire QUARTZ (EA7393). La démarche est basée d’une part sur la décomposition des conditions aux limites libres décrite par Gorman et d’autre part sur la résolution des équations de mouvement en se basant sur les développements en séries de Levy. La matrice de raideur dynamique est ensuite obtenue par projection des déplacements et des efforts de frontières sur des bases fonctionnelles compatibles avec les opérations d’assemblage. Dans un second temps, la formulation des éléments sandwichs et multicouches est décrite par superposition des plaques orthotropes précédemment développées.Les formulations présentées prennent en compte les vibrations de flexion et les vibrations dans le plan, dites vibrations de membrane. La validation de ces éléments est menée par une confrontation systématique de réponses harmoniques non amorties avec celles obtenues par diverses modélisations éléments finis. / This thesis deals with the development of a continuous element for orthotropic, sandwich and multilayer plates. This approach is based essentially on the construction of the dynamic stiffness matrix of orthotropic plates using natural boundary conditions from a reformulation of the isotropic plate elements developed in the QUARTZ laboratory (EA 7393). In order to develop the dynamic stiffness matrix of the studied element we resort on the first hand to the decomposition of free boundary conditions described by Gorman, on the second hand to the resolution of the equations of motion by using Levy series expansions. The dynamic stiffness matrix is then obtained by projecting movements and frontier efforts on functional bases compatible with assembly operations. Finally the continuous sandwich and multilayer plate element is described by superposition of continuous orthotropic plates element previously developed.The formulations presented takes into account the bending vibration and the vibration in the plane, called membrane vibration. The validation of all obtained results is conducted by a systematic comparison of undamped harmonic responses with those obtained by various finite element models.

Uncertainties in load and system properties play a significant role in reliability analysis of vibrating structural systems. The subject of random vibrations has evolved over the last few decades to deal with uncertainties in external loads. A well developed body of literature now exists which documents the status of this subject. Studies on the influ­ence of system property uncertainties on reliability of vibrating structures is, however, of more recent origin. Currently, the problem of dynamic response characterization of sys­tems with parameter uncertainties has emerged as a subject of intensive research. The motivation for this research activity arises from the need for a more accurate assess­ment of the safety of important and high cost structures like nuclear plant installations, satellites and long span bridges. The importance of the problem also lies in understand­ing phenomena like mode localization in nearly periodic structures and deviant system behaviour at high frequencies. It is now well established that these phenomena are strongly influenced by spatial imperfections in the vibrating systems. Design codes, as of now, are unable to systematically address the influence of scatter and uncertainties. Therefore, there is a need to develop robust design algorithms based on the probabilistic description of the uncertainties, leading to safer, better and less over-killed designs. Analysis of structures with parameter uncertainties is wrought with diffi­culties, which primarily arise because the response variables are nonlinearly related to the stochastic system parameters; this being true even when structures are idealized to display linear material and deformation characteristics. The problem is further com­pounded when nonlinear structural behaviour is included in the analysis. The analysis of systems with parameter uncertainties involves modeling of random fields for the system parameters, discretization of these random fields, solutions of stochastic differential and algebraic eigenvalue problems, inversion of random matrices and differential operators, and the characterization of random matrix products. It should be noted that the mathematical nature of many of these problems is substantially different from those which are encountered in the traditional random vibration analysis. The basic problem lies in obtaining the solution of partial differential equations with random coefficients which fluctuate in space. This has necessitated the development of methods and tools to deal with these newer class of problems. An example of this development is the generalization of the finite element methods of structural analysis to encompass problems of stochastic material and geometric characteristics. The present thesis contributes to the development of methods and tools to deal with structural uncertainties in the analysis of vibrating structures. This study is a part of an ongoing research program in the Department, which is aimed at gaining insights into the behaviour of randomly parametered dynamical systems and to evolve computational methods to assess the reliability of large scale engineering structures. Recent studies conducted in the department in this direction, have resulted in the for­mulation of the stochastic dynamic stiffness matrix for straight Euler-Bernoulli beam elements and these results have been used to investigate the transient and the harmonic steady state response of simple built-up structures. In the present study, these earlier formulations are extended to derive the stochastic dynamic stiffness matrix for a more general beam element, namely, the curved Timoshenko beam element. Furthermore, the method has also been extended to study the mean and variance of the stationary response of built-up structures when excited by stationary stochastic forces. This thesis is organized into five chapters and four appendices. The first chapter mainly contains a review of the developments in stochas­tic finite element method (SFEM). Also presented is a brief overview of the dynamics of curved beams and the essence of the dynamic stiffness matrix method. This discussion also covers issues pertaining to modeling rotary inertia and shear deformations in the study of curved beam dynamics. In the context of SFEM, suitability of different methods for modeling system uncertainties, depending on the type of problem, is discussed. The relative merits of several schemes of discretizing random fields, namely, local averaging, series expansions using orthogonal functions, weighted integral approach and the use of system Green functions, are highlighted. Many of the discretization schemes reported in the literature have been developed in the context of static problems. The advantages of using the dynamic stiffness matrix approach in conjunction with discretization schemes based on frequency dependent shape functions, are discussed. The review identifies the dynamic analysis of structures built-up of randomly parametered curved beams, using dynamic stiffness matrix method, as a problem requiring further research. The review also highlights the need for studies on the treatment of non-Gaussian nature of system parameters within the framework of stochastic finite element analysis and simulation methods. The problem of deterministic analysis of curved beam elements is consid­ered first. Chapter 2 reports on the development of the dynamic stiffness matrix for a curved Timoshenko beam element. It is shown that when the beam is uniformly param-etered, the governing field equations can be solved in a closed form. These closed form solutions serve as the basis for the formulation of damping and frequency dependent shape functions which are subsequently employed in the thesis to develop the dynamic stiffness matrix of stochastically inhomogeneous, curved beams. On the other hand, when the beam properties vary spatially, the governing equations have spatially varying coefficients which discount the possibility of closed form solutions. A numerical scheme to deal with this problem is proposed. This consists of converting the governing set of boundary value problems into a larger class of equivalent initial value problems. This set of Initial value problems can be solved using numerical schemes to arrive at the element dynamic stiffness matrix. This algorithm forms the basis for Monte Carlo simulation studies on stochastic beams reported later in this thesis. Numerical results illustrating the formulations developed in this chapter are also presented. A satisfactory agreement of these results has been demonstrated with the corresponding results obtained from independent finite element code using normal mode expansions. The formulation of the dynamic stiffness matrix for a curved, randomly in-homogeneous, Timoshenko beam element is considered in Chapter 3. The displacement fields are discretized using the frequency dependent shape functions derived in the pre­vious chapter. These shape functions are defined with respect to a damped, uniformly parametered beam element and hence are deterministic in nature. Lagrange's equations are used to derive the 6x6 stochastic dynamic stiffness matrix of the beam element. In this formulation, the system property random fields are implicitly discretized as a set of damping and frequency dependent Weighted integrals. The results for a straight Timo- shenko beam are obtained as a special case. Numerical examples on structures made up of single curved/straight beam elements are presented. These examples also illustrate the characterization of the steady state response when excitations are modeled as stationary random processes. Issues related to ton-Gaussian features of the system in-homogeneities are also discussed. The analytical results are shown to agree satisfactorily with corresponding results from Monte Carlo simulations using 500 samples. The dynamics of structures built-up of straight and curved random Tim-oshenko beams is studied in Chapter 4. First, the global stochastic dynamic stiffness matrix is assembled. Subsequently, it is inverted for calculating the mean and variance, of the steady state stochastic response of the structure when subjected to stationary random excitations. Neumann's expansion method is adopted for the inversion of the stochastic dynamic stiffness matrix. Questions on the treatment of the beam characteris­tics as non-Gaussian random fields, are addressed. It is shown that the implementation of Neumann's expansion method and Monte-Carlo simulation method place distinc­tive demands on strategy of modeling system parameters. The Neumann's expansion method, on one hand, requires the knowledge of higher order spectra of beam properties so that the non-Gaussian features of beam parameters are reflected in the analysis. On the other hand, simulation based methods require the knowledge of the range of the stochastic variations and details of the probability density functions. The expediency of implementing Gaussian closure approximation in evaluating contributions from higher order terms in the Neumann expansion is discussed. Illustrative numerical examples comparing analytical and Monte-Carlo simulations are presented and the analytical so­lutions are found to agree favourably with the simulation results. This agreement lends credence to the various approximations involved in discretizing the random fields and inverting the global dynamic stiffness matrix. A few pointers as to how the methods developed in the thesis can be used in assessing the reliability of these structures are also given. A brief summary of contributions made in the thesis together with a few suggestions for further research are presented in Chapter 5. Appendix A describes the models of non-Gaussian random fields employed in the numerical examples considered in this thesis. Detailed expressions for the elements of the covariance matrix of the weighted integrals for the numerical example considered in Chapter 5, are presented in Appendix B; A copy of the paper, which has been ac­cepted for publication in the proceedings of IUTAM symposium on 'Nonlinearity and Stochasticity in Structural Mechanics' has been included as Appendix C.

Two surface wave testing methods, namely, (i) the spectral analysis of surface waves (SASW), and (ii) the multi-channel analysis of surface waves (MASW), form non-destructive and non-intrusive techniques for predicting the shear wave velocity profile of different layers of ground and pavement. These field testing tools are based on the dispersive characteristics of Rayleigh waves, that is, different frequency components of the surface wave travel at different velocities in layered media. The SASW and MASW testing procedure basically comprises of three different components: (i) field measurements by employing geophones/accelerometers, (ii) generating dispersion plots, and (iii) predicting the shear wave velocity profile based on an inversion analysis. For generating the field dispersion plot, the complexities involved while doing the phase unwrapping calculations for the SASW technique, while performing the spectral calculations on the basis of two receivers’ data, makes it difficult to automate since it requires frequent manual judgment. In the present thesis, a new method, based on the sliding Fourier transform, has been introduced. The proposed method has been noted to be quite accurate, computationally economical and it generally overcomes the difficulties associated with the unwrapping of the phase difference between the two sensors’ data. In this approach, the unwrapping of the phase can be carried out without any manual intervention. As a result, an automation of the entire computational process to generate the dispersion plot becomes feasible. The method has been thoroughly validated by including a number of examples on the basis of surface wave field tests as well as synthetic test data. While obtaining the dispersion image by using the MASW method, three different transformation techniques, namely, (i) the Park’s wavefield transform, (ii) the frequency (f) -wavenumber ( ) transform and (iii) the time intercept ( -phase slowness (p) transform have been utilized for generating the multimodal dispersion plots. The performance of these three different methods has been assessed by using synthetic as well as field data records obtained from a ground site by means of 48 geophones. Two-dimensional as well as three-dimensional dispersion plots were generated. The Park’s wavefield transformation method has been found to be especially advantageous since it neither requires a very high sampling rate nor an inclusion of the zero padding of the data in a wavenumber (distance) domain. In the case of an irregular dispersive media, a proper analysis of the higher modes existing in the dispersion plots becomes essential for predicting the shear wave velocity profile of ground on the basis of surface wave tests. In such cases, the establishment of the predominant mode becomes quite significant. In the current investigation for Rayleigh wave propagation, the predominant mode has been computed by maximizing the normalized vertical displacements along the free surface. Eigenvectors computed from the thin layer approach (TLM) approach are analyzed to predict the corresponding predominant mode. It is noted that the establishment of the predominant mode becomes quite important where only two to six sensors are employed and the governing (predominant) modal dispersion curve is usually observed rather than several multiple modes which can otherwise be identified by using around 24 to 48 multiple sensors. By using the TLM, it is, however, not possible to account for the exact contribution of the elastic half space in the dynamic stiffness matrix (DSM) approach. A method is suggested to incorporate the exact contribution of the elastic half space in the TLM. The numerical formulation is finally framed as a quadratic eigenvalue problem which can be easily solved by using the subroutine polyeig in MATLAB. The dispersion plots were generated for several chosen different ground profiles. The numerical results were found to match quite well with the data available from literature. In order to address all the three different aspects of SASW and MASW techniques, a series of field tests were performed on five different ground sites. The ground vibrations were induced by means of (i) a 65 kg mass dropped freely from a height of 5 m, and (ii) by using a 20 pound sledge hammer. It was found that by using a 65 kg mass dropped from a height of 5 m, for stiffer sites, ground exploration becomes feasible even up to a depth of 50-80 m whereas for the softer sites the exploration depth is reduced to about 30 m. By using a 20 lb sledge hammer, the exploration depth is restricted to only 8-10 m due to its low impact energy. Overall, it is expected that the work reported in the thesis will furnish useful guidelines for (i) performing the SASW and MASW field tests, (ii) generating dispersion plots/images, and (iii) predicting the shear wave velocity profile of the site based on an inversion analysis.