多晶体结构应力分析的宏观应力响应.docxVIP

多晶体结构应力分析的宏观应力响应 单晶材料的力学性质 多晶体材料由微晶颗粒组成。矩阵的组成、相结构、尺寸、方向和各向异性是多晶体材料组成的主要因素。因此，整体意义上的多晶材料力学行为实际上是指具有一定数量的晶体晶体的平均力学行为，这是晶体中每个晶体的力学行为合成的结果。然而，在力学分析中，多晶材料的力学行为通常与连续均匀材料的力学行为进行处理。 原则上,多晶体在平均意义上的力学性质可以由单晶体力学性质计算得到.Traylor曾提出利用单晶体力学性质推算多晶体单轴拉伸力学性质的模型和方法.他在模型中假设:单晶体内的塑性变形是均匀的且与多晶体宏观塑性变形相同;单晶体各滑移系的临界分解剪应力相同且服从相同的各向同性硬化规律.于是对于初始无取向占优的多晶体材料可以利用最小变形功率法则求得多晶体应力应变曲线.Budiansky和Wu,Hill和Huchingson也都探讨过由单晶体力学性质推算多晶体力学性质的研究工作.不仅是应力应变曲线,多晶体材料的后继屈服面演化和塑性变形诱导的塑性流动各向异性等也都与多晶材料中单晶晶粒的力学性质及取向演化有关.上述这些研究主要关心多晶材料的“平均化”过程,对材料内部的变形和应力差异并不在意. 由于晶粒的随机取向,多晶体宏观上看来均匀的变形在细观上是不均匀的.实验观测测定多晶体的应变分布证实了这一点.但为了避免分析上的困难,在对材料进行力学分析时人们往往采用均匀(各向同性)和连续的假设.在此基础上关于材料“点”的应力和应变定义实际上抹去了细观尺度上材料微元内的不均匀性,忽略了多晶材料细观尺度上的微结构、各向异性性质和局部细观不均匀应力的影响.产生这种不均匀性的原因主要来自材料的微观结构和材料组成方式,它们的差异对结构的宏观响应也许微不足道,但对材料的微观过程如材料的损伤、断裂和破坏则可能是至关重要的.因为材料的破坏发端于材料局部微小区域的微结构,材料的微观应力及其分布方式对这样微小尺度区域内发生的事件的影响是不能被忽略的.Barbe等利用Voronoi模型研究了多晶体内的应力状态,计算结果表明由于晶粒的相互约束可以导致多晶体内晶粒间非常大的应力差异.但这项研究中多晶体涉及的塑性变形很小. 多晶塑性分析的基础是单晶塑性理论和计算方法,单晶体塑性有限变形的几何学、运动学及本构理论框架是由Hill等和Asaro等学者建立起来的.Peirce等建立了当前构形下单晶塑性有限变形率形式本构计算的改进显式方法.而Kalidindi等,Maniatty等,Sarma和Zacharia在对多晶材料织构演化的研究中发展了在卸载构形下的多晶塑性有限变形计算方法. 本文主要探讨多晶体内的不均匀力学过程,力图用细观尺度上晶体滑移和应力的分析来研究多晶体内由于取向和晶粒几何的随机性造成的不均匀状态(暂不讨论材料的相及相结构的影响).内容涉及:(1)随机生成的Voronoi多晶集合体模型(模型中晶粒是随机生成的多面体且具有随机的取向,假定各晶粒在相同位向上的弹性和滑移特性相同);(2)在当前构形下以Cauchy应力张量为基本变量的单晶黏塑性本构关系迭代计算方法;(3)多晶集合体内非均匀塑性变形的晶体黏塑性分析;(4)不同宏观应力状态下多晶体内的应力分布及统计特征;(5)多晶集合体的宏观响应与晶粒取向的演化. 1 多晶材料的薄观察 1.1 代表性单元的选取 为了考虑多晶体内晶粒尺寸、形状和取向的随机性,本文采用Voronoi多晶集合体作材料的代表性单元(视为宏观的材料“点”).在该集合体内,晶粒尺寸、形状和取向均作无占优趋向的随机安排,对每个单晶晶粒而言材料的弹性和非弹性性质是各向异性的.由于代表性单元包含相当多的晶粒,其整体响应或多晶集合体的宏观性质近似于各向同性. Voronoi多晶集合体模型的计算规模较大,而且非弹性计算比较耗时.出于节省机时的考虑,本文采用的代表性单元有两种:(1)平面二维模型,包含444个晶粒,44104个4节点单元472个3节点单元,53 100个节点;(2)三维立方体模型,包含200个单晶晶粒,50 144个四面体单元,9 865个节点.作为初步分析,模型中每个晶粒的单晶各向异性力学性质是相同的,不考虑晶粒中的微孔洞和微裂纹.用Voronoi方法生成的多晶集合体如图1所示.代表性单元采用宏观均匀变形的假设,在此假设下代表性单元边界表面的变形符合周期条件,保证了物体的宏观连续性. 1.2 单晶体-滑移系分解剪应变的演化 单晶晶粒的弹性和非弹性性质一般都是各向异性的.对满足小弹性应变假设的金属材料而言,扣除塑性后单晶体的弹性本构关系在晶轴坐标系下参照当前构形可以描述为 或 和分别是Cauchy应力、四阶各向异性弹性本构张量和对数弹性应变,△σcrys和是相应的增量.式(1)中全量和增量表述可以对应于不同的求