• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.2
• /
• pp.354-362
• /
• 1995

The finite element modeling is used to study the buckling and postbuckling behavior of composite laminates with an embedded delamination. Degenerated shell element and rigid beam element are utilized for the finite element modeling. In the nonlinear finite element formulation, the updated Lagrangian description method based on the second Piola-Kirchhoff stress tensor and the Green strain tensor is used. The buckling and postbuckling behavior of composite laminates with a delamination are investigated for various delamination sizes, stacking sequences, and boundary conditions. It is shown that the buckling load and postbuckling behavior of composite laminates depend on the buckling model which is determined by the delamination size, stacking sequence and boundary condition. Also, results show that introduction of couplings can reduce greatly the buckling load.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.4
• /
• pp.703-713
• /
• 2001

Mechanical joints such as bolted or riveted joints are widely used in structural components and the reliable determination of the stress intensity factors for corner cracks in mechanical joints is needed to evaluate the safety and fatigue life. This paper analyzes the mixed-mode stress intensity factors of surface and deepest points for quarter elliptical corner cracks in mechanical joints by weight function method and the coefficients included in weight function are determined by finite element analyses for reference loadings. The extended form of the weight function method for two-dimensional mixed-mode to three-dimensional is presented and the number of terms in weight function is determined by comparing the results for the different number of terms. The amount of clearance is an important factor in evaluating the severity of elliptical corner cracks in mechanical joints and even horizontal crack normal to the applied load is under mixed-mode in the case that clearance exists.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.212-217
• /
• 2000

The reliable determination of the stress intensity factors for cracks in bolted Joints is needed to evaluate the safety and fatigue life of them widely used in mechanical components. The weight function method is an efficient technique to calculate the stress intensity factors for various loading conditions using the stresses of an uncracked model. In this paper the mixed-mode stress intensity factors for cracks in bolted joints are obtained by weight function method, in which the coefficients of weight function are determined by finite element analyses far reference loadings. The effects of the magnitude of clearance and factional coefficient on the stress intensity factors are investigated.

• Structural Engineering and Mechanics
• /
• v.58 no.2
• /
• pp.217-230
• /
• 2016

Effects of temperature dependent material properties on mixed mode fracture parameters of functionally graded materials subjected to thermal loading are investigated. A domain form of the $J_k$-integral method including temperature-dependent material properties and its numerical implementation using finite element analysis is presented. Temperature and displacement fields are calculated using finite element analysis and are used to compute mixed mode stress intensity factors using the $J_k$-integral. Numerical results indicate that temperature-dependency of material properties has considerable effect on the mixed-mode stress intensity factors of cracked functionally graded structures.

• Journal of Mechanical Science and Technology
• /
• v.19 no.11
• /
• pp.2016-2024
• /
• 2005

In this work, a mixed beam approach that combines both the stiffness and the flexibility methods has been performed to analyze the coupled composite blades with closed, two-cell cross-sections. The Reissner's semi-complementary energy functional is used to derive the beam force-displacement relations. Only the membrane part of the shell wall is taken into account to make the analysis simple and also to deliver a clear picture of the mixed method. All the cross section stiffness coefficients as well as the distribution of shear across the section are evaluated in a closed-form through the beam formulation. The theory is validated against experimental test data, detailed finite element analysis results, and other analytical results for coupled composite blades with a two-cell airfoil section. Despite the simple kinematic model adopted in the theory, an accuracy comparable to that of two-dimensional finite element analysis has been obtained for cases considered in this study.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.5
• /
• pp.120-127
• /
• 1998

Differential inequalities occurring in problems of obstacle contact problems are recast into variational inequalities and analyzed by finite element methods. A new a posteriori error estimator, which is essential in adaptive finite element method, is introduced to capture the errors in finite element approximations of these variational inequalities. In order to construct a posteriori error estimates, saddle point problems are introduced using Lagrange parameters and upper bounds are provided. The global upper bound is localized by a special mixed formulation, which leads to upper bounds of the element errors. A numerical experiment is performed on an obstacle contact problem to check the effectivity index both in a local and a global sense.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.51-59
• /
• 2000

The Interactions of shock wave with turbulent boundary layers in high-speed flows cause complex flowfields which result in increased adverse pressure gradients, skin friction and temperature. Accurate and reliable prediction of such phenomena is needed in designing high-speed propulsion systems. Such analyses of the complex flowfields require sophisticated numerical scheme that can resolve interactions between shock wave and boundary layers accurately. Therefore the purpose of the present. article is to introduce an accurate and efficient mixed explicit-implicit generalized Galerkin finite element method. To demonstrate the validity of the theory and numerical procedure, several benchmark cases are investigated.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.15 no.1
• /
• pp.57-65
• /
• 2011

In this paper, a good quality mesh generation for the finite element method is investigated for artificial hip joint simulations. In general, bad meshes with a large aspect ratio or mixed elements can give rise to excessively long computational running times and extremely high errors. Typically, hexahedral elements outperform tetrahedral elements during three-dimensional contact analysis using the finite element method. Therefore, it is essential to mesh biologic structures with hexahedral elements. Four meshing schemes for the finite element analysis of an artificial hip joint are presented and compared: (1) tetrahedral elements, (2) wedge and hexahedral elements, (3) open cubic box hexahedral elements, and (4) proposed hexahedral elements. The proposed meshing scheme is to partition a part before seeding so that we have a high quality three-dimensional mesh which consists of only hexahedral elements. The von Mises stress distributions were obtained and analyzed. We also performed mesh refinement convergence tests for all four cases.

• Computational Structural Engineering
• /
• v.9 no.2
• /
• pp.133-142
• /
• 1996

In this paper, two different techniques for mixed-mode type engineering fracture mechanics are investigated to estimate the stress intensity factors by using p-version finite element model. These two techniques are displacement extrapolation with COD and CSD method and J-integral with decomposition method. By decomposing the displacement field obtained from p-version of finite element analysis into symmetric and antisymmetric displacement fields with respect to the crack line, Mode-I and Mode-II stress intensity factors can be determined using aforementioned techniques. The example problems for validating the proposed techniques are centrally and centrally oblique cracked panels under tension. The numerical results associated with the variation of oblique angle and the ratio of crack length and panel width (a /W ratio) are compared with those by theoretical values and empirical solutions in literatures. Very good agreements with the existing solutions are shown.

• Structural Engineering and Mechanics
• /
• v.85 no.2
• /
• pp.147-161
• /
• 2023

Based on the ideas of variational differential quadrature (VDQ) and finite element method (FEM), a numerical approach named as VDQFEM is applied herein to study the large deformations of plate-type structures under static loading with arbitrary shape hole made of functionally graded graphene platelet-reinforced composite (FG-GPLRC) in the context of higher-order shear deformation theory (HSDT). The material properties of composite are approximated based upon the modified Halpin-Tsai model and rule of mixture. Furthermore, various FG distribution patterns are considered along the thickness direction of plate for GPLs. Using novel vector/matrix relations, the governing equations are derived through a variational approach. The matricized formulation can be efficiently employed in the coding process of numerical methods. In VDQFEM, the space domain of structure is first transformed into a number of finite elements. Then, the VDQ discretization technique is implemented within each element. As the last step, the assemblage procedure is performed to derive the set of governing equations which is solved via the pseudo arc-length continuation algorithm. Also, since HSDT is used herein, the mixed formulation approach is proposed to accommodate the continuity of first-order derivatives on the common boundaries of elements. Rectangular and circular plates under various boundary conditions with circular/rectangular/elliptical cutout are selected to generate the numerical results. In the numerical examples, the effects of geometrical properties and reinforcement with GPL on the nonlinear maximum deflection-transverse load amplitude curve are studied.