考虑运动硬化的DD3镍基合金力学行为研究

详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,引入了运动硬化规律,考虑了温度和应变率对晶体塑性变形的影响,通过针对每个滑移系考虑屈服准则和流动规律建立了晶体塑性模型. 对积分过程进行了推导,通过编写ABAQUS材料用户子程序(UMAT), 实现本构模型的有限元积分算法. 在此基础上模拟了DD3镍基单晶合金在单轴拉伸和循环载荷下的响应,并与实验数据进行了对比. 利用该模型可以很好地模拟镍基单晶所具有的各向异性特性,体现了镍基单晶在循环载荷作用下的拉-压不对称性.      

镍基单晶高温合金定向粗化行为及高温蠕变力学性能研究进展

镍基单晶高温合金是一种广泛应用于航空发动机和工业燃气轮机的两相叶片材料,由软的$\gamma $ 基体相和均匀镶嵌在其中的立方状 $\gamma'$ 沉淀强化相组成.它有个显著的特征,即在高温施加应力条件下, $\gamma '$沉淀相会发生定向粗化, 形成筏状.这种筏化行为直接影响了合金的蠕变疲劳寿命,是镍基单晶高温合金强化机制研究的重点. 此外,镍基单晶高温合金无晶界, 不存在高温晶界弱化、纵向晶界裂纹等问题.因此, $\gamma$/$\gamma'$相界面的位错运动、微观结构以及在载荷和温度作用下的演化决定了其蠕变力学性能.本文从镍基单晶高温合金的微观强化机制出发对定向粗化行为及蠕变力学性能进行了综述.重点介绍了定向粗化行为发生的微观机理、驱动力、影响因素和蠕变过程中界面微结构演化、蠕变力学模型以及定向粗化对高温蠕变力学性能的影响,指出了高温蠕变力学性能研究的发展方向和仍待解决的问题.      

含冷却孔镍基合金次级取向效应的应变梯度晶体塑性有限元研究

单晶镍基合金具有优异的耐高温、高强、高韧等性能, 这些力学性能受制造过程引入的次级取向和冷却孔的影响. 已有研究大多关注单孔薄板的变形机理和力学性能, 而工程中应用的往往是多孔薄板, 当前亟需阐明多孔的塑性滑移带变形机理、次级取向效应以及冷却孔引起的应变梯度效应. 文章采用基于位错机制的非局部晶体塑性本构模型对含冷却孔镍基单晶薄板的单拉变形进行了数值模拟. 此模型基于塑性滑移梯度与几何必需位错的关系引入了位错流动项, 因此可有效刻画非均匀变形过程中的应变梯度效应. 为了全面揭示含孔镍基薄板的次级取向效应, 系统研究了[100]和[110]取向(两种次级取向)下镍基薄板的单拉变形行为, 并重点探究了在两种次级取向下冷却孔数量对薄板塑性行为的影响. 此外, 还分析了镍基合金板变形过程中各个滑移系上分切应力变化、主导滑移系开动以及几何必需位错密度的演化过程, 并讨论了塑性滑移量及其分布特征对不同次级取向镍基合金板强度的影响. 研究表明, 单孔和多孔的[110]薄板抗拉强度均低于[100]薄板, 多孔薄板的塑性变形过程比单孔薄板更为复杂且受次级取向影响更大, 并且发生滑移梯度位置主要位于冷却孔附近以及塑性滑移带区域. 研究结果可为工程中镍基合金的设计和服役提供理论指导.      

镁合金晶体塑性本构模型与非均匀孪生变形分析

微结构演化对镁合金材料力学性能有着显著的影响，为了揭示镁合金宏观塑性各向异性特性与非均匀孪生变形的关系，开展了不同路径下的单轴加载试验以及采用含滑移、孪生机制的晶体塑性本构模型对试验条件下的镁合金变形行为进行数值模拟研究。文中本构模型描述了滑移与孪生变形机制以及晶格转动的机制，同时研究采用三维微结构代表性有限元模型，其包含晶粒尺寸、晶向和晶界倾角等微结构参数。研究结果表明，轧制镁合金具有强烈的宏观塑性各向异性行为，并对这种镁合金塑性各向异性行为的模拟结果以及多晶织构的模拟演化结果与试验测量进行对比，结果都基本吻合。对孪生非均匀变形模拟分析表明，镁合金宏观塑性各向异性行为与滑移、孪生变形机制的不同启动组合紧密相关，同时多晶体内应力的非均匀分布受到孪生变形的严重影响。而不同晶粒尺寸的晶粒所发生的孪生变形有比较大的差异，造成孪晶变体在晶粒内的分布极不均匀。本研究可为通过微结构的合理配置来设计和控制材料的力学性能提供理论依据.     

基于微结构演化的V5Cr5Ti合金塑性本构模型

为合理描述V5Cr5Ti合金的塑性变形行为,本文建立了基于微结构演化的塑性本构模型。首先,采用小尺寸试样开展了V5Cr5Ti合金单轴拉伸试验,并对其在不同应变程度下的微结构演化特征进行了分析。研究发现,影响V5Cr5Ti合金塑性变形行为的主要因素是位错密度演化以及团簇状和弥散析出相。据此建立了位错密度演化方程、组分相含量体积分量演化方程,并考虑团簇状和弥散状第二相对V5Cr5Ti合金流动应力的影响,进一步建立了包括非热应力、热激活应力和弥散相强化应力的流动应力关系式。最后,根据隐式应力更新算法对新模型进行了有限元实现,并与实验结果进行比较,验证了新模型的合理性和预测精度。     

基于微结构演化的V5Cr5Ti合金塑性本构模型

为合理描述V5Cr5Ti合金的塑性变形行为，本文建立了基于微结构演化的塑性本构模型。首先，采用小尺寸试样开展了V5Cr5Ti合金单轴拉伸试验，并对其在不同应变程度下的微结构演化特征进行了分析。研究发现，影响V5Cr5Ti合金塑性变形行为的主要因素是位错密度演化以及团簇状和弥散析出相。据此建立了位错密度演化方程、组分相含量体积分量演化方程，并考虑团簇状和弥散状第二相对V5Cr5Ti合金流动应力的影响，进一步建立了包括非热应力、热激活应力和弥散相强化应力的流动应力关系式。最后，根据隐式应力更新算法对新模型进行了有限元实现，并与实验结果进行比较，验证了新模型的合理性和预测精度。     

沉淀相的尺寸与形态对镍基合金蠕变行为的影响

镍基合金具有优良的高温力学性能,广泛应用于涡轮叶片等热端部件。沉淀相的尺寸和形态是影响镍基合金力学性能的重要因素。本文在考虑应变梯度的镍基合金晶体塑性本构模型的基础上,引入了各向异性损伤张量,研究了包含两种不同尺寸和三种不同长细比的沉淀相形态的镍基合金蠕变行为。结果表明,该模型能够很好地反映沉淀相的尺寸对镍基合金蠕变行为的影响,与实验结果符合较好。同时,沉淀相的形态也对镍基合金的力学性能产生重要影响,随着沉淀相长细比的增加,镍基合金的蠕变寿命延长,这体现了粗化和形态对镍基合金蠕变行为影响的一种竞争的机制。     

单晶高温合金中铸造微孔洞的扩长

基于有限变形晶体塑性本构关系及轴对称体胞模型,采用有限元的方法,分析了在不同取向偏角及不同滑移系开动,单晶高温合金中铸造微孔洞扩长的力学行为。分析结果表明:取向偏角对铸造微孔洞的扩长具有重要的意义,但对铸造微孔洞的形状改变影响不大,为改善单晶高温合金热端部件的疲劳、蠕变性能,控制晶体的取向偏角是很必要的;滑移系族开动的类型对铸造微孔洞的扩长有很大的影响,这种影响与Schmid因子、加载的取向相关,为更加准确地分析单晶高温热端部件的寿命,确定滑移系族开动类型至关重要。     

孪生诱发塑性对V-5Cr-5Ti合金应变硬化行为的影响

钒合金(V-Cr-Ti)作为潜在重要的聚变反应堆用结构材料, 近年来受到广泛的关注. 为了研究 V-5Cr-5Ti 合金不同应变率压缩下的应变硬化行为, 特别是孪生对塑性变形的影响, 以位错密度和孪晶演化为基础, 建立了该合金的应变硬化模型. 模型中考虑了孪晶中的位错滑移对材料塑性应变的贡献. 模拟结果表明, 由于孪生诱发塑性, 从而使动态压缩时的位错密度小于准静态加载时的, 这使得 V-5Cr-5Ti 合金在动态压缩时的应变硬化率比准静态加载时的小. 当孪晶形成后, 位错滑移引起的塑性应变率随应变增大而增大, 并逐渐接近加载应变率, 而孪生引起的塑性应变率则随应变增大而减小.      

单晶高温材料中铸造微孔洞的扩长

基于有限变形晶体塑性本构关系及三维体胞模型,采用有限元的方法,分析了在不同应力三维度、不同罗德参数、不同滑移系开动及不同加载取向下,单晶高温合金中铸造微孔洞扩长的力学行为。分析结果表明:累积剪切应变在铸造微孔洞的扩长中起着很重要的作用,大的累积剪切应变对应高含量的铸造微孔洞;开动滑移系族的类型对铸造微孔扩长的影响不容忽视,故准确的确定开动滑移系的类型,对于评估单晶热端部件的寿命至关重要。由于不同的取向具有不同的Schmid因子、弹性模量及开动滑移系,单晶高温合金中的铸造微孔洞的扩长还与取向密切相关,因此根据热端部件工况,合理的选择其取向是有必要的。     

Creep,plasticity, and fatigue of single crystal superalloy

Single crystal components in gas turbine engines are subject to such extreme temperatures and stresses that life prediction becomes highly inaccurate resulting in components that can only be shown to meet their requirements through experience. Reliable life prediction methodologies are required both for design and life management. In order to address this issue we have developed a thermo-viscoplastic constitutive model for single crystal materials. Our incremental large strain formulation additively decomposes the inelastic strain rate into components along the octahedral and cubic slip planes. We have developed a crystallographic-based creep constitutive model able to predict sigmoidal creep behavior of Ni base superalloys. Inelastic shear rate along each slip system is expressed as a sum of a time dependent creep component and a rate independent plastic component. We develop a new robust, computationally efficient rate-independent crystal plasticity approach and combined it with creep flow rule calibrated for Ni-based superalloys. The transient variation of each of the inelastic components includes a back stress for kinematic hardening and latent hardening parameters to account for the stress evolution with inelastic strain as well as the evolution for dislocation densities. The complete formulation accurately predicts both monotonic and cyclic tests at different crystallographic orientations for constant and variable temperature conditions (low cycle fatigue (LCF) and thermo-mechanical fatigue (TMF) tests). Based on the test and modeling results we formulate a new life prediction criterion suitable for both LCF and TMF conditions.     

On the evolution of crystallographic dislocation density in non-homogeneously deforming crystals

A set of evolution equations for dislocation density is developed incorporating the combined evolution of statistically stored and geometrically necessary densities. The statistical density evolves through Burgers vector-conserving reactions based in dislocation mechanics. The geometric density evolves due to the divergence of dislocation fluxes associated with the inhomogeneous nature of plasticity in crystals. Integration of the density-based model requires additional dislocation density/density-flux boundary conditions to complement the standard traction/displacement boundary conditions. The dislocation density evolution equations and the coupling of the dislocation density flux to the slip deformation in a continuum crystal plasticity model are incorporated into a finite element model. Simulations of an idealized crystal with a simplified slip geometry are conducted to demonstrate the length scale-dependence of the mechanical behavior of the constitutive model. The model formulation and simulation results have direct implications on the ability to explicitly model the interaction of dislocation densities with grain boundaries and on the net effect of grain boundaries on the macroscopic mechanical response of polycrystals.     

Dwell fatigue crack nucleation model based on crystal plasticity finite element simulations of polycrystalline titanium alloys

In this paper a crystal plasticity-based crack nucleation model is developed for polycrystalline microstructures undergoing cyclic dwell loading. The fatigue crack nucleation model is developed for dual-phase titanium alloys admitting room temperature creep phenomenon. It is a non-local model that accounts for the cumulative effect of slip on multiple slip systems, and involves evolving mixed-mode stresses in the grain along with dislocation pileups in contiguous grains. Rate dependent, highly anisotropic behavior causes significant localized stress concentration that increases with loading cycles. The crystal plasticity finite element (CPFE) model uses rate and size-dependent anisotropic elasto-crystal plasticity constitutive model to account for these effects. Stress rise in the hard grain is a consequence of time-dependent load shedding in adjacent soft grains, and is the main cause of crack nucleation in the polycrystalline titanium microstructure. CPFE simulation results are post-processed to provide inputs to the crack nucleation model. The nucleation model is calibrated and satisfactorily validated using data available from acoustic microscopy experiments for monitoring crack evolution in dwell fatigue experiments.     

A dislocation density based crystal plasticity finite element model: Application to a two-phase polycrystalline HCP/BCC composites

We present a multiscale model for anisotropic, elasto-plastic, rate- and temperature-sensitive deformation of polycrystalline aggregates to large plastic strains. The model accounts for a dislocation-based hardening law for multiple slip modes and links a single-crystal to a polycrystalline response using a crystal plasticity finite element based homogenization. It is capable of predicting local stress and strain fields based on evolving microstructure including the explicit evolution of dislocation density and crystallographic grain reorientation. We apply the model to simulate monotonic mechanical response of a hexagonal close-packed metal, zirconium (Zr), and a body-centered cubic metal, niobium (Nb), and study the texture evolution and deformation mechanisms in a two-phase Zr/Nb layered composite under severe plastic deformation. The model predicts well the texture in both co-deforming phases to very large plastic strains. In addition, it offers insights into the active slip systems underlying texture evolution, indicating that the observed textures develop by a combination of prismatic, pyramidal, and anomalous basal slip in Zr and primarily {110}〈111〉 slip and secondly {112}〈111〉 slip in Nb.     

On the strain hardening and texture evolution in high manganese steels: Experiments and numerical investigation

We present a systematic investigation on the strain hardening and texture evolution in high manganese steels where twinning induced plasticity (TWIP) plays a significant role for the materials' plastic deformation. Motivated by the stress–strain behavior of typical TWIP steels with compositions of Fe, Mn, and C, we develop a mechanistic model to explain the strain-hardening in crystals where deformation twinning dominates the plastic deformation. The classical single crystal plasticity model accounting for both dislocation slip and deformation twinning are then employed to simulate the plastic deformation in polycrystalline TWIP steels. While only deformation twinning is activated for plasticity, the simulations with samples composed of voronoi grains cannot fully capture the texture evolution of the TWIP steel. By including both twinning deformation and dislocation slip, the model is able to capture both the stress–strain behaviors and the texture evolution in Fe–Mn–C TWIP steel in different boundary-value problems. Further analysis on the strain contributions by both mechanisms suggests that deformation twinning plays the dominant role at the initial stage of plasticity in TWIP steels, and dislocation slip becomes increasingly important at large strains.     

A model for high temperature creep of single crystal superalloys based on nonlocal damage and viscoplastic material behavior

A model for high temperature creep of single crystal superalloys is developed, which includes constitutive laws for nonlocal damage and viscoplasticity. It is based on a variational formulation, employing potentials for free energy, and dissipation originating from plasticity and damage. Evolution equations for plastic strain and damage variables are derived from the well-established minimum principle for the dissipation potential. The model is capable of describing the different stages of creep in a unified way. Plastic deformation in superalloys incorporates the evolution of dislocation densities of the different phases present. It results in a time dependence of the creep rate in primary and secondary creep. Tertiary creep is taken into account by introducing local and nonlocal damage. Herein, the nonlocal one is included in order to model strain localization as well as to remove mesh dependence of finite element calculations. Numerical results and comparisons with experimental data of the single crystal superalloy LEK94 are shown.     

A multi-scale computational model of crystal plasticity at submicron-to-nanometer scales

Plastic flow in crystal at submicron-to-nanometer scales involves many new interesting problems. In this paper, a unified computational model which directly combines 3D discrete dislocation dynamics (DDD) and continuum mechanics is developed to investigate the plastic behaviors at these scales. In this model, the discrete dislocation plasticity in a finite crystal is solved under a completed continuum mechanics framework: (1) an initial internal stress field is introduced to represent the preexisting stationary dislocations in the crystal; (2) the external boundary condition is handled by finite element method spontaneously; and (3) the constitutive relationship is based on the finite deformation theory of crystal plasticity, but the discrete plastic strains induced by the slip of the newly nucleated or propagating dislocations are calculated by dislocation dynamics methodology instead of phenomenological evolution equations used in conventional crystal plasticity. These discrete plastic strains are then localized to the continuum material points by a Burgers vector density function proposed by us. Various processes, such as loop dislocation evolution, dislocation junction formation etc., are simulated to verify the reliability of this computational model. Specifically, a uniaxial compression test for micro-pillars of Cu is simulated by this model to investigate the ‘dislocation starvation hardening’ observed in the recent experiment.