• Advances in nano research
• /
• 제3권3호
• /
• pp.143-168
• /
• 2015

Initiation of crack and its growth simulation requires accurate model of traction - separation law. Accurate modeling of traction-separation law remains always a great challenge. Atomistic simulations based prediction has great potential in arriving at accurate traction-separation law. The present paper is aimed at establishing a method to address the above problem. A method for traction-separation law prediction via utilizing atomistic simulations data has been proposed. In this direction, firstly, a simpler approach of common neighbor analysis (CNA) for the prediction of crack growth has been proposed and results have been compared with previously used approach of threshold potential energy. Next, a scheme for prediction of crack speed has been demonstrated based on the stable crack growth criteria. Also, an algorithm has been proposed that utilizes a variable relaxation time period for the computation of crack growth, accurate stress behavior, and traction-separation atomistic law. An understanding has been established for the generation of smoother traction-separation law (including the effect of free surface) from a huge amount of raw atomistic data. A new curve fit has also been proposed for predicting traction-separation data generated from the molecular dynamics simulations. The proposed traction-separation law has also been compared with the polynomial and exponential model used earlier for the prediction of traction-separation law for the bulk materials.

• Ocean Systems Engineering
• /
• 제6권1호
• /
• pp.23-60
• /
• 2016

The major objective of this study was to develop further understanding of 3D nearshore hydrodynamics under a variety of wave and tidal forcing conditions. The main tool used was a comprehensive 3D numerical model - combining the flow module of Delft3D with the WAVE solver of XBeach - of nearshore hydro- and morphodynamics that can simulate flow, sediment transport, and morphological evolution. Surf-swash zone hydrodynamics were modeled using the 3D Navier-Stokes equations, combined with various turbulence models (${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES). Sediment transport and resulting foreshore profile changes were approximated using different sediment transport relations that consider both bed- and suspended-load transport of non-cohesive sediments. The numerical set-up was tested against field data, with good agreement found. Different numerical experiments under a range of bed characteristics and incident wave and tidal conditions were run to test the model's capability to reproduce 3D flow, wave propagation, sediment transport and morphodynamics in the nearshore at the field scale. The results were interpreted according to existing understanding of surf and swash zone processes. Our numerical experiments confirm that the angle between the crest line of the approaching wave and the shoreline defines the direction and strength of the longshore current, while the longshore current velocity varies across the nearshore zone. The model simulates the undertow, hydraulic cell and rip-current patterns generated by radiation stresses and longshore variability in wave heights. Numerical results show that a non-uniform seabed is crucial for generation of rip currents in the nearshore (when bed slope is uniform, rips are not generated). Increasing the wave height increases the peaks of eddy viscosity and TKE (turbulent kinetic energy), while increasing the tidal amplitude reduces these peaks. Wave and tide interaction has most striking effects on the foreshore profile with the formation of the intertidal bar. High values of eddy viscosity, TKE and wave set-up are spread offshore for coarser grain sizes. Beach profile steepness modifies the nearshore circulation pattern, significantly enhancing the vertical component of the flow. The local recirculation within the longshore current in the inshore region causes a transient offshore shift and strengthening of the longshore current. Overall, the analysis shows that, with reasonable hypotheses, it is possible to simulate the nearshore hydrodynamics subjected to oceanic forcing, consistent with existing understanding of this area. Part II of this work presents 3D nearshore morphodynamics induced by the tides and waves.

• Computers and Concrete
• /
• 제29권 6호
• /
• pp.393-405
• /
• 2022

Lightweight concrete is a superior material due to its light weight and high strength. There however remain significant lacunae in engineering knowledge with regards to shear failure of lightweight fiber reinforced concrete beams. The main aim of the present study is to investigate the optimum usage of steel fibers in lightweight fiber reinforced concrete (LWFRC). Multi-scale finite element model calibrated with experimental results is developed to study the effect of steel fibers on the mechanical properties of LWFRC beams. To decrease the amount of steel fibers, it is preferred to reinforce only the middle section of the LWFRC beams, where the flexural stresses are higher. For numerical simulation, a multi-scale finite element model was developed. The cement matrix was modeled as homogeneous and uniform material and both steel fibers and lightweight coarse aggregates were randomly distributed within the matrix. Considering more realistic assumptions, the bonding between fibers and cement matrix was considered with the Cohesive Zone Model (CZM) and its parameters were determined using the model update method. Furthermore, conformity of Load-Crack Mouth Opening Displacement (CMOD) curves obtained from numerical modeling and experimental test results of notched beams under center-point loading tests were investigated. Validating the finite element model results with experimental tests, the effects of fibers' volume fraction, and the length of the reinforced middle section, on flexural and residual strengths of LWFRC, were studied. Results indicate that using steel fibers in a specified length of the concrete beam with high flexural stresses, and considerable savings can be achieved in using steel fibers. Reducing the length of the reinforced middle section from 50 to 30 cm in specimens containing 10 kg/m3 of steel fibers, resulting in a considerable decrease of the used steel fibers by four times, whereas only a 7% reduction in bearing capacity was observed. Therefore, determining an appropriate length of the reinforced middle section is an essential parameter in reducing fibers, usage leading to more affordable construction costs.

• Steel and Composite Structures
• /
• 제43권4호
• /
• pp.447-456
• /
• 2022

Recent experimental studies showed that deep steel I-shaped profiles classified as high ductility class sections in seismic design international codes exhibit low deformation capacity when subjected to cyclic loading. This paper presents an innovative retrofit solution to increase the rotation capacity of beams using bonded carbon fiber reinforced polymers (CFRP) patches validated with advanced finite element analysis. This investigation focuses on the flexural cyclic behaviour of I-shaped hot rolled steel deep section used as beams in moment-resisting frames (MRF) retrofitted with CFRP patches on the web. The main goal of this CFRP reinforcement is to increase the rotation capacity of the member without increasing the overstrength in order to avoid compromising the strong column-weak beam condition in MRF. A finite element model that simulates the cyclic plasticity behavior of the steel and the damage in the adhesive layer is developed. The damage is modelled using the cohesive zone modelling (CZM) technique that is able to capture the crack initiation and propagation. Details on the modelling techniques including the mesh sensitivity near the fracture zone are presented. The effectiveness of the retrofit solution depends strongly on the selection of the appropriate adhesive. Different adhesive types are investigated where the CZM parameters are calibrated from high fidelity fracture mechanics tests that are thoroughly validated in the literature. This includes a rigid adhesive commonly found in the construction industry and two tough adhesives used in the automotive industry. The results revealed that the CFRP patch can increase the rotation capacity of a steel member considerably when using tough adhesives.

• Steel and Composite Structures
• /
• 제35권4호
• /
• pp.463-474
• /
• 2020

In this paper, the influence of adding multi-walled carbon nanotube (MWCNT) on the pull behavior of steel and GFRP bars in ultra-high-performance concrete (UHPC) was examined experimentally and numerically. For numerical analysis, 3D nonlinear finite element modeling (FEM) with the help of ABAQUS software was used. Mechanical properties of the specimens, including Young's modulus, tensile strength and compressive strength, were extracted from the experimental results of the tests performed on standard cube specimens and for different values of weight percent of MWCNTs. In order to consider more realistic assumptions, the bond between concrete and bar was simulated using adhesive surfaces and Cohesive Zone Model (CZM), whose parameters were obtained by calibrating the results of the finite element model with the experimental results of pullout tests. The accuracy of the results of the finite element model was proved with conducting the pullout experimental test which showed high accuracy of the proposed model. Then, the effect of different parameters such as the material of bar, the diameter of the bar, as well as the weight percent of MWCNT on the bond behavior of bar and UHPC were studied. The results suggest that modifying UHPC with MWCNT improves bond strength between concrete and bar. In MWCNT per 0.01 and 0.3 wt% of MWCNT, the maximum pullout strength of steel bar with a diameter of 16 mm increased by 52.5% and 58.7% compared to the control specimen (UHPC without nanoparticle). Also, this increase in GFRP bars with a diameter of 16 mm was 34.3% and 45%.

• Composites Research
• /
• 제34권6호
• /
• pp.337-344
• /
• 2021

본 연구에서는 복합재료 팬 블레이드 도브테일 요소의 인장하중에 따른 점진적 파손거동을 유한요소 시뮬레이션을 통한 수치적 연구를 수행하고, 인장시험을 통하여 정확도를 검증한다. 도브테일 요소는 터보 팬 엔진의 팬 블레이드를 디스크와 결합시키는 조인트의 하나로, 통상 티타늄 등의 금속 재료로 제작되나 경량화 등의 이유로 복합재료의 적용이 연구되고 있다. 하지만 복합재료를 이용한 팬 블레이드 제조과정에서 드롭오프 플라이(Drop-off ply), 수지 포켓(resin pocket) 등의 제조 결함이 필연적으로 발생한다. 이러한 제조 결함이 복합재료 팬 블레이드 도브테일 요소에 미치는 영향을 확인하기 위해 유한요소모델을 이용한 수치해석을 수행하여 예측 결과와 인장시험 결과를 비교 분석한다. 이때 층간분리(delamination) 거동을 모사 가능한 응집영역 모델을 적용하였다. 결론적으로, 열 잔류응력 및 두께방향 압축하중에 의한 계면 물성 강화 효과를 고려하여 유한요소 해석결과와 시험결과 간의 높은 상사성을 얻을 수 있었다.

• 한국도로학회논문집
• /
• 제8권1호
• /
• pp.139-152
• /
• 2006

Many experimental and numerical approaches have been developed to evaluate paving materials and to predict pavement response and distress. Micromechanical simulation modeling is a technology that can reduce the number of physical tests required in material formulation and design and that can provide more details, e.g., the internal stress and strain state, and energy evolution and dissipation in simulated specimens with realistic microstructural features. A clustered distinct element modeling (DEM) approach was implemented In the two-dimensional particle flow software package (PFC-2D) to study the complex behavior observed in asphalt mixture fracturing. The relationship between continuous and discontinuous material properties was defined based on the potential energy approach. The theoretical relationship was validated with the uniform axial compression and cantilever beam model using two-dimensional plane strain and plane stress models. A bilinear cohesive displacement-softening model was implemented as an intrinsic interface and applied for both homogeneous and heterogeneous fracture modeling in order to simulate behavior in the fracture process zone and to simulate crack propagation. A disk-shaped compact tension test (DC(T)) with heterogeneous microstructure was simulated and compared with the experimental fracture test results to study Mode I fracture. The realistic arbitrary crack propagation including crack deflection, microcracking, crack face sliding, crack branching, and crack tip blunting could be represented in the fracture models. This micromechanical modeling approach represents the early developmental stages towards a 'virtual asphalt laboratory,' where simulations of laboratory tests and eventually field response and distress predictions can be made to enhance our understanding of pavement distress mechanisms, such its thermal fracture, reflective cracking, and fatigue crack growth.

• 토지주택연구
• /
• 제11권2호
• /
• pp.95-104
• /
• 2020

This study aims to explore a purely mechanistic pavement analysis approach where viscoelasticity and fracture of asphalt mixtures are considered to accurately predict deformation and damage behavior of flexible pavements. To do so, the viscoelastic and fracture properties of designated pavement materials are obtained through experiments and a fully mechanistic damage analysis is carried out using a finite element method (FEM). While modeling crack development can be done in various ways, this study uses the cohesive zone approach, which is a well-known fracture mechanics approach to efficiently model crack initiation and propagation. Different pavement configurations and traffic loads are considered based on three main functional classes of roads suggested by FHWA i.e., arterial, collector and local. For each road type, three different material combinations for asphalt concrete (AC) and base layers are considered to study damage behavior of pavement. A concept of the approach is presented and a case study where three different material combinations for AC and base layers are considered is exemplified to investigate progressive damage behavior of pavements when mixture properties and layer configurations were altered. Overall, it can be concluded that mechanistic pavement modeling attempted in this study could differentiate the performance of pavement sections due to varying design inputs. The promising results, although limited yet to be considered a fully practical method, infer that a few mixture tests can be integrated with the finite element modeling of the mixture tests and subsequent structural modeling of pavements to better design mixtures and pavements in a purely mechanistic manner.

• Composites Research
• /
• 제15권4호
• /
• pp.23-31
• /
• 2002

시편 게이지 면적($길이{\;}{\times}{\;}폭$)의 이차원 크기효과가 T300/924 $[45/-45/0/90]_3s$ 탄소섬유/에폭시 적층판의 압축거동에 대해 조사하였다. 개조된 압축시험치구(ICSTM)와 좌굴방지장치가 $30mm{\;}{\times}{\;}30mm,{\;}50mm{\;}{\times}{\;}50mm,{\;}70mm{\;}{\times}{\;}70mm,{\;}90mm{\;}{\times}{\;}90mm$의 게이지 길이와 폭을 가진 시편들의 압축시험에 사용하였다. 모든 경우의 파괴들은 시편 게이지 길이 내에서 주로 갑자기 발생하였다. 파괴 후 분석결과는 $0^{\circ}$층의 섬유의 미소좌굴에 의해 파괴를 시작하여 최종파괴를 일으키는 임계파괴기구일 것으로 생각되었다. 이것은 매트릭스 지배적인 파괴를 의미하며, 초기섬유굴곡에 따라 파괴가 지배적으로 시작된다는 것을 말한다 이것은 또한 제작공정과 품질이 압축강도를 결정하는 중요한 역할을 한다고 볼 수 있다. 좌굴방지장치를 장착하고 시험할 때 장치의 볼트 조임 토크에 따라 시편과의 접촉마찰 등에 의해 실제 압축강도 보다 크게 나타나는 결과를 보였다. 좌굴방지장치의 영향을 유한요소법을 이용하여 해석한 결과 실제 압축강도 보다 7% 정도 크게 나타남을 확인하였다. 부가적으로 홀을 갖는 시편들의 압축시험도 수행되었다. 홀에 의한 국부응력집중이 적층판 강도에 지배적 요인이었다. 파괴강도는 홀 크기와 시편 폭이 증가할수록 감소하였으나 탄성응력집중계수로 예측된 값보다는 일반적으로 크게 나타났다. 이것은 사용된 복합재가 이상적인 취성재질이 아니라는 것을 의미하며 홀 주위에서 다소간의 응력이완이 발생한다고 볼 수 있다. X선 검사 사진분석에서 섬유좌굴과 층간분리형태의 손상이 파괴하중의 약 80%에서 홀 가장자리로부터 시작되었고 임계파괴크랙길이인 2-3mm의 불안정한 상태에 도달하기 전까지는 하중 증가와 더불어 안정되게 파괴가 진전되었다(시편의 기하학적 크기에 의존함). 이 손상과 파괴는 선형 cohesive zone 모델로 해석되었다. 노치없는 시편의 압축강도와 평면 파괴인성의 측정된 적층판 변수들을 사용하여 홀의 크기와 시편 폭의 함수로서 홀을 갖는 적층판의 압축강도를 성공적으로 예측하였다.

• 한국복합재료학회:학술대회논문집
• /
• 한국복합재료학회 2002년도 추계학술발표대회 논문집
• /
• pp.88-91
• /
• 2002

The two dimensional size effect of specimen gauge section (length x width) was investigated on the compressive behavior of a T300/924 [45/-45/0/90]3s, carbon fiber-epoxy laminate. A modified ICSTM compression test fixture was used together with an anti-buckling device to test 3mm thick specimens with a 30$\times$30, 50$\times$50, 70$\times$70, and 90mm$\times$90mm gauge length by width section. In all cases failure was sudden and occurred mainly within the gauge length. Post failure examination suggests that $0^{\circ}$ fiber microbuckling is the critical damage mechanism that causes final failure. This is the matrix dominated failure mode and its triggering depends very much on initial fiber waviness. It is suggested that manufacturing process and quality may play a significant role in determining the compressive strength. When the anti-buckling device was used on specimens, it was showed that the compressive strength with the device was slightly greater than that without the device due to surface friction between the specimen and the device by pretoque in bolts of the device. In the analysis result on influence of the anti-buckling device using the finite element method, it was found that the compressive strength with the anti-buckling device by loaded bolts was about 7% higher than actual compressive strength. Additionally, compressive tests on specimen with an open hole were performed. The local stress concentration arising from the hole dominates the strength of the laminate rather than the stresses in the bulk of the material. It is observed that the remote failure stress decreases with increasing hole size and specimen width but is generally well above the value one might predict from the elastic stress concentration factor. This suggests that the material is not ideally brittle and some stress relief occurs around the hole. X-ray radiography reveals that damage in the form of fiber microbuckling and delamination initiates at the edge of the hole at approximately 80% of the failure load and extends stably under increasing load before becoming unstable at a critical length of 2-3mm (depends on specimen geometry). This damage growth and failure are analysed by a linear cohesive zone model. Using the independently measured laminate parameters of unnotched compressive strength and in-plane fracture toughness the model predicts successfully the notched strength as a function of hole size and width.