• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1994.10a
• /
• pp.229-234
• /
• 1994

The major objective of the present paper is to estabilish finite element simulation technique in order to further analyze the shearing process. For this the ductile fracture criterion and element kill method are used in the present work. It is postulated that a fracture initiation is based on the magnitude of local effective strain. The features of sheared surface are easily observed by the element kill method. The simulation results are compard with existing experimental results. It is found that the results of the present work are close agreement with the existing results.

• Structural Engineering and Mechanics
• /
• v.37 no.6
• /
• pp.633-648
• /
• 2011

This paper investigates the aspect-ratio locking of the isoparametric 8-node quadrilateral (QUAD8) element. An important finding is that, if finite element solution is carried out with in exact arithmetic (i.e., with no truncation and round off errors), the locking tendency of the element is completely avoided even for aspect-ratios as high as 100000. The current finite element codes mostly use floating point arithmetic. Thus, they can only avoid this locking for aspect-ratios up to 100 or 1000. A novel method is proposed in the paper to avoid aspect-ratio locking in floating point computations. In this method, the offending terms of the strain-displacement matrix (i.e., $\mathbf{B}$-matrix) are multiplied by suitable scaling factors to avoid ill-conditioning of stiffness matrix. Numerical examples are presented to demonstrate the efficacy of the method. The examples reveal that aspect-ratio locking is avoided even for aspect-ratios as high as 100000.

• Structural Engineering and Mechanics
• /
• v.51 no.6
• /
• pp.923-937
• /
• 2014

The rational finite element method is different from the standard finite element method, which is constructed using basic solutions of the governing differential equations as interpolation functions in the elements. Therefore, it is superior to the isoparametric approach because of its obvious physical meaning and accuracy; it has successfully been applied to the isotropic elasticity problem. In this paper, the formulation of rational finite elements for plane orthotropic elasticity problems is deduced. This method is formulated directly in the physical domain with full consideration of the requirements of the patch test. Based on the number of element nodes and the interpolation functions, different approaches are applied with complete polynomial interpolation functions. Then, two special stiffness matrixes of elements with four and five nodes are deduced as a representative application. In addition, some typical numerical examples are considered to evaluate the performance of the elements. The numerical results demonstrate that the present method has a high level of accuracy and is an effective technique for solving plane orthotropic elasticity problems.

• Steel and Composite Structures
• /
• v.26 no.5
• /
• pp.595-606
• /
• 2018

A calculation method was presented to calculate the deflection of GFRP-concrete-steel beams with full or partial shear connections. First, the sectional analysis method was improved by considering concrete nonlinearity and shear connection stiffness variation along the beam direction. Then the equivalent slip strain was used to take into consideration of variable cross-sections. Experiments and nonlinear finite element analysis were performed to validate the calculation method. The experimental results showed the deflection of composite beams could be accurately predicted by using the theoretical model or the finite element simulation. Furthermore, more finite element models were established to verify the accuracy of the theoretical model, which included different GFRP plates and different numbers of shear connectors. The theoretical results agreed well with the numerical results. In addition, parametric studies using theoretical method were also performed to find out the effect of parameters on the deflection. Based on the parametric studies, a simplified calculation formula of GFRP-concrete-steel composite beam was exhibited. In general, the calculation method could provide a more accurate theoretical result without complex finite element simulation, and serve for the further study of continuous GFRP-concrete-steel composite beams.

• Structural Engineering and Mechanics
• /
• v.40 no.6
• /
• pp.813-827
• /
• 2011

This paper provides a numerical solution for a finite internally cracked plate using hybrid crack element method (HCE). In the formulation, an inclined crack is placed in any place of a rectangular element and the complex variable method is used. The complex potentials are expressed in a series form, and several undetermined coefficients are involved. The complex potentials for the cracked rectangle are first suggested in this paper. Based on a variational principle, the element stiffness matrix can be evaluated. The next steps are same as in the usual finite element method. Several numerical examples with computed stress intensity factor and T-stress are presented.

• Wind and Structures
• /
• v.26 no.1
• /
• pp.25-34
• /
• 2018

This paper presents a beam finite element model of a vibrate wind blade in large elastic deformation subjected to the aerodynamic, centrifugal, gyroscopic and gravity loads. The gyroscopic loads applied to the blade are induced by her simultaneous vibration and rotation. The proposed beam finite element model is based on a simplex interpolation method and it is mainly intended to the numerical analysis of wind blades vibration in large elastic deformation. For this purpose, the theory of the sheared beams and the finite element method are combined to develop the algebraic equations system governing the three-dimensional motion of blade vibration. The applicability of the theoretical approach is elucidated through an original case study. Also, the static deformation of the used wind blade is assessed by appropriate software using a solid finite element model in order to show the effectiveness of the obtained results. To simulate the nonlinear dynamic response of wind blade, the predictor-corrector Newmark scheme is applied and the stability of numerical process is approved during a large time of blade functioning. Finally, the influence of the modified geometrical stiffness on the amplitudes and frequencies of the wind blade vibration induced by the sinusoidal excitation of gravity is analyzed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.9 s.240
• /
• pp.1225-1234
• /
• 2005

Mode I stress intensity factor $(K_I)$ of Notched Ring Test(NRT) specimen for measuring slow crack growth resistance was found using finite-element method. The theoretical $K_I$ value of NRT was not available in any references and could not be solved analytically. At first, in order to verify the accuracy of the finite-element approach, published $K_I$ values of several cracks were calculated and compared with finite-element results. The results were in good agreement within inherent errors of theoretical $K_I$. Finally the mode I stress intensity factor of NRT was found using 2- and 3-dimensional finite-element methods and expressed as a function of the applied load. This enabled direct comparison of resistance to slow crack growth between NRT and Notched Pipe Test(NPT), which employ different loading regime.

• Journal of Biomedical Engineering Research
• /
• v.27 no.3
• /
• pp.110-116
• /
• 2006

Finite element method (FEM) provides several advantages over other numerical methods such as boundary element method, since it allows truly volumetric analysis and incorporation of realistic electrical conductivity values. Finite element mesh generation is the first requirement in such in FEM to represent the volumetric domain of interest with numerous finite elements accurately. However, conventional mesh generators and approaches offered by commercial packages do not generate meshes that are content-adaptive to the contents of given images. In this paper, we present software that has been implemented to generate content-adaptive finite element meshes (cMESHes) based on the contents of MR images. The software offers various computational tools for cMESH generation from multi-slice MR images. The software named as the Content-adaptive FE Mesh Generation Toolbox runs under the commercially available technical computation software called Matlab. The major routines in the toolbox include anisotropic filtering of MR images, feature map generation, content-adaptive node generation, Delaunay tessellation, and MRI segmentation for the head conductivity modeling. The presented tools should be useful to researchers who wish to generate efficient mesh models from a set of MR images. The toolbox is available upon request made to the Functional and Metabolic Imaging Center or Bio-imaging Laboratory at Kyung Hee University in Korea.

• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.149-160
• /
• 1997

Mechanical structures are composed of substructures connected by joints and boundary elements. While the finite element representation of plain substructures is well developed and reliable, joints have a lot of uncertainties in being accurately modelled and affect dynamic behavior of a total system. In order to improve the accuracy of a finite element model, a new method is proposed, in which reduced finite element model is combined with a system identification technique. After substructures except joints are modelled with finite element method and joint properties are represented by parameter states, non-linear state equation is derived in which parameter states are multiplied by physical states such as displacements and velocities. So the joint parameter identification is transformed into non-linear state estimation problem. The methods are tested and discussed numerically and the feasibility for physical application has been demonstrated through two example structures.

• Structural Engineering and Mechanics
• /
• v.33 no.3
• /
• pp.373-385
• /
• 2009

The purpose of this paper is to study shear locking-free parametric earthquake analysis of thick and thin plates using Mindlin's theory, to determine the effects of the thickness/span ratio, the aspect ratio and the boundary conditions on the linear responses of thick and thin plates subjected to earthquake excitations. In the analysis, finite element method is used for spatial integration and the Newmark-${\beta}$ method is used for the time integration. Finite element formulation of the equations of the thick plate theory is derived by using higher order displacement shape functions. A computer program using finite element method is coded in C++ to analyze the plates clamped or simply supported along all four edges. In the analysis, 17-noded finite element is used. Graphs are presented that should help engineers in the design of thick plates subjected to earthquake excitations. It is concluded that 17-noded finite element can be effectively used in the earthquake analysis of thick and thin plates. It is also concluded that, in general, the changes in the thickness/span ratio are more effective on the maximum responses considered in this study than the changes in the aspect ratio.