Dissertations / Theses on the topic 'Timoshenko beam theory'

In this thesis, transverse vibration of a cracked beam on an elastic foundation and the effect of crack and foundation parameters on transverse vibration natural frequencies are studied. Analytical formulations are derived for a beam with rectangular cross section. The crack is an open type edge crack placed in the medium of the beam and it is uniform along the width of the beam. The cracked beam rests on an elastic foundation. The beam is modeled by two different beam theories, which are Euler-Bernoulli beam theory and Timoshenko beam theory. The effect of the crack is considered by representing the crack by rotational springs. The compliance of the spring that represents the crack is obtained by using fracture mechanics theories. Different foundation models are discussed
these models are Winkler Foundation, Pasternak Foundation, and generalized foundation. The equations of motion are derived by applying Newton'
s 2nd law on an infinitesimal beam element. Non-dimensional parameters are introduced into equations of motion. The beam is separated into pieces at the crack location. By applying the compatibility conditions at the crack location and boundary conditions, characteristic equation whose roots give the non-dimensional natural frequencies is obtained. Numerical solutions are done for a beam with square cross sectional area. The effects of crack ratio, crack location and foundation parameters on transverse vibration natural frequencies are presented. It is observed that existence of crack reduces the natural frequencies. Also the elastic foundation increases the stiffness of the system thus the natural frequencies. The natural frequencies are also affected by the location of the crack.

This dissertation is focused on understanding the performance of a particular fiber-reinforced polymeric composite structural beam, a 91.4 cm (36 inch) deep pultruded double-web beam (DWB) designed for bridge construction. Part 1 focuses on calculating the Timoshenko shear stiffness of the DWB and understanding what factors may introduce error in the experimental measurement of the quantity for this and other sections. Laminated beam theory and finite element analysis (FEA) were used to estimate the shear stiffness. Several references in the literature have hypothesized an increase in the effective measured shear stiffness due to warping. A third order laminated beam theory (TLBT) was derived to explore this concept, and the warping effect was found to be negligible. Furthermore, FEA results actually indicate a decrease in the effective shear stiffness at shorter spans for simple boundary conditions. This effect was attributed to transverse compression at the load points and supports. The higher order sandwich theory of Frostig shows promise for estimating the compression related error in the shear stiffness for thin-walled beams. Part 2 attempts to identify the failure mechanism(s) under quasi-static loading and to develop a strength prediction for the DWB. FEA was utilized to investigate two possible failure modes in the top flange: compression failure of the carbon fiber plies and delamination at the free edges or taper regions. The onset of delamination was predicted using a strength-based approach, and the stress analysis was accomplished using a successive sub-modeling approach in ANSYS. The results of the delamination analyses were inconclusive, but the predicted strengths based on the compression failure mode show excellent agreement with the experimental data. A fatigue life prediction, assuming compression failure, was also developed using the remaining strength and critical element concepts of Reifsnider et al. One DWB fatigued at about 30% of the ultimate capacity showed no signs of damage after 4.9 million cycles, although the predicted number of cycles to failure was 4.4 million. A test on a second beam at about 60% of the ultimate capacity was incomplete at the time of publication. Thus, the success of the fatigue life prediction was not confirmed.
Ph. D.

Se presenta un modelo de elementos finitos que describe el comportamiento de vibración de libre de vigas compuestas laminadas. Se desarrolla el modelo utilizando el principio de Hamilton y la teoría de vigas Timoshenko que incluye deformaciones por corte. Se asume interpolaciones de alto orden para la aproximación de las variables fundamentales. Los laminados compuestos son ortotrópicos con fibras orientadas en diferentes direcciones. Se implementa un programa para materiales compuestos laminado en MATLAB. Se comparan resultados con otros obtenidos en la literatura para validar el modelo. Se realiza un estudio de convergencia y paramétrico con un mismo número de lámina y diferentes direcciones. Se verifica que la formulación que es bastante precisa con resultados satisfactorios en la investigación.
In this work, is presented a finite element model that describes the free vibration behavior of laminated composite beams. The model is developed by the Hamilton principle and the Timoshenko theory that includes shear deformations. Composite laminates are assumed to be orthotropic with fibers oriented in different directions, such as Angle Ply and Cross Ply cases. This investigation works out on a MAPLE program for laminated composites materials that will be completed all in MATLAB program. In order to validate the model, the results are compared with different literatures, also verify the formulation that is quite accurate and obtain quite satisfactory results in the investigation. High order interpolations are assumed to approximate fundamental variables. A convergence study and parametric study will be carried out with the same number of laminas in different directions.
Tesis

A clássica teoria de Euler-Bernoulli para vibrações transversais de vigas elásticas é sabido não ser adequada para vibrações de altas freqüências, como é o caso de vibração de vigas curtas. Esta teoria assume que a deflexão deve-se somente ao momento fletor, uma vez que os efeitos da inércia de rotação e da força cortante são negligenciados. Lord Rayleigh complementou a teoria clássica demonstrando a contribuição da inércia de rotação e Timoshenko estendeu a teoria ao incluir os efeitos da força cortante. A equação resultante é conhecida como sendo a que caracteriza a chamada teoria de viga de Timoshenko. Usando-se a matriz de rigidez dinâmica, as freqüências naturais e a resposta dinâmica de estruturas de barras são determinadas e comparadas de acordo com resultados de quatro modelos de vibração. São estudados o problema de vibração flexional de vigas, pórticos e grelhas, bem como o problema de fundação elástica segundo o modelo de Winkler e também a versão mais avançada que é o modelo de Pasternak.
Classical Euler-Bernoulli theory for transverse vibrations of elastic beams is known to be inadequate to consider high frequency modes which occur for short beams, for example. This theory is derived under the assumption that the deflection is only due to bending. The effects of rotary inertia and shear deformation are ignored. Lord Rayleigh improved the classical theory by considering the effect of rotary inertia. Timoshenko extended the theory to include the effects of shear deformation. The resulting equation is known as Timoshenko beam theory. The natural frequencies and dynamic reponse of framed structures are determined by using the dynamic stiffness matrix and compered according to these theories. The flexional vibration problems of beams, plane frames and grids are analysed, as well problems of elastic foundation according the well known Winkler model and also the more general Pasternak model.

The thesis presents a novel computational method for analysing the static and dynamic behaviour of a multi-damaged beam using local and non-local elasticity theories. Most of the lumped damage beam models proposed to date are based on slender beam theory in classical (local) elasticity and are limited by inaccuracies caused by the implicit assumption of the Euler-Bernoulli beam model and by the spring model itself, which simplifies the real beam behaviour around the crack. In addition, size effects and material heterogeneity cannot be taken into account using the classical elasticity theory due to the absence of any microstructural parameter. The proposed work is based on the inhomogeneous Euler-Bernoulli beam theory in which a Dirac's delta function is added to the bending flexibility at the position of each crack: that is, the severer the damage, the larger is the resulting impulsive term. The crack is assumed to be always open, resulting in a linear system (i.e. nonlinear phenomena associated with breathing cracks are not considered). In order to provide an accurate representation of the structure's behaviour, a new multi-cracked beam element including shear effects and rotatory inertia is developed using the flexibility approach for the concentrated damage. The resulting stiffness matrix and load vector terms are evaluated by the unit-displacement method, employing the closed-form solutions for the multi-cracked beam problem. The same deformed shapes are used to derive the consistent mass matrix, also including the rotatory inertia terms. The two-node multi-damaged beam model has been validated through comparison of the results of static and dynamic analyses for two numerical examples against those provided by a commercial finite element code. The proposed model is shown to improve the computational efficiency as well as the accuracy, thanks to the inclusion of both shear deformations and rotatory inertia. The inaccuracy of the spring model, where for example for a rotational spring a finite jump appears on the rotations' profile, has been tackled by the enrichment of the elastic constitutive law with higher order stress and strain gradients. In particular, a new phenomenological approach based upon a convenient form of non-local elasticity beam theory has been presented. This hybrid non-local beam model is able to take into account the distortion on the stress/strain field around the crack as well as to include the microstructure of the material, without introducing any additional crack related parameters. The Laplace's transform method applied to the differential equation of the problem allowed deriving the static closed-form solution for the multi-cracked Euler-Bernoulli beams with hybrid non-local elasticity. The dynamic analysis has been performed using a new computational meshless method, where the equation of motions are discretised by a Galerkin-type approximation, with convenient shape functions able to ensure the same grade of approximation as the beam element for the classical elasticity. The importance of the inclusion of microstructural parameters is addressed and their effects are quantified also in comparison with those obtained using the classical elasticity theory.

Piston rings are installed in the piston and cylinder wall, which does not have a perfect round shape due to machining tolerances or external loads e.g. head bolts tightening. If the ring cannot follow these deformations, a localized lack of contact will occur and consequently an increase in the engine blow-by and lubricant oil consumption. Current 2D computational methods can not implement such effects – more complex model is necessary. The presented master’s thesis is focused on the developement of a flexible 3D piston ring model able to capture local deformations. It is based on the Timoshenko beam theory in cooperation with MBS software Adams. Model is then compared with FEM using software ANSYS. The validated piston ring model is assembled into the piston/cylinder liner and very basic simulations are run. Finally, future improvements are suggested.

This thesis is about using structural-dynamics based methods to address the existing challenges in the field of Structural Health Monitoring (SHM). Particularly, new structural-dynamics based methods are presented, to model areas of damage, to do damage diagnosis and to estimate and predict the sensitivity of structural vibration properties like natural frequencies to the presence of damage. Towards these objectives, a general analytical procedure, which yields nth-order expressions governing mode shapes and natural frequencies and for damaged elastic structures such as rods, beams, plates and shells of any shape is presented. Features of the procedure include the following: 1. Rather than modeling the damage as a fictitious elastic element or localized or global change in constitutive properties, it is modeled in a mathematically rigorous manner as a geometric discontinuity. 2. The inertia effect (kinetic energy), which, unlike the stiffness effect (strain energy), of the damage has been neglected by researchers, is included in it. 3. The framework is generic and is applicable to wide variety of engineering structures of different shapes with arbitrary boundary conditions which constitute self adjoint systems and also to a wide variety of damage profiles and even multiple areas of damage. To illustrate the ability of the procedure to effectively model the damage, it is applied to beams using Euler-Bernoulli and Timoshenko theories and to plates using Kirchhoff's theory, supported on different types of boundary conditions. Analytical results are compared with experiments using piezoelectric actuators and non-contact Laser-Doppler Vibrometer sensors. Next, the step of damage diagnosis is approached. Damage diagnosis is done using two methodologies. One, the modes and natural frequencies that are determined are used to formulate analytical expressions for a strain energy based damage index. Two, a new damage detection parameter are identified. Assuming the damaged structure to be a linear system, the response is expressed as the summation of the responses of the corresponding undamaged structure and the response (negative response) of the damage alone. If the second part of the response is isolated, it forms what can be regarded as the damage signature. The damage signature gives a clear indication of the damage. In this thesis, the existence of the damage signature is investigated when the damaged structure is excited at one of its natural frequencies and therefore it is called ``partial mode contribution". The second damage detection method is based on this new physical parameter as determined using the partial mode contribution. The physical reasoning is verified analytically, thereupon it is verified using finite element models and experiments. The limits of damage size that can be determined using the method are also investigated. There is no requirement of having a baseline data with this damage detection method. Since the partial mode contribution is a local parameter, it is thus very sensitive to the presence of damage. The parameter is also shown to be not affected by noise in the detection ambience.

Esta investigación se enfoca en el análisis de vibración libre de vigas Timoshenko utilizando el método de elementos finitos. Se desarrolla el modelo utilizando el principio de Hamilton y la teoría de vigas Timoshenko que incluye deformaciones por corte. Se asume interpolaciones de alto orden para la aproximación de las variables fundamentales. Los materiales para emplear son isotrópicos. Se implementa un programa para estos materiales en MATLAB. Se comparan resultados con otros obtenidos en la literatura para validar el modelo. Se realiza un estudio paramétrico con una misma longitud y diferentes esbelteces. Se verifica que la formulación sea bastante precisa con resultados muy satisfactorios.
This research focuses on the free vibration analysis of Timoshkenko beams using the finite element method. The model is developed using the Hamilton principle and the Timoshenko beam theory that includes shear deformations. high order interpolations are assumed for the approximation of the fundamental variables. The materials to be used are isotropic. A program for these materials is implemented in MATLAB. Results are compared with others obtained in the literature to validate the model. A parametric study is carried out with the same length and different slenderness. It is verified that the formulation is quite precise with satisfactory results to the investigation.
Trabajo de investigación

Esta investigación se enfoca en el análisis de vibración libre de vigas Timoshenko utilizando el método de elementos finitos. Se desarrolla el modelo utilizando el principio de Hamilton y la teoría de vigas Timoshenko que incluye deformaciones por corte. Se asume interpolaciones de alto orden para la aproximación de las variables fundamentales. Los materiales para emplear son isotrópicos. Se implementa un programa para estos materiales en MATLAB. Se comparan resultados con otros obtenidos en la literatura para validar el modelo. Se realiza un estudio paramétrico con una misma longitud y diferentes esbelteces. Se verifica que la formulación sea bastante precisa con resultados muy satisfactorios.
This research focuses on the free vibration analysis of Timoshkenko beams using the finite element method. The model is developed using the Hamilton principle and the Timoshenko beam theory that includes shear deformations. high order interpolations are assumed for the approximation of the fundamental variables. The materials to be used are isotropic. A program for these materials is implemented in MATLAB. Results are compared with others obtained in the literature to validate the model. A parametric study is carried out with the same length and different slenderness. It is verified that the formulation is quite precise with satisfactory results to the investigation.
Trabajo de investigación

In wind turbines, flow pressure variations on the air-structure interface cause aerodynamic forces. Consequently the structure deforms and starts to move. The interaction between aerodynamic forces and structural deformations mainly concerns aeroelasticity. Since these two are coupled, they have to be considered simultaneously in cases which the deformations are not negligible in comparison to the other geometric dimensions. The purpose of this work is to improve the simulation model of a vertical axis wind turbine by modifying the structural model from undamped Euler-Bernoulli beam theory with lumped mass matrix to the more advanced Timoshenko beam theory with consistent mass matrix plus an additional damping term. The bending of the beam is then unified with longitudinal and torsional deformations based on a fixed shape cross-section assumption and the Saint-Venant torsion theory. The whole work has been carried out by implementing the finite element method using MATLAB code and implanting it in a previously developed package as a complement. Finally the results have been verified by qualitative comparisons with alternative simulations.

Thesis (M. S.)--Aerospace Engineering, Georgia Institute of Technology, 2009.
Committee Chair: Mavris, Dimitri; Committee Member: Bauchau, Olivier; Committee Member: Schrage, Daniel; Committee Member: Volovoi, Vitali; Committee Member: Yu, Wenbin.

Cette thèse porte sur la sécurité globale des structures en matériaux hétérogènes saturés soumis à des charges extrêmes, et est appliquée à des problèmes d’interaction fluide-structure, tels que l’interaction barrage-réservoir. Un modèle numérique d’interaction est proposé pour prédire les principales tendances et le comportement général d’un barrage en matériau saturé en interaction avec le réservoir dans des analyses de défaillance d’intérêt pratique. Le modèle numérique proposé est d’abord présenté dans un cadre bidimensionnel (2D), puis étendu à un cadre tridimensionnel (3D). La structure est considérée comme un milieu poreux saturé constitué d’un matériau cohésif. On suppose que le fluide externe en interaction avec la structure agit comme une source de saturation des pores. La réponse de la structure en matériau saturé est décrite avec un modèle lattice discrete couplé de type poutre, basé sur la discrétisation du domaine avec la tessellation de Voronoï, où les liens cohésifs sont représentés par des poutres de Timoshenko non linéaires avec un champ de déplacements enrichi en termes de discontinuités fortes. Le couplage entre la phase solide et le fluide dans les pores est traité avec la théorie de Biot et la loi de Darcy décrivant l’écoulement d’un fluide à travers d’un milieu poreux. La prise en compte numérique du couplage interne ajoute un degré de liberté supplémentaire du type pression à chaque nœud de l’élément fini de Timoshenko, qui est ensuite utilisé pour résoudre les problèmes d’interface entre la structure et le fluide. On considère que le fluide externe dans le réservoir est limité à des petits mouvements, ce qui nous permet de le modéliser avec la théorie des ondes acoustiques. Pour cela, la formulation lagrangienne avec l’approximation mixte déplacement-pression est choisie. Le traitement de l’interface fluide-structure dans le modèle numérique d’interaction est résolu d’une manière simple et efficace. Notamment, les éléments finis de la structure et du fluide externe partagent les mêmes degrés de liberté dans les nœuds communs, permettant ainsi la résolution du système d’équations avec une approche de calcul monolithique. Toutes les implémentations et les simulations numériques sont effectués avec la version recherche du code informatique FEAP (Finite Element Analysis Program). Les modèles numériques proposés pour la structure, le fluide externe et le modèle d’interaction sont validés dans le régime élastique linéaire en comparant les résultats calculés avec les valeurs de référence obtenues soit avec des solutions analytiques, soit avec des modèles continus. Les simulations numériques dans le régime non linéaire ont comme but de démontrer les capacités du modèle proposé de capturer la réponse complète à l’échelle macro et les mécanismes de rupture des structures en matériaux saturés. Enfin, la capacité du modèle d’interaction proposé de traiter la défaillance localisée progressive d’un barrage construit en matériau cohésif poreux sous l’interaction barrage-réservoir a été testé pour un programme de chargement spécifique. Pour prendre en compte les effets de la température, le couplage thermique est introduit dans le modèle numérique de la structure
This thesis studies the issue of the overall safety of structures built of heterogeneous and pore-saturated materials under extreme loads in application to fluid-structure interaction problems, such as the dam-reservoir interaction. We propose a numerical model of interaction capable of predicting main tendencies and overall behavior of pore-saturated dam structure interacting with the reservoir in failure analyses of practical interest. The proposed numerical model is first presented in two-dimensional (2D) framework and later extended to three-dimensional (3D) framework. We consider the structure built of porous cohesive material. We assume that the external fluid in interaction with the structure acts as a source of pore saturation. We model the response of the pore-saturated structure with the coupled discrete beam lattice model based on Voronoi cell representation of domain with inelastic Timoshenko beam finite elements enhanced with additional kinematics in terms of embedded strong discontinuities acting as cohesive links. The coupling between the solid phase and the pore fluid is handled with Biot’s porous media theory, and Darcy’s law governing the pore fluid flow. The numerical consideration of internal coupling results with an additional pressure-type degree of freedom placed at each node of the Timoshenko beam finite element, which is later used at the fluidstructure interface. The confined conditions met for external fluid placed in the reservoir enable the modeling of external fluid motion with the acoustic wave theory. For the numerical representation of the external fluid limited to small (irrotational) motion, we choose a Lagrangian formulation and the mixed displacement/pressure based finite element approximation. The end result are the displacement and pressure degrees of freedom per node of external fluid finite elements, which allows for the issue of the fluid-structure interface to be solved in an efficient and straightforward manner by directly connecting the structure and external fluid finite elements at common nodes. As a result, all computations can be performed in a fully monolithic manner. All numerical implementations and computations are performed with the research version of the computer code FEAP (Finite Element Analysis Program). The proposed numerical models of structure, external fluid and ultimately numerical model of interaction are validated in the linear elastic regime of structure response by comparing computed results against reference values obtained either with analytical solutions or continuum models. The numerical simulations in the nonlinear regime of structure response are performed with the aim to demonstrate the proposed coupled discrete beam lattice model capabilities to capture complete macro-scale response and failure mechanisms in pore-saturated structures. Finally, the proposed numerical model of interaction ability to deal with the progressive localized failure of a dam structure built of porous cohesive material under damreservoir interaction for a particular loading program was tested. To account for the temperature effects, the thermal coupling is introduced in the numerical model of the structure

The conventional finite element models (FEM) of problems in structural mechanics are based on the principles of virtual work and the total potential energy. In these models, the secondary variables, such as the bending moment and shear force, are post-computed and do not yield good accuracy. In addition, in the case of the Timoshenko beam theory, the element with lower-order equal interpolation of the variables suffers from shear locking. In both Euler-Bernoulli and Timoshenko beam theories, the elements based on weak form Galerkin formulation also suffer from membrane locking when applied to geometrically nonlinear problems. In order to alleviate these types of locking, often reduced integration techniques are employed. However, this technique has other disadvantages, such as hour-glass modes or spurious rigid body modes. Hence, it is desirable to develop alternative finite element models that overcome the locking problems. Least-squares finite element models are considered to be better alternatives to the weak form Galerkin finite element models and, therefore, are in this study for investigation. The basic idea behind the least-squares finite element model is to compute the residuals due to the approximation of the variables of each equation being modeled, construct integral statement of the sum of the squares of the residuals (called least-squares functional), and minimize the integral with respect to the unknown parameters (i.e., nodal values) of the approximations. The least-squares formulation helps to retain the generalized displacements and forces (or stress resultants) as independent variables, and also allows the use of equal order interpolation functions for all variables. In this thesis comparison is made between the solution accuracy of finite element models of the Euler-Bernoulli and Timoshenko beam theories based on two different least-square models with the conventional weak form Galerkin finite element models. The developed models were applied to beam problems with different boundary conditions. The solutions obtained by the least-squares finite element models found to be very accurate for generalized displacements and forces when compared with the exact solutions, and they are more accurate in predicting the forces when compared to the conventional finite element models.

碩士
國立中興大學
土木工程學系
90
In this study, a composite laminated beam with through-width delamination subjected to axial load along two clamped edges is considered. The fiber of the layers with a composite laminated beam has different angles, and an across-the-width delamination is located between layers. Shear effect and rotary inertia terms, are taken into account in the governing equations of postbuckling deformation and vibration. Based on this model, postbuckling deformations of intact and delaminated beams are found analytically. By using the perturbation method, the frequencies of vibration of postbuckling state are found. The numerical results show that the lengthwise delamination locations, delaminated length and thickness affect the postbuckling deformation and vibration frequency significantly. We found the out-of-phase modes also present in composite beam. It is almost a constant for interface 1 cases, but it decreases to small value or zero as the axial load increase for the others cases. The results of out-of-phase frequencies are different with Jane and Chen (1998) because material properties of composite make it.

博士
國立臺灣科技大學
自動化及控制研究所
99
Roll system is an essential and very important part of a lot of processing machines in textile industry such as processing calender and dye machines. Deflection and vibration, normally occurring in the processing-roll system, have always been large drawbacks of many industrial processes, especially, related to textile industry. Hence, this study introduced a new dynamic control system, including a bending actuator and a proposed controller to suppress vibration, compensate deflection. First, a new mathematic model of a combined-roll is obtained by using eigenfunction expansion method and Timoshenko beam theory. The combined-roll consists of inner steel and outer nylon layers to take full advantages and eliminate drawbacks of materials. Second, the multi-cylinder and oil-roll actuators were introduced for suppressing vibration and deflection. These solutions not only increase significant power of bending actuation system but also reduce transient time of dynamic response and avoid damaging outside surfaces of the rolls. Next, Discrete Model Predictive Controller (DMPC) was designed based on the obtained mathematical model to improve dynamic response performance. The constrained control was also combined with the proposed DMPC controller to obtain optimal control signal and use effective power of actuator system by opening wider operation range of system. Finally, in order to derive high output response performance and robust stability in the control system, the Laguerre function, one of discrete orthonormal function, and the prescribed degree stability were applied to the DMPC. The research results have shown that the response performance is completely higher and the control system is more stable.

Adhesively bonded joints have been investigated for several decades. In most analytical studies, the Bernoulli-Euler beam theory is employed to describe the behaviour of adherends. In the current work, three analytical models are developed for adhesively bonded joints using the Timoshenko beam theory for elastic material and a Bernoulli-Euler beam model for viscoelastic materials. One model is for the peeling test of an adhesively bonded joint, which is described using a Timoshenko beam on an elastic foundation. The adherend is considered as a Timoshenko beam, while the adhesive is taken to be a linearly elastic foundation. Three cases are considered: (1) only the normal stress is acting (mode I); (2) only the transverse shear stress is present (mode II); and (3) the normal and shear stresses co-exist (mode III) in the adhesive. The governing equations are derived in terms of the displacement and rotational angle of the adherend in each case. Analytical solutions are obtained for the displacements, rotational angle, and stresses. Numerical results are presented to show the trends of the displacements and rotational angle changing with geometrical and loading conditions. In the second model, the peeling test of an adhesively bonded joint is represented using a viscoelastic Bernoulli-Euler beam on an elastic foundation. The adherend is considered as a viscoelastic Bernoulli-Euler beam, while the adhesive is taken to be a linearly elastic foundation. Two cases under different stress history are considered: (1) only the normal stress is acting (mode I); and (2) only the transverse shear stress is present (mode II). The governing equations are derived in terms of the displacements. Analytical solutions are obtained for the displacements. The numerical results show that the deflection increases as time and temperature increase. The third model is developed using a symmetric composite adhesively bonded joint. The constitutive and kinematic relations of the adherends are derived based on the Timoshenko beam theory, and the governing equations are obtained for the normal and shear stresses in the adhesive layer. The numerical results are presented to reveal the normal and shear stresses in the adhesive.

碩士
國立成功大學
土木工程學系
104
On the basis of Reissner’s mixed variational theory (RMVT), a nonlocal Timoshenko beam theory (TBT) is developed for the stability analysis of a single-walled carbon nanotube (SWCNT) embedded in an elastic medium, with various boundary conditions and under axial loads. Eringen’s nonlocal elasticity theory is used to account for the small length scale effect. The strong formulations of the RMVT- based nonlocal TBT and its associated possible boundary conditions are presented. The interaction between the SWCNT and its surrounding elastic medium is simulated using the Pasternak foundation models. The critical load parameters of the embedded SWCNT with different boundary conditions are obtained using the differential quadrature (DQ) method, in which the locations of np sampling nodes are selected as the roots of np-order Chebyshev polynomials.

碩士
國立成功大學
土木工程學系
103
A nonlocal Timoshenko beam theory (TBT), based on the Reissner mixed variational theorem (RMVT), is developed for the analysis of a single-walled carbon nanotube (SWCNT) embedded in an elastic medium and with various boundary conditions. The comparisons between the results obtained by using the RMVT-based nonlocal TBT and those of principle of virtual displacement (PVD)-based one. The strong formulations of the RMVT- and PVD-based nonlocal TBTs are derived by using Hamilton’s principle, in which Eringen’s nonlocal constitutive relations are used to account for the small-scale effect. The interaction between the SWCNT and its surrounding elastic medium is simulated using the Winkler and Pasternak foundation models. The static and free vibration of the embedded SWCNT are thus investigated by using these nonlocal TBT combined with the meshless collocation methods, in which the shape functions are constructed by either the differential reproducing kernel (DRK) interpolation method or the differential quadrature (DQ) one. In the implementation of these meshless colocation methods, the results show the performance of RMVT-based nonlocal TBT is superior to that of the PVD-based one. A parametric study with regard to some crucial effects on the static and free vibration characteristics of the embedded SWCNT is undertaken, such as different boundary conditions, nonlocal parameters, aspect ratios, spring constants and shear modulus of the foundation.

Nonlinear finite element models of functionally graded beams with power-law variation of material, accounting for the von-Karman geometric nonlinearity and temperature dependent material properties as well as microstructure dependent length scale have been developed using the Euler-Bernoulli as well as the first-order and third- order beam theories. To capture the size effect, a modified couple stress theory with one length scale parameter is used. Such theories play crucial role in predicting accurate deflections of micro- and nano-beam structures. A general third order beam theory for microstructure dependent beam has been developed for functionally graded beams for the first time using a modified couple stress theory with the von Karman nonlinear strain. Finite element models of the three beam theories have been developed. The thermo-mechanical coupling as well as the bending-stretching coupling play significant role in the deflection response. Numerical results are presented to show the effect of nonlinearity, power-law index, microstructural length scale, and boundary conditions on the bending response of beams under thermo-mechanical loads. In general, the effect of microstructural parameter is to stiffen the beam, while shear deformation has the effect of modeling more realistically as a flexible beam.

Sandwich constructions are the most commonly used structures in aircraft and navy industries, traditionally. These structures are made up of the face sheets and the core, where the face sheets will be taking the load and is connected to other structural members, while the soft core material, will be used to absorb energy during impact like situation. Thus, sandwich constructions are mainly employed in light weight structures where the high energy absorption capability is required. Generally the face sheets will be thin, made up of either metallic or composite material with high stiffness and strength, while the core is light in weight, made up of soft material. Cores generally play very crucial role in achieving the desired properties of sandwich structures, either through geometric arrangement or material properties or both. Foams are in extensive use nowadays as core material due to the ease in manufacturing and their low cost. They are extensively used in automotive and industrial field applications as the desired foam density can be fabricated by adjusting the mixing, curing and heat sink processes. Modeling of sandwich beams play a crucial role in their design with suitable finite elements for face sheets and core, to ensure the compatibility between degrees of freedom at the interfaces. Unless the mathematical model simulates the physics of the model in terms of kinematics, boundary and loading conditions, results predicted will not be accurate. Accurate models helps in obtaining an efficient design of sandwich beams. In Structural Health Monitoring studies, the responses under the impact loading will be captured by carrying out the wave propagation analysis. The loads applied will be for a shorter duration (in the orders of micro seconds), where higher frequency modes will be excited. Wavelengths at such high frequencies are very small and hence, in such cases, very fine mesh generally is employed matching the wavelength requirement of the propagating wave. Traditional Finite element softwares takes enormous time and computational e ort to provide the solution. Various possible models and modeling aspects using the existing Finite element tools for wave propagation analysis are studied in the present work. There exists a huge demand for an accurate, efficient and rapidly convergent finite elements for the analysis of sandwich beams. E orts are made in the present work to address these issues and provide a solution to the sandwich user community. Super convergent and Spectral Finite sandwich Beam Elements with metallic or composite face sheets and soft core are developed. As a philosophy, the sandwich beam finite element is constructed with the combination of two beams representing the face sheets (top and bottom) at their neutral axis. The core effects are captured at the interface boundaries in terms of shear stress and normal transverse stress. In the case of wave propagation analysis, the equations are coupled in time domain and spatial domain and solving them directly is a difficult task. In Spectral Finite Element Method(SFEM), the displacement functions are derived by solving the transformed governing equations in the frequency domain. By transforming them and forces from time domain to frequency domain, the coupled partial differential equations will become coupled ordinary differential equations. These equations in frequency domain, can be solved exactly as they are normally ordinary differential equation with constant coefficients with frequency entering as a parameter. These solutions will be used as interpolating functions for spectral element formulation and in this respect it differs from conventional FE method wherein mostly polynomials are used as interpolating functions. In addition, SFEM solutions are expressed in terms of forward and backward moving waves for all the degrees of freedom involved in the formulations and hence, SFEM provides faster and efficient solutions for wave propagation analysis. In the present work, strong form of the governing differential equations are derived for a given system using Hamilton's principle. Super Convergent elements are developed by solving the static part of the governing differential equations exactly and hence the stiffness matrix derived is exact for point static loads. For wave propagation analysis, as the mass is not exactly represented, these elements are required in the optimal numbers for getting good results. The number of these elements required are generally much lesser than the number of elements required using traditional finite elements since the stiffness distribution is exact. Spectral elements are developed by solving the governing equations exactly in the frequency domain and hence the dynamic stiffness matrix derived is exact for the dynamic loads. Hence, one element between any two joints is enough to solve the whole system under impact loads for simple structures. Developing FE for sandwich beams is quiet challenging. Due to small thickness, the face sheets can be modeled using 1D idealization, while modeling of large core requires 2-D idealization. Hence, most finite or spectral elements requires stitching of these two idealizations into 1-D idealization, which can be accomplished in a variety of ways, some of which are highlighted in this thesis. Variety of finite and spectral finite elements are developed considering Euler and Timoshenko beam theories for modeling the sandwich beams. Simple element models are built with rigid core in both the theories. Models are also developed considering the flexible core with the variation of transverse displacements across depth of the core. This has direct influence on shear stress variation and also transverse normal stress in the core. Simple to higher order models are developed considering different variations in shear stress and transverse normal stress across depth of the core. Development of super convergent finite Euler Bernoulli beam elements Eul4d (4 dof element), Eul10d (10 dof element) are explained along with their results in Chapter 2. Development of different super convergent finite Timoshenko beam elements namely Tim4d (4 dof), Tim7d (7 dof), Tim10d (10 dof) are explained in Chapter 3. Validation of Euler Bernoulli and Timoshenko elements developed in the present work is carried out with test cases available in the open literature for displacements and free vibration frequencies are presented in Chapter 2 and Chapter 3. The results indicates that all developed elements are performing exceedingly well for static loads and free vibration. Super convergence performance for the elements developed is demonstrated with related examples. Spectral elements based on Timoshenko theory STim7d, STim6d, STim6dF are developed and the wave propagation characteristics studies are presented in Chapter 4. Euler spectral elements are derived from Timoshenko spectral elements by enforcing in finite shear rigidity, designated as SEul7d, SEul6d, SEul6dF and are presented. E orts were made in this present work to model the horizontal cracks in top or bottom face sheets using the spectral elements and the methodology is presented in Chapter 4. Wave propagation analysis using general purpose software N AST RAN and the super convergent as well as spectral elements developed in this work, are discussed in detail in Chapter 5. Modeling aspects of sandwich beam in N AST RAN using various combination of elements available and the performance of four possible models simulated were studied. Validation of all four models in N AST RAN, Super convergent Euler, Timoshenko and Spectral Timoshenko finite elements was carried out by simulating a homogenous I beam by comparing the longitudinal and transverse responses. Studies were carried out to find out the response predictions of a sandwich beam with soft core and all the predictions were compared and discussed. The responses in case of cracks in top or bottom face sheets under the longitudinal and transverse loading were studied in this chapter. In Chapter 6, Parametric studies were carried out for bringing out the sensitiveness of the important specific parameters in overall behaviour and performance of a sandwich beam, using Super convergent and Spectral elements developed. This chapter clearly brings out the various aspects of design of sandwich beam such as material selection of core, geometrical configuration of overall beam and core. Effects of shear modulus, mass density on wave propagation characteristics, effects of thick or thin cores with reference to the face sheets and dynamic effects of core are highlighted. Wave propagation characteristics studies includes the study of wave numbers, group speeds, cut off frequencies for a given configuration and identification of frequency zone of operations. The recommendations for improvement in design of sandwich beams based on the parametric studies are made at the end of chapter. The entire thesis, written in seven Chapters, presents a unified treatment of sandwich beam analysis that will be very useful for designers working in the area.

Rotating and non-rotating beams are widely used to model important engineering struc-tures. Hence, the vibration analyses of these beams are an important problem from a structural dynamics point of view. Depending on the beam dimensions, they are mod-eled using different beam theories. In most cases, the governing differential equations of these types of beams do not yield any simple closed-form solutions; hence we look for the inverse problem approach in determining the beam property variations given certain solutions. The long and slender beams are generally modeled using the Euler-Bernoulli beam theory. Under the premise of this theory, we study (i) the second mode tailoring of non-rotating beams having six different boundary conditions, (ii) closed-form solutions for free vibration analysis of free-free beams, (iii) closed-form solutions for free vibration analysis for gravity-loaded cantilever beams, (iv) closed-form solutions for free vibration analysis of rotating cantilever and pinned-free beams and (v) beams with shared eigen-pair. Short and thick beams are generally modeled using the Timoshenko beam theory. Here, we provide analytical closed-form solutions for the free vibration analysis of ro-tating non-homogeneous Timoshenko beams. The Rayleigh beam provides a marginal improvement over the Euler-Bernoulli beam theory without venturing into the math-ematical complexities of the Timoshenko beam theory. Under this theory, we provide closed-form solutions for the free vibration analysis of cantilever Rayleigh beams under three different axial loading conditions - uniform loading, gravity-loading and centrifu-gally loaded. We assume simple polynomial mode shapes which satisfy the different boundary conditions of a particular beam, and derive the corresponding beam property variations. In case of the shared eigenpair, we use the mode shape of a uniform beam which has a closed-form solution and use it to derive the stiffness distribution of a corresponding axially loaded beam having same length, mass variation and boundary condition. For the Timoshenko beam, we assume polynomial functions for the bending displacement and the rotation due to bending. The derived properties are demonstrated as benchmark analytical solutions for approximate and numerical methods used for the free vibration analysis of beams. They can also aid in designing actual beams for a pre-specified frequency or nodal locations in some cases. The effect of different parameters in the derived property variations and the bounds on the pre-specified frequencies and nodal locations are also studied for certain cases. The derived analytical solutions can also serve as a benchmark solution for different statistical simulation tools to find the probabilistic nature of the derived stiffness distri-bution for known probability distributions of the pre-specified frequencies. In presence of uncertainty, this flexural stiffness is treated as a spatial random field. For known probability distributions of the natural frequencies, the corresponding distribution of this field is determined analytically for the rotating cantilever Euler-Bernoulli beams. The derived analytical solutions are also used to derive the coefficient of variation of the stiffness distribution, which is further used to optimize the beam profile to maximize the allowable tolerances during manufacturing.

Facade systems are an integral part of modern day construction, especially with reinforced concrete structures. These facade systems are typically designed to withstand the effects of normal service loads and severe weather conditions. However, these elements are rarely designed to withstand the effects of an external explosion, which is considered as a rising threat to structural safety with the recent escalation of terrorist activities. In addition, these facade systems will act as the first layer of defence against an external explosion limiting the damage to the main structure. The aim of this research is to study the effects of blast loading on reinforced concrete facade systems. More specifically, the behaviour of reinforced concrete facade panels with flexible support conditions will be investigated. The overall aim was pursued by evaluating and utilising experimental studies relevant to this research to undertake 3-D finite element modelling using LS-DYNA and 1-D analytical modelling using a theoretical development. The validated numerical and analytical models were then utilised in a comprehensive parametric study.