A rigorous back analysis of shear strength parameters of landslide slip was presented. Kinematical element method was adopted to determine factor of safety and critical failure surface, which overcomes the disadvantag...A rigorous back analysis of shear strength parameters of landslide slip was presented. Kinematical element method was adopted to determine factor of safety and critical failure surface, which overcomes the disadvantage of limit equilibrium method. The theoretical relationship between the combination of shear strength parameters and stability state was studied. The results show that the location of critical slip surface, F/tan f and F/c depend only on the value of c/tan f. The failure surface moves towards the inside of slope as c/tan f increases. According to the information involving factor of safety and critical failure surface in a specific cross-section, strength parameters can be back calculated based on the above findings. Three examples were given for demonstrating the validity of the present method. The shear strength parameters obtained by back analysis are almost consistent with their correct solutions or test results.展开更多
Currently, virtual assembly technology has attracted increasing attention due to considerations of solving assembly problems in virtual environment before actual assembly in manufactory. Previous studies on kinematic ...Currently, virtual assembly technology has attracted increasing attention due to considerations of solving assembly problems in virtual environment before actual assembly in manufactory. Previous studies on kinematic analysis of mechanism only aim at analyzing motion law of single mechanism, but can not simulate the multi-mechanisms motion process at the same time, let alone simulating the automatic assembly process of products in a whole assembly workshop. In order to simulate the assembly process of products in an assembly workshop and provide effective data for analyzing mechanical performance after finishing assembly simulation in virtual environment, this study investigates the kinematics analysis of mechanisms based on virtual assembly. Firstly, in view of the same function of the kinematic pairs and the assembly constraints on restricting the motion of components （subassembly or part）, the method of identifying kinematic pairs automatically based on assembly constraints is presented. The information of kinematic pairs can be obtained through calculating the constraint degree of the assembly constraints. Secondly, the incidence matrix eliminating element method is proposed in order to search the information and establish the models of mechanisms automatically after finishing assembly simulation in virtual environment. Both methods have important significance for reducing the workload of pretreatment and promoting the level of automation of kinematics analysis. Finally, the method of kinematics analysis of mechanisms is presented. Based on Descartes coordinates, three types of kinematics equations are formed. The parameters, like displacement, velocity, and acceleration, can be obtained by solving these equations. All these data are important to analyze mechanical performance. All the methods are implemented and validated in the prototype system virtual assembly process planning（VAPP）. The mechanism models are established and simulated in the VAPP system, and the result curves are shown accurately. T展开更多
For the purpose of investigating the influence of metastable olivine（MO） phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma...For the purpose of investigating the influence of metastable olivine（MO） phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near～410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab：（i） slab pull force exerts above a depth of～230 km;（ii） MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;（iii） as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;（iv） a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.展开更多
Kinematic models compute the temperature distribution by prescribing a constant convergent velocity for the subducting slab,resulting in an artificial velocity discontinuity,which may accelerate the heating of subduct...Kinematic models compute the temperature distribution by prescribing a constant convergent velocity for the subducting slab,resulting in an artificial velocity discontinuity,which may accelerate the heating of subducting slabs.For the purpose to moderate the influence of such artificial discontinuity, we construct a 2D thermal model for subduction zones with a velocity boundary layer,within which the velocities decrease linearly with the distance from the interfaces of slabs.Temperatures are calculated展开更多
This paper proposes a numerical solution method for upper bound shakedown analysis of perfectly elasto-plastic thin plates by employing the C^(1) natural element method.Based on the Koiter’s theorem and von Mises yie...This paper proposes a numerical solution method for upper bound shakedown analysis of perfectly elasto-plastic thin plates by employing the C^(1) natural element method.Based on the Koiter’s theorem and von Mises yield criterion,the nonlinear mathematical programming formulation for upper bound shakedown analysis of thin plates is established.In this formulation,the trail function of residual displacement increment is approximated by using the C^(1) shape functions,the plastic incompressibility condition is satisfied by introducing a constant matrix in the objective function,and the time integration is resolved by using the Konig’s technique.Meanwhile,the objective function is linearized by distinguishing the non-plastic integral points from the plastic integral points and revising the objective function and associated equality constraints at each iteration.Finally,the upper bound shakedown load multipliers of thin plates are obtained by direct iterative and monotone convergence processes.Several benchmark examples verify the good precision and fast convergence of this proposed method.展开更多
基金Project(51174228)supported by the National Natural Science Foundation of ChinaProject(CX2012B069)supported by Hunan Provincial Innovation Foundation for PostgraduateProject(201003)supported by Transportation Science and Technology Projects of Hunan Province,China
文摘A rigorous back analysis of shear strength parameters of landslide slip was presented. Kinematical element method was adopted to determine factor of safety and critical failure surface, which overcomes the disadvantage of limit equilibrium method. The theoretical relationship between the combination of shear strength parameters and stability state was studied. The results show that the location of critical slip surface, F/tan f and F/c depend only on the value of c/tan f. The failure surface moves towards the inside of slope as c/tan f increases. According to the information involving factor of safety and critical failure surface in a specific cross-section, strength parameters can be back calculated based on the above findings. Three examples were given for demonstrating the validity of the present method. The shear strength parameters obtained by back analysis are almost consistent with their correct solutions or test results.
基金supported by National Natural Science Foundation of China (Grant No. 50805009)National Defense Pre-Research Fund of China during the 11th Five-Year Plan Period (Grant No. 51318010205)
文摘Currently, virtual assembly technology has attracted increasing attention due to considerations of solving assembly problems in virtual environment before actual assembly in manufactory. Previous studies on kinematic analysis of mechanism only aim at analyzing motion law of single mechanism, but can not simulate the multi-mechanisms motion process at the same time, let alone simulating the automatic assembly process of products in a whole assembly workshop. In order to simulate the assembly process of products in an assembly workshop and provide effective data for analyzing mechanical performance after finishing assembly simulation in virtual environment, this study investigates the kinematics analysis of mechanisms based on virtual assembly. Firstly, in view of the same function of the kinematic pairs and the assembly constraints on restricting the motion of components （subassembly or part）, the method of identifying kinematic pairs automatically based on assembly constraints is presented. The information of kinematic pairs can be obtained through calculating the constraint degree of the assembly constraints. Secondly, the incidence matrix eliminating element method is proposed in order to search the information and establish the models of mechanisms automatically after finishing assembly simulation in virtual environment. Both methods have important significance for reducing the workload of pretreatment and promoting the level of automation of kinematics analysis. Finally, the method of kinematics analysis of mechanisms is presented. Based on Descartes coordinates, three types of kinematics equations are formed. The parameters, like displacement, velocity, and acceleration, can be obtained by solving these equations. All these data are important to analyze mechanical performance. All the methods are implemented and validated in the prototype system virtual assembly process planning（VAPP）. The mechanism models are established and simulated in the VAPP system, and the result curves are shown accurately. T
基金supported by the National Natural Science Foundation of China(Nos.40574047 and 40628004)
文摘For the purpose of investigating the influence of metastable olivine（MO） phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near～410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab：（i） slab pull force exerts above a depth of～230 km;（ii） MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;（iii） as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;（iv） a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.
文摘Kinematic models compute the temperature distribution by prescribing a constant convergent velocity for the subducting slab,resulting in an artificial velocity discontinuity,which may accelerate the heating of subducting slabs.For the purpose to moderate the influence of such artificial discontinuity, we construct a 2D thermal model for subduction zones with a velocity boundary layer,within which the velocities decrease linearly with the distance from the interfaces of slabs.Temperatures are calculated
基金supported by the Chinese Postdoctoral Science Foundation(2013M540934)supported by the National Key Research and Development Program of China(2016YFC0801905,2017YFF0210704).
文摘This paper proposes a numerical solution method for upper bound shakedown analysis of perfectly elasto-plastic thin plates by employing the C^(1) natural element method.Based on the Koiter’s theorem and von Mises yield criterion,the nonlinear mathematical programming formulation for upper bound shakedown analysis of thin plates is established.In this formulation,the trail function of residual displacement increment is approximated by using the C^(1) shape functions,the plastic incompressibility condition is satisfied by introducing a constant matrix in the objective function,and the time integration is resolved by using the Konig’s technique.Meanwhile,the objective function is linearized by distinguishing the non-plastic integral points from the plastic integral points and revising the objective function and associated equality constraints at each iteration.Finally,the upper bound shakedown load multipliers of thin plates are obtained by direct iterative and monotone convergence processes.Several benchmark examples verify the good precision and fast convergence of this proposed method.
