金属基复合材料单向和循环变形行为的数值模拟研究

利用晶体细观力学模型，依据组元材料的单晶体变形性质，以Ａｌ－Ａｌ２Ｃｕ自生复合材料为模型材料，用数值法模拟其拉伸和循环拉伸变形过程，得到了与有关实验结果相符的应力－应变曲线，研究了增强相间距和循环加载过程对复合材料变形行为的影响，通过考虑变形过程中组元相及其界面的应力分布规律，分析研究了循环变形过程中基体Ｂａｕｓｃｈｉｎｇｅｒ效应对于复合材料微观变形与损伤机制的作用，结果表明，晶体细观力学模拟计算     

铜双晶体的循环应力-应变响应

比较了铜单晶体和多晶体疲劳行为的异同,提出了研究双晶体疲劳行为的必要性．总结了具有不同晶体取向和晶界的铜双晶体的疲劳行为的最新进展．利用平行晶界铜双晶体的取向因子和晶界影响区,总结了在循环载荷作用下的晶界强化模型．分析了垂直晶界铜双晶体循环塑性变形行为的特点,讨论了组元晶体取向对垂直晶界铜双晶体循环应力－应变曲线的影响．提出了提高单晶体和双晶体疲劳强度的控制因素．     

单晶体和双晶体微观层次变形行为的有限元分析

从微观层次上研究金属材料的变形行为,将位错引入到本构关系中,用硬化函数描述材料的硬化规律,考虑了变形的率相关性,采用三维模型用大变形有限单元法对单晶体在单向拉伸载荷和循环载荷作用下的变形行为、双晶体在单向拉伸作用下滑移系的开动进行了模拟计算,得到了与实验一致的计算结果。     

基于晶体塑性理论的孔洞生长行为研究

采用率相关晶体塑性模型,建立三维胞元计算模型,研究了晶粒取向和晶界对孔洞生长和聚合的影响.比较了不同晶粒取向的单晶和双晶体中孔洞的生长趋势,发现晶粒取向对孔洞生长方向,孔洞形状等有着显著的影响.     

FCC晶体中不均匀变形对空洞长大的影响

通过编制率相关有限元用户子程序,采用一个单胞模型研究了FCC晶体中孔洞在单晶及晶界的长大行为,分析了由于晶体取向及变形失配对孔洞长大和聚合的影响。研究结果表明：孔洞的形状和长大方向与晶体取向密切相关;晶界上孔洞的长大速度大于单晶中孔洞的长大速度;晶粒间的变形失配加速了晶界上孔洞的长大趋势,因而使材料易发生沿晶断裂,随着晶粒间取向因子差异的增加,孔洞越易沿着晶界长大。     

考虑两种晶界的各向异性双晶和三晶体晶界附近弹塑性应力场分析

采用率相关的晶体滑移有限元程序对具有不同晶体取向的双晶体晶界附近及三晶体三晶粒交汇处的弹塑性应力场进行了计算,考虑了几何晶界和物理晶界的影响.计算结果表明:双晶体及三晶体考虑几何晶界和物理晶界时,这两种晶界具有相同的应力分布趋势,只是物理晶界比几何晶界的应力集中程度小,双晶体晶界附近有较大的应力梯度,存在应力集中现象.三晶体三晶粒交汇处可能是应力集中之地也可能不造成应力集中,这主要取决于晶粒晶体取向及加载方向.由此可见,要准确理解金属材料的断裂过程,还需要从细观的角度对晶界的力学响应进行细致和深入的研究.     

含与不含晶界空穴的双晶体蠕变行为研究

基于晶体滑移理论,建立了各向异性镍基合金双晶体的蠕变本构模型和蠕变寿命预测模型,通过MARC用户子程序CRPLAW将上述本构模型进行了有限元实现,并对双晶体蠕变行为进行了计算分析,考虑了:(1)晶体取向的影响;(2)垂直、倾斜和平行于外载方向的三种位向晶界情况;(3)晶界处引进空间空穴的影响。结果表明,双晶体上特别是微空穴和晶界附近区域的蠕变应力应变呈现不同的变化规律,对此晶粒晶体取向和晶界位向有较大的影响;微空穴的存在削弱了双晶体的承载能力,显著地影响了双晶体蠕变持久寿命;相同条件下,垂直晶界对双晶体模型的蠕变损伤影响最为强烈,倾斜晶界次之,平行晶界最小;微空穴的生长与晶界位向和晶体取向有强烈的依赖关系,其中垂直晶界更有利于晶体滑移和微空穴生长。     

双晶体的取向因子

本文通过对双晶体应力分析，提出了双晶体的取向因子ΩＢ，与其组元单晶体Ｇ１，Ｇ２的Ｓｃｈｍｉｄ因子Ω１及Ω２的关系为ΩＢ＝ＶＧ１Ω１＋ＶＧ２Ω２－１．并借助于双晶体的取向因子ΩＢ来比较单晶体与双晶体塑性变形行为的差别，以确定晶界对滑移带开动的晶界阻力．     

表面微空洞长大和相互作用的晶体有限元分析

通过建立空洞长大和相互作用的3D模型,采用晶体塑性有限元模拟研究了FCC晶体表面空洞的长大和相互作用行为,分析了晶体取向和微空洞在表面的深度变化对表面空洞长大和相互作用的影响。模拟结果表明：晶体取向除了影响空洞形状和长大方向外,还会影响空洞长大速度;总体而言,在固定位移边界条件下硬取向晶粒表面的空洞长大和相互作用大于软取向。随着空洞在单晶体表面深度的增加,空洞周围的最大塑性变形增加,变形局部化更加严重,空洞长大速度增加。     

考虑尺寸效应的多晶金属循环塑性分析

本文采用多晶塑性分析方法,设材料点包含一定数量的各向异性单晶晶粒并考虑晶粒尺寸的影响,计算材料点的应力和应变时利用了Taylor假设。模型引入考虑尺寸效应的晶体滑移硬化函数,同时针对晶体滑移引入背应力及其方向性硬化的描述,以反映不同晶粒尺寸材料在循环加载条件下的力学行为。利用该模型,本文第一作者采用显式格式编制了与ABAQUS商用有限元软件接口的用户材料子程序（VUMAT）,实例计算证实该模型可以反映和描述多晶金属材料在材料反复加载条件下的循环塑性行为与尺寸效应。     

Simulation of cyclic deformation of bicrystals: Effects of orientation and deformation properties of constituent grains

The cyclic deformation of aluminium bicrystals was numerically simulated by using a crystalline micromechanical approach to study the effects of crystallographic orientation and deformation properties of constituent grains. Both the macro-and micro-stress-strain responses were analyzed. It was shown that the whole Bauschinger effect and cyclic hardening behavior of bicrystals is primarily controlled by the corresponding properties of the constituent grains, and that the effects of internal stresses induced by the grain boundary and the interaction between grains are secondary. The crystallographic orientation of constituent grains was also shown to have significant influence on the macro-and micro-mechanical behavior of the bicrystals. There were smaller symmetry in the orientation or greater difference in the strength between grains and more remarkable GB effects.     

磁电弹复合材料和刚性导电导磁圆柱压头的二维微动接触分析

本文求解平面应变状态下磁电弹复合材料半平面和刚性导电导磁圆柱压头的二维微动接触问题。假设压头具有良好的导电导磁性,且表面电势和磁势是常数。微动接触依赖载荷的加载历史,所以首先求解单独的法向加载问题,然后在法向加载问题的基础上求解循环变化的切向加载问题。整个接触区可以分为内部的中心粘着区和两个外部的滑移区,其中滑移区满足Coulomb摩擦法则。利用Fourier积分变换,磁电弹半平面的微动接触问题将简化为耦合的Cauchy奇异积分方程组,然后数值离散为线性代数方程组,利用迭代法求解未知的粘着/滑移区尺寸、电荷分布、磁感应强度、法向接触压力和切向接触力。数值算例给出了摩擦系数、总电荷和总磁感应强度对各加载阶段的表面接触应力、电位移和磁感应强度的影响。     

Stress distribution near grain boundary in anisotropic bicrystals and tricrystals

The rate dependent crystallographic finite element program was implemented in ABAQUS as a UMAT for the analysis of the stress distributions near grain boundary in anisotropic bicrystals and tricrystals, taking the different crystallographic orientations into consideration. The numerical results of bicrystals model with the different crystallographic orientations shows that there is a high stress gradient near the grain boundaries. The characteristics of stress structures are dependent on the crystallographic orientations of the two grains. The existing of triple junctions in the tricrystals may result in the stress concentrations, or may not, depending on the crystallographic orientations of the three grains. The conclusion shows that grain boundary with different crystallographic orientations can have different deformation, damage, and faUure behaviors. So it is only on the detail study of the stress distribution can the metal fracture be understood deeply.     

Towards a homogenized plasticity theory which predicts structural and microstructural size effects

We study an idealized bending problem where two types of size effects are present – one induced by the non-uniform (macro) deformation, the other due to the (internal) resistance at grain boundaries. Classical models are not able to capture either of the two types of size dependent behavior. A remedy is to adopt a gradient crystal plasticity formulation which allows one to study the direct influence of different microstructural properties on the material response. However, it is computationally expensive to do so for a typical engineering problem since the discretization has to be done at a sub-granular level. In this paper, a homogenization theory is proposed such that the small deformation gradient crystal plasticity framework by Cermelli and Gurtin [Cermelli, P., Gurtin, M.E., 2002. Geometrically necessary dislocations in viscoplastic single crystals and bicrystals undergoing small deformations. Int. J. Solids Struct. 39, 6281–6309] translates from the micro to macro level consistently. Microstructural properties thus propagate naturally to the macro scale and the homogenized solutions compare well with the fine scale analyses for the two limit cases – microhard and microfree conditions. Three length scale parameters, i.e. the intrinsic length scale, grain size and the foil thickness, manifest themselves in the homogenized solution, thus capturing both types of size effects. We further discuss on the interplay and competition between the two size effects.     

Crystallographic plasticity in fretting of Ti–6AL–4V

The deformation occurring under fretting conditions occurs over length scales of the same order as the grain size. Consequently, the crystallographic orientation of the grains plays a significant role in the deformation response. The cyclic deformation response in the region experiencing fretting predicted by a crystal plasticity model is compared to prediction of an initially isotropic J₂ cyclic plasticity theory with nonlinear kinematic hardening. The crystal plasticity model provides enhanced understanding of the fretting fatigue process, especially with regard to the shakedown and ratchetting limits.     

Microstructure-based crystal plasticity modeling of cyclic deformation of Ti–6Al–4V

Ti–6Al–4V is a dual phase material with range of possible complex microstructures. It is well known that mechanical behavior of Ti–6Al–4V is significantly affected by its texture and microstructure morphology. A three-dimensional microstructure-based constitutive model for monotonic and cyclic deformation of duplex Ti–6Al–4V is developed and implemented. The model includes length scale effects associated with dislocation interactions with different microstructure features, and is calibrated using polycrystalline finite element simulations to fit the measured macroscopic responses (overall stress–strain behavior) of a duplex heat treated Ti–6Al–4V alloy subjected to a complex cyclic loading history. Representative microstructures are simulated using a three-dimensional finite element mesh with periodic boundary conditions imposed in all directions. The measured orientation and misorientation distributions of grains of this duplex Ti–6Al–4V are considered, and similar probability density distributions of the crystallographic orientations are assigned to the finite element mesh. The misorientation distributions are then fit using the simulated annealing method. Effects of microstructural features are examined and compared with the experimental data in terms of their influence on the material yield strength. The results are shown to be in good agreement with the experimental observations.     

Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.