• Structural Engineering and Mechanics
• /
• v.68 no.3
• /
• pp.299-312
• /
• 2018

In this paper, an extended finite element method is proposed to analyze both geometric and material non-linear behavior of general Functionally Graded Material (FGM) plate-shell type structures. A user defined subroutine (UMAT) is developed and implemented in Abaqus/Standard to study the elastoplastic behavior of the ceramic particle-reinforced metal-matrix FGM plates-shells. The standard quadrilateral 4-nodes shell element with three rotational and three translational degrees of freedom per node, S4, is extended in the present study, to deal with elasto-plastic analysis of geometrically non-linear FGM plate-shell structures. The elastoplastic material properties are assumed to vary smoothly through the thickness of the plate-shell type structures. The nonlinear approach is based on Mori-Tanaka model to underline micromechanics and locally determine the effective FGM properties and self-consistent method of Suquet for the homogenization of the stress-field. The elasto-plastic behavior of the ceramic/metal FGM is assumed to follow Ludwik hardening law. An incremental formulation of the elasto-plastic constitutive relation is developed to predict the tangent operator. In order to to highlight the effectiveness and the accuracy of the present finite element procedure, numerical examples of geometrically non-linear elastoplastic functionally graded plates and shells are presented. The effects of the geometrical parameters and the volume fraction index on nonlinear responses are performed.

• Geomechanics and Engineering
• /
• v.6 no.6
• /
• pp.545-560
• /
• 2014

For slope stability analysis, an alternative to the classical limit equilibrium method (LEM) of slices is the shear strength reduction method (SRM), which can be integrated into finite element analysis or finite difference analysis. Recently, probabilistic analysis of earth slopes has been very attractive because it is capable to take the soil uncertainty into account. However, the SRM is less commonly extended to probabilistic framework compared to a variety of probabilistic LEM analysis of earth slopes. To overcome some limitations that hinder the development of probabilistic SRM stability analysis, a new procedure based on recursive algorithm FORM with sensitivity analysis in the space of original variables is proposed. It can be used to deal with correlated non-normal variables subjected to implicit limit state surface. Using the proposed approach, a probabilistic finite element analysis of the stability of an existing earth dam is carried out in this paper.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.209-212
• /
• 2002

Since Berkhoff proposed the mild-slope equation in 1972, it has widely been used for calculation of shallow water wave transformation. Recently, it was extended to give an extended mild-slope equation, which includes the bottom slope squared term and bottom curvature term so as to be capable of modeling wave transformation on rapidly varying topography. These equations were derived by integrating the Laplace equation vertically. In the present study, we develop a finite element model to solve the Laplace equation directly while keeping the same computational efficiency as the mild-slope equation. This model assumes the vertical variation of wave potential as a cosine hyperbolic function as done in the derivation of the mild-slope equation, and the Galerkin method is used to discretize . The computational domain was discretized with proper finite elements, while the radiation condition at infinity was treated by introducing the concept of an infinite element. The upper boundary condition can be either free surface or a solid structure. The applicability of the developed model was verified through example analyses of two-dimensional wave reflection and transmission. .

• Structural Engineering and Mechanics
• /
• v.62 no.1
• /
• pp.33-42
• /
• 2017

We present in this paper a finite element formulation for nonlinear torsional analysis of 3D beams with arbitrary composite cross-sections. Since the proposed formulation employs a continuum mechanics based beam element with kinematics enriched by the extended St. Venant solutions, it can precisely account higher order warping effect and its 3D couplings. We propose a numerical procedure to calculate the extended St. Venant equation and the twisting center of an arbitrary composite cross-section simultaneously. The accuracy and efficiency of the proposed formulation are thoroughly investigated through representative numerical examples.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.4
• /
• pp.174-186
• /
• 1998

A Galerkin's finite element model incorporating infinite elements for modeling of radiation condition at infinity has been developed, which is based on an extended mild-slope equation. To illustrate the validity and applicability of the present model, the example analyses were carried out for a resonance problem in the rectangular harbor of Ippen and Goda (1963) and for wave transformations over circular shoals of Sharp (1968) and Chandrasekera and Cheung (1997). Comparisons with the results obtained by hydraulic experiments and hybrid element method showed that the present model gives very good results in spite of the rapidly varying topography. Numerical experiments were also performed for wave transformations over a circular concave well which may be an alternative to conventional wave barriers.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.4
• /
• pp.335-343
• /
• 2017

The Newman-Raju formula and contour integral-based finite element analyses(FEAs) have been widely used to assess crack growth rates and residual lives at crack locations in ships or offshore structures, but the Newman-Raju formula is known to be less accurate for the complicated weld details and the conventional FEA-based contour integral approach needs concentrated efforts to construct FEA models. Recently, an extended finite element method(XFEM) has been proposed to reduce those modeling efforts with reliable accuracy. Stress intensity factors(SIFs) from the approaches such as the Newman-Raju formula, conventional FEA-based J-integral, and XFEM-based J-integral were compared for an infinitely long plate with a propagating elliptic crack. It was concluded that the XFEM approach was far reliable in terms of prediction ability of SIFs. Assuming a 25 year-aged coast guard patrol ship had the prescribed cracks at the bracket toes attached to longitudinal stiffeners in way of deck and bottom, SIFs were derived based on the three approaches. To obtain axial tension loads acting on the longitudinal stiffeners, long term hull girder bending moments were assumed to obey Weibull distribution of which two parameters were decided from a reference (DNV, 2014). For the complicated weld details, it was concluded that the XFEM approach could cost-effectively and accurately estimate the crack growth rates and residual lives of ship structures.

• Advances in nano research
• /
• v.15 no.2
• /
• pp.141-153
• /
• 2023

The main objective of this paper is to develop the finite element study on the nonlinear free vibration of functionally graded nanocomposite spherical shells reinforced with graphene platelets under the first-order shear deformation shell theory and von Kármán nonlinear kinematic relations. The governing equations are presented by introducing the full asymmetric nonlinear strain-displacement relations followed by the constitutive relations and energy functional. The extended Halpin-Tsai model is utilized to specify the overall Young's modulus of the nanocomposite. Then, the finite element formulation is derived and the quadrilateral 8-node shell element is implemented for finite element discretization. The nonlinear sets of dynamic equations are solved by the use of the harmonic balance technique and iterative method to find the nonlinear frequency response. Several numerical examples are represented to highlight the impact of involved factors on the large-amplitude vibration responses of nanocomposite spherical shells. One of the main findings is that for some geometrical and material parameters, the fundamental vibrational mode shape is asymmetric and the axisymmetric formulation cannot be appropriately employed to model the nonlinear dynamic behavior of nanocomposite spherical shells.

• Structural Engineering and Mechanics
• /
• v.57 no.1
• /
• pp.65-89
• /
• 2016

This paper presents the development of methodologies using Extended Finite Element Method (XFEM) for cracked unstiffened and concentric stiffened panels subjected to constant amplitude tensile fatigue loading. XFEM formulations such as level set representation of crack, element stiffness matrix formulation and numerical integration are presented and implemented in MATLAB software. Stiffeners of the stiffened panels are modelled using truss elements such that nodes of the panel and nodes of the stiffener coincide. Stress Intensity Factor (SIF) is computed from the solutions of XFEM using domain form of interaction integral. Paris's crack growth law is used to compute the number of fatigue cycles up to failure. Numerical investigations are carried out to model the crack growth, estimate the remaining life and generate damage tolerant curves. From the studies, it is observed that (i) there is a considerable increase in fatigue life of stiffened panels compared to unstiffened panels and (ii) as the external applied stress is decreasing number of fatigue life cycles taken by the component is increasing.

• Structural Engineering and Mechanics
• /
• v.87 no.3
• /
• pp.243-254
• /
• 2023

Engineering structures usually suffer from cracks. The crack geometry has an influence on the structural mechanical properties and subsequent crack propagations. However, as an extensively utilized method in fracture analysis, the extended finite element method provided by Abaqus fails to output the specific location and dimensions of fractures. In this study, a technique to capture the crack geometry is proposed. The technique is based on the invariant level set method (I-LSM), which can avoid updating the level set function during crack development. The solution is achieved by an open-source plug-in programmed by Python. Three examples were performed to verify the effectiveness and robustness of the program. The result shows that the developed program can accurately output the crack geometry in both the 2D and 3D models. The open-source plug-in codes are included as supplementary material.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.