20108溷凝土随机损伤力学教育部暑期学校讲座,长沙精选

混凝土随机损伤力学 研究进展,李 杰,同济大学 2010年8月2日,作为一种主导性建筑材料，混凝土在我国重大工程建设中发挥着重要支撑作用。20世纪90年代以来，我国混凝土工程一直保持持续增长势头。,混凝土材料与混凝土结构,至2009年，我国混凝土用量已达 每人每年5.5吨的水平。,混凝土材料与混凝土结构研究,1997至2007年,我国在土木工程领域发表刊物研究论文约66100篇，其中，涉及混凝土材料、混凝土结构和混凝土力学的研究达29600篇，占上述研究论文总数的45%。,如何正确地反映混凝土结构的受力力学行为、尤其是动力荷载作用下的力学行为，仍然是呈现在研究工作者面前的巨大挑战。,1. 混凝土结构分析中的关键科学问题,与混凝土结构分析相关的关键科学问题,混凝土损伤过程中非线性与随机性及其耦合效应 动力加载条件下混凝土损伤的率敏感性 结构滞回耗能的物理机制 损伤扩散与多尺度随机涨落 结构非线性反应的随机演化规律,1)混凝土损伤过程中的非线性与随机性,混凝土损伤过程中 典型的非线性表现,混凝土损伤过程中 典型的随机性表现,在损伤及其演化进程中，非线性与随机性互相耦合! 如何正确认识和科学反映这些特征？,2)动力加载条件下混凝土损伤的率敏感性,细观物理与科学反映方式？,3)结构滞回耗能的物理机制,损伤机制、受力形式与结构滞回耗能机制的关系？,3)损伤扩散与多尺度随机涨落,动力作用下损伤出现弥散！,研究在不同尺度上的损伤扩散、转移、随机涨落及其跨尺度影响规律，是理解结构破坏规律、发现复杂结构整体可靠性形成机制，实现重大工程结构灾变过程控制的重要基础。,4)结构非线性反应的随机演化规律,随机性与非线性的耦合，随机性的输运、演化与涨落，是复杂工程结构性态设计与控制中必须面对的根本问题。,科学问题聚焦,混凝土的本构关系 结构非线性的正确反映 随机性的传播规律,2、混凝土随机损伤本构关系研究,混凝土材料受力非线性特征,经典力学对混凝土本构行为的反映方式,1、弹性力学模型,经典力学对混凝土本构行为的反映方式,3、断裂力学模型,从单个微裂缝及其扩展机制入手研究混凝土的非线性，建立了一套描述混凝土中单个微裂缝及其稳定性的方法； 将裂缝扩展所需能量在给定体积上积分，可以得到本构关系。,经典力学对混凝土本构行为的反映方式,？,断裂力学基本模型,混凝土中的缺陷,研究背景 经典力学（弹性力学、塑性力学、断裂力学等）很难全面反映混凝土的复杂力学性质（例如：刚度退化和强度软化 、单边效应 、拉压软化效应)。 20世纪80年代中期以来，混凝土损伤力学的研究为混凝土本构关系的研究开辟了新的道路。,损伤力学模型的基本原理,在未受损伤部分，虎克定律依然成立,当考虑塑性变形时,一些国际知名的混凝土损伤力学模型,Mazars模型(1984, 1986)：弹性损伤模型，未反映塑性变形的影响； Resende模型(1987) ：弹塑性损伤模型，损伤演化法则完全依赖于试验； Ju模型(1989)：各向同性弹塑性损伤模型，不能反映混凝土材料的单边效应； Lubliner模型(1989)：各向异性塑性-损伤模型，纯经验的损伤演化法则，在Cauchy应力空间考虑塑性应变的演化； Faria 模型(1998)：双标量弹塑性损伤模型，纯经验的损伤演化法则，经验的塑性应变演化法则，未考虑塑性变形对损伤能的影响； Comi 模型 (2001) ：双标量弹性损伤模型，纯经验的损伤演化法则,未考虑塑性变形的影响.,经典理论的局限性 经典连续介质损伤力学对于混凝土本构关系的研究基于唯象学的原理建立混凝土的损伤本构模型，并以确定性损伤变量的形式反映在外部作用下混凝土内部损伤的发展以及损伤对混凝土应力-应变关系的影响。 对于带有根本性的损伤演化法则，经典质损伤力学理论在逻辑上只能采用经验归纳或理论假设的方式给出，而难以说明损伤演化的内在物理机制。,研究切入点 混凝土材料在其形成之初，内部即已存在着随机分布的初始缺陷（如初始裂纹、空隙和微缺陷），从而形成初始随机损伤场。在荷载等外界作用下，结构损伤的演化与初始的损伤状态密切相关，在物理本质上是一个随机的损伤演化过程，具有很强的离散性和不确定性。 概率观点的引入即随机损伤本构关系的研究将有助于更合理地反映混凝土破坏的本质特征。也有助于说明损伤演化的内在物理机制。,1 ） 确定性弹塑性损伤本构关系,从基本的材料损伤物理机制出发选取损伤变量； 同时反映弹性损伤与塑性变形两种机制； 利用连续介质力学和热力学基本原理建模； 全面反映混凝土材料的非线性特性：刚度退化和强度软化 、单边效应 、拉压软化效应、有侧限时材料强度和延性提高等。,李杰，吴建营，2005，“混凝土弹塑性损伤本构模型研究：基本公式”，土木工程学报，38（9）：1420 Wu Jianying, Li Jie * & Faria R., 2006, “An energy release rate-based plastic-damage model for concrete”, International Journal of Solids and Structures, 43(3), 583-612.,基于能量建模的基本原理,建模基本过程（1）,根据材料损伤的物理机制，选取受拉损伤变量和受剪损伤变量描述基本的损伤机制。与之相应，引入有效应力张量的正、负分解,材料的弹性Helmholtz自由能,弹性损伤本构关系,退化：单轴受拉和单轴受压应力应变关系,建模基本过程（2）,在有效应力空间的塑性力学理论基础上，定义塑性Helmholtz自由能势的显式表达式，得到混凝土材料的弹塑性Helmholtz自由能势。,弹塑性自由能,塑性自由能,建模基本过程（3）,根据弹塑性Helmholtz自由能，根据热力学基本原理给出本构关系。在等温绝热条件下， Clausius-Duheim不等式为,由材料的弹塑性Helmholtz自由能势得到损伤能释放率表达式，基于损伤能释放率，建立材料的损伤准则和损伤变量的演化法则.,建模基本过程（4）,损伤能 释放率,损伤准则 与演化法则,与实验结果的对比,单向拉伸,单向压缩,双向压缩,低周反复,强度准则,弹塑性损伤 本构关系 基本框架,2）混凝土弹塑性随机损伤本构关系,经典连续介质损伤力学在本质上属于唯象学的研究框架。因此，对于带有根本性的、具有物理内涵的损伤演化法则，在逻辑上只能采用经验归纳或理论假设（猜想）的方式给出。,损伤变量演化规律,的具体形式?,经验表达,对经典研究的反思,典型损伤演化曲线,在损伤发展过程中，为什么相同的应变增量会导致不同的损伤增量？在不同的损伤阶段，为什么会出现能量的非均匀耗散？,均匀分布时的损伤演化,非均匀分布时的损伤演化,细观物理：随机性导致了损伤发展的非线性,细观弹簧模型：细观断裂应变为随机变量,基于细观随机断裂-滑移模型的混凝土弹塑性损伤本构关系,损伤演化的非线性源于细观层次断裂应变分布的随机性，宏观的损伤演化规律，应该在细观层次的物理分析中寻求其内在机理与建模途径。,李杰，张其云， 2001，“混凝土随机损伤本构关系研究”，同济大学学报，29（10），1135-1141 李杰，卢朝辉，张其云， 2003，“混凝土随机损伤本构关系-单轴受压分析”，同济大学学报，31（5），505-509,在受拉破坏机制中，损伤源于是骨料和水泥砂浆之间界面开裂、或集料及凝胶体中初始微缺陷的扩展，表现为细观受拉弹簧的随机断裂。在以压应力为主的受剪破坏机制中，损伤起源于骨料和水泥砂浆之间界面或初始微缺陷因剪应力作用导致的界面拉开，在模型中表现为受剪弹簧的随机断裂。同时，为了反映水泥砂浆内部以及砂浆、骨料界面间塑性滑移的影响，分别引入细观拉伸塑性变形元件和细观剪切塑性变形元件。基本的细观受拉损伤单元或受剪损伤单元分别由拉伸微弹簧和拉伸塑性变形元件或剪切微弹簧和剪切塑性变形元件串联组成。,随机断裂-滑移模型,损伤变量,细观物理,例：混凝土受拉随机损伤,混凝土劈裂试验示意,声发射传感器位置图,细观微弹簧 应力-应变关系,混凝土声发射试验,能量释放分布图原始数据,能量积分图,单轴受拉声发射能率与时间关系,单轴受拉声发射能率与应变关系,细观损伤单元极限破坏应变随机场的建模,损伤单元破坏概率,损伤单元断裂能量,随机场的建模,随机损伤本构模型与单轴实验结果的对比,单轴受拉,单轴受压,向三维空间损伤演化的扩展: 能量等效应变,受拉损伤能释放率与受剪损伤能释放,上述两式说明：将一维受力损伤能释放率中的弹性应变放大倍，即可得到多维受力状态的损伤能释放率，称 为能量等效应变。根据损伤一致性条件：对于初始损伤相同的两种受力状态，若损伤能释放率相同，则相应损伤相等。因此，通过引入能量等效应变，可以将一维受力状态的损伤演化法则，应用于多维受力状态。,能量等效应变的引入,Jie Li，Xiaodan Ren, 2009，“Stochastic Damage Model for Concrete Based on Energy Equivalent Strain”, International Journal of Structures and Materials, 46（11-12) 2407-2419 .,随机损伤模型对双轴应力强度和本构关系的预测,双向受拉 1：0.55,双向受压 -0.52 : -1,一拉一压 0.1 : -1,双轴强度,确定性损伤模型与随机损伤模型的联系,确定性损伤本构关系与随机损伤本构关系所反映的是同一个物理规律，其间唯一的区别在于损伤变量是取确定性变量还是随机变量。即,若取:,得到的是弹塑性随机损伤本构关系。而当采用统计平均方式、以D 的均值作为损伤演化准则,即取,将得到确定性弹塑性损伤本构关系。,3、结构非线性反应分析原理,结构非线性的基本类型,几何非线性(主要与弹性性质的大变形有关) 物理非线性(主要与材料的损伤、破坏有关),几何非线性,材料非线性,结构非线性全过程分析的必要性,线性分析:不能反映真实的结构内力分布规律!,结构非线性全过程分析原理,结构非线性全过程分析的过程,增量变刚度方法,线性分析与非线性分析根本区别: 线性分析不考虑材料强度的影响!,不同的塑性铰分布方式对应不同的结构承载力和结构延性！,材料强度的影响,问题：怎样在结构层次非线性反应分析中考虑混凝土强度乃至本构关系的随机性？,经典结构非线性全过程分析的困境,4、混凝土结构随机非线性反应分析,系综特征 数字特征值传递关系,样本轨道：确定性物理力学规律,广义概率密度演化方程,反映客观世界的三种基本途径！,概率密度演化,在一般意义上建立了随机函数在系统特征值意义上的近似传递关系。,系综意义上的物理关系,关于随机变量在均值处作级数展开,随机性的传播：数值特征的传递,初始随机源,输出概率分布,动力系统,随机性的传播：密度演化,联系输出与输入概率分布的是概率密度演化方程！,概率守恒原理的随机事件描述,在概率保守系统中，同一样本轨道或一族样本轨道在演化过程中概率守衡。,随机系统中的广义概率密度演化方程,随机动力系统分析,线性系统,非线性系统,概率密度演化方程,线性系统,非线性系统,随机系统中的广义概率密度演化方程,随机静力非线性分析,广义概率密度演化方程,结构非线性反应分析理论的对比,概率密度演化方程的求解, 结合确定性结构分析与有限差分方法求解。,密度演化方法的实施：数值算法基本步骤,I. 随机参数的离散；,II. 进行确定性分析（常规动力分析）；,III. 利用有限差分求解；,IV. 进行数值积分.,Dimension-reduction via mapping,Strategy via tangent spheres,Strategy via number theoretical method