• Proceedings of the Korea Concrete Institute Conference
• /
• 2004.11a
• /
• pp.357-360
• /
• 2004

Estimation of deformation capacity of non-flexural reinforced concrete members is proposed using basic concepts of limit analysis and the virtual work method. This new approach starts with construction of admissible stress field as for an equilibrium set. Failure mechanisms compatible with admissible stress fields are postulated as for displacement set. It is assumed that the ultimate deformations as result of failure mechanisms are controlled by ultimate strain of concrete in compression. The derived formula for deformability of deep beams in shear shows reasonable range of ultimate displacement.

• Structural Engineering and Mechanics
• /
• v.70 no.5
• /
• pp.551-562
• /
• 2019

In the present study, a suitable mathematical model considering parabolic transverse shear strains for dynamic analysis of laminated composite skew plates under different types of impulse and spatial loads was presented for the first time. The proposed mathematical model satisfies zero transverse shear strain at the top and bottom of the plate. On the basis of the cubic variation of thickness coordinate in in-plane displacement fields of the present mathematical model, a 2D finite element (FE) model was developed including skew transformations in the mathematical model. No shear correction factor is required in the present formulation and damping effect was also incorporated. This is the first FE implementation considering a cubic variation of thickness coordinate in in-plane displacement fields including skew transformations to solve the forced vibration problem of composite skew plates. The effect of transverse shear and rotary inertia was incorporated in the present model. The Newmark-${\beta}$ scheme was adapted to perform time integration from step to step. The $C^0$ FE formulation was implemented to overcome the problem of $C^1$ continuity associated with the cubic variation of thickness coordinate in in-plane displacement fields. The numerical studies showed that the present 2D FE model predicts the result close to the analytical results. Many new results varying different parameter such as skew angles, boundary conditions, etc. were presented.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2004.10a
• /
• pp.855-859
• /
• 2004

We develop a new four-node flat shell element which is accurate, efficient, and suitable to be used on general purpose. The new element has a hybrid Trefftz element with drilling degrees of freedom as a membrane part. We define the two independent displacement field: the internal displacement field that satisfies governing equations in the domain a priori and the boundary displacement field that is usually used as a conventional finite element method. The hybrid Trefftz variational formulation connects these two displacement fields on the boundary of the domain. To add drilling degrees of freedom, we introduce the Allman's quadratic displacement field to the boundary displacement field. As a result, our flat shell element has 6 degrees of freedom per a node. We also use the well-known DKMQ plate bending element for the plate part of the proposed element. The DKMQ element satisfies Mindlin-Reissner‘s plate theory along the edge of the element and gives proper behavior regardless of the thickness. A series of numerical experiments shows that the performance of the new element such as accuracy, rate of convergence, robustness to mesh quality, and so on.

• Structural Engineering and Mechanics
• /
• v.31 no.4
• /
• pp.439-451
• /
• 2009

The paper deals with the problem related to the modelling of riveted assemblies for crashworthiness analysis of full-scale complete aircraft structures. Comparisons between experiments and standard FE computations on high-energy accidental situations onto aluminium riveted panels show that macroscopic plastic strains are not sufficiently localised in the FE shells connected to rivet elements. The main reason is related to the structural embrittlement caused by holes, which are currently not modelled. Consequently, standard displacement FE models do not succeed in initialising and propagating the rupture in sheet metal plates and along rivet rows as observed in the experiments. However, the literature survey show that it is possible to formulate super-elements featuring defects that both give accurate singular strain fields and are compatible with standard displacement finite elements. These super-elements can be related to the displacement model of the hybrid-Trefftz principle of the finite element method, which is a kind of domain decomposition method. A feature of hybrid-Trefftz finite elements is that they are mainly used for elastic computations. It is thus proposed to investigate the possibility of formulating a hybrid displacement finite element, including the effects of a hole, dedicated to crashworthiness analysis of full-scale aeronautic structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.5
• /
• pp.1947-1955
• /
• 2013

The purpose of this paper is to present an photoelasticity technique for measuring the displacement and stress distribution in particle assembly subjected to shallow foundation loading. Photoelastic measurement technique was employed to visualize the force transmission of a particle assembly. A model assembly bounded by a steel frame was built by stacking bi-dimensional circular particles made of polycarbonate elastomer. Each particle was coated by a thin photoelastic sheet so that the force transmission represented by bright light stripes can be visualized. In a contacted particle, both magnitude and orientation of principal stress difference can also be measured via the photoelasticity technique. The different distributions of the contact stresses at the initial loading and near the failure were quantitatively compared. The photoelastic patterns and displacement fields observed in the pre-failure state disappears immediately after the buckling of confined force chains.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.11
• /
• pp.1844-1851
• /
• 1997

Displacement fields and interface stresses are obtained by modifying the potential energy functional with a penalty function which enforces the continuity of stresses at the interface of two-materials. Based on the displacement field and the interface stresses, a new methodology to generate a continuous stress field over the entire domain including the interface of the dissimilar materials has been proposed by combining the L$^{2}$ projection method of stress-smoothing and the Loubignac's iterative method of improving the displacement field. Stress analysis was carried out on two examples which are made of highly dissimilar materials. As a result of the analysis, it is found that the proposed method provides improved continuity of the stress field over the entire domain as well as predicting accurate nodal stresses at the interface. In contrast, the conventional displacement-based finite element method provides significant stress discontinuties at the interfaces. In addition, it was found that the total strain energy evaluated from the improved continuous stress field converge to the exact value as increasing the number of iterations in the proposed method.

• Structural Engineering and Mechanics
• /
• v.56 no.1
• /
• pp.137-156
• /
• 2015

In this paper, an analytical solution of displacement, strain and stress field for rotating thick-walled cylinder made of functionally graded material subjected to the uniform external magnetic field and thermal field in plane strain state has been studied. Stress, strain and displacement field as a function of radial coordinates considering magneto-thermo-elasticity are derived analytically. According to the Maxwell electro-dynamic equations, Lorentz force in term of displacement is obtained in cylindrical coordinates. Also, symmetric temperature distribution along the thickness of hollow cylinder is obtained by solving Fourier heat transfer equation in cylindrical coordinates. Using equation of equilibrium and thermo-mechanical constitutive equations associated with Lorentz force, a second-order inhomogeneous differential equation in term of displacement is obtained and will be solved analytically. Except Poisson's ratio, other mechanical properties such as elasticity modulus, density, magnetic permeability coefficient, heat conduction coefficient and thermal expansion coefficient are assumed to vary through the thickness according to a power law. In results analysis, non-homogeneity parameter has been chosen arbitrary and inner and outer surface of cylinder are assumed to be rich metal and rich ceramic, respectively. The effect of rotation, thermal, magnetic field and non-homogeneity parameter of functionally graded material which indicates percentages of cylinder's constituents are studied on displacement, Von Mises equivalent stress and Von Mises equivalent strain fields.

• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.179-188
• /
• 1997

A 4-node element for efficient finite element analysis of plate bending is presented in this paper. This element is formulated based on Mindlin plate theory to take account of shear deformation. To overcome the overestimation of shear stiffness in thin Mindlin plate element, especially in the lower order element, five nonconforming displacement modes are added to the original displacement fields. The proposed nonconforming element does not possess spurious zero-energy mode and does not show shear locking phenomena in very thin plate even for distorted mesh shapes. It was recognized from benchmark numerical tests that the displacement converges to the analytical solutions rapidly and the stress distributions are very smooth. The element also provides good results for the case of high aspect ratio.

• Structural Engineering and Mechanics
• /
• v.51 no.1
• /
• pp.163-180
• /
• 2014

In the displacement based finite element analysis of composite beams that consist of two Euler-Bernoulli beams juxtaposed with a deformable shear connection, the coupling of the displacement fields may cause oscillations in the interlayer slip field and reduction in optimal convergence rate, known as slip-locking. In this study, the B-bar procedure is proposed to alleviate the locking effects. It is also shown that by changing the primary dependent variables in the mathematical model, to be able to interpolate the interlayer slip field directly, oscillations in the slip field can be completely eliminated. Examples are presented to illustrate the performance and the numerical characteristics of the proposed methods.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2003.04a
• /
• pp.143-146
• /
• 2003

A finite element based on the efficient higher order zig-zag theory with multiple delaminations Is developed to refine the predictions of frequency and mode shapes. Displacement field through the thickness are constructed by superimposing linear zig-zag field to the smooth globally cubic varying field. The layer-dependent degrees of freedom of displacement fields are expressed in terms of reference primary degrees of freedom by applying interface continuity conditions including delaminated interfaces as well as free hounding surface conditions of transverse shear stresses. Thus the proposed theory is not only accurate but also efficient. This displacement field can systematically handle the number, shape, size, and locations of delaminations. Throught the dynamic version of variational approach, the dynamic equilibrium equations and variationally consistent boundary conditions are obtained. Through the natural frequency analysis and time response analysis of composite plate with multiple delaminations, the accuracy and efficiency of the present finite element are demonstrated. The present finite element is suitable in the predictions of the dynamic response of the thick composite plate with multiple delaminations.