Search results

The purpose of this study is to simulate the tensile and shear types of fractures using the mixed fracture criteria considering the energy evolution based on the dual bilinear cohesive zone model and investigate the dynamic propagation of tensile and shear fractures induced by an impact load in rock. The propagation of tension and shear at different scales induced by the impact load is also an important aspect of this study.

In this study, based on the well-developed dual bilinear cohesive zone model and combined finite element-discrete element method, the dynamic propagation of tensile and shear fractures induced by the impact load in rock is investigated. Some key technologies, such as the governing partial differential equations, fracture criteria, numerical discretisation and detection and separation, are introduced to form the global algorithm and procedure. By comparing with the tensile and shear fractures induced by the impact load in rock disc in typical experiments, the effectiveness and reliability of the proposed method are well verified.

The dynamic propagation of tensile and shear fractures in the laboratory- and engineering-scale rock disc and rock strata are derived. The influence of mesh sensitivity, impact load velocities and load positions are investigated. The larger load velocities may induce larger fracture width and entire failure. When the impact load is applied near the left support constraint boundary, concentrated shear fractures appear around the loading region, as well as induced shear fracture band, which may induce local instability. The proposed method shows good applicability in studying the propagation of tensile and shear fractures under impact loads.

The proposed method can identify fracture propagation via the stress and energy evolution of rock masses under the impact load, which has potential to be extended into the investigation of the mixed fractures and disturbance of in-situ stresses during dynamic strata mining in deep energy development.

This paper aims to employ a multivariate nonlinear regression analysis to establish a predictive model for the final fracture area, while accounting for the impact of individual parameters.

This analysis is based on the numerical simulation data obtained, using the hybrid finite element–discrete element (FE–DE) method. The forecasting model was compared with the numerical results and the accuracy of the model was evaluated by the root mean square (RMS) and the RMS error, the mean absolute error and the mean absolute percentage error.

The multivariate nonlinear regression model can accurately predict the nonlinear relationships between injection rate, leakoff coefficient, elastic modulus, permeability, Poisson’s ratio, pore pressure and final fracture area. The regression equations obtained from the Newton iteration of the least squares method are strong in terms of the fit to the six sensitive parameters, and the model follow essentially the same trend with the numerical simulation data, with no systematic divergence detected. Least absolutely deviation has a significantly weaker performance than the least squares method. The percentage contribution of sensitive parameters to the final fracture area is available from the simulation results and forecast model. Injection rate, leakoff coefficient, permeability, elastic modulus, pore pressure and Poisson’s ratio contribute 43.4%, −19.4%, 24.8%, −19.2%, −21.3% and 10.1% to the final fracture area, respectively, as they increased gradually. In summary, (1) the fluid injection rate has the greatest influence on the final fracture area. (2)The multivariate nonlinear regression equation was optimally obtained after 59 iterations of the least squares-based Newton method and 27 derivative evaluations, with a decidability coefficient R2 = 0.711 representing the model reliability and the regression equations fit the four parameters of leakoff coefficient, permeability, elastic modulus and pore pressure very satisfactorily. The models follow essentially the identical trend with the numerical simulation data and there is no systematic divergence. The least absolute deviation has a significantly weaker fit than the least squares method. (3)The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field and numerical simulation. Its effectiveness can be trained and optimized by experimental and simulation data, and taking into account more basic data and establishing regression equations, containing more fracturing parameters will be the further research interests.

The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field and numerical simulation. Its effectiveness can be trained and optimized by experimental and simulation data, and taking into account more basic data and establishing regression equations, containing more fracturing parameters will be the further research interests.

The purpose of this paper is to propose a new combined finite-discrete element method (FDEM) to analyze the mechanical properties, failure behavior and slope stability of soil rock mixtures (SRM), in which the rocks within the SRM model have shape randomness, size randomness and spatial distribution randomness.

Based on the modeling method of heterogeneous rocks, the SRM numerical model can be built and by adjusting the boundary between soil and rock, an SRM numerical model with any rock content can be obtained. The reliability and robustness of the new modeling method can be verified by uniaxial compression simulation. In addition, this paper investigates the effects of rock topology, rock content, slope height and slope inclination on the stability of SRM slopes.

Investigations of the influences of rock content, slope height and slope inclination of SRM slopes showed that the slope height had little effect on the failure mode. The influences of rock content and slope inclination on the slope failure mode were significant. With increasing rock content and slope dip angle, SRM slopes gradually transitioned from a single shear failure mode to a multi-shear fracture failure mode, and shear fractures showed irregular and bifurcated characteristics in which the cut-off values of rock content and slope inclination were 20% and 80°, respectively.

This paper proposed a new modeling method for SRMs based on FDEM, with rocks having random shapes, sizes and spatial distributions.

The main aim of the current paper is to find a numerical plan for hydraulic fracturing problem with application in extracting natural gases and oil.

First, time discretization is accomplished via Crank-Nicolson and semi-implicit techniques. At the second step, a high-order finite element method using quadratic triangular elements is proposed to derive the spatial discretization. The efficiency and time consuming of both obtained schemes will be investigated. In addition to the popular uniform mesh refinement strategy, an adaptive mesh refinement strategy will be employed to reduce computational costs.

Numerical results show a good agreement between the two schemes as well as the efficiency of the employed techniques to capture acceptable patterns of the model. In central single-crack mode, the experimental results demonstrate that maximal values of displacements in x- and y- directions are 0.1 and 0.08, respectively. They occur around both ends of the line and sides directly next to the line where pressure takes impact. Moreover, the pressure of injected fluid almost gained its initial value, i.e. 3,000 inside and close to the notch. Further, the results for non-central single-crack mode and bifurcated crack mode are depicted. In central single-crack mode and square computational area with a uniform mesh, computational times corresponding to the numerical schemes based on the high order finite element method for spatial discretization and Crank-Nicolson as well as semi-implicit techniques for temporal discretizations are 207.19s and 97.47s, respectively, with 2,048 elements, final time T = 0.2 and time step size τ = 0.01. Also, the simulations effectively illustrate a further decrease in computational time when the method is equipped with an adaptive mesh refinement strategy. The computational cost is reduced to 4.23s when the governed model is solved with the numerical scheme based on the adaptive high order finite element method and semi-implicit technique for spatial and temporal discretizations, respectively. Similarly, in other samples, the reduction of computational cost has been shown.

This is the first time that the high-order finite element method is employed to solve the model investigated in the current paper.

In the present study, the application of a recent damage plasticity model is presented for nonlinear dynamic analysis of the Koyna gravity dam. This is a single surface isotropic damage plasticity concrete model, which is based on the decomposition of stresses and was proposed in a previous study. The theoretical aspects of the model are initially reviewed, and a few preliminary verification examples are illustrated. Thereafter, the HHT-α (i.e. Hilber–Hughes–Taylor) algorithm is presented for nonlinear dynamic analysis of concrete gravity dams.

Based on the prepared tools, nonlinear behavior of the Koyna Dam is studied by applying the invoked damage plasticity model. For this purpose, three cases are considered for the present study. Case A, which is based on the linear model, is mainly used for comparative purposes. The other two cases (B and C) correspond to the nonlinear (i.e. damage plasticity) model. The basic data for these two cases are similar. However, the employed damping algorithms are different and correspond to constant and variable damping algorithms, respectively. This means that the damping matrix is either kept constant or updated for all iterations of different time increments through the course of analysis.

The time histories of horizontal displacement at the dam crest were initially compared for the three cases: the linear Case A, and two nonlinear Cases B and C. It was observed that nonlinear cases’ responses begin to deviate from the corresponding linear case after the time of about 4.3 s. However, the amount of change for Case C (i.e. variable damping) was much greater than for Case B (i.e. constant damping). This was manifested initially in the peaks of response. It was also noticed that the period of response changed slightly for Case B in comparison with the linear Case A, while this change was significant for Case C. The obtained tensile and compressive damages were subsequently compared for the two nonlinear cases. For constant damping Case B, it was noticed that tensile damage occurred in the D/S face kink and on the U/S face slightly at a lower elevation. Moreover, it had a scattered nature. However, in variable damping Case C, it was noticed that tensile damage was much more localized and acted similar to a discrete crack. Of course, both cases also show tensile damages at the dam’s heel. In regard to compressive damages, it is observed that low values are occurring for both nonlinear cases as expected.

The application of a recent single surface isotropic damage plasticity concrete model is presented for nonlinear dynamic analysis of the Koyna gravity dam. The nonlinear response of the dam is investigated for two different damping algorithms. Moreover, the influence of variable characteristic length is also investigated in the latter part of this study.

1) The OE-MLPG penalty meshfree method is developed to solve cracked structure.(2) Smartening the numerical meshfree method by combining the particle swarm optimization (PSO) optimization algorithms and Voronoi computational geometric algorithm. (3). Selection of base functions, finding optimal penalty factor and distribution of appropriate nodal points to the accuracy of calculation in the meshless local Petrov–Galekrin (MLPG) meshless method.

Using appropriate shape functions and distribution of nodal points in local domains and sub-domains and choosing an approximation or interpolation method has an effective role in the application of meshless methods for the analysis of computational fracture mechanics problems, especially problems with geometric discontinuity and cracks. In this research, computational geometry technique, based on the Voronoi diagram (VD) and Delaunay triangulation and PSO algorithm, are used to distribute nodal points in the sub-domain of analysis (crack line and around it on the crack plane).

By doing this process, the problems caused by too closeness of nodal points in computationally sensitive areas that exist in general methods of nodal point distribution are also solved. Comparing the effect of the number of sentences of basic functions and their order in the definition of shape functions, performing the mono-objective PSO algorithm to find the penalty factor, the coefficient, convergence, arrangement of nodal points during the three stages of VD implementation and the accuracy of the answers found indicates, the efficiency of V-E-MLPG method with Ns = 7 and ß = 0.0037–0.0075 to estimation of 3D-stress intensity factors (3D-SIFs) in computational fracture mechanics.

The present manuscript is a continuation of the studies (Ref. [33]) carried out by the authors, about; feasibility assessment, improvement and solution of challenges, introduction of more capacities and capabilities of the numerical MLPG method have been used. In order to validate the modeling and accuracy of calculations, the results have been compared with the findings of reference article [34] and [35].

Thermal fractures initiated under cooling at the surfaces of a 2-D or 3-D structure propagate, arrest and coalesce, leading to its structural failure and material-property changes, while the same processes can happen in the rock mass between parallel hydraulic fractures filled with cold fluid, leading to enhanced fracture connectivity and permeability.

This study used a 2-D plane strain fracture model for mixed-mode thermal fractures from two parallel cooling surfaces. Fracture propagation was governed by the theory of linear elastic fracture mechanics, while the displacement and temperature fields were discretized using the adaptive finite element method. This model was validated using two numerical benchmarks with strong fracture curvature and then used to simulate the propagation and coalescence of thermal fractures in a long rock mass.

Modeling results show two regimes: (1) thermal fractures from a cooling surface propagate and arrest by following the theoretical solutions of half-plane fractures before the unfractured portion decreases to 20% rock-mass width and (2) some pairs of fractures from the opposite cooling surfaces tend to eventually coalesce. The fracture coalescence time is in a power law with rock-mass width.

These findings are relevant to both subsurface engineering and material engineering: structure failure is a key concern in the latter, while fracture coalescence can enhance the connectivity of thermal and hydraulic fractures and thus reservoir permeability in the former.

Fused Filament Fabrication (FFF) is an extrusion-based manufacturing process using fused thermoplastics. Despite its low cost, the FFF is not extensively used in high-value industrial sectors mainly due to parts' anisotropy (related to the deposition strategy) and residual stresses (caused by successive heating cycles). Thus, this study aims to investigate the process improvement and the optimization of the printed parts.

In this work, a meshless technique – the Radial Point Interpolation Method (RPIM) – is used to numerically simulate the viscoplastic extrusion process – the initial phase of the FFF. Unlike the FEM, in meshless methods, there is no pre-established relationship between the nodes so the nodal mesh will not face mesh distortions and the discretization can easily be modified by adding or removing nodes from the initial nodal mesh. The accuracy of the obtained results highlights the importance of using meshless techniques in this field.

Meshless methods show particular relevance in this topic since the nodes can be distributed to match the layer-by-layer growing condition of the printing process.

Using the flow formulation combined with the heat transfer formulation presented here for the first time within an in-house RPIM code, an algorithm is proposed, implemented and validated for benchmark examples.

The purpose of this study is to develop a new numerical algorithm to simulate the phase-field model.

First, the derivative of the temporal direction is discretized by a second-order linearized finite difference scheme where it conserves the energy stability of the mathematical model. Then, the isogeometric collocation (IGC) method is used to approximate the derivative of spacial direction. The IGC procedure can be applied on irregular physical domains. The IGC method is constructed based upon the nonuniform rational B-splines (NURBS). Each curve and surface can be approximated by the NURBS. Also, a map will be defined to project the physical domain to a simple computational domain. In this procedure, the partial derivatives will be transformed to the new domain by the Jacobian and Hessian matrices. According to the mentioned procedure, the first- and second-order differential matrices are built. Furthermore, the pseudo-spectral algorithm is used to derive the first- and second-order nodal differential matrices. In the end, the Greville Abscissae points are used to the collocation method.

In the numerical experiments, the efficiency and accuracy of the proposed method are assessed through two examples, demonstrating its performance on both rectangular and nonrectangular domains.

This research work introduces the IGC method as a simulation technique for the phase-field crystal model.

This paper develops a neural network surrogate model based on a discrete lattice approach to investigate the influence of complex microstructures on the emergent behavior of biological networks.

The adaptability of network-forming organisms, such as, slime molds, relies on fluid-to-solid state transitions and dynamic behaviors at the level of the discrete microstructure, which continuum modeling methods struggle to capture effectively. To address this challenge, we present an optimized approach that combines lattice spring modeling with machine learning to capture dynamic behavior and develop nonlinear constitutive relationships.

This integrated approach allows us to predict the dynamic response of biological materials with heterogeneous microstructures, overcoming the limitations of conventional trial-and-error lattice design. The study investigates the microstructural behavior of biological materials using a neural network-based surrogate model. The results indicate that our surrogate model is effective in capturing the behavior of discrete lattice microstructures in biological materials.

The combination of numerical simulations and machine learning endows simulations of the slime mold Physarum polycephalum with a more accurate description of its emergent behavior and offers a pathway for the development of more effective lattice structures across a wide range of applications.

The novelty of this research lies in integrating lattice spring modeling and machine learning to explore the dynamic behavior of biological materials. This combined approach surpasses conventional methods, providing a more holistic and accurate representation of emergent behaviors in organisms.