共查询到20条相似文献,搜索用时 46 毫秒
1.
采用离散元方法,利用半径扩展法和重力沉积法分别生成具有初始各向同性和各向异性内结构的试样,并开展三轴不排水压缩和拉伸试验,研究不同制样方法产生的初始各向异性对砂土宏微观力学特性及其临界状态的影响。运用组构张量对砂土的各向异性进行量化,分析不同初始组构各向异性对组构张量演化的影响并确定了组构张量的临界值。试验结果表明:初始组构各向异性对试样的剪胀性有重要影响,由于受重力影响形成初始各向异性,其各向异性程度越大、组构方向与加载方向越一致,剪胀性越显著;初始组构各向异性对试样的临界状态没有影响,砂土的组构张量具有唯一的临界状态值。 相似文献
2.
The strain space multiple mechanism model idealizes the behavior of granular materials based on a multitude of virtual simple shear mechanisms oriented in arbitrary directions. Within this modeling framework, the virtual simple shear stress is defined as a quantity that depends on the contact distribution function as well as the normal and tangential components of inter‐particle contact forces, which evolve independently during the loading process. In other terms, the virtual simple shear stress is an intermediate quantity in the upscaling process from the microscopic level (characterized by the contact distribution and inter‐particle contact forces). The stress space fabric (i.e. the orientation distribution of the virtual simple shear stress) produces macroscopic stress through the tensorial average. Thus, the stress space fabric characterizes the fundamental and higher modes of anisotropy induced in granular materials. Comparing an induced fabric associated with the biaxial shear of plane granular assemblies obtained via a simulation using Discrete Element Method to the strain space multiple mechanism model suggests that the strain space multiple mechanism model has the capability to capture the essential features in the evolution of an induced fabric in granular materials. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
3.
This paper presents a micro-mechanical study on the characteristics of shear-induced anisotropy in granular media. Based on three-dimensional Discrete Element Method (DEM) simulations, the distinct features associated with the evolution of internal granular structure and different anisotropy sources during drained/undrained shearing of granular samples are carefully examined. The study finds that static liquefaction occurs when the geometrical anisotropy in a sample dominates the mechanical anisotropy in the overall shear strength, and the weak force network features an exceptionally high proportion of sliding contacts and develops certain degree of anisotropy. Phase transformation corresponds to a transitional, unstable state associated with a dramatic change in both coordination number and the proportion of sliding contacts in all contacts. The critical state in a granular material is always associated with a highly anisotropic fabric structure wherein both the critical void ratio and critical fabric anisotropy are uniquely related to the mean effective stress. The relations provide a more comprehensive definition for the critical state in granular media with proper reference to the critical fabric anisotropy. 相似文献
4.
This paper presents a micromechanics-based approach to investigate the effects of fabric anisotropy on the behavior of localized failure in granular materials. Based on a micromechanical analysis, the origin of deviatoric stress is decomposed into two components: contact force anisotropy and fabric anisotropy. Using a micro–macro approach, the back stress is interpreted as an contribution to the change of the fabric’s principal direction. The evolution of the back stress is deduced from the stress–fabric relationship and determined with reference to the deviation of the principal directions between the rate of the reduced stress tensor and the actual reduced stress tensor. With this micro–macro framework, a mixed (isotropic–kinematic) hardening model is developed based on the classical isotropic hardening theory. A laboratory simple shear test is first analyzed to validate the proposed model and illustrate the kinematic-hardening mechanism which is usually displayed under non-proportional loading. The analysis further focuses on the anisotropic aspect of localized failure. It has been discovered that the fabric anisotropy can play an important role in the occurrence of shear banding. An increasing degree of fabric anisotropy tends to delay the initiation of the strain localization and result in higher failure strength. The effects of fabric anisotropy have also been illustrated by comparing the theoretical predictions and measured results on the shear band inclination angle, shear strain level and dilatancy at bifurcation. 相似文献
5.
The mechanical response of an assembly of particles depends on the applied boundary conditions. Robust calibration of numerical discrete systems to laboratory results is also a primary step in many studies of granular materials. In this study, a new membrane model was developed for simulating axisymmetric element tests. This membrane model uses a simple algorithm of an array of independently controlled walls and is computationally efficient. The effect of boundary flexibility on the system response was investigated by simulating a series of triaxial tests on dense and loose specimens. At the specimen scale, differences in shear strength and volume change of specimens were observed. It was shown that localization pattern depends on the applied boundary conditions. At the particle scale, particle‐membrane contact forces, coordination number, local void ratio, and anisotropy of fabric were all affected by the boundary flexibility. 相似文献
6.
Natural clays usually show anisotropic stiffness due to their deposition process and anisotropic in situ stress state. The stiffness anisotropy depends on both of the stress anisotropy and fabric anisotropy, while the latter can be quantified by the stiffness anisotropy at isotropic stress states. This paper measures the K0 value (i.e., stress anisotropy) and elastic shear stiffness anisotropy of natural Shanghai clay in a triaxial apparatus with horizontal and vertical bender elements. The results show that the K0 value of Shanghai clay lies in the range of 0.40–0.66, and an empirical equation is proposed to estimate the K0 value based on the plasticity index and initial void ratio. The fabric anisotropy of natural Shanghai clay lies in the range of 1.2–1.4 with a stronger fabric in the horizontal plane. Moreover, the experimental data of the stiffness anisotropy and fabric anisotropy of different clays in the literature are reviewed and analyzed. It reveals that the stiffness anisotropy generally increases, while the fabric anisotropy remains nearly the same during K0 consolidation. For normally consolidated clay, the fabric anisotropy generally lies in the range of 1.1–1.7. For overconsolidated clays, the fabric anisotropy generally increases as the overconsolidation ratio increases. Empirical equations are proposed to approximately estimate the fabric anisotropy of clays based on its stress normalized elastic shear stiffness. 相似文献
7.
This paper presents a fabric tensor-based bounding surface model accounting for anisotropic behaviour (e.g. the dependency of peak strength on loading direction and non-coaxial deformation) of granular materials. This model is developed based on a well-calibrated isotropic bounding surface model. The yield surface is modified by incorporating the back stress which is proportional to a contact normal-based fabric tensor for characterising fabric anisotropy. The evolution law of the fabric tensor, which is dependent on both rates of the stress ratio and the plastic strain, rules that the material fabric tends to align with the loading direction and evolves towards a unique critical state fabric tensor under monotonic shearing. The incorporation of the evolution law leads to a rotational hardening of the yield surface. The anisotropic critical state is assumed to be independent of the initial values of void ratio and fabric tensor. The critical state fabric tensor has the same intermediate stress ratio (i.e. b value) and principal directions as the critical state stress tensor. A non-associated flow rule in the deviatoric plane is adopted, which is able to predict the non-coaxial flow naturally. The stress–strain relation and fabric evolution of model predictions show a satisfactory agreement with DEM simulation results under monotonic shearing with different loading directions. The model is also validated by comparing with laboratory test results of Leighton Buzzard sand and Toyoura sand under various loading paths. The comparison results demonstrate encouraging applicability of the model for predicting the anisotropic behaviour of granular materials. 相似文献
8.
In an effort to study undrained post-liquefaction shear deformation of sand, the discrete element method (DEM) is adopted to conduct undrained cyclic biaxial compression simulations on granular assemblies consisting of 2D circular particles. The simulations are able to successfully reproduce the generation and eventual saturation of shear strain through the series of liquefaction states that the material experiences during cyclic loading after the initial liquefaction. DEM simulations with different deviatoric stress amplitudes and initial mean effective stresses on samples with different void ratios and loading histories are carried out to investigate the relationship between various mechanics- or fabric-related variables and post-liquefaction shear strain development. It is found that well-known metrics such as deviatoric stress amplitude, initial mean effective stress, void ratio, contact normal fabric anisotropy intensity, and coordination number, are not adequately correlated to the observed shear strain development and, therefore, could not possibly be used for its prediction. A new fabric entity, namely the Mean Neighboring Particle Distance (MNPD), is introduced to reflect the space arrangement of particles. It is found that the MNPD has an extremely strong and definitive relationship with the post-liquefaction shear strain development, showing MNPD’s potential role as a parameter governing post-liquefaction behavior of sand. 相似文献
9.
Fabric and its evolution need to be fully considered for effective modeling of the anisotropic behavior of cohesionless granular sand. In this study, a three‐dimensional anisotropic model for granular material is proposed based on the anisotropic critical state theory recently proposed by Li & Dafalias [2012], in which the role of fabric evolution is highlighted. An explicit expression for the yield function is proposed in terms of the invariants and joint invariants of the normalized deviatoric stress ratio tensor and the deviatoric fabric tensor. A void‐based fabric tensor that characterizes the average void size and its orientation of a granular assembly is employed in the model. Upon plastic loading, the material fabric is assumed to evolve continuously with its principal direction tending steadily towards the loading direction. A fabric evolution law is proposed to describe this behavior. With these considerations, a non‐coaxial flow rule is naturally obtained. The model is shown to be capable of characterizing the complex anisotropic behavior of granular materials under monotonic loading conditions and meanwhile retains a relatively simple formulation for numerical implementation. The model predictions of typical behavior of both Toyoura sand and Fraser River sand compare well with experimental data. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
10.
存在临界状态是颗粒材料的一个重要特性。基于孔隙胞元的颗粒离散元方法对二维颗粒体进行双轴加载数值试验,在详细分析数值模拟结果的基础上,从微观几何组构的角度揭示了临界状态的存在机制。基于剪胀性原理,提出了以接触价键表征的微观临界状态理论模型,得到了接触价键与塑性剪切应变的关系表达式,理论模型的结果和二维离散元数值模拟得到的结果吻合较好。通过比较不同情况下数值结果和理论模型中的参数,得到以下结论:表征微观临界状态的参数(临界接触价键和达到临界状态所需要的塑性剪切应变)依赖于颗粒体的微观特性,如颗粒形状、表面摩擦性质、颗粒体的围压和初始孔隙比。 相似文献
11.
通过一系列真三轴离散元数值试验,模拟了不同应力路径下的等b试验中散体材料的强度特征。根据模拟结果详细地分析了三维应力条件下中主应力和应力路径对散体材料峰值强度的影响,研究了峰值摩擦角、峰值应力比的变化规律,并根据真应力的概念和组构张量的演化结果分析了散体材料的强度成因。研究表明,在不同类型的数值试验中峰值偏应力随b参数的变化规律不同,但采用初始围压归一化后的应力-应变曲线规律一致。峰值强度线的斜率只与b值有关而与应力路径无关,且随着b值的增加,峰值应力比qf /pf逐渐减小,数值模拟结果与室内试验结果吻合较好;随着应变的发展,数值试样的组构也随之发生变化,产生了明显的应力诱发各向异性;散体的强度为颗粒摩擦及材料各向异性共同作用的结果;理论上,组构比-应力比坐标系中破坏点位置仅取决于颗粒摩擦角 ,而数值模拟结果与理论值的差异源于颗粒间咬合和滚动摩擦的影响,其影响与颗粒表面摩擦系数有关,也受空间应力状态的影响。 相似文献
12.
In an effort to study the relation of fabrics to the critical states of granular aggregates, the discrete element method (DEM) is used to investigate the evolution of fabrics of virtual granular materials consisting of 2D elongated particles. Specimens with a great variety of initial fabrics in terms of void ratios, preferred particle orientations, and intensities of fabric anisotropy were fabricated and tested with direct shear and biaxial compression tests. During loading of a typical specimen, deformation naturally localizes within shear bands while the remaining of the sample stops deforming. Thus, studying the evolution of fabric requires performing continuous local fabric measurements inside these bands, a suitable task for the proposed DEM methodology. It is found that a common ultimate/critical state is eventually reached by all specimens regardless of their initial states. The ultimate/critical state is characterized by a critical void ratio e which depends on the mean stress p, while the other critical state fabric variables related to particle orientations are largely independent of p. These findings confirm the uniqueness of the critical state line in the e ? p space, and show that the critical state itself is necessarily anisotropic. Additional findings include the following: (1) shear bands are highly heterogeneous and critical states exist only in a statistical sense; (2) critical states can only be reached at very large local shear deformations, which are not always obtained by biaxial compression tests (both physical and numerical); (3) the fabric evolution processes are very complex and highly dependent on the initial fabrics. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
13.
Various factors, such as the volumetric fraction of constituents, mineralogy, and pore fluids, affect heat flow in granular materials. Although the stress applied on granular materials controls the formation of major pathways for heat flow, few studies have focused on a detailed investigation of its significance with regard to the thermal conductivity and anisotropy of the materials. This paper presents a numerical investigation of the stress-induced evolution of anisotropic thermal conductivity of dry granular materials with supplementary experimental results. Granular materials under a variety of stress conditions in element testing are analyzed by the three-dimensional discrete element method, and quantitative variations in their anisotropic effective thermal conductivity are calculated via the network model and conductivity tensor measurements. Results show that the directional development of contact area and fabric under anisotropic stress conditions leads to the evolution of anisotropy in thermal conductivity. The anisotropy induced in thermal conductivity by shear stress is higher than that induced by compressive stress because shear stress causes more significant changes in microstructural configurations and boundary conditions. The shear-stress-induced evolution of anisotropy between principal thermal conductivities depends on dilatancy as well as shearing mode, and the shear-driven discontinuity localizes the conductivity. Factors involved in the stress-induced evolution and their implications on the thermal conductivity characterization are discussed. 相似文献
14.
为了研究粒状材料的各向异性力学行为与细观组构演化之间的关系,采用自主研发的双轴压缩试验系统,以圆形和椭圆形截面的金属棒状材料组成的二维堆积体为试验对象,对不同大主应力方向角?(沉积面与大主应力作用面的夹角)的试样进行了各向等压、常侧向压力、等p剪切3种应力路径试验,并通过分析试样在不同变形阶段的数字照片得到了其细观组构演化规律。发现对于椭圆形截面的试样存在一个卓越剪切方向,随剪应变增大,颗粒长轴呈现出向该方向偏转的趋势,并且在大变形条件下沿该方向形成剪切带;卓越剪切方向与沉积面方向关系不大,而与大主应力作用面方向夹角约为45°+ /2, 为残余内摩擦角;随卓越剪切方向与沉降方向夹角的不同,颗粒偏转程度的不同是导致剪胀特性和峰值强度各向异性的主要原因 相似文献
16.
为了研究颗粒棱角对颗粒材料力学行为的影响,建立了具有不同棱角度的对称多面体颗粒,采用了一种简单并适合任意颗粒形状的接触本构模型,对三维离散元开源程序YADE进行了修改,研究了颗粒棱角度在模拟直剪试验中的影响以及接触力各向异性在剪切过程中的演化规律。研究结果表明,颗粒棱角度越小,颗粒间相互咬合自锁的作用越小,颗粒受剪更易转动,致使颗粒体系的剪切强度和剪胀性下降;竖向加载力越大,颗粒棱角度的影响越明显;法向接触力的各向异性在剪切过程中表现为先增后减最后趋向稳定的趋势;法向接触力的各向异性变化程度随颗粒棱角度的增大而增大。 相似文献
18.
The strength of granular materials during triaxial compression is investigated via a grain scale analysis in this paper. A 3D Discrete Element Method (DEM) program provides the triaxial strength data and helps to validate the micromechanical analysis. Some standard methods in statistics are employed first to quantitatively examine the assumptions made when deriving the stress-force-fabric (SFF) equation. After careful validation, a more concise format for the SFF equation is proposed for triaxial compressions. With this SFF equation, the strength is found to be jointly contributed by the magnitudes of the contact force anisotropy and fabric anisotropy. The influence of the initial void ratio, confining pressure and loading direction on the development of contact force anisotropy and fabric anisotropy is examined and presented. With similar techniques, the “force” term in the SFF equation is further decoupled, and an equation is obtained such that it explicitly links the contact force term with the friction coefficient between grains, a tensor defined as a statistic of the normal contact forces and a tensor defined as a statistic of the mobilisation status of contacts. Based on this equation, another equation regarding the stress ratio of granular assembly is obtained, and it clearly indicates two sources that contribute to the phenomenological friction nature of granular assembly. These two sources are caused by the contact force at the grain scale. The first is the anisotropy of the average normal contact forces, and the second is the mobilisation of contacts. 相似文献
19.
The shear behavior at the interface between the soil and a structure is investigated at the macroscale and particle‐scale levels using a 3‐dimensional discrete element method (DEM). The macroscopic mechanical properties and microscopic quantities affected by the normalized interface roughness and the loading parameters are analyzed. The macro‐response shows that the shear strength of the interface increases as the normalized roughness of the interface increases, and stress softening and dilatancy of the soil material are observed in the tests that feature rough interfaces. The particle‐scale analysis illustrates that a localized band characterized by intense shear deformation emerges from the contact plane and gradually expands as shearing progresses before stabilizing at the residual stress state. The thickness of the localized band is affected by the normalized roughness of the interface and the normal stress, which ranges between 4 and 5 times that of the median grain diameter. A thicker localized band is formed when the soil has a rough shearing interface. After the localized band appears, the granular material structuralizes into 2 regions: the interface zone and the upper zone. The mechanical behavior in the interface zone is representative of the interface according to the local average stress analysis. Certain microscopic quantities in the interface zone are analyzed, including the coordination number and the material fabric. Shear at the interface creates an anisotropic material fabric and leads to the rotation of the major principal stress. 相似文献
20.
It is well known that soil is inherently anisotropic and its mechanical behavior is significantly influenced by its fabric anisotropy. Hypoplasticity is increasingly being accepted in the constitutive modeling for soils, in which many salient features, such as nonlinear stress-strain relations, dilatancy, and critical state failure, can be described by a single tensorial equation. However, within the framework of hypoplasticity, modeling fabric anisotropy remains challenging, as the fabric and its evolution are often vaguely assumed without a sound basis. This paper presents a hypoplastic constitutive model for granular soils based on the newly developed anisotropic critical state theory, in which the conditions of fabric anisotropy are concurrently satisfied along with the traditional conditions at the critical state. A deviatoric fabric tensor is introduced into the Gudehus-Bauer hypoplastic model, and a scalar-valued anisotropic state variable signifying the interplay between the fabric and the stress state is used to characterize its impact on the dilatancy and strength of the soils. In addition, fabric evolution during shearing can explicitly be addressed. Modifications have also been undertaken to improve the performance of the undrained response of the model. The anisotropic hypoplastic model can simulate experimental tests for sand under various combinations of principle stress direction, intermediate principal stress (or mode of shearing), soil densities, and confining pressures, and the associated drastic effect of different principal stress orientations in reference to the material axes of anisotropy can be well captured. 相似文献
|