GREEN FUNCTION MATRICES FOR ELASTICALLY RESTRAINED HETEROGENEOUS CURVED BEAMS

Size: px
Start display at page:

Download "GREEN FUNCTION MATRICES FOR ELASTICALLY RESTRAINED HETEROGENEOUS CURVED BEAMS"

Transcription

1 The th International Conference on Mechanics of Solids, Acoustics and Virations & The 6 th International Conference on "Advanced Composite Materials Engineering" ICMSAV6& COMAT6 Brasov, ROMANIA, -5 Novemer 6 GREEN FUNCTION MATRICES FOR ELASTICALLY RESTRAINED HETEROGENEOUS CURVED BEAMS György Szeidl, László Kiss Institute of Applied Mechanics, University of Miskolc, Miskolc, HUNGARY gyorgy.szeidl@uni-miskolc.hu, mechkiss@uni-miskolc.hu Astract: Virations of heterogeneous curved eams have recently een investigated,,. The thesis and papers cited deal with, among others, the issue what effect the central load has on the viration of pinned-pinned and fixed-fixed curved eams made of heterogeneous material. To find numerical solution, the authors determine the Green function matrices for pre-loaded eams. Then they reduce the eigenvalue prolems, which yield the eigenfrequencies as a function of the load, to eigenvalue prolems governed y a system of Fredholm integral equations. A similar investigation for elastically restrained heterogeneous curved eams is provided in this article. The end-restraints are modeled y linear volute springs. Keywords: Curved eams, functionally graded materials FGM, rotational restraint, natural frequency as a function of the load, Green function matrix. INTRODUCTION Curved eams are widely used in engineering applications let us consider, for instance, arch ridges or roof structures. Research concerning the mechanical ehavior of such structural elements egan in the 9 th century. The free virations of curved eams have een under extensive investigation see, e.g.:, 5 for more details. Considering the virations of pre-loaded circular eams, the numer of the availale articles is much less than for the free virations. Wasserman 6, for example, investigates the load-frequency relationship for spring supported inextensile arches. The load can e dead or follower. Here, similarly as in 7, the Galerkin method was presented as an effective way to get solutions. Chidamparam and Leissa 7 investigate the virations of pinned-pinned and fixed-fixed prestressed homogeneous circular arches under distriuted loads. The extensiility of the centerline is taken into account. We intend to contriute to the literature y improving the mechanical model, extending it for heterogeneous materials and using a numerical technique ased on the Green function matrix.. KINEMATICAL ASSUMPTIONS & GOVERNING EQUATIONS Figure. a Coordinate system D model of the curved eam We have developed a D eam model to investigate the viratory prolem. The curvilinear coordinate system ξ s, η, ζ is attached to the E weighted centerline as shown in Figure. The radius of curvature R is constant there and, moreover, the cross-section geometry and material distriution are uniform. However, the material composition and thus the material parameters can vary over the cross-sectional coordinates η, ζ as long as axis ζ is a symmetry axis oth for the geometry and for the material distriution. Therefore, it is possile to model homogeneous, functionally graded FG and even multi-layered eams, considering each material component to e linearly elastic and isotropic.

2 The centerline intersects the shown cross-section at C, whose location can e found y fulfilling the definition that the E-weighted first moment of the cross-section with respect to the axis η is zero there: Q eη Eη, ζ ζ da. A We consider the validity of the Euler-Bernoulli eam theory for the investigations, i.e., the cross-sections rotate as if they were rigid odies and remain perpendicular to the deformed centerline. Let u o, w o and ϑ e the tangential, radial displacement coordinates and the semi-vertex angle of the eam. Since the radius is constant the coordinate line s and the angle coordinate ϕ are related to each other y s Rϕ. The axial strain ε oξ and the rigid ody rotation ψ oη on the centerline can e expressed in terms of the displacements as ε oξ du o ds w o R, ψ oη u o R dw o ds. The principle of virtual work for the eam shown in Figure yields equilirium equations dn ds dm R ds N M d dm ψ oη f t, R ds ds N M ψ oη N R R f n a which should e fulfilled y the axial force N and the ending moment M. Here f t and f n denote the intensity of the distriuted loads in the tangential and normal directions f f t e t f n e n. Recalling Hooke s law 8, the relation etween the strains and inner forces ecome N I eη R ε oξ M R, A e A M I eη Eη, ζda, I eη d w o ds A w o R, N M R I eη R ε oξ, where Eη, ζζ da, m A er I eη 5 A e is the E-weighted area of the cross-section, I eη is the E-weighted moment of inertia with respect to the ending axis while m is a geometry-heterogeneity parameter the effect of the material distriution is incorporated into the model through the latter one. For simplicity reasons, we introduce dimensionless displacements and a notational convention for the derivatives taken with respect to the angle coordinate: U o u o R, W o w o R ;...n dn... dϕ n, n Z. 6 If we plug equations and into and perform some manipulations these are detailed in we get: Uo m Uo W o mε oξ W o Uo Uo R ft m m W o m ε oξ W o I eη f n. 7 Within the framework of the linear theory, we can freely neglect the effect of the deformations on the equilirium i.e., ε oξ. In the sequel, the increments which occur ecause of the viratory nature of the prolem in the typical quantities are identified y a suscript. Each physical quantity can e given in a form similar to the total tangential displacement which is equal to the sum u o u o. Here u o is the static displacement caused y the pre-load, and u o is the dynamic displacement increment. It turns out that the increments in the axial strain and in the rotation have a similar structure to equations : ε m ε oξ ψ oη ψ oη ε oξ,, ψ oη u o R dw o ds, ε oξ du o ds w o R. 8 The principle of virtual work for the increments yields the equilirium equations see for details : d N M N M ψ oη f t, ds R R R d M ds N R d ds where f t and f n are forces of inertia: N M ψ oη N M ψ oη f n, 9 R R f t ρ a A u o t, f n ρ a A w o t, 9a Thesis is downloadale from the url address

3 in which A as the area of the cross-section and ρ a as the average density on the cross-section. The increments of the inner forces can e given in terms of the displacement increments via Hooke s law: N I eη R mε oξ M d R, M w o I eη ds w o R, N M R I eη R mε oξ. a Sustituting 8 and into 9, we get the equilirium equations in the following form: Uo m Uo W o mε oξ W o m Uo Uo m m ε oξ W o W o R ft I eη f n. As regards the details we refer the reader to. For harmonic virations the amplitudes and Ŵo should satisfy the following differential equations: m m m mε oξ m ε oξ λ ; λ ρ a A R I eη α, where λ and α denote the eigenvalues and eigenfrequencies. The left side of can e rewritten as K y ϕ, ε oξ Py Py Py Py, y T Ŵo. Let rϕ e a dimensionless load: r T r r. The differential equation K y ϕ, ε oξ r descries the ehavior of the pre-loaded eam if the eam is sujected to a further load r given here in a dimensionless form. The free virations ε oξ of heterogeneous circular eams are governed y the differential equations m m m m λ For rotationally restrained eams differential equations, 5 and 6 are associated with the following oundary conditions: ϕ±ϑ ϕ±ϑ Û o,, Ŵ o ± K γŵ o 7 ϕ±ϑ where K γ k γ R/I eη is a dimensionless spring constant given in terms of k γ the spring constant. Equations, 7 and 7 determine an eigenvalue prolem. The i-th eigenfrequency α i depends on the heterogeneity parameters m and ρ a ; and also on the magnitude and the direction of the concentrated force P ζ. The effect of the latter one is accounted through the axial strain it causes: ε oξ ε oξ P, m, ϑ. Here P is a dimensionless load: P P ζ R ϑ/i eη.. NUMERICAL SOLUTION ALGORITHM The definition and the most important properties of the Green function matrix can e found in,. Solution to the inhomogeneous oundary value prolem 5, 7 is sought in the form ϑ G ϕ, ψ G yϕ Gϕ, ψrψdψ, Gϕ, ψ ϕ, ψ, 8 G ϕ, ψ G ϕ, ψ ϑ where G is the Green function matrix and ϕ, ψ are angle coordinates. If we write λy for r in 8, the eigenvalue prolem, 7 is replaced y a homogeneous integral equation system: ϑ yϕ λ Gϕ, ψyψdψ. 9 ϑ Numerical solution to such prolems can e sought e.g., y quadrature methods 9. Consider the integral formula ϑ n Jφ φψ dψ w j φψ j, ψ j ϑ, ϑ, ϑ j 5. 6

4 where ψ j ϕ is a vector and w j are the known weights. Having utilized the latter equation, we otain from 9 that n w j Gϕ, ψ j ỹψ j κỹϕ κ / λ ϑ, ϑ j is the solution, which yields an approximate eigenvalue λ / κ and the corresponding approximate eigenfunction ỹϕ. After setting ϕ to ψ i i,,,..., n, we have n w j Gψ i, ψ j ỹψ j κỹψ i κ / λ ψ i, ψ j ϑ, ϑ, or GDỸ κỹ, j where G Gψ i, ψ j for self-adjoint prolems, while D diagw,..., w... w n,..., w n and ỸT ỹ T ψ ỹ T ψ... ỹ T ψ n. After solving the generalized algeraic eigenvalue prolem, we have the approximate eigenvalues λ r and eigenvectors Y r, while the corresponding eigenfunction is otained via sustituting into : ỹ r ϕ λ r n j w j Gϕ, ψ j ỹ r ψ j r,,,..., n. Divide the interval ϑ, ϑ into equidistant suintervals of length h and apply the integration formula to each suinterval. By repeating the line of thought leading to, the algeraic eigenvalue prolem otained has the same structure as. It is also possile to consider the integral equations 9 as if they were oundary integral equations and apply isoparametric approximation on the suintervals elements. If this is the case, one can approximate the eigenfunction on the e-th element the e-th suinterval which is mapped onto the interval γ, and is denoted y L e y e y N γ e y N γ e y N γ e y, where quadratic local approximation is assumed: N i diagn i, N.5γγ, N γ, N.5γγ e, y i is the value of the eigenfunction yϕ at the left endpoint, the midpoint and the right endpoint of the element, respectively. Upon sustitution of approximation into 9, we have n ỹϕ λ e ey Gx, γn η N γ N γdγ y e e T y, 5 L e e in which, n e is the numer of elements. Using equation 5 as a point of departure, and repeating the line of thought leading to, we get again an algeraic eigenvalue prolem.. THE GREEN FUNCTION MATRICES Based on theses,, the Green function can e given in the form Gϕ, ψ Y j ϕ A j ψ ± B j ψ, 6 }{{} j where a the sign is positivenegative if ϕ ψϕ ψ; the matrices A j and B j have the following structure j j A j A A j j A j A j, Bj B B B j B j j,..., ; 7 j A j A c the coefficients in B j are independent of the oundary conditions. The columns in the matrices Y i are solutions to the homogeneous differential equations K y ϕ, ε oξ. If ε oξ < the concentrated force is compressive, then Y cos ϕ sin ϕ sin ϕ, Y cos ϕ j B j B cos χϕ Mϕ, Y χ sin χϕ However, Y and Y are different when mε oξ > : cosh χϕ Mϕ sinh χϕ Y, Y χ sinh χϕ χ cosh χϕ sin χϕ, Y χ cos χϕ. 8, M m m ε oξ. 9 Consequently we should deal with these two loading cases separately. The Green functions matrix if ε oξ <. Let us now introduce the following notational conventions a B i, B i, c B i, d B i, e B i, f B i.

5 It follows from the structure of the solutions Y l, l,..., that B i B i A i A i. The systems of equations for the unknowns a,..., f can e set up y fulfilling the continuity and discontinuity conditions the Green function matrix, should meet if ϕ ψ. Therefore, if i, we have cos ψ sin ψ cos χψ Mψ sin χψ sin ψ cos ψ χ sin χψ χ cos χψ sin ψ cos ψ χ sin χψ M χ cos χψ cos ψ sin ψ χ cos χψ χ sin χψ sin ψ cos ψ χ sin χψ χ cos χψ cos ψ sin ψ χ cos χψ χ sin χψ from where we get the constants as a c d e f m, a B χ sin ψ χ, χ cos ψ B M m χ, M m c B sin χψ χ M m χ, d B M m, e B cos χψ χ χ, f B M m M ψ m M. χ If i, then cos ψ sin ψ cos χψ Mψ sin χψ sin ψ cos ψ χ sin χψ χ cos χψ sin ψ cos ψ χ sin χψ M χ cos χψ cos ψ sin ψ χ cos χψ χ sin χψ sin ψ cos ψ χ sin χψ χ cos χψ cos ψ sin ψ χ cos χψ χ sin χψ is the equation system, the solution of which assumes the form a c d e f a B cos ψ χ, B sin ψ χ, c B cos χψ χ χ, d B, e B sin χψ χ χ, f B χ. Let α e an aritrary column matrix of size The unknown scalars A i ψ, A i ψ, A i ψ, A i ψ, A i ψ, A i ψ, i, ; ψ ϑ, ϑ in the matrices A j can e determined from the condition that the product Gϕ, ψα should satisfy the oundary conditions 7. This leads to the equation system cos ϑ sin ϑ cos χϑ Mϑ sin χϑ cos ϑ sin ϑ cos χϑ Mϑ sin χϑ sin ϑ cos ϑ χ sin χϑ χ cos χϑ sin ϑ cos ϑ χ sin χϑ χ cos χϑ sin ϑ K γ cos ϑ cos ϑ K γ sin ϑ χ sin χϑ K γχ cos χϑ χ cos χϑ K γχ sin χϑ K γ cos ϑ sin ϑ cos ϑ K γ sin ϑ K γχ cos χϑ χ sin χϑ χ cos χϑ K γχ sin χϑ a cos ϑ sin ϑ c cos χϑ dmϑ e sin χϑ f a cos ϑ sin ϑ c cos χϑ dmϑ e sin χϑ f a sin ϑ cos ϑ cχ sin χϑ d eχ cos χϑ a sin ϑ cos ϑ cχ sin χϑ d eχ cos χϑ a sin ϑ K γ cos ϑ cos ϑ K γ sin ϑ c χ sin χϑ χ cos χϑ e χ cos χϑ K γχ sin χϑ a sin ϑ K γ cos ϑ cos ϑ K γ sin ϑ c χ sin χϑ χ cos χϑ e χ cos χϑ K γχ sin χϑ The closed form solutions are presented in Appendix A. A i A i A i A i A i A i. 5 5

6 Calculation of the Green functions matrix if ε oξ > and mε oξ >. Similarly as aove, we get the following equations if i cos ψ sin ψ cosh χψ Mψ sinh χψ a sin ψ cos ψ χ sinh χψ χ cosh χψ sin ψ cos ψ χ sinh χψ M χ cosh χψ c cos ψ sin ψ χ cosh χψ χ sinh χψ d m. 6 sin ψ cos ψ χ sinh χψ χ cosh χψ e cos ψ sin ψ χ cosh χψ χ sinh χψ f The solutions are as follows χ a B χ M m c B χ χ M m e B sin ψ χ, B χ M m, d B M m, sinh χψ cosh χψ χ χ, f B M m M m Mψ. cos ψ, 7 If i cos ψ sin ψ cosh χψ Mψ sinh χψ a sin ψ cos ψ χ sinh χψ χ cosh χψ sin ψ cos ψ χ sinh χψ M χ cosh χψ c cos ψ sin ψ χ cosh χψ χ sinh χψ d sin ψ cos ψ χ sinh χψ χ cosh χψ e cos ψ sin ψ χ cosh χψ χ sinh χψ f is the equation system to e solved compare it to and the solutions we have otained are: a B cos ψ χ, B sin ψ χ, c B cosh χψ χ χ d B, e B sinh χψ χ χ, f B χ. 8 9 By repeating the steps that resulted in 5 we otain the following equation system for A i,..., A i : cos ϑ sin ϑ cosh χϑ Mϑ sinh χϑ cos ϑ sin ϑ cosh χϑ Mϑ sinh χϑ sin ϑ cos ϑ χ sinh χϑ χ cosh χϑ sin ϑ cos ϑ χ sinh χϑ χ cosh χϑ sin ϑ K γ cos ϑ cos ϑ K γ sin ϑ χ sinh χϑ K γχ cosh χϑ χ cosh χϑ K γχ sinh χϑ K γ cos ϑ sin ϑ cos ϑ K γ sin ϑ χ sinh χϑ K γχ cosh χϑ χ cosh χϑ K γχ sinh χϑ a cos ϑ sin ϑ c cosh χϑ dmϑ e sinh χϑ f a cos ϑ sin ϑ c cosh χϑ dmϑ e sinh χϑ f a sin ϑ cos ϑ cχ sinh χϑ d eχ cosh χϑ a sin ϑ cos ϑ cχ sinh χϑ d eχ cosh χϑ a sin ϑ K γ cos ϑ cos ϑ K γ sin ϑ c χ sinh χϑ K γχ cosh χϑ e χ cosh χϑ K γχ sinh χϑ a sin ϑ K γ cos ϑ cos ϑ K γ sin ϑ c χ sinh χϑ K γχ cosh χϑ e χ cosh χϑ K γχ sinh χϑ The closed form solutions are presented in Appendix B. Assume that K γ. Then the limit lim Kγ Gϕ, ψ yields the Green function matrix for pre-loaded pinned-pinned eams. When K γ, the limit lim Kγ Gϕ, ψ results in the Green function matrix for pre-loaded fixed-fixed eams. 5. COMPUTATIONAL RESULTS In this section we shall present the most important result of the computations only. Let ε oξ crit e the axial strain that elongs to the load P ζ P crit where P crit is the critical load that causes the staility loss of the eam. Further let α i e the i-th eigenfrequency of the loaded eam, while the eigenfrequencies that elong to the free virations then the eam is unloaded are denoted y α i free. Figure represents the quotient α /α free against the quotient ε oξ/ε oξ crit oth for a negative and for a positive P ζ. We remark that this time the suscript always refers to the lowest eigenfrequency for small ϑ the order of the eigenfrequencies will change. The frequencies under compression <tension> decrease <increase> almost linearly A A A A A A 6

7 and independently of m, ϑ and K γ given that m > and ϑ >. These relationships can e approximated with a very good accuracy y the linear functions α α free α α free..98 ε oξ ε oξ crit, if ε oξ <,..986 ε oξ ε oξ crit, if ε oξ >. Note that these results are almost the same as those valid for pinned-pinned or fixed-fixed curved eams. Figure. Results for the two loading cases of elastically restrained curved eams 6. CONCLUDING REMARKS We have presented a new model to clarify the viratory ehavior of circular eams pre-loaded y central load a vertical force at the crown point. The model is ased on the Euler-Bernoulli eam theory and is applicale for heterogeneous materials. The eam-end supports are rotationally restrained pins, which are modeled y linear volute springs having the same spring constant. The effect of the pre-load is incorporated into the model via the strain it causes. An eigenvalue prolem was estalished y using the principle of virtual work. This eigenvalue prolem was transformed to an eigenvalue prolem governed y Fredholm integral equations. The Green function matrix is given in a closed form oth for compressive load and for tensile load. A numerical algorithm is proposed for the solution. Though some computational results are provided, further computations are needed to refine the results and to clarify how the higher eigenfrequencies depend on the load. Acknowledgements y the second author: This research was supported y the National Research, Development and Innovation Office - NKFIH, K57 and y the Zénó Terplán Program. APPENDIX A. SOLUTIONS FOR Let us introduce the following constants: and A i,..., A i IF ε oξ < D sin χϑ sin ϑ sin χϑ K cos ϑ χ sin χϑ sin ϑ Kχ cos χϑ sin ϑ χ sin ϑ sin χϑ K γ cos ϑ sin χϑ χ sin ϑ cos χϑ, D cos ϑ sin χϑ χ sin ϑ cos χϑ Mχϑ χ cos ϑ cos χϑ K γ χ sin ϑ sin χϑ Mχϑ χ cos ϑ sin χϑ sin ϑ cos χϑ. Making use of the constants introduced the solutions sought can e given in the following forms: A i A in /D, A in χ cos ϑ sin χϑ dχ sin χϑ K γ sin ϑ sin χϑ χ cos ϑ cos χϑ d cos χϑ eχ, a A i A in /D, A in aχ cos ϑ cos χϑ amχϑ χ sin ϑ cos χϑ a sin ϑ sin χϑ cχ fχ cos χϑ K γ amχϑ χ sin ϑ sin χϑ cos ϑ cos χϑ a χ cos ϑ sin χϑ cmχ ϑ fχ sin χϑ, A i A in /χd, A in d sin ϑ eχ χ sin ϑ cos χϑ K γ eχ χ sin χϑ sin ϑ cos ϑ cos χϑ d cos ϑ, c A i A in /D, A in χ χ a cos χϑ c cos ϑ f cos ϑ cos χϑ K γ a χ sin χϑ cχ χ sin ϑ fχ sin ϑ cos χϑ χ cos ϑ sin χϑ, d 7

8 A i A in /D, A in a cmϑχ χ cos ϑ sin χϑ c χ sin ϑ sin χϑ cos χϑ cos ϑ f cos ϑ K γ amϑ cmϑχ χ cos ϑ cos χϑ sin ϑ sin χϑ c χ sin ϑ cos χϑ f sin ϑ, e A i A in /χd, A in χ χ sin χϑ dmϑχ χ sin ϑ sin χϑ dχ cos ϑ sin χϑ d sin ϑ cos χϑ eχ χ sin ϑ K γ χ cos χϑ d χ cos ϑ cos χϑ dmϑχ χ sin ϑ cos χϑ cos ϑ sin χϑ eχ χ cos ϑ. f APPENDIX B. SOLUTIONS FOR Let and A i,..., A i IF ε oξ > AND mε oξ >. D χ sin ϑ sinh χϑ K γ χ sin ϑ cosh χϑ cos ϑ sinh χϑ 5 D cos ϑ sinh χϑ χ sin ϑ cosh χϑmϑχ χ cos ϑ cosh χϑk γ χ sin ϑ sinh χϑ Mϑχ χ cos ϑ sinh χϑ sin ϑ cosh χϑ. With the two constants introduced A i A in /D, A in χ cos ϑ sinh χϑ dχ sinh χϑ K γ sin ϑ sinh χϑ χ cos ϑ cosh χϑ dχ cosh χϑ eχ, 6 7a A i A in /D, A in a sin ϑ sinh χϑ aχ cos ϑ cosh χϑ amϑχ χ sin ϑ cosh χϑ fχ cosh χϑ cχ K γ amϑχ χ sin ϑ sinh χϑ cos ϑ cosh χϑ cχ Mϑ a χ cos ϑ sinh χϑ fχ sinh χϑ, 7 A i A in /χd, A in d sin ϑ eχ χ sin ϑ cosh χϑ K γ eχ χ sin ϑ sinh χϑ cos ϑ cosh χϑ d cos ϑ, 7c A i A in /D, A in aχ χ cosh χϑ cχ χ K γ a χ cos ϑ fχ χ sinh χϑ cχ cos ϑ cosh χϑ sin ϑ fχ χ cos ϑ sinh χϑ sin ϑ cosh χϑ, 7d χ A i A in /D, A in a cχ sin ϑ sinh χϑ c cos ϑ cosh χϑ cmϑχ χ cos ϑ sinh χϑ f cos ϑ K γ amϑ c χ sin ϑ cosh χϑ cmϑχ χ cos ϑ cosh χϑ sin ϑ sinh χϑ f sin ϑ, 7e A i A in /χd, are the solutions sought. REFERENCES A in χ χ K γ χ cosh χϑ d sinh χϑdmϑχ χ χ sin ϑ sinh χϑd χ cos ϑ sinh χϑ sin ϑ cosh χϑ eχ χ sin ϑ cos ϑ cosh χϑ dmϑχ χ sin ϑ cosh χϑ cos ϑ sinh χϑ eχ χ cos ϑ. 7f L. P. Kiss. Virations and Staility of Heterogeneous Curved Beams. PhD thesis, Institute of Applied Mechanics, University of Miskolc, Hungary, 5. pp.. L. Kiss, Gy. Szeidl, S. Vlase, B. P. Gálfi, P. Dani, I. R. Munteanu, R. D. Ionescu, and J. Száva. Virations of fixed-fixed heterogeneous curved eams loaded y a central force at the crown point. International Journal for Engineering Modelling, 7-:85,. L. Kiss and Gy. Szeidl. Virations of pinned-pinned heterogeneous circular eams sujected to a radial force at the crown point. Mechanics Based Design of Structures and Machines, : 9, 5. Gy. Szeidl. The Effect of Change in Length on the Natural Frequencies and Staility of Circular Beams. PhD thesis, Department of Mechanics, University of Miskolc, Hungary, 975. in Hungarian. 5 P. Chidamparam and A. W. Leissa. Virations of planar curved eams, rings and arches. Applied Mechanis Review, ASME, 69:67 8, Y. Wasserman. The influence of the ehaviour of the load on the frequencies and critical loads of arches with flexily supported ends. Journal of Sound and Virations, 5:55 56, P. Chidamparam and A. W. Leissa. Influence of centerline extensiility on the in-plane free virations of loaded circular arches. Journal of Sound and Virations, 85: , Gy. Szeidl and L. Kiss. A Nonlinear Mechanical Model For Heterogeneous Curved Beams. In S. Vlase, editor, Proceedings of the th International Conference on Advanced Composite Materials Enginnering, COMAT, volume, pages ,. 8 - Octoer, Braşov, Romania. 9 C. T. H. Baker. The Numerical Treatment of Integral Equations Monographs on Numerical Analysis edited y L. Fox and J. Walsh. Clarendon Press, Oxford,

IN-PLANE VIBRATIONS OF PINNED-FIXED HETEROGENEOUS CURVED BEAMS UNDER A CONCENTRATED FORCE

IN-PLANE VIBRATIONS OF PINNED-FIXED HETEROGENEOUS CURVED BEAMS UNDER A CONCENTRATED FORCE 5th International Conference on "Advanced Composite Materials Engineering" COMAT 4 6-7 OCTOBER 4, Braşov, Romania IN-PLANE VIBRATIONS OF PINNED-FIXED HETEROGENEOUS CURVED BEAMS UNDER A CONCENTRATED FORCE

More information

IN-PLANE BUCKLING OF ROTATIONALLY RESTRAINED HETEROGENEOUS SHALLOW ARCHES SUBJECTED TO A CONCENTRATED FORCE AT THE CROWN POINT

IN-PLANE BUCKLING OF ROTATIONALLY RESTRAINED HETEROGENEOUS SHALLOW ARCHES SUBJECTED TO A CONCENTRATED FORCE AT THE CROWN POINT Journal of Computational and Applied Mechanics, Vol. 9., No. 2., 214, pp. 171 199 IN-PLANE BUCKLING OF ROTATIONALLY RESTRAINED HETEROGENEOUS SHALLOW ARCHES SUBJECTED TO A CONCENTRATED FORCE AT THE CROWN

More information

SIMILARITY METHODS IN ELASTO-PLASTIC BEAM BENDING

SIMILARITY METHODS IN ELASTO-PLASTIC BEAM BENDING Similarity methods in elasto-plastic eam ending XIII International Conference on Computational Plasticity Fundamentals and Applications COMPLAS XIII E Oñate, DRJ Owen, D Peric and M Chiumenti (Eds) SIMILARIT

More information

Exact Shape Functions for Timoshenko Beam Element

Exact Shape Functions for Timoshenko Beam Element IOSR Journal of Computer Engineering (IOSR-JCE) e-iss: 78-66,p-ISS: 78-877, Volume 9, Issue, Ver. IV (May - June 7), PP - www.iosrjournals.org Exact Shape Functions for Timoshenko Beam Element Sri Tudjono,

More information

Lecture 15 Strain and stress in beams

Lecture 15 Strain and stress in beams Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME

More information

Buckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads

Buckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads NASA Technical Paper 3659 Buckling Behavior of Long Symmetrically Laminated Plates Sujected to Shear and Linearly Varying Axial Edge Loads Michael P. Nemeth Langley Research Center Hampton, Virginia National

More information

In-plane free vibration analysis of combined ring-beam structural systems by wave propagation

In-plane free vibration analysis of combined ring-beam structural systems by wave propagation In-plane free viration analysis of comined ring-eam structural systems y wave propagation Benjamin Chouvion, Colin Fox, Stewart Mcwilliam, Atanas Popov To cite this version: Benjamin Chouvion, Colin Fox,

More information

Vibrations and Stability of Heterogeneous Curved Beams

Vibrations and Stability of Heterogeneous Curved Beams University of Miskolc Faculty of Mechanical Engineering and Informatics Institute of Applied Mechanics Vibrations and Stability of Heterogeneous Curved Beams Synopsis of PhD theses by László Péter Kiss

More information

David A. Pape Department of Engineering and Technology Central Michigan University Mt Pleasant, Michigan

David A. Pape Department of Engineering and Technology Central Michigan University Mt Pleasant, Michigan Session: ENG 03-091 Deflection Solutions for Edge Stiffened Plates David A. Pape Department of Engineering and Technology Central Michigan University Mt Pleasant, Michigan david.pape@cmich.edu Angela J.

More information

Chapter 5 Structural Elements: The truss & beam elements

Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations

More information

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)).

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)). Difference Equations to Differential Equations Section 8.5 Applications: Pendulums Mass-Spring Systems In this section we will investigate two applications of our work in Section 8.4. First, we will consider

More information

Symmetric Bending of Beams

Symmetric Bending of Beams Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications

More information

908. Dynamic responses of an elastic beam moving over a simple beam using modal superposition method

908. Dynamic responses of an elastic beam moving over a simple beam using modal superposition method 98. Dynamic responses of an elastic eam moving over a simple eam using modal superposition method Y. J. Wang, J. Shi, Y. Xia 3, School of Civil Engineering, Beijing Jiaotong University, Beijing 44, China,

More information

MECHANICS OF MATERIALS

MECHANICS OF MATERIALS 2009 The McGraw-Hill Companies, Inc. All rights reserved. Fifth SI Edition CHAPTER 3 MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf David F. Mazurek Torsion Lecture Notes:

More information

Homework 6: Energy methods, Implementing FEA.

Homework 6: Energy methods, Implementing FEA. EN75: Advanced Mechanics of Solids Homework 6: Energy methods, Implementing FEA. School of Engineering Brown University. The figure shows a eam with clamped ends sujected to a point force at its center.

More information

1 Nonlinear deformation

1 Nonlinear deformation NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected

More information

Mechanical Design in Optical Engineering

Mechanical Design in Optical Engineering OPTI Buckling Buckling and Stability: As we learned in the previous lectures, structures may fail in a variety of ways, depending on the materials, load and support conditions. We had two primary concerns:

More information

Comb resonator design (2)

Comb resonator design (2) Lecture 6: Comb resonator design () -Intro Intro. to Mechanics of Materials School of Electrical l Engineering i and Computer Science, Seoul National University Nano/Micro Systems & Controls Laboratory

More information

A SHEAR LOCKING-FREE BEAM FINITE ELEMENT BASED ON THE MODIFIED TIMOSHENKO BEAM THEORY

A SHEAR LOCKING-FREE BEAM FINITE ELEMENT BASED ON THE MODIFIED TIMOSHENKO BEAM THEORY Ivo Senjanović Nikola Vladimir Dae Seung Cho ISSN 333-4 eissn 849-39 A SHEAR LOCKING-FREE BEAM FINITE ELEMENT BASED ON THE MODIFIED TIMOSHENKO BEAM THEORY Summary UDC 534-6 An outline of the Timoshenko

More information

Solving Homogeneous Trees of Sturm-Liouville Equations using an Infinite Order Determinant Method

Solving Homogeneous Trees of Sturm-Liouville Equations using an Infinite Order Determinant Method Paper Civil-Comp Press, Proceedings of the Eleventh International Conference on Computational Structures Technology,.H.V. Topping, Editor), Civil-Comp Press, Stirlingshire, Scotland Solving Homogeneous

More information

Stresses in Curved Beam

Stresses in Curved Beam Stresses in Curved Beam Consider a curved beam subjected to bending moment M b as shown in the figure. The distribution of stress in curved flexural member is determined by using the following assumptions:

More information

MODELLING OF ACOUSTIC POWER RADIATION FROM MOBILE SCREW COMPRESSORS

MODELLING OF ACOUSTIC POWER RADIATION FROM MOBILE SCREW COMPRESSORS 11 th International Conference on Viration Prolems Z. Dimitrovová et.al. (eds.) Lison, Portugal, 9 12 Septemer 213 MODELLING OF ACOUSTIC POWER RADIATION FROM MOBILE SCREW COMPRESSORS Jan Dupal* 1, Jan

More information

1 Systems of Differential Equations

1 Systems of Differential Equations March, 20 7- Systems of Differential Equations Let U e an open suset of R n, I e an open interval in R and : I R n R n e a function from I R n to R n The equation ẋ = ft, x is called a first order ordinary

More information

P = ρ{ g a } + µ 2 V II. FLUID STATICS

P = ρ{ g a } + µ 2 V II. FLUID STATICS II. FLUID STATICS From a force analysis on a triangular fluid element at rest, the following three concepts are easily developed: For a continuous, hydrostatic, shear free fluid: 1. Pressure is constant

More information

2 Basic Equations in Generalized Plane Strain

2 Basic Equations in Generalized Plane Strain Boundary integral equations for plane orthotropic bodies and exterior regions G. Szeidl and J. Dudra University of Miskolc, Department of Mechanics 3515 Miskolc-Egyetemváros, Hungary Abstract Assuming

More information

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric

More information

ME 323 Examination #2

ME 323 Examination #2 ME 33 Eamination # SOUTION Novemer 14, 17 ROEM NO. 1 3 points ma. The cantilever eam D of the ending stiffness is sujected to a concentrated moment M at C. The eam is also supported y a roller at. Using

More information

Available online at ScienceDirect. Procedia IUTAM 13 (2015 ) 82 89

Available online at   ScienceDirect. Procedia IUTAM 13 (2015 ) 82 89 Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 13 (215 ) 82 89 IUTAM Symposium on Dynamical Analysis of Multibody Systems with Design Uncertainties The importance of imperfections

More information

Exact Free Vibration of Webs Moving Axially at High Speed

Exact Free Vibration of Webs Moving Axially at High Speed Eact Free Viration of Wes Moving Aially at High Speed S. HATAMI *, M. AZHARI, MM. SAADATPOUR, P. MEMARZADEH *Department of Engineering, Yasouj University, Yasouj Department of Civil Engineering, Isfahan

More information

Advanced Vibrations. Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian

Advanced Vibrations. Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian Advanced Vibrations Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Exact Solutions Relation between Discrete and Distributed

More information

Luis Manuel Santana Gallego 100 Investigation and simulation of the clock skew in modern integrated circuits. Clock Skew Model

Luis Manuel Santana Gallego 100 Investigation and simulation of the clock skew in modern integrated circuits. Clock Skew Model Luis Manuel Santana Gallego 100 Appendix 3 Clock Skew Model Xiaohong Jiang and Susumu Horiguchi [JIA-01] 1. Introduction The evolution of VLSI chips toward larger die sizes and faster clock speeds makes

More information

Riveted Joints and Linear Buckling in the Steel Load-bearing Structure

Riveted Joints and Linear Buckling in the Steel Load-bearing Structure American Journal of Mechanical Engineering, 017, Vol. 5, No. 6, 39-333 Availale online at http://pus.sciepu.com/ajme/5/6/0 Science and Education Pulishing DOI:10.1691/ajme-5-6-0 Riveted Joints and Linear

More information

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012 MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to

More information

Design Parameter Sensitivity Analysis of High-Speed Motorized Spindle Systems Considering High-Speed Effects

Design Parameter Sensitivity Analysis of High-Speed Motorized Spindle Systems Considering High-Speed Effects Proceedings of the 2007 IEEE International Conference on Mechatronics and Automation August 5-8, 2007, Harin, China Design Parameter Sensitivity Analysis of High-Speed Motorized Spindle Systems Considering

More information

Math 216 Second Midterm 28 March, 2013

Math 216 Second Midterm 28 March, 2013 Math 26 Second Midterm 28 March, 23 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that

More information

Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian

Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:

More information

BEM Formulation plane orthotropic bodies a modification for exterior regions and its proof

BEM Formulation plane orthotropic bodies a modification for exterior regions and its proof Ŕ periodica polytechnica Civil Engineering 5/ 7 3 web: http:// www.pp.bme.hu/ ci c Periodica Polytechnica 7 BEM Formulation plane orthotropic bodies a modification for exterior regions its proof RESEARCH

More information

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002 student personal identification (ID) number on each sheet. Do not write your name on any sheet. #1. A homogeneous, isotropic, linear elastic bar has rectangular cross sectional area A, modulus of elasticity

More information

Concepts of Stress and Strain

Concepts of Stress and Strain Concepts of Stress and Strain One of our principal concerns in this course is material ehavior (Strength). But strength models are often intimately related to stress. Thus, we need to e ale to compute

More information

BALKANTRIB O5 5 th INTERNATIONAL CONFERENCE ON TRIBOLOGY JUNE Kragujevac, Serbia and Montenegro

BALKANTRIB O5 5 th INTERNATIONAL CONFERENCE ON TRIBOLOGY JUNE Kragujevac, Serbia and Montenegro BALKANTRIB O5 5 th INTERNATIONAL CONFERENCE ON TRIBOLOGY JUNE5-8 5 Kragujevac, Seria and Montenegro THE TEMPERATURE DISTRIBUTION IN THE WHEEL/BRAKE SHOES CONTACT C Popescu Wagon Depot of Craiova Romania

More information

4. BEAMS: CURVED, COMPOSITE, UNSYMMETRICAL

4. BEAMS: CURVED, COMPOSITE, UNSYMMETRICAL 4. BEMS: CURVED, COMPOSITE, UNSYMMETRICL Discussions of beams in bending are usually limited to beams with at least one longitudinal plane of symmetry with the load applied in the plane of symmetry or

More information

Quadratic non-condon effect in optical spectra of impurity paramagnetic centers in dielectric crystals

Quadratic non-condon effect in optical spectra of impurity paramagnetic centers in dielectric crystals Journal of Physics: Conference Series Quadratic non-condon effect in optical spectra of impurity paramagnetic centers in dielectric crystals To cite this article: R Yu Yunusov and O V Solovyev J Phys:

More information

Dynamic response of foundations on threedimensional layered soil using the scaled boundary finite element method

Dynamic response of foundations on threedimensional layered soil using the scaled boundary finite element method IOP Conference Series: Materials Science and ngineering Dynamic response of foundations on threedimensional layered soil using the scaled oundary finite element method To cite this article: Carolin irk

More information

Advanced Structural Analysis EGF Section Properties and Bending

Advanced Structural Analysis EGF Section Properties and Bending Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear

More information

Consider an elastic spring as shown in the Fig.2.4. When the spring is slowly

Consider an elastic spring as shown in the Fig.2.4. When the spring is slowly .3 Strain Energy Consider an elastic spring as shown in the Fig..4. When the spring is slowly pulled, it deflects by a small amount u 1. When the load is removed from the spring, it goes back to the original

More information

Stress Analysis Lecture 3 ME 276 Spring Dr./ Ahmed Mohamed Nagib Elmekawy

Stress Analysis Lecture 3 ME 276 Spring Dr./ Ahmed Mohamed Nagib Elmekawy Stress Analysis Lecture 3 ME 276 Spring 2017-2018 Dr./ Ahmed Mohamed Nagib Elmekawy Axial Stress 2 Beam under the action of two tensile forces 3 Beam under the action of two tensile forces 4 Shear Stress

More information

UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION

UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume, Numer 4/0, pp. 9 95 UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION

More information

202 Index. failure, 26 field equation, 122 force, 1

202 Index. failure, 26 field equation, 122 force, 1 Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic

More information

Lecture 3: The Concept of Stress, Generalized Stresses and Equilibrium

Lecture 3: The Concept of Stress, Generalized Stresses and Equilibrium Lecture 3: The Concept of Stress, Generalized Stresses and Equilibrium 3.1 Stress Tensor We start with the presentation of simple concepts in one and two dimensions before introducing a general concept

More information

AN IMPROVED EFFECTIVE WIDTH METHOD BASED ON THE THEORY OF PLASTICITY

AN IMPROVED EFFECTIVE WIDTH METHOD BASED ON THE THEORY OF PLASTICITY Advanced Steel Construction Vol. 6, No., pp. 55-547 () 55 AN IMPROVED EFFECTIVE WIDTH METHOD BASED ON THE THEORY OF PLASTICITY Thomas Hansen *, Jesper Gath and M.P. Nielsen ALECTIA A/S, Teknikeryen 34,

More information

Comb Resonator Design (2)

Comb Resonator Design (2) Lecture 6: Comb Resonator Design () -Intro. to Mechanics of Materials Sh School of felectrical ti lengineering i and dcomputer Science, Si Seoul National University Nano/Micro Systems & Controls Laboratory

More information

The Non Linear Behavior of the Microplane Model in COMSOL

The Non Linear Behavior of the Microplane Model in COMSOL The on Linear Behavior of the Microplane Model in COMSOL A. Frigerio 1 1 RSE S.p.A. *via Ruattino, 54 20134 Milan (Italy), antonella.frigerio@rse-we.it Astract: umerical models ased on the Finite Element

More information

More Examples Of Generalized Coordinates

More Examples Of Generalized Coordinates Slides of ecture 8 Today s Class: Review Of Homework From ecture 7 Hamilton s Principle More Examples Of Generalized Coordinates Calculating Generalized Forces Via Virtual Work /3/98 /home/djsegal/unm/vibcourse/slides/ecture8.frm

More information

STATICALLY INDETERMINATE STRUCTURES

STATICALLY INDETERMINATE STRUCTURES STATICALLY INDETERMINATE STRUCTURES INTRODUCTION Generally the trusses are supported on (i) a hinged support and (ii) a roller support. The reaction components of a hinged support are two (in horizontal

More information

BENDING VIBRATIONS OF ROTATING NON-UNIFORM COMPOSITE TIMOSHENKO BEAMS WITH AN ELASTICALLY RESTRAINED ROOT

BENDING VIBRATIONS OF ROTATING NON-UNIFORM COMPOSITE TIMOSHENKO BEAMS WITH AN ELASTICALLY RESTRAINED ROOT BENDING VIBRATIONS OF ROTATING NON-UNIFORM COMPOSITE TIMOSHENKO BEAMS WITH AN EASTICAY RESTRAINED ROOT Sen Yung ee and Jin Tang Yang Department of Mechanical Engineering, National Cheng Kung University,

More information

IN-PLANE BUCKLING BEHAVIOR OF PITCHED ROOF STEEL FRAMES WITH SEMI-RIGID CONNECTIONS N. Silvestre 1, A. Mesquita 2 D. Camotim 1, L.

IN-PLANE BUCKLING BEHAVIOR OF PITCHED ROOF STEEL FRAMES WITH SEMI-RIGID CONNECTIONS N. Silvestre 1, A. Mesquita 2 D. Camotim 1, L. IN-PLANE BUCKLING BEHAVIOR OF PITCHED ROOF STEEL FRAMES WITH SEMI-RIGID CONNECTIONS N. Silvestre, A. Mesquita D. Camotim, L. Silva INTRODUCTION Traditionally, the design and analysis of steel frames assumes

More information

Simulation of inelastic cyclic buckling behavior of steel box sections

Simulation of inelastic cyclic buckling behavior of steel box sections Computers and Structures 85 (27) 446 457 www.elsevier.com/locate/compstruc Simulation of inelastic cyclic uckling ehavior of steel ox sections Murat Dicleli a, *, Anshu Mehta a Department of Engineering

More information

IN-PLANE VIBRATIONS OF CIRCULAR CURVED BEAMS WITH A TRANSVERSE OPEN CRACK 1. INTRODUCTION

IN-PLANE VIBRATIONS OF CIRCULAR CURVED BEAMS WITH A TRANSVERSE OPEN CRACK 1. INTRODUCTION Mathematical and Computational Applications, Vol. 11, No. 1, pp. 1-10, 006. Association for Scientific Research IN-PLANE VIBRATIONS OF CIRCULAR CURVED BEAMS WITH A TRANSVERSE OPEN CRACK Department of Mechanical

More information

A REFINED SLIP MODEL FOR PREDICTING THE VARIABLE BENDING STIFFNESS IN A SINGLE LAYERED CABLE

A REFINED SLIP MODEL FOR PREDICTING THE VARIABLE BENDING STIFFNESS IN A SINGLE LAYERED CABLE International Journal of Mechanical Engineering and Technology (IJMET) Volume 9, Issue 3, Decemer 8, pp. 3, Article ID: IJMET_9_3_ Availale online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=9&itype=3

More information

1.571 Structural Analysis and Control Prof. Connor Section 3: Analysis of Cable Supported Structures

1.571 Structural Analysis and Control Prof. Connor Section 3: Analysis of Cable Supported Structures .57 Structural Analysis and Control Prof. Connor Section 3: Analysis of Cale Supported Structures Many high performance structures include cales and cale systems. Analysis and design of cale systems is

More information

Slender Structures Load carrying principles

Slender Structures Load carrying principles Slender Structures Load carrying principles Basic cases: Extension, Shear, Torsion, Cable Bending (Euler) v017-1 Hans Welleman 1 Content (preliminary schedule) Basic cases Extension, shear, torsion, cable

More information

Condition number of the BEM matrix arising from the Stokes equations in 2D

Condition number of the BEM matrix arising from the Stokes equations in 2D Condition numer of the BEM matrix arising from the Stokes equations in D W. Dijkstra a, R.M.M. Mattheij Department of Mathematics and Computing Science, Eindhoven University of Technology, P.O. Box 53,

More information

AA242B: MECHANICAL VIBRATIONS

AA242B: MECHANICAL VIBRATIONS AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:

More information

M. Vable Mechanics of Materials: Chapter 5. Torsion of Shafts

M. Vable Mechanics of Materials: Chapter 5. Torsion of Shafts Torsion of Shafts Shafts are structural members with length significantly greater than the largest cross-sectional dimension used in transmitting torque from one plane to another. Learning objectives Understand

More information

Dynamic Optimization of Geneva Mechanisms

Dynamic Optimization of Geneva Mechanisms Dynamic Optimization of Geneva Mechanisms Vivek A. Suan and Marco A. Meggiolaro Department of Mechanical Engineering Massachusetts Institute of Technology Camridge, MA 09 ABSTRACT The Geneva wheel is the

More information

CHAPTER 5. Beam Theory

CHAPTER 5. Beam Theory CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions

More information

Free Vibrations of Timoshenko Beam with End Mass in the Field of Centrifugal Forces

Free Vibrations of Timoshenko Beam with End Mass in the Field of Centrifugal Forces Mechanics and Mechanical Engineering Vol. 18, No. 1 2014) 37 51 c Lodz University of Technology Free Vibrations of Timoshenko Beam with End Mass in the Field of Centrifugal Forces A. I. Manevich V. Yu.

More information

Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section

Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section Applied and Computational Mechanics 2 (2008) 397 408 Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section M. Zajíček a, a Faculty of Applied Sciences, UWB in Pilsen,

More information

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,

More information

dw 2 3(w) x 4 x 4 =L x 3 x 1 =0 x 2 El #1 El #2 El #3 Potential energy of element 3: Total potential energy Potential energy of element 1:

dw 2 3(w) x 4 x 4 =L x 3 x 1 =0 x 2 El #1 El #2 El #3 Potential energy of element 3: Total potential energy Potential energy of element 1: MAE 44 & CIV 44 Introduction to Finite Elements Reading assignment: ecture notes, ogan. Summary: Pro. Suvranu De Finite element ormulation or D elasticity using the Rayleigh-Ritz Principle Stiness matri

More information

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 2, 2011

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 2, 2011 Volume, No, 11 Copyright 1 All rights reserved Integrated Pulishing Association REVIEW ARTICLE ISSN 976 459 Analysis of free virations of VISCO elastic square plate of variale thickness with temperature

More information

The Theory of a Scroll Profile

The Theory of a Scroll Profile Purdue Universit Purdue e-pus International Compressor Engineering Conference School of Mechanical Engineering 006 The Theor of a Scroll Profile Huert Bukac Little Dnamics Follow this and additional works

More information

Modeling of Water Flows around a Circular Cylinder with the SPH Method

Modeling of Water Flows around a Circular Cylinder with the SPH Method Archives of Hydro-Engineering and Environmental Mechanics Vol. 61 (2014), No. 1 2, pp. 39 60 DOI: 10.1515/heem-2015-0003 IBW PAN, ISSN 1231 3726 Modeling of Water Flows around a Circular Cylinder with

More information

On the asymptotic derivation of Winkler-type energies from 3D elasticity

On the asymptotic derivation of Winkler-type energies from 3D elasticity On the asymptotic derivation of Winkler-type energies from 3D elasticity Andrés A. León Baldelli 1 and B. Bourdin 1,2 1 Center for Computation & Technology, Louisiana State University, Baton Rouge LA 70803,

More information

CHAPTER -6- BENDING Part -1-

CHAPTER -6- BENDING Part -1- Ishik University / Sulaimani Civil Engineering Department Mechanics of Materials CE 211 CHAPTER -6- BENDING Part -1-1 CHAPTER -6- Bending Outlines of this chapter: 6.1. Chapter Objectives 6.2. Shear and

More information

Mechanical Design in Optical Engineering

Mechanical Design in Optical Engineering Torsion Torsion: Torsion refers to the twisting of a structural member that is loaded by couples (torque) that produce rotation about the member s longitudinal axis. In other words, the member is loaded

More information

An internal crack problem in an infinite transversely isotropic elastic layer

An internal crack problem in an infinite transversely isotropic elastic layer Int. J. Adv. Appl. Math. and Mech. 3() (05) 6 70 (ISSN: 347-59) Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics An internal crack prolem in an infinite

More information

A consistent dynamic finite element formulation for a pipe using Euler parameters

A consistent dynamic finite element formulation for a pipe using Euler parameters 111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,

More information

Mathematical model of static deformation of micropolar elastic circular thin bar

Mathematical model of static deformation of micropolar elastic circular thin bar Mathematical model of static deformation of micropolar elastic circular thin bar Mathematical model of static deformation of micropolar elastic circular thin bar Samvel H. Sargsyan, Meline V. Khachatryan

More information

COPYRIGHTED MATERIAL. Index

COPYRIGHTED MATERIAL. Index Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,

More information

RESEARCH REPORT IN MECHANICS

RESEARCH REPORT IN MECHANICS No. XX ISSN XXXX XXXX Oct. 8 Semi-analytical postuckling analysis of stiffened plates with a free edge y ars Bruak RESEARCH REPORT IN MECHANICS UNIVERSITY OF OSO DEPARTMENT OF MATHEMATICS MECHANICS DIVISION

More information

3. BEAMS: STRAIN, STRESS, DEFLECTIONS

3. BEAMS: STRAIN, STRESS, DEFLECTIONS 3. BEAMS: STRAIN, STRESS, DEFLECTIONS The beam, or flexural member, is frequently encountered in structures and machines, and its elementary stress analysis constitutes one of the more interesting facets

More information

LINEAR AND NONLINEAR SHELL THEORY. Contents

LINEAR AND NONLINEAR SHELL THEORY. Contents LINEAR AND NONLINEAR SHELL THEORY Contents Strain-displacement relations for nonlinear shell theory Approximate strain-displacement relations: Linear theory Small strain theory Small strains & moderate

More information

LATERAL STABILITY OF DEEP BEAMS WITH SHEAR-BEAM SUPPORT

LATERAL STABILITY OF DEEP BEAMS WITH SHEAR-BEAM SUPPORT U. FOREST SERVICE RESEARCH PAPER FPL 43 OCTOBER U. S. DEPARTMENT OF AGRICULTURE FOREST SERVICE FOREST PRODUCTS LABORATORY MADISON, WIS. LATERAL STABILITY OF DEEP BEAMS WITH SHEAR-BEAM SUPPORT The FOREST

More information

726. The analysis of a composite beam with piezoelectric actuator based on the approximate method

726. The analysis of a composite beam with piezoelectric actuator based on the approximate method 76. The analysis of a comosite eam with iezoelectric actuator ased on the aroximate method Andrzej Buchacz, Marek Płaczek The Silesian University of Technology, Faculty of Mechanical ngineering Institute

More information

The Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density

The Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density Applied Mathematics & Information Sciences 23 2008, 237 257 An International Journal c 2008 Dixie W Publishing Corporation, U. S. A. The Rotating Inhomogeneous Elastic Cylinders of Variable-Thickness and

More information

Appendix C. Modal Analysis of a Uniform Cantilever with a Tip Mass. C.1 Transverse Vibrations. Boundary-Value Problem

Appendix C. Modal Analysis of a Uniform Cantilever with a Tip Mass. C.1 Transverse Vibrations. Boundary-Value Problem Appendix C Modal Analysis of a Uniform Cantilever with a Tip Mass C.1 Transverse Vibrations The following analytical modal analysis is given for the linear transverse vibrations of an undamped Euler Bernoulli

More information

6. Bending CHAPTER OBJECTIVES

6. Bending CHAPTER OBJECTIVES CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where

More information

ENG2000 Chapter 7 Beams. ENG2000: R.I. Hornsey Beam: 1

ENG2000 Chapter 7 Beams. ENG2000: R.I. Hornsey Beam: 1 ENG2000 Chapter 7 Beams ENG2000: R.I. Hornsey Beam: 1 Overview In this chapter, we consider the stresses and moments present in loaded beams shear stress and bending moment diagrams We will also look at

More information

Handout 6: Rotational motion and moment of inertia. Angular velocity and angular acceleration

Handout 6: Rotational motion and moment of inertia. Angular velocity and angular acceleration 1 Handout 6: Rotational motion and moment of inertia Angular velocity and angular acceleration In Figure 1, a particle b is rotating about an axis along a circular path with radius r. The radius sweeps

More information

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is

More information

Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.

Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.

More information

Analytical Research of a Tension the Warp Yarns for the Cycle of Work of the Weaving Loom

Analytical Research of a Tension the Warp Yarns for the Cycle of Work of the Weaving Loom ISSN: 350-038 International ournal of AdvancedResearch in Science, Vol. 5, Issue 0, Octoer 08 Analytical Research of a Tension the Warp Yarns for the Cycle of Work of the Weaving Loom RasulovKh.Yu, Rakhimkhodjaev

More information

arxiv:math/ v1 [math.sp] 27 Mar 2006

arxiv:math/ v1 [math.sp] 27 Mar 2006 SHARP BOUNDS FOR EIGENVALUES OF TRIANGLES arxiv:math/0603630v1 [math.sp] 27 Mar 2006 BART LOMIEJ SIUDEJA Astract. We prove that the first eigenvalue of the Dirichlet Laplacian for a triangle in the plane

More information

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the

More information

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e

More information

Basic Energy Principles in Stiffness Analysis

Basic Energy Principles in Stiffness Analysis Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting

More information

Chapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx

Chapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx Chapter 1 Lecture Notes Chapter 1 Oscillatory Motion Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx When the mass is released, the spring will pull

More information

A Continuous Vibration Theory for Beams with a Vertical Edge Crack

A Continuous Vibration Theory for Beams with a Vertical Edge Crack Transaction B: Mechanical Engineering Vol. 17, No. 3, pp. 194{4 c Sharif University of Technology, June 1 A Continuous Viration Theory for Beams with a Vertical Edge Crac Astract. M. Behzad 1;, A. Erahimi

More information

Introduction to Finite Element Method. Dr. Aamer Haque

Introduction to Finite Element Method. Dr. Aamer Haque Introduction to Finite Element Method 4 th Order Beam Equation Dr. Aamer Haque http://math.iit.edu/~ahaque6 ahaque7@iit.edu Illinois Institute of Technology July 1, 009 Outline Euler-Bernoulli Beams Assumptions

More information