Impact dynamics of flexible solids is important in engineering practice. Obtaining its dynamic response is a challenging task and usually achieved by numerical methods. The objectives of the study are twofold. Firstly...Impact dynamics of flexible solids is important in engineering practice. Obtaining its dynamic response is a challenging task and usually achieved by numerical methods. The objectives of the study are twofold. Firstly, the discrete singular convolution （DSC） is used for the first time to analyze the impact dynamics. Secondly, the efficiency of various numerical methods for dynamic analysis is explored via an example of a flexible rod hit by a rigid ball. Three numerical methods, including the conventional finite element （FE） method, the DSC algorithm, and the spectral finite element （SFE） method, and one proposed modeling strategy, the improved spectral finite element （ISFE） method, are involved. Numerical results are compared with the known analytical solutions to show their efficiency. It is demonstrated that the proposed ISFE modeling strategy with a proper length of con- ventional FE yields the most accurate contact stress among the four investigated models. It is also found that the DSC algorithm is an alternative method for collision problems.展开更多
A hybrid finite-discrete element method was implemented to study the fracture process of rough rock joints under direct shearing. The hybrid method reproduced the joint shear resistance evolution process from asperity...A hybrid finite-discrete element method was implemented to study the fracture process of rough rock joints under direct shearing. The hybrid method reproduced the joint shear resistance evolution process from asperity sliding to degradation and from gouge formation to grinding. It is found that, in the direct shear test of rough rock joints under constant normal displacement loading conditions, higher shearing rate promotes the asperity degradation but constraints the volume dilation, which then results in higher peak shear resistance, more gouge formation and grinding, and smoother new joint surfaces. Moreover, it is found that the joint roughness affects the joint shear resistance evolution through influencing the joint fracture micro mechanism. The asperity degradation and gouge grinding are the main failure micro-mechanism in shearing rougher rock joints with deeper asperities while the asperity sliding is the main failure micro-mechanism in shearing smoother rock joints with shallower asperities. It is concluded that the hybrid finite-discrete element method is a valuable numerical tool better than traditional finite element method and discrete element method for modelling the joint sliding, asperity degradation, gouge formation, and gouge grinding occurred in the direct shear tests of rough rock joints.展开更多
On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate an...On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.展开更多
The discrete element method (DEM) was used to simulate the flow characteristic and strength characteristic of the conditioned sands in the earth pressure balance (EPB) tunneling. In the laboratory the conditioned sand...The discrete element method (DEM) was used to simulate the flow characteristic and strength characteristic of the conditioned sands in the earth pressure balance (EPB) tunneling. In the laboratory the conditioned sands were reproduced and the slump test and the direct shear test of the conditioned sands were implemented. A DEM equivalent model that can simulate the macro mechanical characteristic of the conditioned sands was proposed,and the corresponding numerical models of the slump test and the shear test were established. By selecting proper DEM model parameters,the errors of the slump values between the simulation results and the test results are in the range of 10.3%-14.3%,and the error of the curves between the shear displacement and the shear stress calculated with the DEM simulation is 4.68%-16.5% compared with that of the laboratory direct shear test. This illustrates that the proposed DEM equivalent model can approximately simulate the mechanical characteristics of the conditioned sands,which provides the basis for further simulation of the interaction between the conditioned soil and the chamber pressure system of the EPB machine.展开更多
The petroleum industry has shown great interest in the study of drilling optimization on pre-salt formations given the low rates of penetration observed so far. Rate of penetration is the key to economically drill the...The petroleum industry has shown great interest in the study of drilling optimization on pre-salt formations given the low rates of penetration observed so far. Rate of penetration is the key to economically drill the pre-salt carbonate rock. This work presents the results of numerical modeling through finite element method and discrete element method for single cutter drilling in carbonate samples. The work is relevant to understand the mechanics of drill bit-rock interaction while drilling deep wells and the results were validated with experimental data raised under simulated downhole conditions. The numerical models were carried out under different geometrical configurations, varying the cutter chamfer size and back-rake angles. The forces generated on the cutter are translated into mechanical specific energy as this parameter is often used to measure drilling efficiency. Results indicate that the chamfer size does not change significantly the mechanical specific energy values, characteristic. Results also show there is a significant increase although the cutter aggressiveness is influenced by this geometrical in drilling resistance for larger values of back-rake angle.展开更多
In this work the authors simulate a contaminant transport problem in three dimensions that takes place in the soil of waste disposals. Such problem is modeled by a diffusion-dominated equation. The solution of this eq...In this work the authors simulate a contaminant transport problem in three dimensions that takes place in the soil of waste disposals. Such problem is modeled by a diffusion-dominated equation. The solution of this equation is addressed by using mixed finite element method for the spatial discretization of the equation. The resulting linear algebraic system is handled by an iterative domain decomposition procedure. This procedure is naturally parallelizable, and permits to implement computational codes in distributed memory machines in order to save on CPU time. Numerical results of the serial and parallel codes were compared with experimental results, and their performance measures were evaluated. The results indicate that the parallelizable procedure is an efficient tool for performing simulations of the problem.展开更多
采用耦合的有限元-离散元法(finite-discrete element method, FDEM)研究了圆形隧洞底板大变形灾变机制,并研究了地应力量值、侧压系数、底板位置和隧洞形状对底鼓灾害的影响。结果表明:1)在静水压力状态下,圆形隧洞底板大变形灾变力学...采用耦合的有限元-离散元法(finite-discrete element method, FDEM)研究了圆形隧洞底板大变形灾变机制,并研究了地应力量值、侧压系数、底板位置和隧洞形状对底鼓灾害的影响。结果表明:1)在静水压力状态下,圆形隧洞底板大变形灾变力学机制为隧洞开挖造成围岩径向应力降低和切向应力升高,当超过岩体强度后产生共轭剪切破裂,并伴随拉伸断裂;浅部破裂岩块在深部围岩弹性变形恢复和剪胀效应下发生向隧洞内的翻转大运动,进而引发底板大变形灾害;2)随着地应力量值、侧压系数和底板位置的变化,围岩破坏率、损伤破裂区半径和隧洞表面围岩最大位移量等也发生变化,但仍为破裂碎胀性的大变形;3)传统的滑移线场理论难以解释直墙拱形底板围岩的共轭剪切破裂和深部岩体破裂,采用FDEM和双轴压缩力学模型可成功解释上述现象,为隧洞底鼓大变形灾变机制提供新的研究思路。展开更多
Based on the flux equivalent principle of a single fracture, the discrete fracture concept was developed, in which the macroscopic fractures are explicitly described as （n-l） dimensional geometry element. On the fun...Based on the flux equivalent principle of a single fracture, the discrete fracture concept was developed, in which the macroscopic fractures are explicitly described as （n-l） dimensional geometry element. On the fundamental of this simplification, the discrete-fractured model was developed which is suitable for all types of fractured porous media. The principle of discrete-fractured model was introduced in detail, and the general mathematical model was expressed subsequently. The fully coupling discrete-fractured mathematical model was implemented using Galerkin finite element method. The validity and accuracy of the model were shown through the Buckley-Leverett problem in a single fracture. Then the discrete-fractured model was applied to the two different type fractured porous media. The effect of fractures on the water flooding in fractured reservoirs was investigated. The numerical results showed that the fractures made the porous media more heterogeneous and anisotropic, and that the orientation, size, type of fracture and the connectivity of fractures network have important impacts on the two-phase flow.展开更多
A combination of the classical Newton Method and the multigrid method, i.e., a Newton multigrid method is given for solving quasilinear parabolic equations discretized by finite elements. The convergence of the algori...A combination of the classical Newton Method and the multigrid method, i.e., a Newton multigrid method is given for solving quasilinear parabolic equations discretized by finite elements. The convergence of the algorithm is obtained for only one step Newton iteration per level. The asymptotically computational cost for quasilinear parabolic problems is O（NNk） similar to multigrid method for linear parabolic problems.展开更多
基金Supported by the National Natural Science Foundation of China(50830201)the Priority Academic Program Development of Jiangsu Higher Education Institutions~~
文摘Impact dynamics of flexible solids is important in engineering practice. Obtaining its dynamic response is a challenging task and usually achieved by numerical methods. The objectives of the study are twofold. Firstly, the discrete singular convolution （DSC） is used for the first time to analyze the impact dynamics. Secondly, the efficiency of various numerical methods for dynamic analysis is explored via an example of a flexible rod hit by a rigid ball. Three numerical methods, including the conventional finite element （FE） method, the DSC algorithm, and the spectral finite element （SFE） method, and one proposed modeling strategy, the improved spectral finite element （ISFE） method, are involved. Numerical results are compared with the known analytical solutions to show their efficiency. It is demonstrated that the proposed ISFE modeling strategy with a proper length of con- ventional FE yields the most accurate contact stress among the four investigated models. It is also found that the DSC algorithm is an alternative method for collision problems.
基金The first author would like to thank the supports of the NARGS, IRGS and AAS grants of Australia, and the National Science Foundation grants （No. 51574060 and No. 51079017） of China, in which the first author is the intemational collaborator. The academic visits of the third and fourth authors to the University of Tasmania are partly supported by a PhD visiting scholarship and an academic visiting scholarship, respectively, provided by the China Scholarship Council, which are greatly appreciated.
文摘A hybrid finite-discrete element method was implemented to study the fracture process of rough rock joints under direct shearing. The hybrid method reproduced the joint shear resistance evolution process from asperity sliding to degradation and from gouge formation to grinding. It is found that, in the direct shear test of rough rock joints under constant normal displacement loading conditions, higher shearing rate promotes the asperity degradation but constraints the volume dilation, which then results in higher peak shear resistance, more gouge formation and grinding, and smoother new joint surfaces. Moreover, it is found that the joint roughness affects the joint shear resistance evolution through influencing the joint fracture micro mechanism. The asperity degradation and gouge grinding are the main failure micro-mechanism in shearing rougher rock joints with deeper asperities while the asperity sliding is the main failure micro-mechanism in shearing smoother rock joints with shallower asperities. It is concluded that the hybrid finite-discrete element method is a valuable numerical tool better than traditional finite element method and discrete element method for modelling the joint sliding, asperity degradation, gouge formation, and gouge grinding occurred in the direct shear tests of rough rock joints.
基金supported by the Key Projects of the National Science Foundation of China(Nos.52178190,52078250,11832013)
文摘On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.
基金Project (2007CB714006) supported by the National Basic Research Program of China
文摘The discrete element method (DEM) was used to simulate the flow characteristic and strength characteristic of the conditioned sands in the earth pressure balance (EPB) tunneling. In the laboratory the conditioned sands were reproduced and the slump test and the direct shear test of the conditioned sands were implemented. A DEM equivalent model that can simulate the macro mechanical characteristic of the conditioned sands was proposed,and the corresponding numerical models of the slump test and the shear test were established. By selecting proper DEM model parameters,the errors of the slump values between the simulation results and the test results are in the range of 10.3%-14.3%,and the error of the curves between the shear displacement and the shear stress calculated with the DEM simulation is 4.68%-16.5% compared with that of the laboratory direct shear test. This illustrates that the proposed DEM equivalent model can approximately simulate the mechanical characteristics of the conditioned sands,which provides the basis for further simulation of the interaction between the conditioned soil and the chamber pressure system of the EPB machine.
文摘The petroleum industry has shown great interest in the study of drilling optimization on pre-salt formations given the low rates of penetration observed so far. Rate of penetration is the key to economically drill the pre-salt carbonate rock. This work presents the results of numerical modeling through finite element method and discrete element method for single cutter drilling in carbonate samples. The work is relevant to understand the mechanics of drill bit-rock interaction while drilling deep wells and the results were validated with experimental data raised under simulated downhole conditions. The numerical models were carried out under different geometrical configurations, varying the cutter chamfer size and back-rake angles. The forces generated on the cutter are translated into mechanical specific energy as this parameter is often used to measure drilling efficiency. Results indicate that the chamfer size does not change significantly the mechanical specific energy values, characteristic. Results also show there is a significant increase although the cutter aggressiveness is influenced by this geometrical in drilling resistance for larger values of back-rake angle.
文摘In this work the authors simulate a contaminant transport problem in three dimensions that takes place in the soil of waste disposals. Such problem is modeled by a diffusion-dominated equation. The solution of this equation is addressed by using mixed finite element method for the spatial discretization of the equation. The resulting linear algebraic system is handled by an iterative domain decomposition procedure. This procedure is naturally parallelizable, and permits to implement computational codes in distributed memory machines in order to save on CPU time. Numerical results of the serial and parallel codes were compared with experimental results, and their performance measures were evaluated. The results indicate that the parallelizable procedure is an efficient tool for performing simulations of the problem.
文摘采用耦合的有限元-离散元法(finite-discrete element method, FDEM)研究了圆形隧洞底板大变形灾变机制,并研究了地应力量值、侧压系数、底板位置和隧洞形状对底鼓灾害的影响。结果表明:1)在静水压力状态下,圆形隧洞底板大变形灾变力学机制为隧洞开挖造成围岩径向应力降低和切向应力升高,当超过岩体强度后产生共轭剪切破裂,并伴随拉伸断裂;浅部破裂岩块在深部围岩弹性变形恢复和剪胀效应下发生向隧洞内的翻转大运动,进而引发底板大变形灾害;2)随着地应力量值、侧压系数和底板位置的变化,围岩破坏率、损伤破裂区半径和隧洞表面围岩最大位移量等也发生变化,但仍为破裂碎胀性的大变形;3)传统的滑移线场理论难以解释直墙拱形底板围岩的共轭剪切破裂和深部岩体破裂,采用FDEM和双轴压缩力学模型可成功解释上述现象,为隧洞底鼓大变形灾变机制提供新的研究思路。
基金supported by the National Basic Research Program of China("973"Program)(Grant No.2011CB20100)the Important National Science and Technology Project of China(Grant No.2011ZX05014- 005-003HZ)+1 种基金the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20090133110006)the Fundamental Research Funds for the Central Universities(Grant No. 09CX04005A)
文摘Based on the flux equivalent principle of a single fracture, the discrete fracture concept was developed, in which the macroscopic fractures are explicitly described as （n-l） dimensional geometry element. On the fundamental of this simplification, the discrete-fractured model was developed which is suitable for all types of fractured porous media. The principle of discrete-fractured model was introduced in detail, and the general mathematical model was expressed subsequently. The fully coupling discrete-fractured mathematical model was implemented using Galerkin finite element method. The validity and accuracy of the model were shown through the Buckley-Leverett problem in a single fracture. Then the discrete-fractured model was applied to the two different type fractured porous media. The effect of fractures on the water flooding in fractured reservoirs was investigated. The numerical results showed that the fractures made the porous media more heterogeneous and anisotropic, and that the orientation, size, type of fracture and the connectivity of fractures network have important impacts on the two-phase flow.
基金This research is supported by the National Natural Science Foundation of China(10471011).
文摘A combination of the classical Newton Method and the multigrid method, i.e., a Newton multigrid method is given for solving quasilinear parabolic equations discretized by finite elements. The convergence of the algorithm is obtained for only one step Newton iteration per level. The asymptotically computational cost for quasilinear parabolic problems is O（NNk） similar to multigrid method for linear parabolic problems.
