This paper focuses on eliminating the unphysical negative susceptibility which appears when magnetic field is at unsaturated excitation level and reduces from extremity of the hysteresis loop in one-dimension coupled ...This paper focuses on eliminating the unphysical negative susceptibility which appears when magnetic field is at unsaturated excitation level and reduces from extremity of the hysteresis loop in one-dimension coupled hysteresis model. The domain flexing function c(H) is used to replace the domain flexing constant c in one-dimension coupled hysteresis model. The feasibility and rationality of proposed modification are convinced by comparing the magnetization and magnetostriction curves with experimental data and another typical modification results. The effects of pre-stress and temperature on magnetic-elastic-thermal coupling property and hysteresis behavior are investigated.展开更多
The frequency-dependent dynamic effective properties (phase velocity, attenuation and elastic modulus) of porous materials are studied numerically. The coherent plane longitudinal and shear wave equations, which are o...The frequency-dependent dynamic effective properties (phase velocity, attenuation and elastic modulus) of porous materials are studied numerically. The coherent plane longitudinal and shear wave equations, which are obtained by averaging on the multiple scattering fields, are used to evaluate the frequency-dependent dynamic effective properties of a porous material. It is found that the prediction of the dynamic effective properties includes the size effects of voids which are not included in most prediction of the traditional static effective properties. The prediction of the dynamic effective elastic modulus at a relatively low frequency range is compared with that of the traditional static effective elastic modulus, and the dynamic effective elastic modulus is found to be very close to the Hashin-Shtrikman upper bound.展开更多
On the basis of ANSYS finite element model(FEM) software, the deep-level rockburst in Fuxin coalfield was simulated numerically. Based on Haizhou Mine and Wulong Mine as two typical deep-level rockburst examples in Fu...On the basis of ANSYS finite element model(FEM) software, the deep-level rockburst in Fuxin coalfield was simulated numerically. Based on Haizhou Mine and Wulong Mine as two typical deep-level rockburst examples in Fuxin coalfield, the rules and characteristics of the deep-level rockburst were analyzed. And the models were established. For Haizhou mine, the relationship between mining distance and rockburst was presented when 100, 300, 600 m were mined in 3313 working face. When 300 m were mined, the rockburst began to emerge. When 600 m were mined, the rockburstwas the most possible to happen and the compression stress of the working face reached to the maximum value. The effect of tectonic stress on synclinal axis is also a key factor to rockburst occurrence. This was verified by the rockburst happened when 496 m were mined. For Wulong mine, based on the 311 working face as an example,the contours of Y stress in the roof and floor were obtained when the mining distance were 100, 200, 300 and 400 m. When 100 and 400 m were mined, the high stress concentration regions occurred in the front of working face. This shows the rockburst is easy to happen. It is confirmed by the rockburst when 91m were mined in 311 working plane.The above indicates that the numerical simulation has instructive rule to study the deep-Ievel rockburst in Fuxin coalfield.展开更多
A level set method, the TVD scheme of second-order upwind procedure coupled with flux-limiter, and SIMPLE algorithm were incorporated to simulate the flow and interfacial motion of immiscible two-fluids with large den...A level set method, the TVD scheme of second-order upwind procedure coupled with flux-limiter, and SIMPLE algorithm were incorporated to simulate the flow and interfacial motion of immiscible two-fluids with large density ratio and viscosity ratios, large topology distortion and surface tension. As a numerical example axi-symmetric rising bubbles were investigated. It is found that the method is numerically stable and has good convergence property and the results are in good agreement with other works.展开更多
A mechanical model for strain softening pillar is proposed considering the characteristics of progressive shear failure and strain localization. The pillar undergoes elastic, strain softening and slabbing stages. In t...A mechanical model for strain softening pillar is proposed considering the characteristics of progressive shear failure and strain localization. The pillar undergoes elastic, strain softening and slabbing stages. In the elastic stage, vertical compressive stress and deformation at upper end of pillar are uniform, while in the strain softening stage there appears nonuniform due to occurrence of shear bands, leading to the decrease of load-carrying capacity.In addition, the size of failure zone increases in the strain softening stage and reaches its maximum value when slabbing begins. In the latter two stages, the size of elastic core always decreases. In the slabbing stage, the size of failure zone remains a constant and the pillar becomes thinner. Total deformation of the pillar is derived by linearly elastic Hooke's law and gradient-dependent plasticity where thickness of localization band is determined according to the characteristic length. Post-peak stiffness is proposed according to analytical solution of averaged compressive stressaverage deformation curve. Instability criterion of the pillar and roof strata system is proposed analytically using instability condition given by Salamon. It is found that the constitutive parameters of material of pillar, the geometrical size of pillar and the number of shear bands influence the stability of the system; stress gradient controls the starting time of slabbing, however it has no influence on the post-peak stiffness of the pillar.展开更多
A method for calculation of temperature distribution in adiabatic shear band is proposed in terms of gradient-dependent plasticity where the characteristic length describes the interactions and interplaying among micr...A method for calculation of temperature distribution in adiabatic shear band is proposed in terms of gradient-dependent plasticity where the characteristic length describes the interactions and interplaying among microstructures. First, the increment of the plastic shear strain distribution in adiabatic shear band is obtained based on gradient-dependent plasticity. Then, the plastic work distribution is derived according to the current flow shear stress and the obtained increment of plastic shear strain distribution. In the light of the well-known assumption that 90% of plastic work is converted into the heat resulting in increase in temperature in adiabatic shear band, the increment of the temperature distribution is presented. Next, the average temperature increment in the shear band is calculated to compute the change in flow shear stress due to the thermal softening effect. After the actual flow shear stress considering the thermal softening effect is obtained according to the Johnson-Cook constitutive relation, the increment of the plastic shear strain distribution, the plastic work and the temperature in the next time step are recalculated until the total time is consumed. Summing the temperature distribution leads to rise in the total temperature distribution. The present calculated maximum temperature in adiabatic shear band in titanium agrees with the experimental observations. Moreover, the temperature profiles for different flow shear stresses are qualitatively consistent with experimental and numerical results. Effects of some related parameters on the temperature distribution are also predicted.展开更多
Fault rockburst is treated as a strain localization problem under dynamic loading condition considering strain gradient and strain rate. As a kind of dynamic fracture phenomena, rockburst has characteristics of strain...Fault rockburst is treated as a strain localization problem under dynamic loading condition considering strain gradient and strain rate. As a kind of dynamic fracture phenomena, rockburst has characteristics of strain localization, which is considered as a one-dimensional shear problem subjected to normal compressive stress and tangential shear stress. The constitutive relation of rock material is bilinear (elastic and strain softening) and sensitive to shear strain rate. The solutions proposed based on gradientdependent plasticity show that intense plastic strain is concentrated in fault band and the thickness of the band depends on the characteristic length of rock material. The post-peak stiffness of the fault band was determined according to the constitutive parameters of rock material and shear strain rate. Fault band undergoing strain softening and elastic rock mass outside the band constitute a system and the instability criterion of the system was proposed based on energy theory. The criterion depends on the constitutive relation of rock material, the structural size and the strain rate. The static result regardless of the strain rate is the special case of the present analytical solution. High strain rate can lead to instability of the system.展开更多
The inclination angle ofsbear band is analyzed considering heterogeneity of rock material when a single shear hand is formed in the center of specimen under triaxial compression. The analytical solution of post-peak a...The inclination angle ofsbear band is analyzed considering heterogeneity of rock material when a single shear hand is formed in the center of specimen under triaxial compression. The analytical solution of post-peak axial stress-axial strain curve is deduced using the assumption that the total post-peak deformation is composed of entire uniform elastic deformation and localized shear plastic deformation dependent on the thickness of shear band. The obtained solution shows that the post-peak stiffness is related to the inclination angle of shear band, confining pressure, thickness of shear hand and elastic modulus, etc. Using the solution, the expression for the inclination angle of shear hand can be presented easily and it is dependent on constitutive parameters of rock material and geometry parameters of rock specimen. Larger dilation angle or loading rate leads to increment of the inclination angle. In addition, the inclination angle increases with the thickness of the shear band, which cannot be explained or forecasted by other existing solutions, such as Coulomb inclination, Roscoe inclination and Arthur inclination, etc.展开更多
The paper concerns the issue of size law, localized deformation and dilation or compaction due to shearlocalization. It is assumed that the shear localization initiates at the peak shear stress in the form of single s...The paper concerns the issue of size law, localized deformation and dilation or compaction due to shearlocalization. It is assumed that the shear localization initiates at the peak shear stress in the form of single shear band, andbased on gradient-dependent plasticity, an analytical solution on size effect or snap-back is obtained. The results show thatthe post peak response becomes steeper and even exhibits snap-back with increasing of length. For small specimen, therelative shear displacement when specimen failure occurs is lower than that of larger specimen and the shearstress-relative displacement curve becomes steeper. The theoretical solution on non-uniformity of strains in shear band isobtained and evolution of the relative shear displacement is represented. By resorting to the linear relation between localplastic shear strain and local plastic volumetric strain, the dilation and compaction within shear band are analyzed.Relation between apparent shear strain and apparent normal strain and relation between shear displacement and verticaldisplacement are established.展开更多
The paper was numerically focused on investigation of deformation, failure and instability of shear band-surrounding elastic reck system in plane strain direct shear test considering shear dilatancy according to fast ...The paper was numerically focused on investigation of deformation, failure and instability of shear band-surrounding elastic reck system in plane strain direct shear test considering shear dilatancy according to fast lagrangian analysis of continua (FLAC). The adopted failure criterion was a composite Mohr-Coulomb criterion with tension cut-off and post-peak constitutive relation of rock, i.e. linear strain-softening. Numerical results show that dilation angle affects the responses of elements, the number and the position of yielded elements. Increasing dilation angle results in higher load-carrying capability of elements, higher deformation or strain corresponding to peak stress, less brittle post-peak stress-deformation curve. Strain-hardening behavior can occur if dilation angle is high.Therefore, shear band-elastic reck body system tends to be stable and rock burst does not occur as dilation angle is increased. Moreover, the number of yielded elements is increased with dilation angle increase and two' parallel plastic zones initially generated in the middle of two loading ends of specimen no longer remain collinear, reflecting increase in deformation resistant of the system. Evolution of volumetric strain rate was investigated based on three-dimensional curved surface diagrams. Approximately, volumetric strain rate concentration regions coincide with plastic zones. Volumetric strain rate in yielded elements is very evident. However, in elastic zones volumetric strain has not been nearly changed throughout the numerical tests.展开更多
The capacity of energy absorption by fault bands after rock burst was calculated quantitatively according to shear stressshear deformation curves considering the interactions and interplaying among microstructures due...The capacity of energy absorption by fault bands after rock burst was calculated quantitatively according to shear stressshear deformation curves considering the interactions and interplaying among microstructures due to the heterogeneity of strain softening rock materials. The post-peak stiffness of rock specimens subjected to direct shear was derived strictly based on gradientdependent plasticity, which can not be obtained from the classical elastoplastic theory. Analytical solutions for the dissipated energy of rock burst were proposed whether the slope of the post-peak shear stress-shear deformation curve is positive or not. The analytical solutions show that shear stress level, confining pressure, shear strength, brittleness, strain rate and heterogeneity of rock materials have important influence on the dissipated energy. The larger value of the dissipated energy means that the capacity of energy dissipation in the form of shear bands is superior and a lower magnitude of rock burst is expected under the condition of the same work done by external shear force. The possibility of rock burst is reduced for a lower softening modulus or a larger thickness of shear bands.展开更多
The problem of adaptive generalized predictive control which consists of output prediction errors for a class of switched systems is studied. The switching law is determined by the output predictive errors of a finite...The problem of adaptive generalized predictive control which consists of output prediction errors for a class of switched systems is studied. The switching law is determined by the output predictive errors of a finite number of subsystems. For the single subsystem and multiple subsystems cases, it is proved that the given direct algorithm of generalized predictive control guarantees the global convergence of the system. This algorithm overcomes the inherent drawbacks of the slow convergence and large transient errors for the conventional adaptive control.展开更多
The suboptimal control program via memoryless state feedback strategies for LQ differential games with multiple players is studied in this paper. Sufficient conditions for the existence of the suboptimal strategies fo...The suboptimal control program via memoryless state feedback strategies for LQ differential games with multiple players is studied in this paper. Sufficient conditions for the existence of the suboptimal strategies for LQ differential games are presented. It is shown that the suboptimal strategies of LQ differential games are associated with a coupled algebraic Riccati inequality. Furthermore, the problem of designing suboptimal strategies is considered. A non-convex optimization problem with BMI constrains is formulated to design the suboptimal strategies which minimizes the performance indices of the closed-loop LQ differential games and can be solved by using LMI Toolbox of MATLAB. An example is given to illustrate the proposed results.展开更多
Stability of liquid-gas coexistence equilibrium in a relaxation model for isothermal phase transition in a sealed one-dimensional tube was discussed. With matched asymptotic expansion, a linear system for first order ...Stability of liquid-gas coexistence equilibrium in a relaxation model for isothermal phase transition in a sealed one-dimensional tube was discussed. With matched asymptotic expansion, a linear system for first order perturbations was derived formally. By solving this system analytically, it is shown that small initial perturbations are damped out in general; yet they may maintain at certain level for special cases.Numerical evidence is presented. This manifests the regularization effects of relaxation.展开更多
The fluid flow induced by light-density, low-stiffness structures was treated as inviscid, incompressible irrotational and steady plane flow. On the basis of the dipole configuration method, a singularity distribution...The fluid flow induced by light-density, low-stiffness structures was treated as inviscid, incompressible irrotational and steady plane flow. On the basis of the dipole configuration method, a singularity distribution method of distributing sources/sinks and dipoles on interfaces of the structure and fluid was developed to solve the problem of fluid flow induced by the vibration of common structures, such as columns and columns with fins,deduce the expression of kinetic energy of the fluid flow, and obtain the added mass finally.The calculational instances with analytical solutions prove the reliability of this method.展开更多
Lyapunov's second method was used to study the nonlinear stability of parallel shear flows for stress-free boundaries. By introducing an energy functional, it was shown that the plane Couette and plane Poiseuille ...Lyapunov's second method was used to study the nonlinear stability of parallel shear flows for stress-free boundaries. By introducing an energy functional, it was shown that the plane Couette and plane Poiseuille flows are conditionally and asymptotically stable for all Reynolds numbers. In particular, to two-dimensional perturbations, by defining new energy functionals the unconditional stability of the basic flows was proved.展开更多
The hemodynamic mechanism of rolling manipulation (RM) of traditional Chinese medical massage (TCMM) is investigated. An axisymmetrical nonlinear model and an arbitrary Lagrangian-Eulerian finite element method (ALE-F...The hemodynamic mechanism of rolling manipulation (RM) of traditional Chinese medical massage (TCMM) is investigated. An axisymmetrical nonlinear model and an arbitrary Lagrangian-Eulerian finite element method (ALE-FEM) with rezoning algorithm were introduced to study the viscous flow through an axisymmetrical rigid tube with axially moving stenosis to simulate the rolling manipulation. Flow rate and wall shear stress were obtained by solving complete Navier-Stokes equations numerically. The numerical results show that the stenosis moving frequency, namely the frequency of rolling manipulation, has great effect on the disturbance of flow and the wall shear stress. The stenosis coefficient, which characterizes the severity of the stenosis, another adjustable parameter in rolling manipulation, also shows the significant effect on flow rate and wall shear stress. These numerical results may provide some data that can be taken into consideration when massage is used in clinic.展开更多
We propose a simple swarm model to study collective behaviour of a group of mobile autonomous agents interacting through a long range attraction and short range repulsion function. It is shown that the individuals (ag...We propose a simple swarm model to study collective behaviour of a group of mobile autonomous agents interacting through a long range attraction and short range repulsion function. It is shown that the individuals (agents) will aggregate and eventually form a cohesive cluster of finite size around the swarm centre in a finite time, and the size depends only on the parameters of the swarm model. Furthermore, it is also shown that all the individuals will converge to equilibrium positions of the swarm model, and thus the configuration of the swarm converges to a constant constellation. Numerical simulations are also worked out to illustrate the analytical results.展开更多
Although the coefficient of restitution was originally thought to be only a material property, the coefficient of restitution also depends upon initial conditions as well as on the frictional effect for oblique collis...Although the coefficient of restitution was originally thought to be only a material property, the coefficient of restitution also depends upon initial conditions as well as on the frictional effect for oblique collisions. The objective of this paper is to demonstrate a method for obtaining the coefficient of restitution for oblique collisions and thereby to provide a theoretical guide for collision experiments. In this paper, we derive expressions for the energetic coefficient of restitution ( e* ) based on general normal contact deformation law, by which the value of e* can be obtained according to the initial conditions. An example shows that the results calculated by the derived expressions are reasonable.展开更多
This paper presents the applications of digital image correlation technique to the mesoscopic damage and fracture study of some granular based composite materials including steelfiber reinforced concrete, sandstone an...This paper presents the applications of digital image correlation technique to the mesoscopic damage and fracture study of some granular based composite materials including steelfiber reinforced concrete, sandstone and crystal-polymer composite. The deformation fields of the composite materials resulted from stress localization were obtained by the correlation computation of the surface images with loading steps and thus the related damage prediction and fracture parameters were evaluated. The correlation searching could be performed either directly based on the gray levels of the digital images or from the wavelet transform (WT) coefficients of the transform spectrum. The latter was developed by the authors and showed higher resolution and sensitivity to the singularity detection. Because the displacement components came from the rough surfaces of the composite materials without any coats of gratings or fringes of optical interferometry, both surface profiles and the deformation fields of the composites were visualized which was helpful to compare each other to analyze the damage of those heterogeneous materials.展开更多
基金Project supported by the Fund of Natural Science Foundation of China （Nos. 10972094, 11032006, 11121202 and 11202087） and the Fundamental Research Funds for the Central Universities （No. lzujbky-2011-6） and Specialized Research Fund for the Doctoral Program of Higher Education （No. 20110211120027）.
文摘This paper focuses on eliminating the unphysical negative susceptibility which appears when magnetic field is at unsaturated excitation level and reduces from extremity of the hysteresis loop in one-dimension coupled hysteresis model. The domain flexing function c(H) is used to replace the domain flexing constant c in one-dimension coupled hysteresis model. The feasibility and rationality of proposed modification are convinced by comparing the magnetization and magnetostriction curves with experimental data and another typical modification results. The effects of pre-stress and temperature on magnetic-elastic-thermal coupling property and hysteresis behavior are investigated.
基金国家自然科学基金，the Talent Foundation of the University of Sciences and Technology Beijing
文摘The frequency-dependent dynamic effective properties (phase velocity, attenuation and elastic modulus) of porous materials are studied numerically. The coherent plane longitudinal and shear wave equations, which are obtained by averaging on the multiple scattering fields, are used to evaluate the frequency-dependent dynamic effective properties of a porous material. It is found that the prediction of the dynamic effective properties includes the size effects of voids which are not included in most prediction of the traditional static effective properties. The prediction of the dynamic effective elastic modulus at a relatively low frequency range is compared with that of the traditional static effective elastic modulus, and the dynamic effective elastic modulus is found to be very close to the Hashin-Shtrikman upper bound.
基金Supported by National Nature Science Foundation of China (50490275)
文摘On the basis of ANSYS finite element model(FEM) software, the deep-level rockburst in Fuxin coalfield was simulated numerically. Based on Haizhou Mine and Wulong Mine as two typical deep-level rockburst examples in Fuxin coalfield, the rules and characteristics of the deep-level rockburst were analyzed. And the models were established. For Haizhou mine, the relationship between mining distance and rockburst was presented when 100, 300, 600 m were mined in 3313 working face. When 300 m were mined, the rockburst began to emerge. When 600 m were mined, the rockburstwas the most possible to happen and the compression stress of the working face reached to the maximum value. The effect of tectonic stress on synclinal axis is also a key factor to rockburst occurrence. This was verified by the rockburst happened when 496 m were mined. For Wulong mine, based on the 311 working face as an example,the contours of Y stress in the roof and floor were obtained when the mining distance were 100, 200, 300 and 400 m. When 100 and 400 m were mined, the high stress concentration regions occurred in the front of working face. This shows the rockburst is easy to happen. It is confirmed by the rockburst when 91m were mined in 311 working plane.The above indicates that the numerical simulation has instructive rule to study the deep-Ievel rockburst in Fuxin coalfield.
文摘A level set method, the TVD scheme of second-order upwind procedure coupled with flux-limiter, and SIMPLE algorithm were incorporated to simulate the flow and interfacial motion of immiscible two-fluids with large density ratio and viscosity ratios, large topology distortion and surface tension. As a numerical example axi-symmetric rising bubbles were investigated. It is found that the method is numerically stable and has good convergence property and the results are in good agreement with other works.
文摘A mechanical model for strain softening pillar is proposed considering the characteristics of progressive shear failure and strain localization. The pillar undergoes elastic, strain softening and slabbing stages. In the elastic stage, vertical compressive stress and deformation at upper end of pillar are uniform, while in the strain softening stage there appears nonuniform due to occurrence of shear bands, leading to the decrease of load-carrying capacity.In addition, the size of failure zone increases in the strain softening stage and reaches its maximum value when slabbing begins. In the latter two stages, the size of elastic core always decreases. In the slabbing stage, the size of failure zone remains a constant and the pillar becomes thinner. Total deformation of the pillar is derived by linearly elastic Hooke's law and gradient-dependent plasticity where thickness of localization band is determined according to the characteristic length. Post-peak stiffness is proposed according to analytical solution of averaged compressive stressaverage deformation curve. Instability criterion of the pillar and roof strata system is proposed analytically using instability condition given by Salamon. It is found that the constitutive parameters of material of pillar, the geometrical size of pillar and the number of shear bands influence the stability of the system; stress gradient controls the starting time of slabbing, however it has no influence on the post-peak stiffness of the pillar.
文摘A method for calculation of temperature distribution in adiabatic shear band is proposed in terms of gradient-dependent plasticity where the characteristic length describes the interactions and interplaying among microstructures. First, the increment of the plastic shear strain distribution in adiabatic shear band is obtained based on gradient-dependent plasticity. Then, the plastic work distribution is derived according to the current flow shear stress and the obtained increment of plastic shear strain distribution. In the light of the well-known assumption that 90% of plastic work is converted into the heat resulting in increase in temperature in adiabatic shear band, the increment of the temperature distribution is presented. Next, the average temperature increment in the shear band is calculated to compute the change in flow shear stress due to the thermal softening effect. After the actual flow shear stress considering the thermal softening effect is obtained according to the Johnson-Cook constitutive relation, the increment of the plastic shear strain distribution, the plastic work and the temperature in the next time step are recalculated until the total time is consumed. Summing the temperature distribution leads to rise in the total temperature distribution. The present calculated maximum temperature in adiabatic shear band in titanium agrees with the experimental observations. Moreover, the temperature profiles for different flow shear stresses are qualitatively consistent with experimental and numerical results. Effects of some related parameters on the temperature distribution are also predicted.
文摘Fault rockburst is treated as a strain localization problem under dynamic loading condition considering strain gradient and strain rate. As a kind of dynamic fracture phenomena, rockburst has characteristics of strain localization, which is considered as a one-dimensional shear problem subjected to normal compressive stress and tangential shear stress. The constitutive relation of rock material is bilinear (elastic and strain softening) and sensitive to shear strain rate. The solutions proposed based on gradientdependent plasticity show that intense plastic strain is concentrated in fault band and the thickness of the band depends on the characteristic length of rock material. The post-peak stiffness of the fault band was determined according to the constitutive parameters of rock material and shear strain rate. Fault band undergoing strain softening and elastic rock mass outside the band constitute a system and the instability criterion of the system was proposed based on energy theory. The criterion depends on the constitutive relation of rock material, the structural size and the strain rate. The static result regardless of the strain rate is the special case of the present analytical solution. High strain rate can lead to instability of the system.
文摘The inclination angle ofsbear band is analyzed considering heterogeneity of rock material when a single shear hand is formed in the center of specimen under triaxial compression. The analytical solution of post-peak axial stress-axial strain curve is deduced using the assumption that the total post-peak deformation is composed of entire uniform elastic deformation and localized shear plastic deformation dependent on the thickness of shear band. The obtained solution shows that the post-peak stiffness is related to the inclination angle of shear band, confining pressure, thickness of shear hand and elastic modulus, etc. Using the solution, the expression for the inclination angle of shear hand can be presented easily and it is dependent on constitutive parameters of rock material and geometry parameters of rock specimen. Larger dilation angle or loading rate leads to increment of the inclination angle. In addition, the inclination angle increases with the thickness of the shear band, which cannot be explained or forecasted by other existing solutions, such as Coulomb inclination, Roscoe inclination and Arthur inclination, etc.
基金Supported by the National Natural Science Foundation of China(50309004)。
文摘The paper concerns the issue of size law, localized deformation and dilation or compaction due to shearlocalization. It is assumed that the shear localization initiates at the peak shear stress in the form of single shear band, andbased on gradient-dependent plasticity, an analytical solution on size effect or snap-back is obtained. The results show thatthe post peak response becomes steeper and even exhibits snap-back with increasing of length. For small specimen, therelative shear displacement when specimen failure occurs is lower than that of larger specimen and the shearstress-relative displacement curve becomes steeper. The theoretical solution on non-uniformity of strains in shear band isobtained and evolution of the relative shear displacement is represented. By resorting to the linear relation between localplastic shear strain and local plastic volumetric strain, the dilation and compaction within shear band are analyzed.Relation between apparent shear strain and apparent normal strain and relation between shear displacement and verticaldisplacement are established.
基金Supported by the National Natural Science Foundation of China(50309004)
文摘The paper was numerically focused on investigation of deformation, failure and instability of shear band-surrounding elastic reck system in plane strain direct shear test considering shear dilatancy according to fast lagrangian analysis of continua (FLAC). The adopted failure criterion was a composite Mohr-Coulomb criterion with tension cut-off and post-peak constitutive relation of rock, i.e. linear strain-softening. Numerical results show that dilation angle affects the responses of elements, the number and the position of yielded elements. Increasing dilation angle results in higher load-carrying capability of elements, higher deformation or strain corresponding to peak stress, less brittle post-peak stress-deformation curve. Strain-hardening behavior can occur if dilation angle is high.Therefore, shear band-elastic reck body system tends to be stable and rock burst does not occur as dilation angle is increased. Moreover, the number of yielded elements is increased with dilation angle increase and two' parallel plastic zones initially generated in the middle of two loading ends of specimen no longer remain collinear, reflecting increase in deformation resistant of the system. Evolution of volumetric strain rate was investigated based on three-dimensional curved surface diagrams. Approximately, volumetric strain rate concentration regions coincide with plastic zones. Volumetric strain rate in yielded elements is very evident. However, in elastic zones volumetric strain has not been nearly changed throughout the numerical tests.
文摘The capacity of energy absorption by fault bands after rock burst was calculated quantitatively according to shear stressshear deformation curves considering the interactions and interplaying among microstructures due to the heterogeneity of strain softening rock materials. The post-peak stiffness of rock specimens subjected to direct shear was derived strictly based on gradientdependent plasticity, which can not be obtained from the classical elastoplastic theory. Analytical solutions for the dissipated energy of rock burst were proposed whether the slope of the post-peak shear stress-shear deformation curve is positive or not. The analytical solutions show that shear stress level, confining pressure, shear strength, brittleness, strain rate and heterogeneity of rock materials have important influence on the dissipated energy. The larger value of the dissipated energy means that the capacity of energy dissipation in the form of shear bands is superior and a lower magnitude of rock burst is expected under the condition of the same work done by external shear force. The possibility of rock burst is reduced for a lower softening modulus or a larger thickness of shear bands.
文摘The problem of adaptive generalized predictive control which consists of output prediction errors for a class of switched systems is studied. The switching law is determined by the output predictive errors of a finite number of subsystems. For the single subsystem and multiple subsystems cases, it is proved that the given direct algorithm of generalized predictive control guarantees the global convergence of the system. This algorithm overcomes the inherent drawbacks of the slow convergence and large transient errors for the conventional adaptive control.
文摘The suboptimal control program via memoryless state feedback strategies for LQ differential games with multiple players is studied in this paper. Sufficient conditions for the existence of the suboptimal strategies for LQ differential games are presented. It is shown that the suboptimal strategies of LQ differential games are associated with a coupled algebraic Riccati inequality. Furthermore, the problem of designing suboptimal strategies is considered. A non-convex optimization problem with BMI constrains is formulated to design the suboptimal strategies which minimizes the performance indices of the closed-loop LQ differential games and can be solved by using LMI Toolbox of MATLAB. An example is given to illustrate the proposed results.
文摘Stability of liquid-gas coexistence equilibrium in a relaxation model for isothermal phase transition in a sealed one-dimensional tube was discussed. With matched asymptotic expansion, a linear system for first order perturbations was derived formally. By solving this system analytically, it is shown that small initial perturbations are damped out in general; yet they may maintain at certain level for special cases.Numerical evidence is presented. This manifests the regularization effects of relaxation.
文摘The fluid flow induced by light-density, low-stiffness structures was treated as inviscid, incompressible irrotational and steady plane flow. On the basis of the dipole configuration method, a singularity distribution method of distributing sources/sinks and dipoles on interfaces of the structure and fluid was developed to solve the problem of fluid flow induced by the vibration of common structures, such as columns and columns with fins,deduce the expression of kinetic energy of the fluid flow, and obtain the added mass finally.The calculational instances with analytical solutions prove the reliability of this method.
文摘Lyapunov's second method was used to study the nonlinear stability of parallel shear flows for stress-free boundaries. By introducing an energy functional, it was shown that the plane Couette and plane Poiseuille flows are conditionally and asymptotically stable for all Reynolds numbers. In particular, to two-dimensional perturbations, by defining new energy functionals the unconditional stability of the basic flows was proved.
文摘The hemodynamic mechanism of rolling manipulation (RM) of traditional Chinese medical massage (TCMM) is investigated. An axisymmetrical nonlinear model and an arbitrary Lagrangian-Eulerian finite element method (ALE-FEM) with rezoning algorithm were introduced to study the viscous flow through an axisymmetrical rigid tube with axially moving stenosis to simulate the rolling manipulation. Flow rate and wall shear stress were obtained by solving complete Navier-Stokes equations numerically. The numerical results show that the stenosis moving frequency, namely the frequency of rolling manipulation, has great effect on the disturbance of flow and the wall shear stress. The stenosis coefficient, which characterizes the severity of the stenosis, another adjustable parameter in rolling manipulation, also shows the significant effect on flow rate and wall shear stress. These numerical results may provide some data that can be taken into consideration when massage is used in clinic.
文摘We propose a simple swarm model to study collective behaviour of a group of mobile autonomous agents interacting through a long range attraction and short range repulsion function. It is shown that the individuals (agents) will aggregate and eventually form a cohesive cluster of finite size around the swarm centre in a finite time, and the size depends only on the parameters of the swarm model. Furthermore, it is also shown that all the individuals will converge to equilibrium positions of the swarm model, and thus the configuration of the swarm converges to a constant constellation. Numerical simulations are also worked out to illustrate the analytical results.
文摘Although the coefficient of restitution was originally thought to be only a material property, the coefficient of restitution also depends upon initial conditions as well as on the frictional effect for oblique collisions. The objective of this paper is to demonstrate a method for obtaining the coefficient of restitution for oblique collisions and thereby to provide a theoretical guide for collision experiments. In this paper, we derive expressions for the energetic coefficient of restitution ( e* ) based on general normal contact deformation law, by which the value of e* can be obtained according to the initial conditions. An example shows that the results calculated by the derived expressions are reasonable.
基金The project supported by the National Natural Science Foundation of China (10125211 and 10072002),the Scientific Committee of Yunnan Province for the Program of Steel Fiber Reinforced Concrete,and the Institute of Chemical Materials,CAEP at Mianyang
文摘This paper presents the applications of digital image correlation technique to the mesoscopic damage and fracture study of some granular based composite materials including steelfiber reinforced concrete, sandstone and crystal-polymer composite. The deformation fields of the composite materials resulted from stress localization were obtained by the correlation computation of the surface images with loading steps and thus the related damage prediction and fracture parameters were evaluated. The correlation searching could be performed either directly based on the gray levels of the digital images or from the wavelet transform (WT) coefficients of the transform spectrum. The latter was developed by the authors and showed higher resolution and sensitivity to the singularity detection. Because the displacement components came from the rough surfaces of the composite materials without any coats of gratings or fringes of optical interferometry, both surface profiles and the deformation fields of the composites were visualized which was helpful to compare each other to analyze the damage of those heterogeneous materials.