• International Journal of Fuzzy Logic and Intelligent Systems
• /
• v.5 no.2
• /
• pp.119-123
• /
• 2005

This paper describes a fuzzy-based system for analyzing the stress intensity factors (SIFs) of three-dimensional (3D) cracks. 3D finite element method(FEM) was used to obtain the SIF for subsurface cracks and surface cracks existing in inhomogeneous materials. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model by using the fuzzy theory. Nodes are generated by the bucketing method, and ten-noded quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. The singular elements such that the mid-point nodes near crack front are shifted at the quarter-points, and these are automatically placed along the 3D crack front. The complete FE model is generated, and a stress analysis is performed. The SIFs are calculated using the displacement extrapolation method. The results were compared with those surface cracks in homogeneous materials. Also, this system is applied to analyze cladding effect of surface cracks in inhomogeneous materials.

• Coupled systems mechanics
• /
• v.3 no.3
• /
• pp.305-318
• /
• 2014

The study deals with the physical modeling of a typical building frame resting on pile foundation and embedded in cohesive soil mass using complete three-dimensional finite element analysis. Two different pile groups comprising four piles ($2{\times}2$) and nine piles ($3{\times}3$) are considered. Further, three different pile diameters along with the various pile spacings are considered. The elements of the superstructure frame and those of the pile foundation are descretized using twenty-node isoparametric continuum elements. The interface between the pile and pile and soil is idealized using sixteen-node isoparametric surface elements. The current study is an improved version of finite element modeling for the soil elements compared to the one reported in the literature (Chore and Ingle 2008). The soil elements are discretized using eight-, nine- and twelve-node continuum elements. Both the elements of superstructure and substructure (i.e., foundation) including soil are assumed to remain in the elastic state at all the time. The interaction analysis is carried out using sub-structure approach in the parametric study. The total stress analysis is carried out considering the immediate behaviour of the soil. The effect of various parameters of the pile foundation such as spacing in a group and number piles in a group, along with pile diameter, is evaluated on the response of superstructure. The response includes the displacement at the top of the frame and bending moment in columns. The soil-structure interaction effect is found to increase displacement in the range of 58 -152% and increase the absolute maximum positive and negative moments in the column in the range of 14-15% and 26-28%, respectively. The effect of the soil- structure interaction is observed to be significant for the configuration of the pile groups and the soil considered in the present study.

• Structural Engineering and Mechanics
• /
• v.35 no.3
• /
• pp.283-299
• /
• 2010

This paper addresses the numerical simulation of fatigue crack growth in arbitrary 2D geometries under constant amplitude loading by the using a new finite element software. The purpose of this software is on the determination of 2D crack paths and surfaces as well as on the evaluation of components Lifetimes as a part of the damage tolerant assessment. Throughout the simulation of fatigue crack propagation an automatic adaptive mesh is carried out in the vicinity of the crack front nodes and in the elements which represent the higher stresses distribution. The fatigue crack direction and the corresponding stress-intensity factors are estimated at each small crack increment by employing the displacement extrapolation technique under facilitation of singular crack tip elements. The propagation is modeled by successive linear extensions, which are determined by the stress intensity factors under linear elastic fracture mechanics (LEFM) assumption. The stress intensity factors range history must be recorded along the small crack increments. Upon completion of the stress intensity factors range history recording, fatigue crack propagation life of the examined specimen is predicted. A consistent transfer algorithm and a crack relaxation method are proposed and implemented for this purpose. Verification of the predicted fatigue life is validated with relevant experimental data and numerical results obtained by other researchers. The comparisons show that the program is capable of demonstrating the fatigue life prediction results as well as the fatigue crack path satisfactorily.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1479-1486
• /
• 2019

The purpose of the present study is to develop the 3D static and noise simulator based on Galerkin Finite Element Method (GFEM) using the unstructured hexahedral elements. The 3D, 2G neutron diffusion and noise equations are discretized using the unstructured hexahedral by considering the linear approximation of the shape function in each element. The validation of the static calculation is performed via comparison between calculated results and reported data for the VVER-1000 benchmark problem. A sensitivity analysis of the calculation to the element type (unstructured hexahedral or tetrahedron elements) is done. Finally, the neutron noise calculation is performed for the neutron noise source of type of variable strength using the Green function technique. It is shown that the error reduction in the static calculation is considerable when the unstructured tetrahedron elements are replaced with the hexahedral ones. Since the neutron flux distribution and neutron multiplication factor are appeared in the neutron noise equation, the more accurate calculation of these parameters leads to obtaining the neutron noise distribution with high accuracy. The investigation of the changes of the neutron noise distribution in axial direction of the reactor core shows that the 3D neutron noise analysis is required instead of 2D.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.289-294
• /
• 2010

Numerical stud of an oscillating body in incompressible fluid is performed. Stabilized finite element method comprising of Streamline-Upwind/Petrov-Galerkin (SUPG) and Pressure-Stabilizing/Petrov-Galerkin (PSPG) formulations of linear triangular elements was employed to solve 2D incompressible Navier-Stokes equations whereas the motion of the body was considered by incorporating the arbitrary Langrangian-Eulerian(ALE) formulation. An algebraic moving mesh strategy is utilized for obtaining body conforming mesh deformation at each time step. Two tests cases, namely motion of a circular cylinder and of an airfoil in incompressible flow were analyzed. The model is first validated against the stationary cases and then the capability to handle moving boundaries is demonstrated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.2
• /
• pp.231-244
• /
• 2000

This paper presents a dynamic analysis technique for a 2-D soil-structure interaction problem in the frequency domain, which can directly be applied as an analysis tool for seismic response analyses of underground structures, tunnels, embankments, and so on. In this method, the structure and near-field soil is modeled by the standard finite elements, while the unbounded far-field soil is represented using the dynamic infinite elements in the frequency domain. The earthquake-input motion is regarded as traveling P and SV waves which are incident vertically from the far-field of underlying half-space to the near-field of layered medium. The equivalent earthquake forces are then calculated utilizing so-called fixed-exterior-boundary-method and the free-field responses including displacements and tractions. For the verification of the present study, seismic response analyses are carried out for a multi-layered half-space free-field soil medium and a cylindrical cavity embedded in a homogeneous half-space. Comparisons of the present results with solutions by other approaches indicate that the proposed methodology gives accurate estimates. Finally, an application example of seismic response analysis for a subway station is presented, which demonstrates the applicability of the present study.

• Coupled systems mechanics
• /
• v.6 no.4
• /
• pp.501-522
• /
• 2017

This paper introduces a numerical procedure to incorporate elasto-plastic local deformation effects in the dynamic analysis of beams. The appealing feature is that simple beam type finite elements can be used for the global model which needs not to be altered by the localized elasto-plastic deformations. An overlapping local sophisticated 2D membrane model replaces the internal forces of the beam elements in the predefined region where the localized deformations take place. An iterative coupling technique is used to perform this replacement. Comparisons with full membrane analysis are provided in order to illustrate the accuracy and efficiency of the method developed herein. In this study, the membrane formulation is able to capture the elasto-plastic material behaviour based on the von Misses yield criterion and the associated flow rule for plane stress. The Newmark time integration method is adopted for the step-by-step dynamic analysis.

• Transactions of Materials Processing
• /
• v.11 no.2
• /
• pp.179-186
• /
• 2002

In the finite element analysis of metal forming Processes using general Lagrangian formulation, element nodes in the mesh move and elements are distorted as the material is deformed. The excessive degeneracy of mesh interrupts finite element analysis and thus increases the error of plastic deformation energy, In this study, a remeshing scheme using so-called mesh compression method is proposed to effectively analyze the flash which is generated usually in hot forging processes. In order to verify the effectiveness of the method, several examples are tested in two-dimensional and three-dimensional problems.

• Interaction and multiscale mechanics
• /
• v.3 no.2
• /
• pp.192-211
• /
• 2010

In most framed structures the nonlinearities and the damages are localized, extending over a limited length of the structural member. In order to capture the details of the local damage, the segments of a member that have entered the nonlinear range may need to be analyzed using the three-dimensional element (3D) model whereas the rest of the member can be analyzed using the simpler one-dimensional (1D) element model with fewer degrees of freedom. An Element-Coupling model was proposed to couple the small scale solid 3D elements with the large scale 1D beam elements. The mixed dimensional coupling is performed imposing the kinematic coupling hypothesis of the 1D model on the interfaces of the 3D model. The analysis results are compared with test results of a reinforced concrete pipe column and a structure consisting of reinforced concrete columns and a steel space truss subjected to static and dynamic loading. This structure is a reduced scale model of a direct air-cooled condenser support platform built in a thermal power plant. The reduction scale for the column as well as for the structure was 1:8. The same structures are also analyzed using 3D solid elements for the entire structure to demonstrate the validity of the Element-Coupling model. A comparison of the accuracy and the computational effort indicates that by the proposed Element-Coupling method the accuracy is almost the same but the computational effort is significantly reduced.

• Structural Engineering and Mechanics
• /
• v.53 no.4
• /
• pp.703-723
• /
• 2015

In this paper optimization of volume fraction distribution in a thick hollow cylinder with finite length made of two-dimensional functionally graded material (2D-FGM) and subjected to steady state thermal and mechanical loadings is considered. The finite element method with graded material properties within each element (graded finite elements) is used to model the structure. Volume fractions of constituent materials on a finite number of design points are taken as design variables and the volume fractions at any arbitrary point in the cylinder are obtained via cubic spline interpolation functions. The objective function selected as having the normalized effective stress equal to one at all points that leads to a uniform stress distribution in the structure. Genetic Algorithm jointed with interior penalty-function method for implementing constraints is effectively employed to find the global solution of the optimization problem. Obtained results indicates that by using the uniform distribution of normalized effective stress as objective function, considerably more efficient usage of materials can be achieved compared with the power law volume fraction distribution. Also considering uniform distribution of safety factor as design criteria instead of minimizing peak effective stress affects remarkably the optimum volume fractions.