• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회A
• /
• 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
• /
• 제72권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
• /
• 제16권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.

• 한국해양공학회지
• /
• 제30권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.

• 자원환경지질
• /
• 제28권5호
• /
• pp.503-512
• /
• 1995

2차원 불연속 절리 모델 (Discrete Fracture Model)을 사용하여 절리망 내에서 지하수 흐름을 시뮬레이션하였다. 불연속 절리 모델에서는 지하수가 오직 절리망을 통해서 흐른다고 가정한다. 이와 같은 분석은 결정질암 같이 지질 매체의 투수율이 매우 낮은 경우에 유용하다. 하지만 불연속 절리망을 완벽하게 구현하는 것이 불가능하므로, 이에 접근하는 방법으로 확률 통계적 모델이 제안되었다. 확률-통계적 모텔은 특성인자(밀도, 방향, 길이, 틈새두께 등)가 특별한 분포 유형을 갖는다고 가정한다. 확률-통계적 모델은 가정된 분포를 따르도록 특성인자를 생성한다. 이 후 본 모델을 통해 분석된 몇몇 특성인자를 가지고 절리망을 생성한다. 절리망을 생성한 이 후 지하수의 유통을 계산하기 위해 유한요소법을 적용하였다. 이 때 일차원 선요소가 유한요소망의 주요 요소이다. 시뮬레이션 결과는 절리망 내의 주요 흐름 경로를 통해 보여진다. 절리망 내의 지하수 속도를 비교하여 주요 흐름 경로를 찾아낸다. 본 연구에서 개발된 모델은 절리망 내의 지하수 흐름에 특성인자들이 미치는 영향을 평가할 수 있는 방법을 제공한다. 이를 위하여 30번의 생성을 하는 몬테카를로 시뮬레이션을 통해서 여러 특성 인자들이 지하수 흐름에 미치는 영향을 평가하였다.

• Structural Engineering and Mechanics
• /
• 제18권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.

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 1998년도 가을 학술발표논문집(II)
• /
• 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
• /
• 제24권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
• /
• 제16권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.

• 대한기계학회논문집A
• /
• 제21권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.