• Computers and Concrete
• /
• 제3권2_3호
• /
• pp.79-90
• /
• 2006

This paper describes the results of an investigation into high mass-low velocity impact behaviour of reinforced concrete beams. Tests have been conducted on fifteen 2.7 m or 1.5 m span beams under drop-weight loads. A high-speed video camera has been used at rates of up to 4,500 frames per second in order to record the crack formation, propagation, particle spallation and scabbing. In some tests the strain in the reinforcement has been recorded using "Durham" strain gauged bars, a technique developed by Scott and Marchand (2000) in which the strain gauges are embedded in the bars, so that the strains in the reinforcement can be recorded without affecting the bond between the concrete and the reinforcement. The impact force acting on the beams has been measured using a load cell placed within the impactor. A high-speed data logging system has been used to record the impact load, strains, accelerations, etc., so that time histories can be obtained. This research has led to the development of computational techniques based on combined continuum/discontinuum methods (finite/discrete element methods) to permit the simulation of impact loaded reinforced concrete beams. The implementation has been within the software package ELFEN (2004). Beams, similar to those tested, have been analysed using ELFEN a good agreement has been obtained for both the load-time histories and the crack patterns.

• Smart Structures and Systems
• /
• 제22권4호
• /
• pp.459-468
• /
• 2018

In this paper, the effect of non-persistent joints was determined on the behavior of concrete specimens subjected to biaxial loading through numerical modeling using particle flow code in two dimensions (PFC2D). Firstly, a numerical model was calibrated by uniaxial, Brazilian and triaxial experimental results to ensure the conformity of the simulated numerical model's response. Secondly, sixteen rectangular models with dimension of 100 mm by 100 mm were developed. Each model contains two non-persistent joints with lengths of 40 mm and 20 mm, respectively. The angularity of the larger joint changes from $30^{\circ}$ to $90^{\circ}$. In each configuration, the small joint angularity changes from $0^{\circ}$ to $90^{\circ}$ in $30^{\circ}$ increments. All of the models were under confining stress of 1 MPa. By using of the biaxial test configuration, the failure process was visually observed. Discrete element simulations demonstrated that macro shear fractures in models are because of microscopic tensile breakage of a large number of bonded discs. The failure pattern in Rock Bridge is mostly affected by joint overlapping whereas the biaxial strength is closely related to the failure pattern.

• Structural Engineering and Mechanics
• /
• 제68권5호
• /
• pp.537-547
• /
• 2018

This research presents the effect of anisotropy of the hollow disc mode under Brazilian test using PFC3D. The Brazilian tensile strength test was performed on the hollow disc specimens containing the bedding layers and then these specimens were numerically modeled by using the two dimensional discrete element code (PFC3D) to calibrate this computer code for the simulation of the cracks propagation and cracks coalescence in the anisotropic bedded rocks. The thickness of each layer within the specimens varied as 5 mm, 10 mm and 20 mm and the layers angles were changed as $0^{\circ}$, $25^{\circ}$, $50^{\circ}$, $75^{\circ}$ and $90^{\circ}$. The diameter of internal hole was taken as 15 mm and the loading rate during the testing process kept as 0.016 mm/s. It has been shown that for layers angles below $25^{\circ}$ the tensile cracks produce in between the layers and extend toward the model boundary till interact and break the specimen. The failure process of the specimen may enhance as the layer angle increases so that the Brazilian tensile strength reaches to its minimum value when the bedding layers is between $50^{\circ}$ and $75^{\circ}$ but its value reaches to maximum at a layer angle of $90^{\circ}$. The number of tensile cracks decreases as the layers thickness increases and with increasing the layers angle, less layer mobilize in the failure process.

• Computers and Concrete
• /
• 제33권5호
• /
• pp.545-557
• /
• 2024

A real-case incident occurred where a 9-meter-high segment of a pre-fabricated concrete separation wall unexpectedly collapsed. This collapse was triggered by improperly depositing excavated soil against the wall's back, a condition for which the wall segments were not designed to withstand lateral earth pressure, leading to a flexural failure. The event's analysis, integrating technical data and observational insights, revealed that internal forces at the time of failure significantly exceeded the wall's capacity per standard design. The Lattice Discrete Particle Model (LDPM) further replicates the collapse mechanism. Our approach involved defining various parameter sets to replicate the concrete's mechanical response, consistent with the tested compressive strength. Subsequent stages included calibrating these parameters across different scales and conducting full-scale simulations. These simulations carried out with various parameter sets, were thoroughly analyzed to identify the most representative failure mechanism. We developed an equation from this analysis that quickly correlates the parameters to the wall's load-carry capacity, aligned with the simulation. Additionally, our study examined the wall's post-peak behavior, extending up to the point of collapse. This aspect of the analysis was essential for preventing failure, providing crucial time for intervention, and potentially averting a disaster. However, the reinforced concrete residual state is far from being fully understood. While it's impractical for engineers to depend on the residual state of structural elements during the design phase, comprehending this state is essential for effective response and mitigation strategies after initial failure occurs.

• 화약ㆍ발파
• /
• 제22권4호
• /
• pp.7-15
• /
• 2004

본 연구에서는 PFC3D를 사용한 폭원모델링 기법을 제안하고, 제안된 기법을 시멘트 모르타르와 같은 연약재료의 발파에 적용하여 그 적용성을 시험해 보았다. PFC3D는 개별요소법(DEM)을 기반으로 하고 있어 응력파의 전파와 재료의 동적 파괴현상을 모사하는데 적합한 코드로 분류된다. 폭원모델링 과정에서는 공내입자들의 반경을 팽창/수축시키는 기법을 통해 공벽입자들에 접촉력의 형태로 폭발압력을 부여하는 방법을 사용하였으며, 입력하중에 따라 공벽에서 유발되는 접촉력을 계산단계마다 측정 및 보정함으로써 폭발압력의 크기를 제어할 수 있도록 하였다. 시멘트 모르타르 블록의 발파모델링 과정에서는 기존의 외력을 이용하는 방법과 본 연구에서 제안하고 있는 접촉력을 이용하는 기법을 각기 적용함으로써 연약재료의 파괴과정을 정성적으로 비교하여 보았다. 해석결과, 제안된 폭원모델링 기법을 적용한다면 암석이나 콘크리트와 같은 공학재료들이 발파과정에서 보이는 파괴거동을 수치적으로 보다 유사하게 모사 할 수 있을 것으로 판단된다.

• Advances in concrete construction
• /
• 제13권 2호
• /
• pp.199-213
• /
• 2022

Experimental and discrete element approaches were used to examine the effects of F shape non-persistent joints on the failure behaviour of concrete under uniaxial compressive test. concrete specimens with dimensions of 200 cm×200 cm×50 cm were provided. Within the specimen, F shape non-persistent joint consisting three joints were provided. The large joint length was 6 cm, and the length of two small joints were 2 cm. Vertical distance between two small joints change from 1.5 cm to 4.5 cm with increment of 1.5 cm. In constant joint lengths, the angle of large joint change from 0° to 90° with increments of 30°. Totally 12 different models were tested under compression test. The axial load rate on the model was 0.05 mm/min. Concurrent with experimental tests, numerical simulation (Particle flow code in two dimension) were performed on the models containing F shape non-persistent joint. Distance between small joints and joint angles were similar to experimental one. the results indicated that the failure process was mostly governed by both of the Distance between small joints and joint angles. The axial loading rate on the model was 0.05 mm/min. The compressive strengths of the samples were related to the fracture pattern and failure mechanism of the discontinuities. Furthermore, it was shown that the compressive behaviour of discontinuities is related to the number of the induced tensile cracks which are increased by increasing the joint angle. In the first, there were only a few acoustic emission (AE) hits in the initial stage of loading, and then AE hits rapidly grow before the applied stress reached its peak. Furthermore, a large number of AE hits accompanied every stress drop. Finally, the failure pattern and failure strength are similar in both approaches i.e., the experimental testing and the numerical simulation approaches.

• Geomechanics and Engineering
• /
• 제33권2호
• /
• pp.133-140
• /
• 2023

Soil erosion can cause scouring and failures of underwater structures, therefore, various soil improvement techniques are used to increase the soil erosion resistance. The microbially induced calcium carbonate precipitation (MICP) method is proposed to increase the erosion resistance, however, there are only limited experimental and numerical studies on the use of MICP treatment for improvement of surface erosion resistance. Therefore, this study investigates the improvement in surface erosion resistance of sands by MICP through laboratory experiments and numerical modeling. The surface erosion behaviors of coarse sands with various calcium carbonate contents were first investigated via the erosion function apparatus (EFA). The test results showed that MICP treatment increased the overall erosion resistance, and the contribution of the precipitated calcium carbonate to the erosion resistance and critical shear stress was quantified in relation to the calcium carbonate contents. Further, these surface erosion processes occurring in the EFA test were simulated through the coupled computational fluid dynamics (CFD) and discrete element method (DEM) with the cohesion bonding model to reflect the mineral precipitation effect. The simulation results were compared with the experimental results, and the developed CFD-DEM model with the cohesion bonding model well predicted the critical shear stress of MICP-treated sand. This work demonstrates that the MICP treatment is effective in improving soil erosion resistance, and the coupled CFD-DEM with a bonding model is a useful and promising tool to analyze the soil erosion behavior for MICP-treated sand at a particle scale.

• Structural Engineering and Mechanics
• /
• 제81권1호
• /
• pp.81-91
• /
• 2022

The effects of interaction between concrete-gypsum interface and edge crack on the failure behavior of the specimens in senicircular bend (SCB) test were studied in the laboratory and also simulated numerically using the discrete element method. Some quarter circular specimens of gypsum and concrete with 5 cm radii and hieghts were separately prepared. Then the semicircular testing specimens were made by attaching one gypsum and one concrete sample to one another using a special glue and one edge crack is produced (in the interface) by do not using the glue in that part of the interface. The tensile strengths of concrete and gypsum samples were separately measured as 2.2 MPa and 1.3 MPa, respectively. during all testing performances a constant loading rate of 0.005 mm/s were stablished. The proposed testing method showed that the mechanism of failure and fracture in the brittle materials were mostly governed by the dimensions and number of discontinuities. The fracture toughnesses of the SCB samples were related to the fracture patterns during the failure processes of these specimens. The tensile behaviour of edge notch was related to the number of induced tensile cracks which were increased by decreasing the joint length. The fracture toughness of samples was constant by increasing the joint length. The failure process and fracture pattern in the notched semi-circular bending specimens were similar for both methods used in this study (i.e., the laboratory tests and the simulation procedure using the particle flow code (PFC2D)).

• Geomechanics and Engineering
• /
• 제28권2호
• /
• pp.117-133
• /
• 2022

Determination of the mechanical behaviour of jointed rock masses has been a challenge for rock engineers for decades. This problem is more pronounced for non-persistent jointed rock masses due to complicated interaction of rock bridges on the overall behaviour. This paper aims to study the effect of a non-persistent joint set configuration on the mechanical behaviour of rock materials under both uniaxial and biaxial compression tests using a discrete element code. The numerical simulation of biaxial compressive strength of rock masses has been challenging in the past due to shortcomings of bonded particle models in reproducing the failure envelope of rock materials. This problem was resolved in this study by employing the flat-joint contact model. The validity of the numerical model was investigated through a comprehensive comparative study against physical uniaxial and biaxial compression experiments. Good agreement was found between numerical and experimental tests in terms of the recorded peak strength and the failure mode in both loading conditions. Studies on the effect of joint orientation on the failure mode showed that four zones of intact, transition to block rotation, block rotation and transition to intact failure occurs when the joint dip angle varies from 0° to 90°. It was found that the applied confining stress can significantly alter the range of these zones. It was observed that the minimum strength occurs at the joint dip angle of around 45 degrees under different confining stresses. It was also found that the joint orientation can alter the post peak behaviour and the lowest brittleness was observed at the block rotation zone.

• Wind and Structures
• /
• 제30권4호
• /
• pp.433-450
• /
• 2020

The strong turbulence characteristic of typhoon not only will significantly change flow field characteristics surrounding the large-scale wind turbine and aerodynamic force distribution on surface, but also may cause morphological evolution of coast dune and thereby form sand storms. A 5MW horizontal-axis wind turbine in a wind power plant of southeastern coastal areas in China was chosen to investigate the distribution law of additional loads caused by wind-sand coupling movement of coast dune at landing of strong typhoons. Firstly, a mesoscale Weather Research and Forecasting (WRF) mode was introduced in for high spatial resolution simulation of typhoon "Megi". Wind speed profile on the boundary layer of typhoon was gained through fitting based on nonlinear least squares and then it was integrated into the user-defined function (UDF) as an entry condition of small-scaled CFD numerical simulation. On this basis, a synchronous iterative modeling of wind field and sand particle combination was carried out by using a continuous phase and discrete phase. Influencing laws of typhoon and normal wind on moving characteristics of sand particles, equivalent pressure distribution mode of structural surface and characteristics of lift resistance coefficient were compared. Results demonstrated that: Compared with normal wind, mesoscale typhoon intensifies the 3D aerodynamic distribution mode on structural surface of wind turbine significantly. Different from wind loads, sand loads mainly impact on 30° ranges at two sides of the lower windward region on the tower. The ratio between sand loads and wind load reaches 3.937% and the maximum sand pressure coefficient is 0.09. The coupling impact effect of strong typhoon and large sand particles is more significant, in which the resistance coefficient of tower is increased by 9.80% to the maximum extent. The maximum resistance coefficient in typhoon field is 13.79% higher than that in the normal wind field.