Grouting has been the most effective approach to mitigate water inrush disasters in underground engineering due to its ability to plug groundwater and enhance rock strength.Nevertheless,there is a lack of potent numer...Grouting has been the most effective approach to mitigate water inrush disasters in underground engineering due to its ability to plug groundwater and enhance rock strength.Nevertheless,there is a lack of potent numerical tools for assessing the grouting effectiveness in water-rich fractured strata.In this study,the hydro-mechanical coupled discontinuous deformation analysis(HM-DDA)is inaugurally extended to simulate the grouting process in a water-rich discrete fracture network(DFN),including the slurry migration,fracture dilation,water plugging in a seepage field,and joint reinforcement after coagulation.To validate the capabilities of the developed method,several numerical examples are conducted incorporating the Newtonian fluid and Bingham slurry.The simulation results closely align with the analytical solutions.Additionally,a set of compression tests is conducted on the fresh and grouted rock specimens to verify the reinforcement method and calibrate the rational properties of reinforced joints.An engineering-scale model based on a real water inrush case of the Yonglian tunnel in a water-rich fractured zone has been established.The model demonstrates the effectiveness of grouting reinforcement in mitigating water inrush disaster.The results indicate that increased grouting pressure greatly affects the regulation of water outflow from the tunnel face and the prevention of rock detachment face after excavation.展开更多
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 land...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.展开更多
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 ...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.展开更多
Since its proposal in 1985,the Discontinuous Deformation Analysis(DDA)has developed into a powerful tool for simulating the mechanical behaviour of rock masses and other discontinuous systems and has been widely appli...Since its proposal in 1985,the Discontinuous Deformation Analysis(DDA)has developed into a powerful tool for simulating the mechanical behaviour of rock masses and other discontinuous systems and has been widely applied in the field of geotechnical engineering.This paper systematically reviews the theoretical development of DDA,including the basic framework for solving two-dimensional and three-dimensional problems,improvements in contact detection algorithms,coupling with other numerical methods,and high-performance computing.In addition,the application of DDA in engineering problems such as wave propagation,rock blasting,slope instability,underground structures,and hydraulic structures is discussed in detail,demonstrating its wide applicability in solving complex geological and engineering problems.In the future,DDA deeply integrated with artificial intelligence(AI)may enhance its performance in simulation accuracy,computational efficiency,and engineering applications.Automated modelling,non-parametric simulation,and physical-numerical fusion in AI-based DDA may become new research directions.展开更多
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 m...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.展开更多
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...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.展开更多
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 tha...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.展开更多
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 ...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.展开更多
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 r...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.展开更多
Originated in the early 1990s, SCGM(1, m ) c model has enjoyed widespread application in the fields of urban planning, society economy prediction and modal control in recent years. However, none of these applications ...Originated in the early 1990s, SCGM(1, m ) c model has enjoyed widespread application in the fields of urban planning, society economy prediction and modal control in recent years. However, none of these applications have taken account of time lag effects in the modeling process. Aiming at overcoming the defect, the authors introduced time lag items into SCGM(1, m ) c model and developed a SCGM(1, m ) c model with time lag, then discusses in detail some principal problems in the model, such as parameters estimation, model verifying, model prediction, etc. The model was used on a real slope monitoring project and compared with the conventional SCGM(1, m ) c model. The results show an improvement of average models precision from 1.321 to 0.238 and total average of relative prediction errors from 12.41% to 7.98% when the modeling data length ranges from 29 to 48 in the slope monitoring case.展开更多
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 horizont...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.展开更多
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 synthesi...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.展开更多
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 f...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.展开更多
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...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.展开更多
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 ...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.展开更多
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...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.展开更多
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 p...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 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 displa...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.展开更多
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...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.展开更多
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 comp...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.展开更多
基金supported by the China Scholarship Council(CSC,Grant No.202108050072)JSPS KAKENHI(Grant No.JP19KK0121)。
文摘Grouting has been the most effective approach to mitigate water inrush disasters in underground engineering due to its ability to plug groundwater and enhance rock strength.Nevertheless,there is a lack of potent numerical tools for assessing the grouting effectiveness in water-rich fractured strata.In this study,the hydro-mechanical coupled discontinuous deformation analysis(HM-DDA)is inaugurally extended to simulate the grouting process in a water-rich discrete fracture network(DFN),including the slurry migration,fracture dilation,water plugging in a seepage field,and joint reinforcement after coagulation.To validate the capabilities of the developed method,several numerical examples are conducted incorporating the Newtonian fluid and Bingham slurry.The simulation results closely align with the analytical solutions.Additionally,a set of compression tests is conducted on the fresh and grouted rock specimens to verify the reinforcement method and calibrate the rational properties of reinforced joints.An engineering-scale model based on a real water inrush case of the Yonglian tunnel in a water-rich fractured zone has been established.The model demonstrates the effectiveness of grouting reinforcement in mitigating water inrush disaster.The results indicate that increased grouting pressure greatly affects the regulation of water outflow from the tunnel face and the prevention of rock detachment face after excavation.
基金supported by the National Natural Science Foundations of China(grant numbers U22A20601 and 52209142)the Opening fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(grant number SKLGP2022K018)+1 种基金the Science&Technology Department of Sichuan Province(grant number 2023NSFSC0284)the Science and Technology Major Project of Tibetan Autonomous Region of China(grant number XZ202201ZD0003G)。
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No.42372310).
文摘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.
文摘Since its proposal in 1985,the Discontinuous Deformation Analysis(DDA)has developed into a powerful tool for simulating the mechanical behaviour of rock masses and other discontinuous systems and has been widely applied in the field of geotechnical engineering.This paper systematically reviews the theoretical development of DDA,including the basic framework for solving two-dimensional and three-dimensional problems,improvements in contact detection algorithms,coupling with other numerical methods,and high-performance computing.In addition,the application of DDA in engineering problems such as wave propagation,rock blasting,slope instability,underground structures,and hydraulic structures is discussed in detail,demonstrating its wide applicability in solving complex geological and engineering problems.In the future,DDA deeply integrated with artificial intelligence(AI)may enhance its performance in simulation accuracy,computational efficiency,and engineering applications.Automated modelling,non-parametric simulation,and physical-numerical fusion in AI-based DDA may become new research directions.
基金The Scientific Innovation Research of College Graduates in Jiangsu Province(No.CXLX11_0143)
文摘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.
基金The project supported by the National Natural Science Foundation of China (50139010)
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.4167230041902290)Opening fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(No.SKLGP2017K015)。
文摘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.
基金supported by the Key Project of Chinese National Programs for Fundamental Research and Development(2010CB731502)the National Natural Science Foundation of China(50978745)
文摘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.
基金The authors would like to thank the National Natural Science Foundation of China(Grant Nos.51879184 and 52079091)for funding this work.
文摘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.
基金TheNationalNaturalScienceFoundationofChina! (No .497742 0 9)
文摘Originated in the early 1990s, SCGM(1, m ) c model has enjoyed widespread application in the fields of urban planning, society economy prediction and modal control in recent years. However, none of these applications have taken account of time lag effects in the modeling process. Aiming at overcoming the defect, the authors introduced time lag items into SCGM(1, m ) c model and developed a SCGM(1, m ) c model with time lag, then discusses in detail some principal problems in the model, such as parameters estimation, model verifying, model prediction, etc. The model was used on a real slope monitoring project and compared with the conventional SCGM(1, m ) c model. The results show an improvement of average models precision from 1.321 to 0.238 and total average of relative prediction errors from 12.41% to 7.98% when the modeling data length ranges from 29 to 48 in the slope monitoring case.
基金This research was funded by the China Scholarship Council(CSC)and partially supported by the Project 911(Vietnam).The data analysis was carried out as a part of the second author’s PhD studies at the School of Geodesy and Geomatics,Wuhan University,People’s Republic of China[grant number 2011GXZN02].
文摘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.
文摘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.
基金support from the Australian Research Council(LP180100048)Baosteel Company。
文摘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.
基金supported by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (11421091)。
文摘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.
文摘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.
文摘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.
文摘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 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.
基金supported by the Shandong Province Scientific Research Fund Project of Binzhou University under Grant No. BZXYLG200704
文摘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.
基金Supported by National Natural Science Foundation of China(Grant No.51875033)Fundamental Research Funds for the Central Universities of China(Grant No.2021YJS137).
文摘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.
