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.

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.

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.

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.

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.

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.

Dynamic process analysis of the Niumiangou landslide triggered by the Wenchuan earthquake using the finite difference method and a modified discontinuous deformation analysis

The Niumiangou landslide was the largest landslide triggered by the 2008 Wenchuan earthquake,which was significantly affected by the amplification effect of seismic acceleration.The ringshear experiments indicated that the materials in the source area of the Niumiangou landslide were subjected to friction degradation under a big shear displacement,which may result in rapid movement of the landslide.In order to better understand the landslide movement and study the effect of the friction degradation on movement mechanisms,the dynamic process of Niumiangou landslide was simulated with a new numerical method,which combines the finite difference method(FDM)and the discontinuous deformation analysis(DDA).First,the FDM was used to study the initiation time,amplification effect and velocity of the landslide.Afterwards,these initiation velocities were applied to the blocks in the DDA model by corresponding coordination in the FDM model.A displacementdependent friction model of the sliding surface was incorporated into DDA code to further understand the kinetic behavior of the landslide.The results show that the displacement-dependent friction strongly decreases the friction coefficient of sliding surface under a big displacement,which can obviously promote the run-out and velocity of landslide.The model output well matches the topographic map formed by the landslide.This implies that the proposed model can be applied to the simulation of earthquake-induced landslides with amplification effect,and the friction degradation model is important to clarify the movement mechanism of high-speed and long-distance landslides.

Deformation Analysis and Fixture Design of Thin-walled Cylinder in Drilling Process Based on TRIZ Theory

Thin-walled cylindrical workpiece is easy to deform during machining and clamping processes due to the insufficient rigidi.Moreover,it’s also difficult to ensure the perpendicularity of flange holes during drilling process.In this paper,the element birth and death technique is used to obtain the axial deformation of the hole through finite element simulation.The measured value of the perpendicularity of the hole was compared with the simulated value to verify then the rationality of the simulation model.To reduce the perpendicularity error of the hole in the drilling process,the theory of inventive principle solution(TRIZ)was used to analyze the drilling process of thin-walled cylinder,and the corresponding fixture was developed to adjust the supporting surface height adaptively.Three different fixture supporting layout schemes were used for numerical simulation of drilling process,and the maximum,average and standard deviation of the axial deformation of the flange holes and their maximum hole perpendicularity errors were comparatively analyzed,and the optimal arrangement was optimized.The results show that the proposed deformation control strategy can effectively improve the drilling deformation of thin-walled cylindrical workpiece,thereby significantly improving the machining quality of the parts.

A new nano-scale surface marking technique for the deformation analysis of Mg-based alloys

In this work a new nanoscale surface marking technique,namely electron beam damage induced surface marking(EBDISM),is developed and tested for the first time on a fine-grained pure Mg surface.This technique utilizes focused high-energy electron beam of a scanning electron microscope to“burn”dense arrays of nano-scale grid patterns on the sample surface,and it is proved to be very effective for identifying and measuring localised deformation behaviours.However,the surface marking deposited by EBDISM is not permanent and it tends to deteriorate overtime.Cheap,easy to use and versatile,the EBDISM technique has a huge potential for quantitative measurement of strain field and nano-scale deformation analysis.

Methods of configuration test and deformation analysis for large airship

In recent years, high-altitude aerostats have been increasingly developed in the direction of multi-functionality and large size. Due to the large size and the high flexibility, new challenges for large aerostats have appeared in the configuration test and the deformation analysis. The methods of the configuration test and the deformation analysis for large airship have been researched and discussed. A tested method of the configuration,named internal scanning, is established to quickly obtain the spatial information of all surfaces for the large airship by the three-dimensional(3D) laser scanning technology. By using the surface wrap method, the configuration parameters of the large airship are calculated. According to the test data of the configuration, the structural dimensions such as the distances between the characteristic sections are measured. The method of the deformation analysis for the airship contains the algorithm of nonuniform rational B-splines(NURBS) and the finite element(FE)method. The algorithm of NURBS is used to obtain the reconfiguration model of the large airship. The seams are considered and the seam areas are divided. The FE model of the middle part of the large airship is established. The distributions of the stress and the strain for the large airship are obtained by the FE method. The position of the larger deformation for the airship is found.

A MIXED VARIATIONAL FORMULATION FOR LARGE DEFORMATION ANALYSIS OF PLATES

A mixed vuriational formulation for large deformation analysis of plates is introduced. In this formulation the equilibrium ami compatibility equations are satisfied identically by means of stress functions and displacement components, respectively, and the constitu,lye equations are satisfied in a least square sense. An example is solved and the results are compared with those available in the literature.Further, the functional is particularized for buckling analysis of plates and a simple example is solved to illustrate the theory.

Large Deformation Analysis and Synthesis of Elastic Closed-loop Mechanism Made of a Certain Spring Wire Described by Free Curves

Recently novel mechanisms with compact size and without many mechanical elements such as bearing are strongly required for medical devices such as surgical operation devices. This paper describes analysis and synthesis of elastic link mechanisms of a single spring beam which can be manufactured by NC coiling machines. These mechanisms are expected as disposable micro forceps. Smooth Curvature Model（SCM） with 3rd order Legendre polynomial curvature functions is applied to calculate large deformation of a curved cantilever beam by taking account of the balance between external and internal elastic forces and moments. SCM is then extended to analyze large deformation of a closed-loop curved elastic beam which is composed of multiple free curved beams. A closed-loop elastic link is divided into two free curved cantilever beams each of which is assumed as serially connected free curved cantilever beams described with SCM. The sets of coefficients of Legendre polynomials of SCM in all free curved cantilever beams are determined by taking account of the force and moment balance at connecting point where external input force is applied. The sets of coefficients of Legendre polynomials of a nonleaded closed-loop elastic link are optimized to design a link mechanism which can generate specified output motion due to input force applied at the assumed dividing point. For example, two planar micro grippers with a single pulling input force are analyzed and designed. The elastic deformation analyzed with proposed method agrees very well with that calculated with FEM. The designed micro gripper can generate the desired pinching motion. The proposed method can contribute to design compact and simple elastic mechanisms without high calculation costs.

Deformation Analysis and Optimization on Hollow Spherical Rotor in Electrostactically Suspended Gyroscope

Non-spherical of rotor was described with solution method. Electrostatically suspended gyroscope＇s hollow spherical rotor＇s structure was presented. The simulative analysis of static deformation, dynamic deformation and synthesize deformation of rotor under different working conditions using the finite element software were carried out. Its deformation law and volume were obtained. The structural parameters of the rotor were optimized. The value of pressure required when the rotor was machined was calculated. The analysis has important theoretical reference value to the design for hollow spherical rotor in electrostatically suspended gyroscope.

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.

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.

A continuous/discontinuous deformation analysis(CDDA)method based on deformable blocks for fracture modeling

In the framework of finite element meshes,a novel continuous/discontinuous deformation analysis(CDDA)method is proposed in this paper for modeling of crack problems.In the present CDDA,simple polynomial interpolations are defined at the deformable block elements,and a link element is employed to connect the adjacent block elements.The CDDA is particularly suitable for modeling the fracture propagation because the switch from continuous deformation analysis to discontinuous deformation analysis is natural and convenient without additional procedures.The SIFs(stress intensity factors)for various types of cracks,such as kinked cracks or curved cracks,can be easily computed in the CDDA by using the virtual crack extension technique(VCET).Both the formulation and implementation of the VCET in CDDA are simple and straightforward.Numerical examples indicate that the present CDDA can obtain high accuracy in SIF results with simple polynomial interpolations and insensitive to mesh sizes,and can automatically simulate the crack propagation without degrading accuracy.

Simulation of blast induced crater in jointed rock mass by discontinuous deformation analysis method

Rock blasting is a dynamic process accom panied with the propagations of shock waves and the dispersion of the explosion gas.This paper adopts the discontinuous deformation analysis(DDA)method to simulate the rock blasting process.A dynamic parameter adjustment and the non-reflecting boundary condition are implemented in the DDA method.The sub-block DDA method to simulate fracture problems is used.The blasting process in jointed rock mass is simulated by application of the explosion gas pressure on the expanding borehole walls and induced connected fracture surfaces around the boreholes.The blast craters with different overburdens are derived.The whole process including the explosion gas dispersion,borehole expansion,rock mass failure and cast,and the formation of the final blasting piles in rock blasting are well reproduced numerically.Parametric study for different overburdens is carried out,and the results are analyzed and discussed.

Parallelization of spherical discontinuous deformation analysis(SDDA) for geotechnical problems based on cloud computing environment

The computational efficiency of the traditional serial spherical discontinuous deformation analysis(SDDA) program has limited its application in geotechnical engineering problems that need a large number of spheres. The cloud computing technology is used to parallel the SDDA program for the first time in this research. The most computationally intensive portions of the SDDA program, i.e., contact detection and matrix solution, are parallelized with proposed algorithms. The accuracy of the cloud-based parallel SDDA program(CB-PSDDA) is verified first. Further efficiency tests show that significant speedup can be obtained with an 8-server configuration and the computing scale can be up to several tens of thousands of particles. The cloud-based parallelized SDDA program increases its potential in applications of deformation and failure analysis of large-scale and realistic geotechnical engineering problems.

Dam deformation analysis based on BPNN merging models

Hydropower has made a significant contribution to the economic development of Vietnam,thus it is important to monitor the safety of hydropower dams for the good of the country and the people.In this paper,dam horizontal displacement is analyzed and then forecasted using three methods:the multi-regression model,the seasonal integrated auto-regressive moving average(SARIMA)model and the back-propagation neural network(BPNN)merging models.The monitoring data of the Hoa Binh Dam in Vietnam,including horizontal displacement,time,reservoir water level,and air temperature,are used for the experiments.The results indicate that all of these three methods can approximately describe the trend of dam deformation despite their different forecast accuracies.Hence,their short-term forecasts can provide valuable references for the dam safety.

Contact-Stress-Based Stress Recovery Methods for Discontinuous Deformation Analysis

Discontinuous deformation analysis(DDA)has been widely applied for the simulation of block systems that have many discontinuous surfaces.The penalty method is utilized to ensure that there are no penetrations between blocks.A linear polynomial function for displacement leads to a constant stress for a block,which cannot precisely describe the stress field within the block.Therefore,a high-order polynomial displacement function and a flue mesh are always used to improve the precision of the stress field.However,these means are not practical for simulating block systems that have many contacts.In this paper,the contact-stress-based stress recovery methods are proposed for DDA.High-precision solutions for the contact stresses on the boundaries of the blocks are utilized.The first-order Gaussian point of a block is the block's centroid,where the constant stress obtained via DDA is of higher precision.The high-precision solutions for the stresses are utilized in the least squares method to recover a single block's inner stress field.The proposed methods enhance the resolution of the stress field inside a single block without increasing the computational effort in the main iterative process for displacement in DDA.Numerical examples are simulated using both the finite element method(FEM)with a fine mesh and the proposed DDA program.The recovered DDA results can accurately describe the distribution of the stresses in a single block and,in some areas,have the same precision as the FEM results.Moreover,the precision of the proposed methods improves as the gradient of the contact stress on the boundary decreases.