• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2001.05a
• /
• pp.224-227
• /
• 2001

Adaptive analysis of multilayered composite and sandwich plates is carried out. The adaptive analysis is based on a finite element error form, which measures the difference between the through-the-thickness distribution of finite element displacement and the actual displacement. The region where the error-measure exceeds the prescribed admitted error value, the finite element mesh locally refined in the thickness direction using the mesh superposition technique. Several numerical tests are conducted to validate the effectiveness of the current approach for adaptive analysis of laminated plates.

• Coupled systems mechanics
• /
• v.8 no.5
• /
• pp.393-414
• /
• 2019

This paper presents seismic analysis of concrete gravity dams considering soil-structure-fluid interaction. Displacement based plane strain finite element formulation is considered for the dam and foundation domain whereas pressure based finite element formulation is considered for the reservoir domain. A direct coupling method has been adopted to obtain the interaction effects among the dam, foundation and reservoir domain to obtain the dynamic responses of the dam. An efficient absorbing boundary condition has been implemented at the truncation surfaces of the foundation and reservoir domains. A parametric study has been carried out considering each domain separately and collectively based on natural frequencies, crest displacement and stress at the neck level of the dam body. The combined frequency of the entire coupled system is very less than that of the each individual sub-system. The crest displacement and neck level stresses of the dam shows prominent enhancement when coupling effect is taken into consideration. These outcomes suggest that a complete coupled analysis is necessary to obtain the actual responses of the concrete gravity dam. The developed methodology can easily be implemented in finite element code for analyzing the coupled problem to obtain the desired responses of the individual subdomains.

• Advances in biomechanics and applications
• /
• v.1 no.1
• /
• pp.1-14
• /
• 2014

In this paper, a simple and accurate finite element model coupled to quasi-brittle damage law able to describe the multiple cracks initiation and their progressive propagation is developed in order to predict the complete force-displacement curve and the fracture pattern of human proximal femur under quasi-static load. The motivation of this work was to propose a simple and practical FE model with a good compromise between complexity and accuracy of the simulation considering a limited number of model parameters that can predict proximal femur fracture more accurately and physically than the fracture criteria based models. Different damage laws for cortical and trabecular bone are proposed based on experimental results to describe the inelastic damage accumulation under the excessive load. When the damage parameter reaches its critical value inside an element of the mesh, its stiffness matrix is set to zero leading to the redistribution of the stress state in the vicinity of the fractured zone (crack initiation). Once a crack is initiated, the propagation direction is simulated by the propagation of the broken elements of the mesh. To illustrate the potential of the proposed approach, the left femur of a male (age 61) previously investigated by Keyak and Falkinstein, 2003 (Model B: male, age 61) was simulated till complete fracture under one-legged stance quasi-static load. The proposed finite element model leads to more realistic and precise results concerning the shape of the force-displacement curve (yielding and fracturing) and the profile of the fractured edge.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.9 no.1
• /
• pp.9-13
• /
• 2009

A simple outline for teaching adaptive scheme based finite element method for planar problems as a part of computer aided teaching of structural engineering curriculum is presented. Displacement based finite element formulation for planar problems and representative strain value based adaptive scheme for mesh generation are considered. As examples, a cantilever beam with a concentrated load treated as a planar problem and stretching of a plate with a circular hole are analyzed with displacement based finite element method with adaptive meshes. The examples and outlines show how adaptive based finite element method may become an essential part of computer aided teaching of structural engineering.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.04a
• /
• pp.17-27
• /
• 1995

For the post-buckling and elasto-plastic analysis of shell structures, the total Lagrangian formulation is presented based upon the degenerated shell element. Geometrically correct formulation is developed by updating the direction of normal vectors in the iteration process and evaluating the total Green-Lagrange stain corresponding U total displacements. In the calculation of the stiffness matrix, the element formulation takes into account the effect of finite rotation increments by retaining second order rotation terms in the incremental displacement field. The selective or reduced integration scheme using the heterosis element is applied in order to overcome both shear locking phenomena and the zero energy mode. The load/displacement incremental scheme is adopted for geometric non-linear F .E. analysis. Based on such methodology, the computer program is developed and numerical examples to demonstrate the accuracy and the effectiveness of the proposed shell element are presented and compared with references's results.

• Computational Structural Engineering
• /
• v.1 no.2
• /
• pp.83-90
• /
• 1988

The present work is concerned with the improvement of finite element for the analysis of plate bending structures. The element formulation is based upon Mindlin plate concept. The displacement field of this element is formed by adding nonconforming modes to two rotational displacement components of a 'heterosis plate element. The element has the requisite numbers of zero eigenvalues associated with rigid body modes to avoid the spurious zero energy mode. It is shown that the results obtained by the element converged to the exact solutions very rapidly as the mesh is refined and exhibited reliable solutions through numerical studies for standard benchmark problems. This element is shown to overcome the shear locking problem completely in very thin plate situation even for irregular meshes.

• Computational Structural Engineering
• /
• v.5 no.2
• /
• pp.105-118
• /
• 1992

Static performance was compared for the triangular plate elements through some numerical experiments. Four Kirchhoff elements and six Mindlin elements were selected for the comparison. Numerical tests were executed for the problems of rectangular plates with regular and distorted meshes, rhombic plates, circular plates and cantilever plates. Among the Kirchhoff 9 DOF elements, the discrete Kirchhoff theory element was the best. Element distortion and the aspect ratio were shown to have negligible effects on the displacement behaviour. The Specht's element resulted in better results than the Bergan's but it was sensitive to the aspect ratio. The element based on the hybrid stress method also resulted in good results but it assumed to be less reliable. Among the linear Mindlin elements, the discrete shear triangle was the best in view of reliability, accuracy and convergence. Since the thin plate behaviour of it was as good as the DKT element, it can be used effectively in the finite element code regardless of the thickness. As a quadratic Mindlin element, the MITC7 element resulted in best results in almost all cases considered. The results were at least as good as those of doubly refined meshes of linear elements.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.6 no.3
• /
• pp.11-21
• /
• 2002

The evaluation of displacement ductility is performed by direct method through tracking the inelastic hysteretic behavior of RC bridge columns subject to cyclic loading using a flexibility-based fiber element mode. To reasonably track the inelastic behavior until the RC bridge column reaches its ultimate state, the average stress-average strain relations and joint elements, which agree well with experiments, are modified and applied considering the tension stiffening behavior and discontinuous displacement between the column and its base. In addition the evaluation of displacement ductility is performed by a direct method easily applicable to numerical analysis. Locations for the integration points, values for the post-crushing concrete strength and low-cycle fatigue failure of longitudinal reinforcement that affect the calculation of yielding and ultimate displacements are proposed for the application to flexibility-based fiber element model. Since less than 10% of error occurs during the displacement ductility analysis, the yielding and ultimate displacements evaluated by the applied analysis method and model appear to be valid.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.20 no.6
• /
• pp.703-713
• /
• 2011

In this study, to minimize trial and error in the design and manufacturing processes of a high-precision large-surface micro-grooving machine which is able to fabricate the molds for 42 inch LCD light guide panels, the effects of the structural deformation of the micro-grooving machine according to the positions of the X-axis, Y-axis and Z-axis feed systems were examined on the tool tip displacement errors associated with the machining accuracy. The virtual prototype (finite element model) of the micro-grooving machine was constructed to include the joint stiffnesses of the hydrostatic bearings, hydrostatic guideways and linear motors, and then the tool tip displacement errors were measured from the virtual prototype. Especially, to establish the prediction model of the tool tip displacement errors, which was constructed using the positions of the X-axis, Y-axis and Z-axis feed systems as independent variables, the response surface method based on the central composite design was introduced. The reliability of the prediction model was verified by the fact that the tool tip displacement errors obtained from the prediction model coincided well those measured from the virtual prototype. And the causes of the tool tip displacement errors were identified through the analysis of interactions between the positions of the X-axis, Y-axis and Z-axis feed systems.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.2982-2994
• /
• 1993

Low-velocity impact responses of composite laminates are investigated using the finite element method based on various theories. In two-dimensional nonlinear analysis, a displacement field considering higher order shear deformation and large deflection of the laminate is assumed and a finite element formulation is developed using a C$^{o}$-continuous 9-node plate element. Also, three-dimensional linear analysis based on the infinitesimal strain-displacement assumptions is performed using 8-node brick elements with incompatible modes. A modified Hertzian contact law is incorporated into the finite element program to evaluate the impact force. In the time integration, the Newmark constant acceleration algorithm is used in conjuction with successive iterations within each time step. Numerical results from static analysis as well as the impact response analysis are presented including impact force histories, deflections, strains in the laminate. Impact responses according to two typical low-velocity impact conditions are compared each other.