Based on the fact that the shear stress along anchorage segment is neither linearly nor uniformly distributed, the load transfer mechanism of the tension type anchor was studied and the mechanical characteristic of an...Based on the fact that the shear stress along anchorage segment is neither linearly nor uniformly distributed, the load transfer mechanism of the tension type anchor was studied and the mechanical characteristic of anchorage segment was analyzed. Shear stress?strain relationship of soil surrounding anchorage body was simplified into three-folding-lines model consisting of elastic phase, elasto-plastic phase and residual phase considering its softening characteristic. Meanwhile, shear displacement method that has been extensively used in the analysis of pile foundation was introduced. Based on elasto-plastic theory, the distributions of displacement, shear stress and axial force along the anchorage segment of tension type anchor were obtained, and the formula for calculating the elastic limit load was also developed accordingly. Finally, an example was given to discuss the variation of stress and displacement in the anchorage segment with the loads exerted on the anchor, and a program was worked out to calculate the anchor maximum bearing capacity. The influence of some parameters on the anchor bearing capacity was discussed, and effective anchorage length was obtained simultaneously. The results show that the shear stress first increases and then decreases and finally trends to the residual strength with increase of distance from bottom of the anchorage body, the displacement increases all the time with the increase of distance from bottom of the anchorage body, and the increase of velocity gradually becomes greater.展开更多
A new method for array calibration of array gain and phase uncertainties, which severely degrade the performance of spatial spectrum estimation, is presented. The method is based on the idea of the instrumental sensor...A new method for array calibration of array gain and phase uncertainties, which severely degrade the performance of spatial spectrum estimation, is presented. The method is based on the idea of the instrumental sensors method （ISM）, two well-calibrated sensors are added into the original array. By applying the principle of estimation of signal parameters via rotational invariance techniques （ESPRIT）, the direction-of-arrivals （DOAs） and uncertainties can be estimated simultaneously through eigen-decomposition. Compared with the conventional ones, this new method has less computational complexity while has higher estimation precision, what＇s more, it can overcome the problem of ambiguity. Both theoretical analysis and computer simulations show the effectiveness of the proposed method.展开更多
A revised displacement discontinuity method(DDM) program is developed for the simulation of rock joint propagation and dilatancy analysis. The non-linear joint model used in the program adopts Barton-Bandis normal def...A revised displacement discontinuity method(DDM) program is developed for the simulation of rock joint propagation and dilatancy analysis. The non-linear joint model used in the program adopts Barton-Bandis normal deformation model, Kulhaway shear deformation model and Mohr-Coulomb criterion. The joint propagation criterion is based on the equivalent stress intensity factor which can be obtained by regression analysis. The simulated rock joint propagation accords well with the existing knowledge. The closure and opening of joint is investigated by DDM, and it is shown that if the opening volume of propagated joint is larger than closure volume of the old joint, the joint dilatancy occurs. The dilatancy condition is mainly controlled by the normal stiffness of the rock joint. When the normal stiffness is larger than the critical value, joint dilatancy occurs. The critical normal stiffness of rock joint changes with the joint-load angle, and joint dilatancy is most possible to occur at 30°.展开更多
The higher order displacement discontinuity method(HODDM) utilizing special crack tip elements has been used in the solution of linear elastic fracture mechanics(LEFM) problems. The paper has selected several example ...The higher order displacement discontinuity method(HODDM) utilizing special crack tip elements has been used in the solution of linear elastic fracture mechanics(LEFM) problems. The paper has selected several example problems from the fracture mechanics literature(with available analytical solutions) including center slant crack in an infinite and finite body, single and double edge cracks, cracks emanating from a circular hole. The numerical values of Mode Ⅰ and Mode Ⅱ SIFs for these problems using HODDM are in excellent agreement with analytical results(reaching up to 0.001% deviation from their analytical results). The HODDM is also compared with the XFEM and a modified XFEM results. The results show that the HODDM needs a considerably lower computational effort(with less than 400 nodes) than the XFEM and the modified XFEM(which needs more than 10000 nodes) to reach a much higher accuracy. The proposed HODDM offers higher accuracy and lower computation effort for a wide range of problems in LEFM.展开更多
The present research is focused on the numerical crack coalescence analysis of the micro-cracks and cracks produced during the cutting action of TBM disc cutters. The linear elastic fracture mechanics(LEFM) concepts a...The present research is focused on the numerical crack coalescence analysis of the micro-cracks and cracks produced during the cutting action of TBM disc cutters. The linear elastic fracture mechanics(LEFM) concepts and the maximum tangential stress criterion are used to investigate the micro crack propagation and its direction underneath the excavating discs. A higher order displacement discontinuity method with quadratic displacement discontinuity elements is used to estimate the stress intensity factors near the crack tips. Rock cutting mechanisms under single and double type discs are simulated by the proposed numerical method.The main purposes of the present modeling are to simulate the chip formation process of indented rocks by single and double discs.The effects of specific disc parameters(except speed) on the thrust force Ft, the rolling force Fr, and the specific energy ES are investigated. It has been shown that the specific energy(energy required to cut through a unit volume of rock) of the double disc is less than that of the single disc. Crack propagation in rocks under disc cutters is numerically modeled and the optimum ratio of disc spacing S to penetration depth Pd(i.e. S/Pd ratio) of about 10 is obtained, which is in good agreement with the theoretical and experimental results cited in the literature.展开更多
In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat condu...In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat conduction techniques, this calibration approach does not require explicit input of the probe locations, thermophysical properties of the host material and temperature sensor parameters related to thermal contact resistance, sensor capacitance and conductive lead losses. All those parameters and properties are inherently contained in the calibration framework in terms of Volterra integral equation of the first kind. The Laplace transform technique is applied and the frequency domain manipulations of the heat equation are performed for deriving the calibration integral equation. Due to the ill-posed nature, regularization is required for the inverse heat conduction problem, a future-time method or singular value decomposition (SVD) can be used for stabilizing the ill-posed Volterra integral equation of the first kind.展开更多
The calculation precision and convergence speed of streamline strip element method are increased by (using) the method whose initial value of the exit lateral displacement is determined with strip element variation me...The calculation precision and convergence speed of streamline strip element method are increased by (using) the method whose initial value of the exit lateral displacement is determined with strip element variation method, and the accurate tension lateral distribution model is adopted based on the original third power spline function streamline strip element method. The basic theory of the strip element method is developed. The calculated results by the improved streamline strip element method and the original streamline strip element method are compared with the measured results, showing that the calculated results of the improved method are in good agreement with the measured results.展开更多
In engineering practice simplified methods are essential to the seismic design of embedded earth retaining walls,as fullydynamic numerical analyses are costly,time-consuming and require specific expertise.Recently dev...In engineering practice simplified methods are essential to the seismic design of embedded earth retaining walls,as fullydynamic numerical analyses are costly,time-consuming and require specific expertise.Recently developed pseudostatic methods provide earth stresses and internal forces,even in those cases in which the strength of the soil surrounding the structure is not entirely mobilised.Semiempirical correlations or Newmark sliding block method provide an estimate of earthquake-induced permanent displacements.However,the use of these methods is hindered by uncertainties in the evaluation of a few input parameters,affecting the reliability of the methods.This study uses 1 D site response analyses and 2 D fully-dynamic finite element analyses to show that simplified methods can provide a reasonable estimate of the maximum bending moment and permanent displacements for stiff cantilever walls embedded in uniform sand,providing that a few input parameters are evaluated through semiempirical correlations and a simple 1 D site response analysis.展开更多
为研究农业机械与水田壤土间的相互作用,需获取水田壤土的物理及接触参数。结合物理堆积试验,以休止角作为响应值,采用离散元法(DEM)并选取Hertz-Mindlin with JKR(Johnson-Kendall-Roberts)接触模型对长江中游地区水田壤土展开参数标...为研究农业机械与水田壤土间的相互作用,需获取水田壤土的物理及接触参数。结合物理堆积试验,以休止角作为响应值,采用离散元法(DEM)并选取Hertz-Mindlin with JKR(Johnson-Kendall-Roberts)接触模型对长江中游地区水田壤土展开参数标定研究。首先,通过物理堆积试验获取了壤土休止角(AoR)与含水率间的定量关系,由不同含水率土壤的堆积结果筛分出4种代表性堆积形态,由于水田壤土堆积体轮廓外形比较独特,因此仅对其左右两侧轮廓采用三次多项式进行局部拟合,计算其休止角。以长江中游地区水田壤土成因和预试验为依据来确定其离散元模型中9个参数的高低水平值,通过Plackett-Burman试验设计进行方差分析,发现壤土剪切模量、壤土间动摩擦因数、壤土与不锈钢间静摩擦因数和JKR表面能对AoR影响明显。然后,采用基于响应面法(RSM)原理的Box-Behnken试验设计(BBD)建立了AoR与4个显著性参数间的二次多项式回归模型。依据二次多项式回归模型对目标响应进行预测,得到最优参数组合。以此为基础对壤土AoR进行离散元仿真,AoR数值计算结果(45.4°)与试验结果(44.6°)相对误差为1.79%。最后,选取含水率分别为44.4%、48.7%的壤土进行堆积角仿真模拟,计算结果与堆积试验相对误差分别为2.8%、7.14%。研究表明:回归模型可以根据壤土含水率或AoR预测长江中游地区水田壤土的相关本征参数和接触参数。展开更多
基金Project(20050532021) supported by the Research Fund for the Doctoral Program of Higher Education
文摘Based on the fact that the shear stress along anchorage segment is neither linearly nor uniformly distributed, the load transfer mechanism of the tension type anchor was studied and the mechanical characteristic of anchorage segment was analyzed. Shear stress?strain relationship of soil surrounding anchorage body was simplified into three-folding-lines model consisting of elastic phase, elasto-plastic phase and residual phase considering its softening characteristic. Meanwhile, shear displacement method that has been extensively used in the analysis of pile foundation was introduced. Based on elasto-plastic theory, the distributions of displacement, shear stress and axial force along the anchorage segment of tension type anchor were obtained, and the formula for calculating the elastic limit load was also developed accordingly. Finally, an example was given to discuss the variation of stress and displacement in the anchorage segment with the loads exerted on the anchor, and a program was worked out to calculate the anchor maximum bearing capacity. The influence of some parameters on the anchor bearing capacity was discussed, and effective anchorage length was obtained simultaneously. The results show that the shear stress first increases and then decreases and finally trends to the residual strength with increase of distance from bottom of the anchorage body, the displacement increases all the time with the increase of distance from bottom of the anchorage body, and the increase of velocity gradually becomes greater.
文摘A new method for array calibration of array gain and phase uncertainties, which severely degrade the performance of spatial spectrum estimation, is presented. The method is based on the idea of the instrumental sensors method （ISM）, two well-calibrated sensors are added into the original array. By applying the principle of estimation of signal parameters via rotational invariance techniques （ESPRIT）, the direction-of-arrivals （DOAs） and uncertainties can be estimated simultaneously through eigen-decomposition. Compared with the conventional ones, this new method has less computational complexity while has higher estimation precision, what＇s more, it can overcome the problem of ambiguity. Both theoretical analysis and computer simulations show the effectiveness of the proposed method.
基金Project(2009318000046) supported by the Western Transport Technical Program of the Ministry of Transport,China
文摘A revised displacement discontinuity method(DDM) program is developed for the simulation of rock joint propagation and dilatancy analysis. The non-linear joint model used in the program adopts Barton-Bandis normal deformation model, Kulhaway shear deformation model and Mohr-Coulomb criterion. The joint propagation criterion is based on the equivalent stress intensity factor which can be obtained by regression analysis. The simulated rock joint propagation accords well with the existing knowledge. The closure and opening of joint is investigated by DDM, and it is shown that if the opening volume of propagated joint is larger than closure volume of the old joint, the joint dilatancy occurs. The dilatancy condition is mainly controlled by the normal stiffness of the rock joint. When the normal stiffness is larger than the critical value, joint dilatancy occurs. The critical normal stiffness of rock joint changes with the joint-load angle, and joint dilatancy is most possible to occur at 30°.
文摘The higher order displacement discontinuity method(HODDM) utilizing special crack tip elements has been used in the solution of linear elastic fracture mechanics(LEFM) problems. The paper has selected several example problems from the fracture mechanics literature(with available analytical solutions) including center slant crack in an infinite and finite body, single and double edge cracks, cracks emanating from a circular hole. The numerical values of Mode Ⅰ and Mode Ⅱ SIFs for these problems using HODDM are in excellent agreement with analytical results(reaching up to 0.001% deviation from their analytical results). The HODDM is also compared with the XFEM and a modified XFEM results. The results show that the HODDM needs a considerably lower computational effort(with less than 400 nodes) than the XFEM and the modified XFEM(which needs more than 10000 nodes) to reach a much higher accuracy. The proposed HODDM offers higher accuracy and lower computation effort for a wide range of problems in LEFM.
文摘The present research is focused on the numerical crack coalescence analysis of the micro-cracks and cracks produced during the cutting action of TBM disc cutters. The linear elastic fracture mechanics(LEFM) concepts and the maximum tangential stress criterion are used to investigate the micro crack propagation and its direction underneath the excavating discs. A higher order displacement discontinuity method with quadratic displacement discontinuity elements is used to estimate the stress intensity factors near the crack tips. Rock cutting mechanisms under single and double type discs are simulated by the proposed numerical method.The main purposes of the present modeling are to simulate the chip formation process of indented rocks by single and double discs.The effects of specific disc parameters(except speed) on the thrust force Ft, the rolling force Fr, and the specific energy ES are investigated. It has been shown that the specific energy(energy required to cut through a unit volume of rock) of the double disc is less than that of the single disc. Crack propagation in rocks under disc cutters is numerically modeled and the optimum ratio of disc spacing S to penetration depth Pd(i.e. S/Pd ratio) of about 10 is obtained, which is in good agreement with the theoretical and experimental results cited in the literature.
文摘In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat conduction techniques, this calibration approach does not require explicit input of the probe locations, thermophysical properties of the host material and temperature sensor parameters related to thermal contact resistance, sensor capacitance and conductive lead losses. All those parameters and properties are inherently contained in the calibration framework in terms of Volterra integral equation of the first kind. The Laplace transform technique is applied and the frequency domain manipulations of the heat equation are performed for deriving the calibration integral equation. Due to the ill-posed nature, regularization is required for the inverse heat conduction problem, a future-time method or singular value decomposition (SVD) can be used for stabilizing the ill-posed Volterra integral equation of the first kind.
文摘The calculation precision and convergence speed of streamline strip element method are increased by (using) the method whose initial value of the exit lateral displacement is determined with strip element variation method, and the accurate tension lateral distribution model is adopted based on the original third power spline function streamline strip element method. The basic theory of the strip element method is developed. The calculated results by the improved streamline strip element method and the original streamline strip element method are compared with the measured results, showing that the calculated results of the improved method are in good agreement with the measured results.
文摘In engineering practice simplified methods are essential to the seismic design of embedded earth retaining walls,as fullydynamic numerical analyses are costly,time-consuming and require specific expertise.Recently developed pseudostatic methods provide earth stresses and internal forces,even in those cases in which the strength of the soil surrounding the structure is not entirely mobilised.Semiempirical correlations or Newmark sliding block method provide an estimate of earthquake-induced permanent displacements.However,the use of these methods is hindered by uncertainties in the evaluation of a few input parameters,affecting the reliability of the methods.This study uses 1 D site response analyses and 2 D fully-dynamic finite element analyses to show that simplified methods can provide a reasonable estimate of the maximum bending moment and permanent displacements for stiff cantilever walls embedded in uniform sand,providing that a few input parameters are evaluated through semiempirical correlations and a simple 1 D site response analysis.
文摘为研究农业机械与水田壤土间的相互作用,需获取水田壤土的物理及接触参数。结合物理堆积试验,以休止角作为响应值,采用离散元法(DEM)并选取Hertz-Mindlin with JKR(Johnson-Kendall-Roberts)接触模型对长江中游地区水田壤土展开参数标定研究。首先,通过物理堆积试验获取了壤土休止角(AoR)与含水率间的定量关系,由不同含水率土壤的堆积结果筛分出4种代表性堆积形态,由于水田壤土堆积体轮廓外形比较独特,因此仅对其左右两侧轮廓采用三次多项式进行局部拟合,计算其休止角。以长江中游地区水田壤土成因和预试验为依据来确定其离散元模型中9个参数的高低水平值,通过Plackett-Burman试验设计进行方差分析,发现壤土剪切模量、壤土间动摩擦因数、壤土与不锈钢间静摩擦因数和JKR表面能对AoR影响明显。然后,采用基于响应面法(RSM)原理的Box-Behnken试验设计(BBD)建立了AoR与4个显著性参数间的二次多项式回归模型。依据二次多项式回归模型对目标响应进行预测,得到最优参数组合。以此为基础对壤土AoR进行离散元仿真,AoR数值计算结果(45.4°)与试验结果(44.6°)相对误差为1.79%。最后,选取含水率分别为44.4%、48.7%的壤土进行堆积角仿真模拟,计算结果与堆积试验相对误差分别为2.8%、7.14%。研究表明:回归模型可以根据壤土含水率或AoR预测长江中游地区水田壤土的相关本征参数和接触参数。
