• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.1301-1302
• /
• 2008

In this study, the fracture strength for fracture mechanism porous calcium phosphate made from sintered with ${\beta}$-tricalcium phosphate obtained by wet precipitation procedure is analyzed using finite element method and experiment measurement. First, three $3{\times}3{\times}3mm^3$ and $5{\times}5{\times}5mm^3$ specimens are prepared and tomographic images of one $5{\times}5{\times}5mm^3$ specimen are obtained by micro focus X-ray CT. The compression tests using the specimens are carried out to measure the elastic modulus and fracture strength to analyze the fracture mechanism of porous calcium phosphate specimen. The tomographic images are reconstructed by 3D reconstruction program. The finite elements are directly built up in the reconstructed specimen. The numerical simulation for the compression tests is performed using the element. The mechanism of calcium phosphate of simulation are obtained by the compression tests using there cylindric specimen of height 19.5 mm and diameter 10 mm. From the results, the applicability of porous calcium phosphate is evaluated to care fracture and vacant bone of a patient as the reinforcement material.

• Structural Engineering and Mechanics
• /
• v.72 no.2
• /
• pp.191-201
• /
• 2019

Finite element analysis is one of the most used tools for studying femoral neck fracture. Nerveless, consensus concerning either the choice of material characteristics, damage law and /or geometric models (linear on nonlinear) remains unreached. In this work, we propose a numerical quasi-brittle damage model to describe the behavior of the proximal femur associated with two methods to evaluate the Young modulus. Eight proximal femur finite elements models were constructed from CT scan data (4 donors: 3 women; 1 man). The numerical computations showed a good agreement between the numerical curves (load - displacement) and the experimental ones. A very encouraging result is obtained when a comparison is made between the computed fracture loads and the experimental ones ($R^2=0.825$, Relative error =6.49%). All specific numerical computation provided very fair qualitative matches with the fracture patterns for the sideway fall simulation. Finally, the comparative study based on 32 simulations adopting linear and nonlinear meshing led to the conclusion that the quantitatively results are improved when a nonlinear mesh is used.

• Computers and Concrete
• /
• v.16 no.6
• /
• pp.933-961
• /
• 2015

Damage by high-speed impact fracture is a dominant mode of failure in several applications of concrete structures. Numerical modelling can play a crucial role in understanding and predicting complex fracture processes. The commonly used mesh-based Finite Element Method has difficulties in accurately modelling the high deformation and disintegration associated with fracture, as this often distorts the mesh. Even with careful re-meshing FEM often fails to handle extreme deformations and results in poor accuracy. Moreover, simulating the mechanism of fragmentation requires detachment of elements along their boundaries, and this needs a fine mesh to allow the natural propagation of damage/cracks. Smoothed Particle Hydrodynamics (SPH) is an alternative particle based (mesh-less) Lagrangian method that is particularly suitable for analysing fracture because of its capability to model large deformation and to track free surfaces generated due to fracturing. Here we demonstrate the capabilities of SPH for predicting brittle fracture by studying a slender concrete structure (column) under the impact of a high-speed projectile. To explore the effect of the projectile material behaviour on the fracture process, the projectile is assumed to be either perfectly-elastic or elastoplastic in two separate cases. The transient stress field and the resulting evolution of damage under impact are investigated. The nature of the collision and the constitutive behaviour are found to considerably affect the fracture process for the structure including the crack propagation rates, and the size and motion of the fragments. The progress of fracture is tracked by measuring the average damage level of the structure and the extent of energy dissipation, which depend strongly on the type of collision. The effect of fracture property (failure strain) of the concrete due to its various compositions is found to have a profound effect on the damage and fragmentation pattern of the structure.

• Journal of Ocean Engineering and Technology
• /
• v.30 no.6
• /
• pp.474-483
• /
• 2016

Accidental events such as collisions, groundings, and hydrocarbon explosions in marine structures can cause catastrophic damage. Thus, it is extremely important to predict the extent of such damage, which determines the total amount of oil spills and the residual hull girder strength. Punching fracture tests were conducted by Choung (2009b), where various sizes of indenters and circular unstiffened steel plates with different thicknesses were used to quasi-statically realize damage extents. A three-dimensional fracture strain surface was developed based on a reference (Choung et al., 2015b), where the average stress triaxiality and average normalized Lode angle were used as the parameters governing the fracture of ductile steels. In this study, new numerical analyses were performed using very fine axisymmetric elements in combination with an Abaqus user-subroutine to implement the three-dimensional fracture strain surface. Conventional numerical analyses were also conducted for the tests to identify the best fit fracture strain values by changing the fracture strains. Based on the phenomenon of the average normalized Lode angle starting out positive and then becoming slightly negative, it was inferred that the shear stress primarily dominates in determining the fractures locations, with a partial contribution from the compressive stress. It should be stated that the three-dimensional fracture surface effectively predicted at least the shear stress-dominant fracture behavior of a mild steel.

• Economic and Environmental Geology
• /
• v.28 no.5
• /
• pp.503-512
• /
• 1995

Groundwater flow in fracture networks is simulated using a discrete fracture (DF) model which assume that groundwater flows only through the fracture network. This assumption is available if the permeability of rock matrix is very low. It is almost impossible to describe fracture networks perfectly, so a stochastic approach is used. The stochastic approach assumes that the characteristic parameters in fracture network have special distribution patterns. The stochastic model generates fracture networks with some characteristic parameters. The finite element method is used to compute fracture flows. One-dimensional line element is the element type of the finite elements. The simulation results are shown by dominant flow paths in the fracture network. The dominant flow path can be found from the simulated groundwater flow field. The model developed in this study provides the tool to estimate the influences of characteristic parameters on groundwater flow in fracture networks. The influences of some characteristic parameters on the frcture flow are estimated by the Monte Carlo simulation based on 30 realizations.

• Structural Engineering and Mechanics
• /
• v.18 no.6
• /
• pp.731-744
• /
• 2004

Modified virtual crack closure integral (MVCCI) technique has become very popular for computation of strain energy release rate (SERR) and stress intensity factor (SIF) for 2-D crack problems. The objective of this paper is to propose a numerical integration procedure for MVCCI so as to generalize the technique and make its application much wider. This new procedure called as numerically integrated MVCCI (NI-MVCCI) will remove the dependence of MVCCI equations on the type of finite element employed in the basic stress analysis. Numerical studies on fracture analysis of 2-D crack (mode I and II) problems have been conducted by employing 4-noded, 8-noded (regular & quarter-point), 9-noded and 12-noded finite elements. For non-singular (regular) elements at crack tip, NI-MVCCI technique generates the same results as MVCCI, but the advantage for higher order regular and singular elements is that complex equations for MVCCI need not be derived. Gauss numerical integration rule to be employed for 8-noded singular (quarter-point) element for accurate computation of SERR and SIF has been recommended based on the numerical studies.

• Proceedings of the Korea Concrete Institute Conference
• /
• 1998.10a
• /
• pp.450-455
• /
• 1998

In this paper, finite element with embedded discontinuous line is introduced in order to avoid the difficulties of adding new nodal points along with crack growth in discrete crack model. With the discontinuous element using discontinuous shape function, stiffness matrix of finite element is derived and dual mapping technique for numerical integration is employed. Using the finite element program made with employed algorithms, algorithm is verified and fracture analysis of simple concrete beam is performed.

• Computers and Concrete
• /
• v.24 no.3
• /
• pp.207-222
• /
• 2019

This paper has presented an effective and accurate meso-scale finite element model for simulating the fracture process of concrete under compression-shear loading. In the proposed model, concrete is parted into four important phases: aggregates, cement matrix, interfacial transition zone (ITZ), and the initial defects. Aggregate particles were modelled as randomly distributed polygons with a varying size according to the sieve curve developed by Fuller and Thompson. With regard to initial defects, only voids are considered. Cohesive elements with zero thickness are inserted into the initial mesh of cement matrix and along the interface between aggregate and cement matrix to simulate the cracking process of concrete. The constitutive model provided by ABAQUS is modified based on Wang's experiment and used to describe the failure behaviour of cohesive elements. User defined programs for aggregate delivery, cohesive element insertion and modified facture constitutive model are developed based on Python language, and embedded into the commercial FEM package ABAQUS. The effectiveness and accuracy of the proposed model are firstly identified by comparing the numerical results with the experimental ones, and then it is used to investigate the effect of meso-structure on the macro behavior of concrete. The shear strength of concrete under different pressures is also involved in this study, which could provide a reference for the macroscopic simulation of concrete component under shear force.

• Structural Engineering and Mechanics
• /
• v.16 no.6
• /
• pp.713-730
• /
• 2003

This paper contains the results of the study on the development of fracture and crack propagation in quasi-brittle materials, such as concrete or rocks, using the Discrete Element Method (DEM). A new discrete element numerical model is proposed as the basis for analyzing the inelastic evolution and growth of cracks up to the point of gross material failure. The model is expected to predict the fracture behavior for the quasi-brittle material structure using the elementary aggregate level, the interaction between aggregate materials, and bond cementation. The algorithms generate normal and shear forces between two interfacing blocks and contains two kinds of contact logic, one for connected blocks and the other one for blocks that are not directly connected. The Mohr-Coulomb theory has been used for the fracture limit. In this algorithm the particles are moving based on the connected block logic until the forces increase up to the fracture limit. After passing the limit, the particles are governed by the discrete block logic. In setting up a discrete polygon element model, two dimensional polygons are used to investigate the response of an assembly of different shapes, sizes, and orientations with blocks subjected to simple applied loads. Several examples involving assemblies of particles are presented to show the behavior of the fracture and the failure process.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.3
• /
• pp.519-527
• /
• 1997

AISI 316 steel has been used extensively for heater and boiler tube of the structural plants such as power, chemical and petroleum plants under severe operating conditions. Usually, material degradation due to microcrack or precipitation of carbides and segregation of impurity elements, is occured by damage accumulated for long-term service at high temperature in this material. In this study, the effect of aging time on fracture toughness was investigated to evaluate the measurement of material degradation. The elastic-plastic fracture toughness behaviour of AISI 316 steel pipe aged at $550^{\circ}C$for 1h-10000h (the aged material) was characterized using the single specimen J-R curve technique and eletric potential drop method at normal loading rate(load-line displacement speed of 0.2mm/min) in room temperature and air environment. The fracture toughness data from above experiments is compared with the $J_{in}$ obtained from predicted values of crack initiation point using potential drop method.