Measurement of length-scale and solution of cantilever beam in couple stress elasto-plasticity

Owing to the absence of proper analytical solution of cantilever beams for couple stress/strain gradient elasto-plastic theory, experimental studies of the cantilever beam in the micro-scale are not suitable for the determination of material length-scale. Based on the couple stress elasto-plasticity, an analytical solution of thin cantilever beams is firstly presented, and the solution can be regarded as an extension of the elastic and rigid-plastic solutions of pure bending beam. A comparison with numerical results shows that the current analytical solution is reliable for the case of σ0 〈〈 H 〈〈 E, where σ0 is the initial yield strength, H is the hardening modulus and E is the elastic modulus. Fortunately, the above mentioned condition can be satisfied for many metal materials, and thus the solution can be used to determine the material length-scale of micro-structures in conjunction with the experiment of cantilever beams in the micro-scale.     

STUDY ON TENSILE FAILURE MECHANISM AND HYBRID EFFECT OF INTRAPLY HYBRID COMPOSITES

The random critical-core model is adapted to investigate the tensile failure mech-anism and hybrid effect of unidirectionally arrayed hybrid composites with alternating low andhigh elongation fibers.By utilizing the model in conjunction with the results of the stress con-centration analysis in which the interfacial damage between fiber and matrix is considered,amicroscopic statistical analysis of both the first failure and ultimate failure of hybrids is per-formed.The variations of the first failure strain,the ultimate failure strain and the hybrid effectas the interfacial shear strength are obtained quantitatively.The concept of the hybrid effect forstrains has been clarified.The present results are compared with available experiment data anda reasonable agreement is found between the analytical predictions and the experimental results.     

THE LAYOUT OPTIMIZATION OF STIFFENERS FOR PLATE-SHELL STRUCTURES

The plate-shell structures with stiffeners are widely used in a broad range of engineering structures. This study presents the layout optimization of stiffeners. The minimum weight of stiffeners is taken as the objective function with the global stiffness constraint. In the layout optimization, the stiffeners should be placed at the locations with high strain energy/or stress. Conversely, elements of stiffeners with a small strain energy/or stress are considered to be used inefficiently and can be removed. Thus, to identify the element efficiency so that most inefficiently used elements of stiffeners can be removed, the element sensitivity of the strain energy of stiffeners is introduced, and a search criterion for locations of stiffeners is presented. The layout optimization approach is given for determining which elements of the stiffeners need to be kept or removed. In each iterative design, a high efficiency reanalysis approach is used to reduce the computational effort. The present approach is implemented for the layout optimization of stiffeners for a bunker loaded by the hydrostatic pressure. The numerical results show that the present approach is effective for dealing with layout optimization of stiffeners for plate-shell structures.     

THE APPLICATION OF DISCRETE ELEMENT METHOD IN SOLVING THREE-DIMENTIONAL IMPACT DYNAMICS PROBLEMS

A three-dimensional discrete element model of the connective type is presented. Moreover, a three- dimensional numerical analysis code, which can carry out the transitional process from connective model (for continuum) to contact model (for non-continuum), is developed for simulating the mechanical process from continuum to non-continuum. The wave propagation process in a concrete block (as continuum) made of cement grout under impact loading is numerically simulated with this code. By comparing its numerical results with those by LS-DYNA, the calculation accuracy of the model and algorithm is proved. Furthermore, the failure process of the concrete block under quasi-static loading is demonstrated, showing the basic dynamic transitional process from continuum to non-continuum. The results of calculation can be displayed by animation. The damage modes are similar to the experimental results. The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum. It also shows that the discrete element method (DEM) will have broad prospects for development and application.     

Theoretical study of void closure in nonlinear plastic materials

Void closing from a spherical shape to a crack is investigated quantitatively in the present study. The constitutive relation of the Void-free matrix is assumed to obey the Norton power law. A representative volume element （RVE） which includes matrix and void is employed and a Rayleigh-Ritz procedure is developed to study the deformation-rates of a spherical void and a penny-shaped crack. Based on an approximate interpolation scheme, an analytical model for void closure in nonlinear plastic materials is established. It is found that the local plastic flows of the matrix material are the main mechanism of void deformation. It is also shown that the relative void volume during the deformation depends on the Norton exponent, on the far-field stress triaxiality, as well as on the far-field effective strain. The predictions of void closure using the present model are compared with the corresponding results in the literature, showing good agreement. The model for void closure provides a novel way for process design and optimization in terms of elimination of voids in billets because the model for void closure can easily be applied in the CAE analysis.     

THREE-DIMENSIONAL MODELING FOR THIN PLATE-LIKE STRUCTURES INCLUDING SURFACE EFFECTS BY USING STATE SPACE METHOD

<正>A three-dimensional(3-D)approach based on the state space method is proposed to study size-dependent mechanical properties of ultra-thin plate-like elastic structures considering surface effects.The structure is modeled as a laminate composed of a bulk bounded with upper and bottom surface layers,which are allowed to have different material properties from the bulk layer.State equations,including the surface properties of the structure,can be established on the basis of 3-D fundamental elasticity to analyze the size-dependent static characteristics of the thin plate-like structure.Compared with two-dimensional plate theories based size-dependent models for thin film structures in literature,the present 3-D approach is exact,which can provide benchmark results to assess the accuracy of 2-D plate theories and various numerical approaches. To show the feasibility of the proposed approach,a 3-D analytical solution for a simply supported plate-like thin structure including surface layers is derived.An algorithm is proposed for the calculation of the state equations obtained to ensure that the numerical results can reveal the surface effects clearly even for extremely thin surface layers.Numerical examples are carried out to exhibit the surface effects and some discussions are provided based on the results obtained.     

Mechanism of unsteady aerodynamic heating with sudden change in surface temperature

The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is presented with both analytical and numerical approaches. With the analytical method, the unsteady compressible boundary layer equation is solved. In the neighborhood of the initial and final steady states, the transient responses can be expressed with a steady-state solution plus a perturbation series. By combining these two solutions, a complete solution in the entire time domain is achieved. In the region in which the analytical approach is applicable, numerical results are in good agreement with the analytical results, showing reliability of the methods. The result shows two distinct features of the unsteady response. In a short period just after a sudden increase in the wall temperature, the direction of the wall heat flux is reverted, and a new inflexion near the wall occurs in the profile of the thermal boundary layer. This is a typical unsteady characteristic. However, these unsteady responses only exist in a very short period in hypersonic flows, meaning that, in a long-term aerodynamic heating process considering only unsteady surface temperature, the unsteady characteristics of the flow can be ignored, and the traditional quasi-steady aerodynamic heating prediction methods are still valid.     

Analytical solution for the stress field of hierarchical defects:multiscale framework and applications

Hierarchical defects are defined as adjacent defects at different length scales.Involved are the two scales where the stress field distribution is interrelated.Based on the complex variable method and conformal mapping,a multiscale framework for solving the problems of hierarchical defects is formulated.The separated representations of mapping function,the governing equations of potentials,and the stress field are subsequently obtained.The proposed multiscale framework can be used to solve a variety of simplified engineering problems.The case in point is the analytical solution of a macroscopic elliptic hole with a microscopic circular edge defect.The results indicate that the microscopic defect aggregates the stress concentration on the macroscopic defect and likely leads to global propagation and rupture.Multiple micro-defects have interactive effects on the distribution of the stress field.The level of stress concentration may be reduced by the coalescence of micro-defects.This work provides a unified method to analytically investigate the influence of edge micro-defects within the scope of multiscale hierarchy.The formulated multiscale approach can also be potentially applied to materials with hierarchical defects,such as additive manufacturing and bio-inspired materials.     

DISCRETE AND CONTINUUM MODELLING OF GRANULAR FLOW

This paper analyses three popular methods simulating granular flow at different time and length scales： discrete element method （DEM）, averaging method and viscous, elastic-plastic continuum model. The theoretical models of these methods and their applications to hopper flows are discussed. It is shown that DEM is an effective method to study the fundamentals of granular flow at a particle or microscopic scale. By use of the continuum approach, granular flow can also be described at a continuum or macroscopic scale. Macroscopic quantities such as velocity and stress can be obtained by use of such computational method as FEM. However, this approach depends on the constitutive relationship of materials and ignores the effect of microscopic structure of granular flow. The combined approach of DEM and averaging method can overcome this problem. The approach takes into account the discrete nature of granular materials and does not require any global assumption and thus allows a better understanding of the fundamental mechanisms of granular flow. However, it is difficult to adapt this approach to process modelling because of the limited number of particles which can be handled with the present computational capacity, and the difficulty in handling non-spherical particles. Further work is needed to develoo an aoorooriate aooroach to overcome these problems.     

New discrete element models for elastoplastic problems

The discrete element method （DEM） has attractive features for problems with severe damages, but lack of theoretical basis for continua behavior especially for nonlinear behavior has seriously restricted its application. The present study proposes a new approach to developing the DEM as a general and robust technique for modeling the elastoplastic behavior of solid materials. New types of connective links between elements are proposed, the interelement parameters are theoretically determined based on the principle of energy equivalence and a yield criterion and a flow rule for DEM are given for describing nonlinear behavior of materials. Moreover, a numerical scheme, which can be applied to modeling the behavior of a continuum as well as the transformation from a continuum to a discontinuum, is obtained by introducing a fracture criterion and a contact model into the DEM. The elastoplastic stress wave propagations and the tensile failure process of a steel plate are simulated, and the numerical results agree well with those obtained from the finite element method （FEM） and corresponding experiment, and thus the accuracy and efficiency of the DEM scheme are demonstrated.     

3-D numerical simulation of fracture processes in heterogeneous brittle materials

By using the lattice model combined with finite element methods and statistical techniques, a numerical approach is developed to establish mechanical models of three-dimensional heterogeneous brittle materials. A special numerical code is introduced, in which a lattice model and statistical approaches are used to simulate the initial heterogeneity of material properties. The size of displacement-load step is adaptively determined so that only few elements would fail in each load step. When the tensile principal strain in an element exceeds the ultimate strain of this element, the element is considered broken and its Young's modulus is set to be very low. Some important behaviors of heterogeneous brittle materials are indicated using this code. Load-displacement curves and figures of three-dimensional fracture patterns are also numerically obtained, which are similar to those observed in laboratory tests.     

三维非均匀脆性材料破坏过程的数值模拟

采用有限元方法模拟了三维均匀固体材料在宏观等效力学性质和破坏过程。首先采用格形（lattice）方法把试件离散成三维均匀网格,在每个单元格中将材料按照均匀处理,根据给定的统计规则来确定不同单元格中的材料常数以反映材料的非均匀性。然后对非均匀脆性材料选用简单的本构关系与断裂准则,采用自适应选取载荷步长对试件进行加载,通过非平衡迭代技术对刚度矩阵进行不断修正,实现了非均匀脆性材料的弹性行为及破坏过程的数值模拟。在此基础上,通过数值计算研究了材料的非均匀分布对宏观等效力学性质和破坏过程的影响,给出了破坏全过程的非线性载荷－位移曲线以及不同载荷阶段的三维损伤破坏的演化图。     

The application of the Eshelby equivalent inclusion method for unifying modulus and transformation toughening

When a crack is lodged in an inclusion, both difference between the modulus of the inclusion and matrix material and stress-free transformation strain of the inclusion will cause the near-tip stress intensity factor to be greater (amplification effect) or less (shielding or toughening effect) than that prevailing in a homogeneous material. In this paper, the inclusion may represent a second phase particle in composites and a transformation or microcracked process zone in brittle materials, which may undergo a stress-free transformation strain induced by phase transformation, microcracking, thermal expansion mismatch and so forth. A close form of solution is derived for predicting the toughening (or amplification) effect. The derivation is based on Eshelby equivalent inclusion approach that provides rigorous theoretical basis to unify the modulus and transformation contributions to the near-tip field. As validated by numerical examples, the developed formula has excellent accuracy for different application cases.     

热环境中旋转运动功能梯度圆板的强非线性固有振动

研究热环境中旋转运动功能梯度圆板的非线性固有振动问题．针对金属－陶瓷功能梯度圆板,考虑几何非线性、材料物理属性参数随温度变化以及材料组分沿厚度方向按幂律分布的情况,应用哈密顿原理推得热环境中旋转运动功能梯度圆板的非线性振动微分方程．考虑周边夹支边界条件,利用伽辽金法得到了横向非线性固有振动方程,并确定了静载荷引起的静挠度．用改进的多尺度法求解强非线性方程,得出非线性固有频率表达式．通过算例,分析了旋转运动功能梯度圆板固有频率随转速、温度等参量的变化情况．结果表明,非线性固有频率随金属含量的增加而降低;随转速和圆板厚度的增大而升高;随功能梯度圆板表面温度的升高而降低．     

2D编织陶瓷基复合材料应力-应变行为的试验研究和模拟

本文对2D编织陶瓷基复合材料拉伸应力-应变行为进行了试验研究和理论模拟。将2D编织结构简化为：正交铺层结构和纤维束波动结构。基于基体随机开裂、纤维随机断裂的统计分布理论,得到正交铺层结构的应力-应变关系;基于体积平均方法,将纤维束波动部分分割为若干子单元;由于纤维束的波动使各子单元材料方向与加载方向不一致,因此考虑了各子单元的线性行为和非线性行为对材料响应的影响,同时引入强度分析模型,得到纤维束波动部分的应力-应变关系。结合正交铺层部分和纤维束波动部分的应力-应变关系,得到2D编织结构的应力-应变行为,理论与试验吻合较好。     

岩石紧凑拉伸断裂过程及其尺寸效应的三维数值模拟

采用最近开发的三维岩石破裂过程分析软件RFPA3D模拟单边裂纹紧凑拉伸断裂过程。试验中五个不同尺寸的岩样具有相同的力学性质参数分布,模拟结果得到了裂纹扩展中的应力场、位移场和声发射的空间分布以及单边裂纹扩展贯通的过程。单边裂纹拉伸断裂的路径是一个复杂的空间三维曲面,三维裂纹比二维裂纹更为复杂。分析了岩石试样的峰值强度和试样尺寸之间的关系。随着岩样尺寸的增加,峰值强度逐渐减小,并且延性破坏特征更加明显,模拟结果满足岩石的尺寸效应规律。最后模拟了三组不同均匀性的试样拉伸破坏过程,结果表明细观上的非均匀性对岩石尺寸效应有很大影响,随着非均匀性的增加,岩石宏观强度随之提高,即使在均匀材料中一样存在尺寸效应。     

岩石单轴压缩作用下变形局部化的梯度塑性解

采用梯度塑性理论研究单轴压缩作用下岩石变形局部化,得到了单轴压缩作用下岩石变形局部化带宽度的一维、二维解析解,为实验测定内部材料长度参数提供了理论依据.     

含水率和围压对安家岭泥岩峰后力学特性影响的试验研究

为分析含水率和围压对泥岩峰后力学特性的影响,从山西安家岭矿取泥岩,制成不同含水率的5组试样,在YAW2000电液伺服试验机上开展了三轴试验,获得了不同含水率泥岩试样三轴全程应力应变曲线,使用激光共聚焦显微镜观测了泥岩增水过程中微观结构变化。利用试验结果,分析了围压和含水率对泥岩峰值强度、残余强度、弹性模量、破坏应变和脆性模量的影响规律和泥岩增水过程中微观结构的变化规律。引入脆性模量系数和强度退化指数来描述围压对泥岩峰后强度退化过程和残余强度的影响,与FLAC中的SS模型结合,建立了考虑围压影响的泥岩应变软化力学模型,模拟了围压对泥岩应变软化行为的影响。研究结果表明:(1)随着围压增加,泥岩的峰值强度、偏应力峰值、破坏应变和残余强度都增长,峰后强度降低速率趋缓,强度退化指数和脆性模量系数可以较好地描述围压对泥岩残余强度和峰后强度退化过程的影响。(2)泥岩增水过程中,岩样内微裂隙及尺寸增长,泥岩的力学性质劣化。随着含水率增加,泥岩的弹性模量、峰值强度和残余强度降低,破坏应变增长。含水率与泥岩的弹性模量、峰值强度和破坏应变之间近似服从线性关系。(3)本文基于脆性模量系数的岩石应变软化模型能较好地描述三轴压缩泥岩的全程变形行为。     

基于裂纹扩展脆性岩石动力压缩力学特性研究

动态压缩荷载作用下，脆性岩石内部动态细观裂纹扩展特性，对岩石宏观动态力学特性有着重要的影响。然而，对岩石内部动态细观裂纹扩展与宏观动态力学特性的关系研究较少。基于准静态裂纹扩展作用下的应力-应变本构模型、准静态与动态裂纹扩展断裂韧度关系、裂纹速率与应变率关系模型及应变率与动态断裂韧度关系，提出了一种基于细观力学的动态应力-应变本构模型。其中裂纹速率与应变率关系，是根据裂纹长度与应变关系的时间导数推出；应变率与动态断裂韧度关系，是根据推出的裂纹速率及应变率关系，与裂纹速率及断裂韧度关系相结合而得到。研究了应变率对应力-应变本构关系及动态压缩强度影响。并通过试验结果验证了模型的合理性。讨论了岩石初始损伤、围压、模型中参数m、ε₀和R对应力-应变关系、动态压缩强度和动态弹性模量的影响。研究结果可为动态压缩荷载作用下深部地下工程脆性围岩稳定性分析提供了一定的理论支持。