• Structural Engineering and Mechanics
• /
• v.64 no.6
• /
• pp.751-763
• /
• 2017

This article deals with the buckling behaviour of multilayered magneto-electro-elastic (MEE) plate subjected to uniaxial and biaxial compressive (in-plane) loads. The constitutive equations of MEE material are used to derive a finite element (FE) formulation involving the coupling between electric, magnetic and elastic fields. The displacement field corresponding to first order shear deformation theory (FSDT) has been employed. The in-plane stress distribution within the MEE plate existing due to the enacted force is considered to be equivalent to the applied in-plane compressive load in the pre-buckling range. The same stress distribution is used to derive the potential energy functional. The non-dimensional critical buckling load is accomplished from the solution of allied linear eigenvalue problem. Influence of stacking sequence, span to thickness ratio, aspect ratio, load factor and boundary condition on critical buckling load and their corresponding mode shape is investigated. In addition, static deflection of MEE plate under the sinusoidal and the uniformly distributed load has been studied for different stacking sequences and boundary conditions.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.37 no.3
• /
• pp.345-352
• /
• 2013

In this research, a study on stress shape functions was conducted to analyze the contact stress problem by using a hybrid photoelasticity. Because the contact stress problem is generally solved as a half-plane problem, the relationship between two analytical stress functions, which are compositions of the Airy stress function, was similar to one of the crack problem. However, this relationship in itself could not be used to solve the contact stress problem (especially one with singular points). Therefore, to analyze the contact stress problem more correctly, stress shape functions based on the condition of two contact end points had to be considered in the form of these two analytical stress functions. The four types of stress shape functions were related to the stress singularities at the two contact end points. Among them, the primary two types used for the analysis of an O-ring were selected, and their validities were verified in this work.

• Structural Engineering and Mechanics
• /
• v.78 no.2
• /
• pp.135-143
• /
• 2021

In this study, continuous contact problem in the functionally graded (FG) layer loaded with a FG flat punch resting on the elastic semi-infinite plane was analyzed by the finite element method (FEM). It was assumed that the shear modulus and density of the layer and punch varied according to exponentially throughout their depth. FG layer's weight was included to the problem and additionally all surfaces were considered as frictionless. Analysis of FG materials was performed with a special macro which was added to the ANSYS program. Firstly, the shear modulus of the punch was considered to be very rigid and the results of initial separation load (λcr) and distance (xcr) were compared with the analytical solution. Afterwards, results obtained from the contact analysis made according to the inhomogeneity parameters (β, γ) between FG punch-FG layer which had been unprecedented in the literature were discussed. As a result, FG punch's stress values at the punch edges where stress accumulations occurred were found to be smaller than the rigid punch. The security of the structure, longer life of the material and ease of production are directly related to the reduction of the stress values. The results obtained in this study are important in this respect. Also this work is the first study that investigates the effect of FG punch on the FG layer.

• Advances in materials Research
• /
• v.3 no.1
• /
• pp.237-257
• /
• 2014

In this paper, we adopted a two-dimensional analytical electro-elastic model to predict the stress distributions of the piezoelectric actuator in 3D case. The actuator was embedded in an elastic host structure under electrical loadings. The problem is reduced to the solution of singular integral equations of the first kind. The interfacial stresses and the axial normal stress in both plane stress state and plane strain state were obtained to study the actuation effects being transferred from the actuator to the host. The stress distributions of the PZT actuator in different length and different thickness were analyzed to guarantee the generality. The validity of the present model has been demonstrated by application of specific examples and comparisons with the corresponding results obtained from the Finite Element Method.

• Structural Engineering and Mechanics
• /
• v.20 no.2
• /
• pp.225-239
• /
• 2005

Applications of the Trefftz boundary element method (BEM) to anti-plane electroelastic problems are presented in this paper. Both direct and indirect methods with domain decomposition are discussed in details. Each crack is treated as semi-infinite thin slit defined in a subregion, for which a particular solution of the anti-plane problem, satisfying exactly the crack-face condition, is derived. The stress intensity factors defined at each crack tip can be directly computed from the coefficients of the particular solution. The performance of the proposed formulation is assessed by two examples and comparison is made with results obtained by other approaches. The Trefftz boundary element approach is demonstrated to be suitable for the analysis of the anti-plane problem of piezoelectric materials.

• Structural Engineering and Mechanics
• /
• v.25 no.4
• /
• pp.383-403
• /
• 2007

The contact problem for an elastic layer resting on an elastic half plane is considered according to the theory of elasticity with integral transformation technique. External loads P and Q are transmitted to the layer by means of two dissimilar rigid flat punches. Widths of punches are different and the thickness of the layer is h. All surfaces are frictionless and it is assumed that the layer is subjected to uniform vertical body force due to effect of gravity. The contact along the interface between elastic layer and half plane will be continuous, if the value of load factor, ${\lambda}$, is less than a critical value, ${\lambda}_{cr}$. However, if tensile tractions are not allowed on the interface, for ${\lambda}$ > ${\lambda}_{cr}$ the layer separates from the interface along a certain finite region. First the continuous contact problem is reduced to singular integral equations and solved numerically using appropriate Gauss-Chebyshev integration formulas. Initial separation loads, ${\lambda}_{cr}$, initial separation points, $x_{cr}$, are determined. Also the required distance between the punches to avoid any separation between the punches and the layer is studied and the limit distance between punches that ends interaction of punches, is investigated. Then discontinuous contact problem is formulated in terms of singular integral equations. The numerical results for initial and end points of the separation region, displacements of the region and the contact stress distribution along the interface between elastic layer and half plane is determined for various dimensionless quantities.

• Structural Engineering and Mechanics
• /
• v.75 no.3
• /
• pp.357-367
• /
• 2020

The objective of this article is investigation of dynamic response of thick multilayer functionally graded (FG) beam under generalized dynamic forces. The plane stress problem is exploited to describe the constitutive equation of thick FG beam to get realistic and accurate response. Applied dynamic forces are assumed to be sinusoidal harmonic, sinusoidal pulse or triangle in time domain and point load. Equations of motion of deep FG beam are derived based on the Hamilton principle from kinematic relations and constitutive equations of plane stress problem. The numerical finite element procedure is adopted to discretize the space domain of structure and transform partial differential equations of motion to ordinary differential equations in time domain. Numerical time integration method is used to solve the system of equations in time domain and find the time responses. Numerical parametric studies are performed to illustrate effects of force type, graduation parameter, geometrical and stacking sequence of layers on the time response of deep multilayer FG beams.

• Proceedings of the KSME Conference
• /
• 2001.06a
• /
• pp.144-150
• /
• 2001

Considered is non-symmetric contact traction induced by the tilting of a contact body and/or by a far field bulk tensile load to the other body. The problem is under the regime of plane strain. General profile of the contact end is incorporated and partial slip condition is supposed. As an example contact configuration, an indentation of a punch with rounded corners onto a half plane is studied. The variation of the internal stress field due to the tilting and the bulk tension is investigated. An edge crack problem is analyzed to examine the influence of the non-symmetric traction. It is shown that the tilting of a punch does not influence the behaviour of the crack. Rather, the effect of the bulk tension on the cracking behaviour is found considerable.

• Structural Engineering and Mechanics
• /
• v.37 no.3
• /
• pp.285-296
• /
• 2011

The plane crack-contact problem for an infinite elastic layer with two symmetric rectangular rigid stamps on its upper and lower surfaces is considered. The elastic layer having an internal crack parallel to its surfaces is subjected to two concentrated loads p on its upper and lower surfaces trough the rigid rectangular stamps and a pair of uniform compressive stress $p_0$ along the crack surface. It is assumed that the contact between the elastic layer and the rigid stamps is frictionless and the effect of the gravity force is neglected. The problem is reduced to a system of singular integral equations in which the derivative of the crack surface displacement and the contact pressures are unknown functions. The system of singular integral equations is solved numerically by making use of an appropriate Gauss-Chebyshev integration formula. Numerical results for stress-intensity factor, critical load factor, $\mathcal{Q}_c$, causing initial closure of the crack tip, the crack surface displacements and the contact stress distribution are presented and shown graphically for various dimensionless quantities.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.04a
• /
• pp.182-189
• /
• 1996

The focus of this study is on the problem of the design of structure of undetermined topology. This problem has been regarded as being the most challenging of structural optimization problems, because of the difficulty of allowing topology to change. Conventional approaches break down when element sizes approach to zero, due to stiffness matrix singularity. In this study, a novel nonlinear Programming formulation of the topology Problem is developed and examined. Its main feature is the ability to account for topology variation through zero element sizes. Stiffness matrix singularity is avoided by embedding the equilibrium equations as equality constraints in the optimization problem. Although the formulation is general, two dimensional plane elasticity examples are presented. The design problem is to find minimum weight of a plane structure of fixed geometry but variable topology, subject to constraints on stress and displacement. Variables are thicknesses of finite elements, and are permitted to assume zero sizes. The examples demonstrate that the formulation is effective for finding at least a locally minimal weight.