The dissertation examines advancements in the methods and techniques used in the field of econometrics. These advancements include: (i) a re-examination of the underlying statistical foundations of statistical models with binary dependent variables. (ii) using feed-forward backpropagation artificial neural networks for modeling dichotomous choice processes, and (iii) the estimation of unconditional demand elasticities using the flexible multistage demand system with asymmetric partitions and fixed effects across time. The first paper re-examines the underlying statistical foundations of statistical models with binary dependent variables using the probabilistic reduction approach. This re-examination leads to the development of the Bernoulli Regression Model, a family of statistical models arising from conditional Bernoulli distributions. The paper provides guidelines for specifying and estimating a Bernoulli Regression Model, as well as, methods for generating and simulating conditional binary choice processes. Finally, the Multinomial Regression Model is presented as a direct extension. The second paper empirically compares the out-of-sample predictive capabilities of artificial neural networks to binary logit and probit models. To facilitate this comparison, the statistical foundations of dichotomous choice models and feed-forward backpropagation artificial neural networks (FFBANNs) are re-evaluated. Using contingent valuation survey data, the paper shows that FFBANNs provide an alternative to the binary logit and probit models with linear index functions. Direct comparisons between the models showed that the FFBANNs performed marginally better than the logit and probit models for a number of within-sample and out-of-sample performance measures, but in the majority of cases these differences were not statistically significant. In addition, guidelines for modeling contingent valuation survey data and techniques for estimating median WTP measures using FFBANNs are examined. The third paper estimates a set of unconditional price and expenditure elasticities for 49 different processed food categories using scanner data and the flexible and symmetric translog (FAST) multistage demand system. Due to the use of panel data and the presence of heterogeneity across time, temporal fixed effects were incorporated into the model. Overall, estimated price elasticities are larger, in absolute terms, than previous estimates. The use of disaggregated product groupings, scanner data, and the estimation of unconditional elasticities likely accounts for these differences. / Ph. D.

Structural Health Monitoring (SHM) is the implementation of damage detection and characterization algorithms using in vitro sensing and actuation for rapidly determining faults in structural systems before the damage leads to catastrophic failure. SHM systems provide near real time information on the state of the integrity of civil, mechanical and aerospace structures. A roadblock in implementing SHM systems in practice is the possibility of false positives introduced by environmental changes. In particular, temperature changes can cause many SHM algorithms to indicate damage when no damage exists. While several experimentally based efforts have been attempted to alleviate temperature effects on SHM algorithms, fundamental research on the effects of temperature on SHM has not been investigated. The work presented in this dissertation composes of two main parts: the first part focuses on the experimental studies of different mechanical structures of aluminum beams, lug samples and railroad switch bolts. The experimental study of the aluminum lug samples and beams is done to propose and examine methods and models for in situ interrogation and detection of damage (in the form of a fatigue crack) in these specimen and to quantify the smallest detectable crack size in aluminum structures. This is done by applying the electrical impedance-based SHM method and using piezoceramic sensors and actuators. Moreover, in order to better extract the damage features from the measured electrical impedance, the ARX non-linear feature extraction is employed. This non-linear feature extraction, compared to the linear one, results in detection of damages in the micro-level size and improves the early detection of fatigue cracks in structures. Experimental results also show that the temperature variation is an important factor in the structural health monitoring applications and its effect on the impedance-based monitoring of the initiation and growth of fatigue cracks in the lug samples is experimentally investigated. The electrical impedance-based SHM technique is also applied in monitoring the loosening of bolted joints in a full-scale railroad switch and the sensitivity of this technique to different levels of loosening of the bolts is investigated. The second part of the work presented here focuses on the analytical study and better understanding of the effect of temperature on the vibration-based SHM. This is done by analytical modeling of the vibratory response of an Euler-Bernoulli beam with two different support conditions of simply supported and clamped-clamped and with a single, non-breathing fatigue crack at different locations along the length of the beam. The effect of temperature variations on the vibratory response of the beam structure is modeled by considering the two effects of temperature-dependent material properties and thermal stress formations inside the structure. The inclusion of thermal effects from both of these points of view (i.e. material properties variations and generation of thermal stresses) as independent factors is investigated and justified by studying the formulations of Helmholtz free energy and stresses inside a body. The effect of temperature variations on the vibratory response of the cracked beam are then studied by integrating these two temperature-related effects into the analytical modeling. The effect of a growing fatigue crack as well as temperature variations and thermal loadings is then numerically studied on the deflection of the beam and the output voltage of a surface-bonded piezoceramic sensor. / Ph. D.

Distributed parameter systems (DPS) are systems defined on infinite dimensional spaces. This includes problems governed by partial differential equations (PDEs) and delay differential equations. In order to numerically implement a controller for a physical system we often first approximate the PDE and the PDE controller using some finite dimensional scheme. However, control design at this level will typically give rise to controllers that are inherently large-scale. This presents a challenge since we are interested in the design of robust, real-time controllers for physical systems. Therefore, a reduction in the size of the model and/or controller must take place at some point. Traditional methods to obtain lower order controllers involve reducing the model from that for the PDE, and then applying a standard control design technique. One such model reduction technique is balanced truncation. However, it has been argued that this type of method may have an inherent weakness since there is a loss of physical information from the high order, PDE approximating model prior to control design. In an attempt to capture characteristics of the PDE controller before the reduction step, alternative techniques have been introduced that can be thought of as controller reduction methods as opposed to model reduction methods. One such technique is LQG balanced truncation. Only recently has theory for LQG balanced truncation been developed in the infinite dimensional setting. In this work, we numerically investigate the viability of LQG balanced truncation as a suitable means for designing low order, robust controllers for distributed parameter systems. We accomplish this by applying both balanced reduction techniques, coupled with LQG, MinMax and central control designs for the low order controllers, to the cable mass, Klein-Gordon, and Euler-Bernoulli beam PDE systems. All numerical results include a comparison of controller performance and robustness properties of the closed loop systems. / Ph. D.

This thesis develops an immersed finite element (IFE) space for numerical simulations arising from beam design with multiple materials. This IFE space is based upon meshes that can be independent of interface of the materials used to form a beam. Both the forward and inverse problems associated with the beam equation are considered. The order of accuracy of this IFE space is numerically investigated from the point of view of both the interpolation and finite element solution of the interface boundary value problems. Both single and multiple interfaces are considered in our numerical simulation. The results demonstrate that this IFE space has the optimal order of approximation capability. / Master of Science

Dans la première partie de la thèse, les schémas d’intégration conservatifs sont appliqués aux poutres co-rotationnelles 2D. Les cinématiques d'Euler-Bernoulli et de Timoshenko sont abordées. Ces formulations produisent des expressions de l'énergie interne et l'énergie cinétique complexe et fortement non-linéaires. L’idée centrale de l’algorithme consiste à définir, par intégration, le champ des déformations en fin de pas à partir du champ de vitesses de déformations et non à partir du champ des déplacements au travers de la relation déplacement-déformation. La même technique est appliquée aux termes d’inerties. Ensuite, une poutre co-rotationnelle plane avec rotules généralisées élasto-(visco)-plastiques aux extrémités est développée et comparée au modèle fibre avec le même comportement pour des problèmes d'impact. Des exemples numériques montrent que les effets de la vitesse de déformation influencent sensiblement la réponse de la structure. Dans la seconde partie de cette thèse, une théorie de poutre spatiale d’Euler-Bernoulli géométriquement exacte est développée. Le principal défi dans la construction d’une telle théorie réside dans le fait qu’il n’existe aucun moyen naturel de définir un trièdre orthonormé dans la configuration déformée. Une nouvelle méthodologie permettant de définir ce trièdre et par conséquent de développer une théorie de poutre spatiale en incorporant l'hypothèse d'Euler- Bernoulli est fournie. Cette approche utilise le processus d'orthogonalisation de Gram-Schmidt couplé avec un paramètre rotation qui complète la description cinématique et décrit la rotation associée à la torsion. Ce processus permet de surmonter le caractère non-unique de la procédure de Gram-Schmidt. La formulation est étendue au cas dynamique et un schéma intégration temporelle conservant l'énergie est également développé. De nombreux exemples démontrent l’efficacité de cette formulation. / In the first part of the thesis, energymomentum conserving algorithms are designed for planar co-rotational beams. Both Euler-Bernoulli and Timoshenko kinematics are addressed. These formulations provide us with highly complex nonlinear expressions for the internal energy as well as for the kinetic energy which involve second derivatives of the displacement field. The main idea of the algorithm is to circumvent the complexities of the geometric non-linearities by resorting to strain velocities to provide, by means of integration, the expressions for the strain measures themselves. Similarly, the same strategy is applied to the highly nonlinear inertia terms. Next, 2D elasto-(visco)-plastic fiber co-rotational beams element and a planar co-rotational beam with generalized elasto-(visco)-plastic hinges at beam ends have been developed and compared against each other for impact problems. In the second part of this thesis, a geometrically exact 3D Euler-Bernoulli beam theory is developed.The main challenge in defining a three-dimensional Euler-Bernoulli beam theory lies in the fact that there is no natural way of defining a base system at the deformed configuration. A novel methodology to do so leading to the development of a spatial rod formulation which incorporates the Euler-Bernoulli assumption is provided. The approach makes use of Gram-Schmidt orthogonalisation process coupled to a one-parametric rotation to complete the description of the torsional cross sectional rotation and overcomes the non-uniqueness of the Gram-Schmidt procedure. Furthermore, the formulation is extended to the dynamical case and a stable, energy conserving time-stepping algorithm is developed as well. Many examples confirm the power of the formulation and the integration method presented.

Les robots manipulateurs sont généralement des machines rigides, conçues pour que leurflexibilité ne perturbe pas leurs mouvements. En effet, des flexibilités mécaniques importantesdans la structure d’un système introduisent des degrés de liberté supplémentaires dont le comportementest complexe et difficile à maîtriser. Cependant, la réduction de la masse d’un systèmeest bénéfique du point de vue des coûts, de la performance énergétique, de la sécurité et des performancesdynamiques. Afin de faciliter l’accès aux nombreux avantages d’une structure légèremalgré la présence de fortes flexibilités, cette thèse porte sur la modélisation, la conception et lacommande de robots manipulateurs flexibles. Elle est motivée par le projet YAKA, dont l’applicationest le lancement et la récupération de drones à voilure fixe depuis un navire faisant route.Cette application nécessite une importante dynamique sur un vaste espace de travail, bien au-delàdes spécifications des robots rigides classiques. Les outils de modélisation, de conception et decommande proposés prennent en compte la flexibilité des segments et des articulations, pour unnombre quelconque de degrés de liberté et de segments flexibles. Le modèle dynamique flexibleest obtenu par le formalisme de Lagrange, les poutres flexibles sont représentées par le modèled’Euler-Bernoulli. Le schéma de commande proposé se décompose en une inversion de modèledynamique rigide et un bloc de précommande par Input Shaping adapté aux robots manipulateursflexibles. Les outils de conception proposés permettent de baser le processus de conceptionsur des performances prédites du système complet muni de ses actionneurs et de son contrôleuravec une simulation réaliste. Les validations expérimentales effectuées sur le robot YAKA permettentde valider la pertinence de la démarche suivie. Les résultats du projet YAKA confirment lafaisabilité de la mise en oeuvre d’un robot flexible de grande envergure et à forte dynamique dansun contexte industriel, en particulier pour le lancement et la récupération d’un drone à voilurefixe depuis un navire faisant route. / Robot manipulators are generally stiff machines, designed in a way that flexibility does not affecttheir movements. Indeed, significant flexibility introduces additional degrees of freedom witha complex behavior. However, reducing the mass of a system allows for costs, performance, andsafety improvements. In order to allow those benefits despite important flexibility, this thesis focuseson modeling, design and control of flexible robot manipulators. It is motivated by the YAKAproject, which aims at developing a robot to launch and recover fixed wing UAVs from a movingship. It implies reaching very high dynamics on a large workspace, way beyond the specificationsof common rigid robots. The proposed tools for modeling, design and control allow for taking intoaccount both joint and link flexibility, for any number of degrees of freedom and flexible links.The elastodynamic model is obtained with Lagrange principle, each flexible link being representedwith one ormany Euler-Bernouilli beams. The proposed control scheme uses a nonlinear rigiddynamic inversion and extends classical Input Shaping techniques to flexible robot manipulators.The proposed design tools allow for performance prediction of the system including its actuatorsand controllers thanks to a realistic simulation. Experiments conducted with the YAKA robot validatedthe proposed approach. The results of the YAKA project confirmed the feasibility of usinga large scale, highly dynamic flexible robot in an industrial context, in particular for UAVs launchand recovery operations from amoving ship.

Les Structures atomiques Quasi-périodiques (QP) possèdent des propriétés particulières, notamment dans le domaine vibrationnel. Il pourrait être intéressant de pouvoir transférer ces propriétés à des méta-matériaux macroscopiques. Des réseaux de poutres quasi-périodiques 2D sont étudiés dans cette thèse dans le cadre du modèle élément finis (EF) poutre Euler Bernoulli. Ces réseaux de poutres peuvent facilement être produits par fabrication additive ou par découpe laser. Il est possible de faire varier l'élancement des poutres (le ratio hauteur sur longueur) qui est un paramètre intéressant pour modifier la réponse mécanique du réseau. En utilisant la méthode EF, l'influence de l'élancement des poutres sur la réponse vibratoire des réseaux de poutres QP est étudiée. La méthode numérique Kernel Polynomial (KPM) est adaptée avec succès de la dynamique moléculaire aux réseaux de poutres pour étudier leurs modes de vibration sans avoir à diagonaliser complètement la matrice dynamique. Les réseaux de poutres QP présentent des propriétés similaires à leur compère atomique: en particulier la localisation de modes sur des sous-structures et une relation de dispersion hiérarchisée. Le comportement à la fracture est aussi étudié étant donné que les symétries présentes dans les QP pourraient permettre des réseaux de poutres ne présentant pas de plans faibles pour la propagation de fissures. Cela a été démontré d'après des calculs EF statiques avec un critère de fracture fragile sur l'énergie de déformation. Les simulations statiques ne suffisent pas car elles ne peuvent pas capturer les phénomènes dynamiques complexes qui apparaissent lors de la fissuration fragile. Les propriétés de vibration du QP pourraient aussi avoir un impact sur la propagation dynamique de fissure. Un modèle dynamique de fissuration est développé afin d'étudier l'impact de l'élancement sur la capacité des réseaux de poutres QP à dissiper de l'énergie par fissuration. Finalement une méthode Coarse Graining est développée pour identifier un milieu Cosserat continu équivalant au réseau de poutres QP pour différentes échelles. Cette méthode permet d'identifier la densité, les déformations, les contraintes et donc les modules d'élasticité du milieu Cosserat équivalent, permettant ainsi une meilleure compréhension du rôle des sous structures précédemment identifiées. / Quasi periodic (QP) structures have shown peculiar properties in the atomistic domain, especially the vibrational one. It could be interesting to be able to transpose these properties in macroscopic meta-materials. Quasi periodic 2D beam lattices are studied in this thesis due to the simplicity of the Euler Bernoulli finite element (FE) model. These beam lattices can easily be produced by additive manufacturing or by laser cutting. It is possible to vary the beam slenderness (i.e the ratio of height over length) that is a interesting parameter to modify the mechanical response of the lattice. Using finite element method, the influence of the beam slenderness over the vibration behavior of the QP beam lattices will be studied. The Kernel Polynomial numerical Method (KPM) is successfully adapted from molecular dynamics simulations in order to study vibrational modes of FE beam lattices without having to fully diagonalize the dynamical matrix. The QP lattices show similar properties as their atomic counterpart e.g mode localization over sub-stuctures and hierarchical dispersion relation. The fracture behavior is also studied, as the special symmetries allowed by the quasi periodicity could result in beam lattices without weak planes for crack propagation. It was proved to be true from static FE simulations with a brittle strain energy breaking criterion. Static simulations were not enough and do not grasp the complex dynamical phenomena taking place in brittle fracture. A dynamic crack propagation model was thus developed. The vibrational properties of quasi periodic structures could also have an impact on the dynamic crack propagation. Several simulations are run in order to study the impact of the slenderness on the energy dissipated by fracture of QP lattices. Finally, a coarse graining method (CG) was developed to identify a continuous Cosserat medium at different scales from the FE beam model. This CG method allows to identify, density, strain, stress and elastic moduli of an equivalent continuous Cosserat. This allows a better understanding of the role of previously identified characteristic sub structures.

La détection et le pistage de cibles de surface, maritimes ou terrestres, constituent l’un des champs d’application de la surveillance par radar aéroporté. Dans ce contexte spécifique, il s’agit d’estimer les trajectoires d’un ou de plusieurs objets mobiles au cours du temps à partir de mesures radar bruitées. Cependant, plusieurs contraintes s’additionnent au problème d’estimation des trajectoires :1. le nombre d’objets présents dans la région d’intérêt est inconnu et peut évoluer au cours du temps,2. les mesures fournies par le radar ne correspondent pas toutes à des objets mobiles car certaines sont dues à l’environnement ; il s’agit de fausses alarmes,3. une mesure n’est pas toujours disponible pour chaque objet à chaque instant ; il s’agit de non-détections,4. les cibles de surface peuvent être très diverses en termes de capacité de manoeuvre.Pour tenir compte des trois premières exigences, les modèles d’ensembles finis aléatoires peuvent être envisagés pour procéder aux estimations simultanées du nombre d’objets et de leur trajectoire dans un formalisme bayésien. Pour répondre à la quatrième contrainte, une classification des objets à pister peut s’avérer utile. Aussi, dans le cadre de cette thèse, nous nous intéressons à deux traitements adaptatifs qui intègrent ces deux principes.Tout d’abord, nous proposons une approche conjointe de pistage et de classification dédiée au cas d’un objet évoluant en présence de fausses alarmes. Notre contribution réside dans le développement d’un algorithme incorporant un filtre fondé sur un ensemble fini aléatoire de Bernoulli. L’algorithme résultant combine robustesse aux fausses alarmes et capacité à classer l’objet. Cette classification peut être renforcée grâce à l’estimation d’un paramètre discriminant comme la longueur, qui est déduite d’une mesure d’étalement distance.Le second traitement adaptatif présenté dans cette thèse est une technique de pistage de groupes de cibles dont les mouvements sont coordonnés. Chaque groupe est caractérisé par un paramètre commun définissant la coordination des mouvements de ses cibles. Cependant, ces dernières conservent une capacité de manoeuvre propre par rapport à la dynamique de groupe. S’appuyant sur le formalisme des ensembles finis aléatoires, la solution proposée modélise hiérarchiquement la configuration multi-groupes multi-cibles. Au niveau supérieur, la situation globale est représentée par un ensemble fini aléatoire dont les éléments correspondent aux groupes de cibles. Ils sont constitués du paramètredu groupe et d’un ensemble fini aléatoire multi-cibles. Ce dernier contient les vecteurs d’état des cibles du groupe dont le nombre peut évoluer au cours du temps. L’algorithme d’estimation développé est lui-aussi organisé de manière hiérarchique. Un filtre multi-Bernoulli labélisé (LMB) permet d’estimer le nombre de groupes, puis pour chacun d’entre eux, leur probabilité d’existence ainsi que leur paramètre commun. Pour ce faire, le filtre LMB interagit avec un banc de filtres multi-cibles qui opèrent conditionnellement à une hypothèse de groupe. Chaque filtre multi-cibles estime le nombre et les vecteurs d’état des objets du groupe. Cette approche permet de fournir à l’opérationnel des informations sur la situation tactique. / Detecting and tracking maritime or ground targets is one of the application fields for surveillance by airborne radar systems. In this specific context, the goal is to estimate the trajectories of one or more moving objects over time by using noisy radar measurements. However, several constraints have to be considered in addition to the problem of estimating trajectories:1. the number of objects inside the region of interest is unknown and may change over time,2. the measurements provided by the radar can arise from the environment and do not necessarily correspond to a mobile object; the phenomenon is called false detection,3. a measurement is not always available for each object; the phenomenon is called non-detection,4. the maneuverability depends on the surface targets.Concerning the three first points, random finite set models can be considered to simultaneously estimate the number of objects and their trajectories in a Bayesian formalism. To deal with the fourth constraint, a classification of the objects to be tracked can be useful. During this PhD thesis, we developped two adaptive approaches that take into account both principles.First of all, we propose a joint target tracking and classification method dedicated to an object with the presence of false detections. Our contribution is to incorporate a filter based on a Bernoulli random finite set. The resulting algorithm combines robustness to the false detections and the ability to classify the object. This classification can exploit the estimation of a discriminating parameter such as the target length that can be deduced from a target length extent measurement.The second adaptive approach presented in this PhD dissertation aims at tracking target groups whose movements are coordinated. Each group is characterized by a common parameter defining the coordination of the movements of its targets. However, the targets keep their own capabilities of maneuvering relatively to the group dynamics. Based on the random finite sets formalism, the proposed solution represents the multi-target multi-group configuration hierarchically. At the top level, the overall situation is modeled by a random finite set whose elements correspond to the target groups. They consist of the common parameter of the group and a multi-target random finite set. The latter contains the state vectors of the targets of the group whose number may change over time. The estimation algorithm developed is also organized hierarchically. A labeled multi-Bernoulli filter (LMB) makes it possible to estimate the number of groups, and for each of them, to obtain their probability of existence as well as their common parameter. For this purpose, the LMB filter interacts with a bank of multi-target filters working conditionally to a group hypothesis. Each multi-target filter estimates the number and state vectors of the objects in the group. This approach provides operational information on the tactical situation.

D'Alembert est l'auteur d'un quatrième et dernier traité d'hydrodynamique méconnu des historiens des sciences : le Mémoire 57 du t. VIII de ses "Opuscules mathématiques" (1780). Il y répond, pour l'essentiel, à un écrit de Borda de 1766 contenant une sévère remise en question des principes et concepts fondateurs de son "Traité des fluides" (1744). Cette polémique entre les deux géomètres s'inscrit dans une période de crise de l'hydrodynamique touchant à la concordance entre théorie et expérience. Les équations aux dérivées partielles obtenues par D'Alembert dans son "Essai d'une nouvelle théorie de la résistance des fluides" (1752), puis par Euler dans ses célèbres mémoires de 1755, n'ont pu être résolues. L'approche unidimensionnelle initiée par Daniel Bernoulli dans l'"Hydrodynamica" (1738) puis adoptée par Jean Bernoulli dans l'"Hydraulica" (1743) et par D'Alembert dans le "Traité des fluides" (1744), repose, quant à elle, sur une approximation, l'hypothèse du parallélisme des tranches, trop éloignée des résultats expérimentaux disponibles. C'est cet écart que Borda tente de pallier en 1766 afin de faire sortir la discipline de l'impasse dans laquelle elle se trouve.<br />Dans le Mémoire 57, D'Alembert, directement visé par ce dernier, donne son point de vue sur le sujet, répond implicitement aux critiques de Borda, revient sur ses écrits antérieurs, et propose de nouvelles pistes de recherche théoriques. L'étude de ce traité, tant du point de vue physique que mathématique, engage ici un réexamen des principes et concepts fondateurs de sa théorie des écoulements, la plupart n'ayant été que partiellement abordés par l'historiographie. Il nous permet également de situer ses recherches vis-à-vis de celles des principaux artisans du développement de l'hydrodynamique à cette époque, Daniel Bernoulli, Jean Bernoulli, Euler et Lagrange, et de donner une idée de leur réception et de leur pérennité.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / This work researches on free and forced vibrations of a double beam coupled system. The system is regarded as two Euler-Bernoulli beams which are parallel, have the same length, are simply supported and are connected through a viscoelastic layer. Natural frequencies and their mode shapes, also called eigenfunctions, of the coupled system are obtained through a uniform beam methodology which uses the free dynamical basis to represent the solution of the the modal equation. This study uses modal analysis and block matrix formulation, while the dynamical basis used to represent the modal solution is obtained from the dynamical solution of a fourth order differential equation whose coefficients are just those of the original problem. The natural frequencies and mode shapes of the undamped system are determined for several values of beam parameters. For the damped case, damping ratios of each beam and also of the viscoelastic layer (which characterizes the coupling the system) are considered. The forced response is represented using matrix impulse response, which is the solution of an initial value problem with impulsive initial conditions. / Neste trabalho é realizado um estudo sobre vibrações livres e forçadas de um sistema de dupla viga acoplado. O sistema é composto por duas vigas do tipo Euler-Bernoulli, paralelas, de mesmo comprimento, simplesmente apoiadas e conectadas por uma camada viscoelástica. São obtidas as frequências naturais e os modos de vibração ou autofunções do sistema acoplado utilizando uma metodologia para vigas uniformes, que usa a base dinâmica para escrever a solução da equação modal. O estudo é realizado através da análise modal e de uma formulação matricial em blocos, e a base dinâmica usada para escrever a solução da equação modal é gerada pela solução dinâmica de uma equação diferencial de quarta ordem cujos coeficientes são os mesmos do problema considerado. As frequências naturais e os modos de vibração para o sistema não amortecido são determinados para vários valores dos parâmetros da viga. Para o caso amortecido, consideramos o amortecimento individual em cada viga e o amortecimento que compõe a camada viscoelástica o qual caracteriza o acoplamento no sistema. A resposta forçada do sistema é escrita em função da resposta impulso matricial que é solução de um problema de valor inicial com condições iniciais impulsivas.