共查询到20条相似文献,搜索用时 11 毫秒
1.
2.
In this study, a homogenization theory based on the Gurtin strain gradient formulation and its finite element discretization are developed for investigating the size effects on macroscopic responses of periodic materials. To derive the homogenization equations consisting of the relation of macroscopic stress, the weak form of stress balance, and the weak form of microforce balance, the Y-periodicity is used as additional, as well as standard, boundary conditions at the boundary of a unit cell. Then, by applying a tangent modulus method, a set of finite element equations is obtained from the homogenization equations. The computational stability and efficiency of this finite element discretization are verified by analyzing a model composite. Furthermore, a model polycrystal is analyzed for investigating the grain size dependence of polycrystal plasticity. In this analysis, the micro-clamped, micro-free, and defect-free conditions are considered as the additional boundary conditions at grain boundaries, and their effects are discussed. 相似文献
3.
K.C. HwangH. Jiang Y. Huang H. GaoN. Hu 《Journal of the mechanics and physics of solids》2002,50(1):81-99
Plastic deformation exhibits strong size dependence at the micron scale, as observed in micro-torsion, bending, and indentation experiments. Classical plasticity theories, which possess no internal material lengths, cannot explain this size dependence. Based on dislocation mechanics, strain gradient plasticity theories have been developed for micron-scale applications. These theories, however, have been limited to infinitesimal deformation, even though the micro-scale experiments involve rather large strains and rotations. In this paper, we propose a finite deformation theory of strain gradient plasticity. The kinematics relations (including strain gradients), equilibrium equations, and constitutive laws are expressed in the reference configuration. The finite deformation strain gradient theory is used to model micro-indentation with results agreeing very well with the experimental data. We show that the finite deformation effect is not very significant for modeling micro-indentation experiments. 相似文献
4.
《Comptes Rendus de l'Académie des Sciences》2001,329(11):797-802
Extending the previous work by Chambon et al. [2] to the finite deformation regime, a local second gradient theory of plasticity for isotropic materials with microstructure is developed based on the multiplicative decomposition of the deformation gradient, the additive decomposition of the second deformation gradient and the principle of maximum dissipation. 相似文献
5.
借助于Cosserat连续介质模型,探讨了应力函数和位移对避免有限元C$^{1}$
连续性困难的互补性作用. 通过对应力函数对偶理论的深入分析,为将应力函数列式得到的
余能单元转化为具有一般位移自由度的势能单元提供了严格的理论基础,在此基础上,
给出应用应力函数构造有限元的一般方法. 相似文献
6.
Mitsutoshi Kuroda Viggo Tvergaard 《Journal of the mechanics and physics of solids》2008,56(8):2573-2584
For higher-order gradient crystal plasticity, a finite deformation formulation is presented. The theory does not deviate much from the conventional crystal plasticity theory. Only a back stress effect and additional differential equations for evolution of the geometrically necessary dislocation (GND) densities supplement the conventional theory within a non-work-conjugate framework in which there is no need to introduce higher-order microscopic stresses that would be work-conjugate to slip rate gradients. We discuss its connection to a work-conjugate type of finite deformation gradient crystal plasticity that is based on an assumption of the existence of higher-order stresses. Furthermore, a boundary-value problem for simple shear of a constrained thin strip is studied numerically, and some characteristic features of finite deformation are demonstrated through a comparison to a solution for the small deformation theory. As in a previous formulation for small deformation, the present formulation applies to the context of multiple and three-dimensional slip deformations. 相似文献
7.
Strain gradient plasticity for finite deformations is addressed within the framework of nonlocal continuum thermodynamics, featured by the concepts of (nonlocality) energy residual and globally simple material. The plastic strain gradient is assumed to be physically meaningful in the domain of particle isoclinic configurations (with the director vector triad constant both in space and time), whereas the objective notion of corotational gradient makes it possible to compute the plastic strain gradient in any domain of particle intermediate configurations. A phenomenological elastic–plastic constitutive model is presented, with mixed kinematic/isotropic hardening laws in the form of PDEs and related higher order boundary conditions (including those associated with the moving elastic/plastic boundary). Two fourth-order projection tensor operators, functions of the elastic and plastic strain states, are shown to relate the skew-symmetric parts of the Mandel stress and back stress to the related symmetric parts. Consistent with the thermodynamic restrictions therein derived, the flow laws for rate-independent associative plasticity are formulated in a six-dimensional tensor space in terms of symmetric parts of Mandel stresses and related work-conjugate generalized plastic strain rates. A simple shear problem application is presented for illustrative purposes. 相似文献
8.
循环硬化材料本构模型的隐式应力积分和有限元实现 总被引:1,自引:0,他引:1
针对新发展的、能够描述循环硬化行为应变幅值依赖性的粘塑性本构模型,讨论了它的数值实现方法。首先,为了能够对材料的循环棘轮行为(Ratcheting)和循环应力松弛现象进行描述,对已有的本构模型进行了改进;然后,在改进模型的基础上,建立了一个新的、全隐式应力积分算法,进而推导了相应的一致切线刚度(Consistent Tangent Modulus)矩阵的表达式;最后,通过ABAQUS用户材料子程序UMAT将上述本构模型进行了有限元实现,并通过一些算例对一些构件的循环变形行为进行了有限元数值模拟,讨论了该类本构模型有限元实现的必要性和合理性。 相似文献
9.
黄晓凡 《应用数学和力学(英文版)》1986,7(8):785-794
We consider the first boundary value problem of the second order elliptic equation and serendipity rectangular elements. Papers [2,3,9] proved that the gradients of finite element solution possess superconvergence at Gaussianpoint. In this paper, we extend the results in papers [2,3,9] in the sense that the coefficients of the elliptic equations are discontinuous on a curve S which lies in the bounded domain . 相似文献
10.
挠曲电效应是一种存在于所有电介质材料中的特殊的力电耦合效应,本质上是应变梯度与电极化之间的线性耦合。然而,应变梯度会引入位移的高阶偏量,常给挠曲电问题的理论求解带来困难。且已有研究表明应变梯度弹性项会影响纳米结构中的力电耦合响应,但是现有的挠曲电研究大多忽略了应变梯度弹性的影响。因此,本文提出了一种既考虑应变梯度弹性,又考虑挠曲电效应的有效数值方法。基于全应变梯度弹性理论,建立了包含3个独立材料尺度参数的纳米欧拉梁的理论模型和有限元模型,提出了满足C2弱连续的两节点六自由度单元。基于本文的有限单元法,以简支欧拉梁为例,通过分析讨论挠度、电势和能量效率,得到了挠曲电效应和应变梯度弹性项对梁的力电响应的影响。结果表明,挠曲电效应存在尺寸依赖性,且应变梯度弹性项在纳米电介质结构的挠曲电研究中的影响不可忽略。 相似文献
11.
12.
13.
The formulation and implementation of a singular finite element for analyzing homogeneous anistropic materials is presented in this paper. Lekhnitskii's stress function method is used to formulate the boundary value problem with the stress function expressed as a Laurent series. The development of the element stiffness matrix and the method of integrating the element to conventional displacement based finite element programs is shown. The stiffness matrix generation is based on a least squates collocation technique to satisfy displacement continuity boundary conditions at the element interface. Implementation of the element is demonstrated for cracked anisotropic materials subjected to inplane loading. Center cracked, on and off-axis coupons under tensile loading are analyzed using the element. It is shown that the stress distributions and intensity factors compare well with those obtained using other methods. 相似文献
14.
Mitsutoshi Kuroda 《International Journal of Solids and Structures》2011,48(24):3382-3394
A finite-strain higher-order gradient crystal plasticity model accounting for the backstress effect originating from the existence of geometrically necessary dislocations (GNDs) is applied to plane strain finite element analysis. Different element types are tested to seek out an element formulation that is reliable and useful for solving problems involving severe plastic deformation. In the present finite element formulation, the GND density rates are chosen to be additional nodal degrees of freedom. Different orders of shape functions are employed for the interpolation of displacement rates and GND density rates. Their effects on solutions are examined in detail by considering three boundary value problems: a simple shear of a constrained layer (a film), a compression problem with loading surfaces impenetrable to dislocations, and a tension problem involving shear band formation. In all the cases, the formulation in which eight-node elements with reduced integration and four-node elements with full integration are used respectively for displacement rates and the GND density rates gives reasonable solutions. In addition to the discussion on the choice of finite elements, detailed behavior in gradient-dependent solids, such as the accumulation of GND density and the distribution of backstress on each slip system, is investigated by utilizing the reliable computational results obtained. 相似文献
15.
A new high-precision finite element for analysis of shell structures is presented. It is derived from a slightly generalized equilibrium principle. Accordingly both stresses and displacements are obtained as primary result of analysis. At the assembly level the element has 45 degrees of freedom, all of them generalized displacements. For the price of some additional computational effort on the elemental level of analysis the proposed element is believed to gain certain advantages over the recently developed high-precision displacement elements. Thin as well as thick shell structures of arbitrary shape and loading can be equally analyzed. Engineering accuracy is attained with only very few elements. A variety of numerical examples demonstrates the applicability of the new element to all kinds of situations occuring in practice. A review of the existing high-precision shell elements is also included. 相似文献
16.
影响函数与有限元应力计算 总被引:2,自引:1,他引:2
用有限元法得到位移场后,总要计算应力场。通常的做法是对位移进行微商计算应变,再根据应力-应变关系计算应力。有限元位移计算的精度比较高,但通过用位移微商来计算应力,精度会大大降低。本文利用Hamilton对偶体系的已有成果,解析求解位移和应力的影响函数,利用有限元法计算得到的位移和节点力,通过功的互等定理,可以求得一点的应力值。因影响函数是分析解,而且计算应力时不必进行微商,应力精度大幅提高。数值结果表明该方法是可行的和有效的。由该方法编制成的计算程序,可作为有限元通用程序应力计算的一个模块,将较大地提高有限元应力计算的精度和稳定性。 相似文献
17.
A new finite element method for strain gradient theories and applications to fracture analyses 总被引:2,自引:0,他引:2
A new compatible finite element method for strain gradient theories is presented. In the new finite element method, pure displacement derivatives are taken as the fundamental variables. The new numerical method is successfully used to analyze the simple strain gradient problems – the fundamental fracture problems. Through comparing the numerical solutions with the existed exact solutions, the effectiveness of the new finite element method is tested and confirmed. Additionally, an application of the Zienkiewicz–Taylor C1 finite element method to the strain gradient problem is discussed. By using the new finite element method, plane-strain mode I and mode II crack tip fields are calculated based on a constitutive law which is a simple generalization of the conventional J2 deformation plasticity theory to include strain gradient effects. Three new constitutive parameters enter to characterize the scale over which strain gradient effects become important. During the analysis the general compressible version of Fleck–Hutchinson strain gradient plasticity is adopted. Crack tip solutions, the traction distributions along the plane ahead of the crack tip are calculated. The solutions display the considerable elevation of traction within the zone near the crack tip. 相似文献
18.
A numerical procedure to calculate the pre-buckling and postbuckling response of general structures is presented. This procedure
is based on the pseudo-arclength algorithm suggested by E. Riks et al., which has some numerical difficulties during implementation
of large applied analysis programs. To overcome these difficulties, a scheme based on rank-1 modification of the matrix is
proposed. Some example show this procedure behaves well in passing through the limit point and is rather efficient.
recommended by Prof. Chen Yaosong and Prof. Wu Jike 相似文献
19.
《Acta Mechanica Solida Sinica》2017,(4)
In this paper, a 3-node triangular element for couple stress theory is proposed based on the assumed stress quasi-conforming method. The formulation starts from polynomial approximation of stresses. Then the stress-function matrix is treated as the weighted function to weaken the strain-displacement equations. Finally, the string-net functions are introduced to calculate strain integration and the stress smooth technique is adopted to improve the stress accuracy. Numerical results show that the proposed new model can pass the C~(0-1) patch test with excellent precision, does not exhibit extra zero energy modes and can capture the scale effects of microstructure. 相似文献
20.
《International Journal of Solids and Structures》2002,39(10):2731-2743
The deformation behavior of materials in the micron scale has been experimentally shown to be size dependent. In the absence of stretch and dilatation gradients, the size dependence can be explained using classical couple stress theory in which the full curvature tensor is used as deformation measures in addition to the conventional strain measures. In the couple stress theory formulation, only conventional equilibrium relations of forces and moments of forces are used. The couple's association with position is arbitrary. In this paper, an additional equilibrium relation is developed to govern the behavior of the couples. The relation constrained the couple stress tensor to be symmetric, and the symmetric curvature tensor became the only properly conjugated high order strain measures in the theory to have a real contribution to the total strain energy of the system. On the basis of this modification, a linear elastic model for isotropic materials is developed. The torsion of a cylindrical bar and the pure bending of a flat plate of infinite width are analyzed to illustrate the effect of the modification. 相似文献