• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.3
• /
• pp.11-20
• /
• 1988

A general finite element method is developed to analyze reinforced concrete plates under dead loads and monotonically increasing live loads. This method can be used to trace the load-deformation response and crack propagation through elastic, inelastic and ultimate ranges. The internal concrete and steel stresses can also be determined for any stage of the response history. A layered 8 node isoparametric element taking account of coupling effect between the membrane and the bending action is developed. An incremental tangent stiffness method is used to obtain a numerical solution. Validity of the method is studied by comparing the numerical solutions with other results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.7 s.94
• /
• pp.1719-1730
• /
• 1993

This paper presents the theoretical natural frequencies of out-of-plane deformable ring based on the variables such as out-of-plane deflection, torsional rotation and shear rotation. Based on the same variables, a finite element eigen analysis is carried out by using the $C^0$-continuous, isoparametric element which has three nodes per element and three degrees-of-freedom at each node. Numerical experiments are peformed to find the integration scheme which produces accurate natural frequencies, natural modes and correct rigid body motion. The uniformly reduced integration and the selective reduced integration give more accurate numerical frequencies than the uniformly full integration, but the uniformly reduced integration produces incorrect rigid body motion while selective reduced integration does correct one. Therefore, the ring element based on the three variables which employes selective reduced integration is recommended to avoid spurious modes, to alleviate the error due to shear locking and to produce correct rigid body motion, simultaneously.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.351-356
• /
• 2000

Various transition elements are generally used for the effective analysis of a complicated mechanical structure. In this paper, a solid-to-beam transition finite element which connects a continuum element and a $c^1-continuity$ beam element each other is proposed. The shape functions of the transition finite elements, which a 8-noded hexahedral solid element fur 3D analysis and a 4-noded quadrilateral plane element fur 2D analysis are connected to a Euler's beam element, are explicitely formulated. In order to show the effectiveness and convergence characteristics of the proposed transition elements. numerical tests are performed for various examples and their results are compared with those obtained by other methods. As the result of this study. following conclusions are obtained: (1)The proposed transition finite elements show the monotonic convergence characteristics because of having used the compatible displacement folds. (2)As being used the transition element in the finite element analysis, the finite element modelings are more convenient and the analysis results are more accurate because of the formulation characteristies of the Euler's beam element.

• Steel and Composite Structures
• /
• v.23 no.1
• /
• pp.41-51
• /
• 2017

A $C^0$ FE model developed based on an efficient higher order zigzag theory is used for hygrothermal analysis of laminated composite plates. The $C^0$ FE model satisfies the inter-laminar shear stress continuity at the interfaces and zero transverse shear stress conditions at plate top and bottom. In this model the first derivatives of transverse displacement have been treated as independent variables to circumvent the problem of $C^1$ continuity associated with the above plate theory. In the present theory the above mentioned $C^0$ continuity of the present element is compensated in the stiffness matrix formulation by using penalty parameter approach. In order to avoid stress oscillations observed in the displacement based finite element, the stress field derived from temperature/moisture fields (initial strains) must be consistent with total strain field. Special steps are introduced by field consistent approach (e.g., sampling at gauss points) to compensate this problem. A nine noded $C^0$ continuous isoparametric element is used in the proposed FE model. Comparison of present numerical results with other existing solutions shows that the proposed FE model is efficient, accurate and free of locking.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.65-73
• /
• 1989

Continuum formulations for the expressions of dynamic energy release rates and computational methods for dynamic stress intensity factors are developed for the analysis of dynamic fracture problems subjected to stress wave loading. Explicit volume integral expressions for instantaneous dynamic energy release rates are derived by modeling virtual crack extensions with the dynamic Eulerian-Lagrangian kinematic description. In the finite element applications a finite region around a crack-tip is modeled by using quarter-point singular isoparametric elements, and the volume integrals are evaluated for each crack-tip element during virtual crack extensions while the singularity is maintained. It is shown that the use of the present method is more reliable and accurate for the dynamic fracture analysis than that of other path-independent integral methods when the effects of stress waves are significant.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.5
• /
• pp.787-797
• /
• 1986

An implicit approach is employed to obtain a general boundary integral formulation of axisymmetric elastic problems in terms of a pair of singular integral equations. The corresponding kernel functions from the solutions of Navier's equation are derived by applying a three dimensional integral and a direct axisymmetrical approach. A numerical discretization schem including the evaluation of Cauchy principal values of the singular integral is described. Finally the typical axisymmetric elastic models are analyzed, i.e. the hollow sphere, the constant thickness and the V-notched round bar.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.2
• /
• pp.77-86
• /
• 1999

In the previous $paper^{(1),(2)}$, a geometrically non-linear finite element formulation of spatial cables subjected to self-weights and support motions was presented using multiple noded cable elements and how to determine the initial equililbrium state of cables was addressed. In this paper, in order to perform dynamic non-linear analysis of ocean cables subjected to support motions and earthquakes, a numerical method to calculate Morison forces and incorporate effects of earthquake motions is presented based on the Newmark method. Challenging example problems are presented in order to investigate dynamic non-linear behaviors of ocean cables subjected to support motions and earthquake loadings.

• Structural Engineering and Mechanics
• /
• v.22 no.4
• /
• pp.483-510
• /
• 2006

The dynamic instability characteristics of laminated composite stiffened shell panels subjected to in-plane harmonic edge loading are investigated in this paper. The eight-noded isoparametric degenerated shell element and a compatible three-noded curved beam element are used to model the shell panels and the stiffeners respectively. As the usual formulation of degenerated beam element is found to overestimate the torsional rigidity, an attempt has been made to reformulate it in an efficient manner. Moreover the new formulation for the beam element requires five degrees of freedom per node as that of shell element. The method of Hill's infinite determinant is applied to analyze the dynamic instability regions. Numerical results are presented to demonstrate the effects of various parameters like shell geometry, lamination scheme, stiffening scheme, static and dynamic load factors and boundary conditions, on the dynamic instability behaviour of laminated composite stiffened panels subjected to in-plane harmonic loads along the boundaries. The results of free vibration and buckling of the laminated composite stiffened curved panels are also presented.

• International Journal of Aerospace System Engineering
• /
• v.3 no.1
• /
• pp.10-16
• /
• 2016

A trigonometric Layerwise higher order shear deformation theory (TLHSDT) is developed and implemented for free vibration and buckling analysis of laminated composite and sandwich plates by analytical and finite element formulation. The present model assumes parabolic variation of out-plane stresses through the depth of the plate and also accomplish the zero transverse shear stresses over the surface of the plate. Thus a need of shear correction factor is obviated. The present zigzag model able to meet the transverse shear stress continuity and zigzag form of in-plane displacement continuity at the plate interfaces. Hence, botheration of shear correction coefficient is neglected. In the case of analytical method, the governing differential equation and boundary conditions are obtained from the principle of virtual work. For the finite element formulation, an efficient eight noded $C^0$ continuous isoparametric serendipity element is established and employed to examine the dynamic analysis. Like FSDT, the considered mathematical model possesses similar number of variables and which decides the present models computationally more effective. Several numerical predictions are carried out and results are compared with those of other existing numerical approaches.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.5
• /
• pp.1069-1078
• /
• 1990

It is shown that the performance of finite element based on energy orthogonal functions may be superior to conventional formulation for plane stress problem. Using this finite element, it is then attempted to show the distribution of stress concentration effect for subsurface under loading point. It turned out that the stress concentration effect for subsurface is not dependent on the width of the member but the loading area. And then it is shown that the solution attained by taking the stress function as a Fourier series is not satisfactory in y<0.1B.