Particle Discontinuous Deformation Analysis of Static and Dynamic Crack Propagation in Brittle Material

Crack propagation in brittle material is not only crucial for structural safety evaluation,but also has a wideranging impact on material design,damage assessment,resource extraction,and scientific research.A thorough investigation into the behavior of crack propagation contributes to a better understanding and control of the properties of brittle materials,thereby enhancing the reliability and safety of both materials and structures.As an implicit discrete elementmethod,the Discontinuous Deformation Analysis(DDA)has gained significant attention for its developments and applications in recent years.Among these developments,the particle DDA equipped with the bonded particle model is a powerful tool for predicting the whole process of material from continuity to failure.The primary objective of this research is to develop and utilize the particle DDAtomodel and understand the complex behavior of cracks in brittle materials under both static and dynamic loadings.The particle DDA is applied to several classical crack propagation problems,including the crack branching,compact tensile test,Kalthoff impact experiment,and tensile test of a rectangular plate with a hole.The evolutions of cracks under various stress or geometrical conditions are carefully investigated.The simulated results are compared with the experiments and other numerical results.It is found that the crack propagation patterns,including crack branching and the formation of secondary cracks,can be well reproduced.The results show that the particle DDA is a qualified method for crack propagation problems,providing valuable insights into the fracture mechanism of brittle materials.

Application of modified discontinuous deformation analysis to dynamic modelling of the Baige landslide in the Jinsha River

Accurate dynamic modeling of landslides could help understand the movement mechanisms and guide disaster mitigation and prevention.Discontinuous deformation analysis(DDA)is an effective approach for investigating landslides.However,DDA fails to accurately capture the degradation in shear strength of rock joints commonly observed in high-speed landslides.In this study,DDA is modified by incorporating simplified joint shear strength degradation.Based on the modified DDA,the kinematics of the Baige landslide that occurred along the Jinsha River in China on 10 October 2018 are reproduced.The violent starting velocity of the landslide is considered explicitly.Three cases with different violent starting velocities are investigated to show their effect on the landslide movement process.Subsequently,the landslide movement process and the final accumulation characteristics are analyzed from multiple perspectives.The results show that the violent starting velocity affects the landslide motion characteristics,which is found to be about 4 m/s in the Baige landslide.The movement process of the Baige landslide involves four stages:initiation,high-speed sliding,impact-climbing,low-speed motion and accumulation.The accumulation states of sliding masses in different zones are different,which essentially corresponds to reality.The research results suggest that the modified DDA is applicable to similar high-level rock landslides.

General Kinetostatic Modeling and Deformation Analysis of a Two-Module Rod-Driven Continuum Robot with Friction Considered

Continuum robots actuated by flexible rods have large potential applications,such as detection and operation tasks in confined environments,since the push and pull actuation of flexible rods withstand tension and compressive force,and increase the structure's rigidity.In this paper,a generalized kinetostatics model for multi-module and multi-segment continuum robots considering the effect of friction based on the Cosserat rod theory is established.Then,the model is applied to a two-module rod-driven continuum robot with winding ropes to analyze its deformation and load characteristics.Four different in-plane configurations under the external load term as S1,S2,C1,and C2 are defined.Taking a bending plane as an example,the tip deformation along thex-axis of these shapes is simulated and compared,which shows that the load capacity of C1 and C2 is generally larger than that of S1 and S2.Furthermore,the deformation experiments and simulations show that the maximum error ratio without external loads relative to the total length is no more than 3%,and it is no more than 4.7%under the external load.The established kinetostatics model is proven sufficient to accurately analyze the rod-driven continuum robot with the consideration of internal friction.

On the calibration and verification of Voronoi-based discontinuous deformation analysis for modeling rock fracture

Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to rock fracture simulation.However,parameter calibration,a fundamental issue in discontinuum methods,has not received enough attention in DDA.In this study,the parameter calibration of DDA for intact rock is carefully studied.To this end,a subblock DDA with Voronoi tessellation is presented first.Then,a modified contact constitutive law is introduced,in which the tensile and shear meso-strengths are modified to be independent of the bond lengths.This improvement can prevent the unjustified preferential failure of short edges.A method for imposing confining pressure is also introduced.Thereafter,sensitivity analysis is performed to investigate the influence of the calculated parameters and meso-parameters on the mechanical properties of modeled rock.Based on the sensitivity analysis,a unified calibration procedure is suggested for both cases with and without confining pressure.Finally,the calibration procedure is applied to two examples,including a biaxial compression test.The results show that the proposed Voronoi-based DDA can simulate rock fracture with and without confining pressure very well after careful parameter calibration.

Slope deformation partitioning and monitoring points optimization based on cluster analysis

The scientific and fair positioning of monitoring locations for surface displacement on slopes is a prerequisite for early warning and forecasting.However,there is no specific provision on how to effectively determine the number and location of monitoring points according to the actual deformation characteristics of the slope.There are still some defects in the layout of monitoring points.To this end,based on displacement data series and spatial location information of surface displacement monitoring points,by combining displacement series correlation and spatial distance influence factors,a spatial deformation correlation calculation model of slope based on clustering analysis was proposed to calculate the correlation between different monitoring points,based on which the deformation area of the slope was divided.The redundant monitoring points in each partition were eliminated based on the partition's outcome,and the overall optimal arrangement of slope monitoring points was then achieved.This method scientifically addresses the issues of slope deformation zoning and data gathering overlap.It not only eliminates human subjectivity from slope deformation zoning but also increases the efficiency and accuracy of slope monitoring.In order to verify the effectiveness of the method,a sand-mudstone interbedded CounterTilt excavation slope in the Chongqing city of China was used as the research object.Twenty-four monitoring points deployed on this slope were monitored for surface displacement for 13 months.The spatial location of the monitoring points was discussed.The results show that the proposed method of slope deformation zoning and the optimized placement of monitoring points are feasible.

Correlation analysis of longitudinal cracks and vertical deformation within asphalt pavement of cold regions

The asphalt pavement longitudinal crack is a common distress in cold regions,resulting from uneven deformation of the subgrade.Analysis of the correlation law between uneven deformation and crack distress is of positive significance for understanding the mechanism of crack initiation,and putting forward treatment measures.In view of the complexity of longitudinal crack inducement and road surface deformation,the grey relational method was used to analyze this relationship.Through long-term monitoring of the vertical deformation data of typical road sections,the vertical deformation law of the pavement surface and its deformation characteristics under the action of temperature field are analyzed.Parameters such as vertical relative deformation,vertical relative deformation rate and vertical differential deformation VDSr were constructed to describe vertical deformation characteristics.Typical distribution characteristics of longitudinal fractures and their length and distribution characteristics are also described.The grey correlation analysis theory was utilized to analyze the relationship between deformation characteristics of sections,cross sections and monitoring points and longitudinal crack characteristics(length and location).The analysis reveals a linear positive correlation or a high correlation between several indicators.This study can provide a deeper understanding of the occurrence and development mechanism of longitudinal cracks in asphalt pavement of cold areas,and give references for the research of road engineering structure,materials and distress prevention.

The Deformation Analysis of the 3D Alignment Control Network Based on the Multiple Congruence Models

In the construction and maintenance of particle accelerators,all the accelerator elements should be installed in the same coordinate system,only in this way could the devices in the actual world be consistent with the design drawings.However,with the occurrence of the movements of the reinforced concrete cover plates at short notice or building deformations in the long term,the control points upon the engineering structure will be displaced,and the fitness between the subnetwork and the global control network may be irresponsible.Therefore,it is necessary to evaluate the deformations of the 3D alignment control network.Different from the extant investigations,in this paper,to characterize the deformations of the control network,all of the congruent models between the points measured in different epochs have been identified,and the congruence model with the most control points is considered as the primary or fundamental model,the remaining models are recognized as the additional ones.Furthermore,the discrepancies between the primary S-transformation parameters and the additional S-transformation parameters can reflect the relative movements of the additional congruence models.Both the iterative GCT method and the iterative combinatorial theory are proposed to detect multiple congruence models in the control network.Considering the actual work of the alignment,it is essential to identify the competitive models in the monitoring network,which can provide us a hint that,even the fitness between the subnetwork and the global control network is good,there are still deformations which may be ignored.The numerical experiments show that the suggested approaches can describe the deformation of the 3D alignment control network roundly.

Application of neural network merging model in dam deformation analysis

In order to improve the prediction accuracy and test the generalization ability of the dam deformation analysis model, the back-propagation（BP） neural network model for dam deformation analysis is studied, and the merging model is built based on the neural network BP algorithm and the traditional statistical model. The three models mentioned above are calculated and analyzed according to the long-term deformation observation data in Chencun Dam. The analytical results show that the average prediction accuracies of the statistical model and the BP neural network model are ~ 0.477 and ＋- 0.390 mm, respectively, while the prediction accuracy of the merging model is ~0. 318 mm, which is improved by 33% and 18% compared to the other two models, respectively. And the merging model has a better generalization ability and broad applicability.

Numerical Simulation of Surrounding Rock Deformation and Grouting Reinforcement of Cross-Fault Tunnel under Different Excavation Methods

Tunnel construction is susceptible to accidents such as loosening, deformation, collapse, and water inrush, especiallyunder complex geological conditions like dense fault areas. These accidents can cause instability and damageto the tunnel. As a result, it is essential to conduct research on tunnel construction and grouting reinforcementtechnology in fault fracture zones to address these issues and ensure the safety of tunnel excavation projects. Thisstudy utilized the Xianglushan cross-fault tunnel to conduct a comprehensive analysis on the construction, support,and reinforcement of a tunnel crossing a fault fracture zone using the three-dimensional finite element numericalmethod. The study yielded the following research conclusions: The excavation conditions of the cross-fault tunnelarray were analyzed to determine the optimal construction method for excavation while controlling deformationand stress in the surrounding rock. The middle partition method (CD method) was found to be the most suitable.Additionally, the effects of advanced reinforcement grouting on the cross-fault fracture zone tunnel were studied,and the optimal combination of grouting reinforcement range (140°) and grouting thickness (1m) was determined.The stress and deformation data obtained fromon-site monitoring of the surrounding rock was slightly lower thanthe numerical simulation results. However, the change trend of both sets of data was found to be consistent. Theseresearch findings provide technical analysis and data support for the construction and design of cross-fault tunnels.

Identification and hazard analysis of landslides triggered by earthquakes and rainfall

This study aims to utilize the Small Baseline Subset Interferometric Synthetic Aperture Radar(SBAS-In SAR)technique and Google Earth optical remote sensing images to analyze the area within 20 km around the epicenter of a M 3.9, earthquake that occurred in Tanchang County, Gansu Province, on December 28, 2020. The objective is to identify potential earthquake-induced landslides, assess their scale, and determine their impact range. The study results reveal the successful identification of two potential landslides in the 20 km radius around the epicenter. Through time-series deformation analysis, it was observed that these potential landslides were significantly influenced by both the earthquake and rainfall. Further estimation of these potential landslides indicates maximum depths of 7.4 m and 14.1 m for the failure surfaces, with volumes of 9.02 × 10~4m~3and 25.5 ×10~4m~3, respectively. Finally, based on the simulation analysis of Massflow software, the maximum thickness of soil accumulation in the final accumulation area after sliding of the potential landslide in Shangyaai is 12 m, the area of the final accumulation area is 1.75 × 10~4m~2, and the farthest movement distance is 1124 m. The maximum thickness of soil accumulation in the final accumulation area after sliding of the potential landslide in Wangshancun is 8 m, the area of the final accumulation area is 7.89 × 10~4m~2, and the farthest movement distance is 742 m.

Static and Thermal Analysis of Aluminium (413,390,384 and 332) Piston Using Finite Element Method

The main objective of this research was to examine the suitability of aluminium alloy to design a piston of an internal combustion engine for improvement in weight and cost reduction. The piston was modelled using Autodesk Inventor 2017 software. The modelled piston was then imported into Ansys for further analysis. Static structural and thermal analysis were carried out on the pistons of the four different materials namely: Al 413 alloy, Al 384 alloy, Al 390 alloy and Al332 alloy to determine the total deformation, equivalent Von Mises stress, maximum shear stress, and the safety factor. The results of the study revealed that, aluminium 332 alloy piston deformed less compared to the deformations of aluminium 390 alloy piston, aluminium 384 alloy piston and aluminium 413 alloy piston. The induced Von Mises stresses in the pistons of the four different materials were found to be far lower than the yield strengths of all the materials. Hence, all the selected materials including the implementing material have equal properties to withstand the maximum gas load. All the selected materials were observed to have high thermal conductivity enough to be able to withstand the operating temperature in the engine cylinders.

Numerical stress-deformation analysis of cut-off wall in clay-core rockfill dam on thick overburden

The cut-off wall in a clay-core rockfill dam built on a thick overburden layer is subjected to a large compressive pressure under the action of the loads such as the dead weight of both the dam and the overburden layer, the frictional force induced by the differential settlement between the cut-off wall and surrounding soils, and the water pressure. Thus, reduction of the stress of the cut-off wall has become one of the main problems for consideration in engineering design. In this paper, numerical analysis of a core rockfill dam built on a thick overburden layer was conducted and some factors influencing the stress-strain behaviors of the cut-off wall were investigated. The factors include the improvement of the overburden layer, the modeling approach for interfacial contact between the cut-off wall and surrounding soils, the modulus of the cut-off wall concrete, and the connected pattern between the cut-off wall and the clay core. The result shows that improving the overburden layer,selecting plastic concrete with a low modulus and high strength, and optimizing the connection between the cut-off wall and the clay core of the dam are effective measures of reducing the deformations and compressive stresses of the cut-off wall. In addition, both the Goodman element and the mud-layer element are suitable for simulating the interfacial contact between the cut-off wall and surrounding soils.

Long-term deformation analysis of Shuibuya concrete face rockfill dam based on response surface method and improved genetic algorithm

Due to the size effects of rockfill materials, the settlement difference between numerical simulation and in situ monitoring of rockfill dams is a topic of general concern.The constitutive model parameters obtained from laboratory triaxial tests often underestimate the deformation of high rockfill dams.Therefore, constitutive model parameters obtained by back analysis were used to calculate and predict the long-term deformation of rockfill dams.Instead of using artificial neural networks (ANNs), the response surface method (RSM) was employed to replace the finite element simulation used in the optimization iteration.Only 27 training samples were required for RSM, improving computational efficiency compared with ANN, which required 300 training samples.RSM can be used to describe the relationship between the constitutive model parameters and dam settlements.The inversion results of the Shuibuya concrete face rockfill dam (CFRD) show that the calculated settlements agree with the measured data, indicating the accuracy and efficiency of RSM.

Analysis and application of automatic deformation monitoring data for buildings and structures of mining area

The buildings and structures of mines were monitored automatically using modern surveying technology. Through the analysis of the monitoring data, the deformation characteristics were found out from three aspects containing points, lines and regions, which play an important role in understanding the stable state of buildings and structures. The stability and deformation of monitoring points were analysed, and time-series data of monitoring points were denoised with wavelet analysis and Kalman filtering, and exponent function and periodic function were used to get the ideal deformation trend model of monitoring points. Through calculating the monitoring data obtained, analyzing the deformation trend, and cognizing the deformation regularity, it can better service mine safety production and decision-making.

FORMULATIONS OF THE THREE-DIMENSIONAL DISCONTINUOUS DEFORMATION ANALYSIS METHOD

This paper extends the original 2D discontinuous deformation analysis(DDA)method proposed by Shi to 3D cases,and presents the formulations of the 3D DDA.The formulations maintain the characteristics of the original 2D DDA approach.Contacts between the blocks are detected by using Common-Plane (C-P) approach and the non-smooth contact,such as of vertex-to-vertex,vertex- to-edge and edge-to-edge types,can be handled easily based on the C-P method.The matrices of equilibrium equations have been given in detail for programming purposes.The C program codes for the 3D DDA are developed.The ability and accuracy of the formulations and the program are verified by the analytical solutions of several dynamic examples.The robustness and versatility of the algorithms presented in this paper are demonstrated with the aid of an example of scattering of densely packed cubes.Finally,implications and future extensions are discussed.

Analysis of deformation and internal flow patterns for rising single bubbles in different liquids

Gas–liquid multiphase flow is a significant phenomenon in chemical processes. The rising behaviors of single bubbles in the quiescent liquids have been investigated but the internal flow patterns and deformation rules of bubbles, which influence the mass transfer efficiency to a large extent, have received much less attention. In this paper, the volume of fluid method was used to calculate the bubble shapes, pressure, velocity distributions,and the flow patterns inside the bubbles. The rising behavior of the bubbles with four different initial diameters,i.e., 3 mm, 5 mm, 7 mm and 9 mm was investigated in four various liquids including water, 61.23% glycerol,86.73% glycerol and 100% glycerol. The results show that the liquid properties and bubble initial diameters have great impacts on bubble shapes. Moreover, flow patterns inside the bubbles with different initial diameters were analyzed and classified into three types under the condition of different bubble shapes. Three correlations for predicting the maximum internal circulation inside the bubbles in 86.73% glycerol were presented and the R-square values were all bigger than 0.98. Through analyzing the pressure and velocity distributions around the bubbles, four rules of bubble deformation were also obtained to explain and predict the shapes.

Isogeometric analysis of free-form Timoshenko curved beams including the nonlinear effects of large deformations

In the present paper, the isogeometric analysis（IGA） of free-form planar curved beams is formulated based on the nonlinear Timoshenko beam theory to investigate the large deformation of beams with variable curvature. Based on the isoparametric concept, the shape functions of the field variables（displacement and rotation） in a finite element analysis are considered to be the same as the non-uniform rational basis spline（NURBS） basis functions defining the geometry. The validity of the presented formulation is tested in five case studies covering a wide range of engineering curved structures including from straight and constant curvature to variable curvature beams. The nonlinear deformation results obtained by the presented method are compared to well-established benchmark examples and also compared to the results of linear and nonlinear finite element analyses. As the nonlinear load-deflection behavior of Timoshenko beams is the main topic of this article, the results strongly show the applicability of the IGA method to the large deformation analysis of free-form curved beams. Finally, it is interesting to notice that, until very recently, the large deformations analysis of free-form Timoshenko curved beams has not been considered in IGA by researchers.

Elastoplastic 3D Deformation and Stress Analysis of Strip Rolling

An elastoplastic method for analyzing the 3D deformation, stress and transverse distribution of tension stress during cold strip rolling is developed. The analysis is based on the elastoplastic variational principle in which a kinematically admissible velocity field is constructed with the lateral flow function as an unknown function. The stress distribution and volume strain distribution are obtained by solving the simultaneous equations formed by the longitudinal differential equation of equilibrium and constitutive equations. The lateral flow function is determined by minimizing the total energy dissipation rate. Experimental investigation was carried out on a reversible cold mill. The front tension stress distributions in cold rolled strips were measured by a multi roll segmented tension sensing shapemeter. The calculated results are in good agreement with the measured ones.

Thermal deformation of aluminum alloy casting materials for tire mold by numerical analysis

Thermal deformation of aluminum alloy casting materials for manufacturing the tire mold was numerically investigated.The AC7A and AC4C casting material was selected as casting material and the metal casting device was used in order to manufacture the mold product of automobile tire in the actual industrial field.The temperature distribution and the cooling time of casting materials were numerically calculated by finite element analysis (FEA).Also,the thermal deformation such as displacement and stress distribution was calculated from the temperature results.The thermal deformation was closely related to the temperature difference between the surface and inside of the casting.The numerical analysis results reveal that the thermal deformation of AC7A casting material is higher than that of AC4C casting material.Also,the thermal deformation results at the central part are larger than that on the side of casting because of the shrinkage caused by the cooling speed difference.

Validation and application of three-dimensional discontinuous deformation analysis with tetrahedron finite element meshed block

In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.