微平面模型与经典混凝土理论模型的比较.docx

《微平面模型与经典混凝土理论模型的比较.docx》由会员分享，可在线阅读，更多相关《微平面模型与经典混凝土理论模型的比较.docx（14页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、微平面模型与经典混凝土理论模型的比较 摘要： 针对在混凝土非线性分析时所采纳的宏观本构模型中，没有考虑混凝土软化阶段或考虑不当，以及在多轴受力状况下宏观模型不能很好地模拟混凝土受力状态的问题，基于微平面理论，利用MATLAB编写混凝土材料的非线性本构关系程序，从理论和计算结果两方面与经典的基于应力张量及其不变量的本构模型进行比较. 比较结果表明，在弹性变形阶段，经典理论模型和微平面模型计算结果与试验结果无明显区分；在非弹性变形阶段，经典理论模型与微平面模型有较大的差异，后者与试验结果吻合较好. 微平面理论在混凝土三向受力和软化阶段更能精确地描述混凝土材料本构关系，且本构关系更简洁、清楚. 关键

2、词： 混凝土； 非线性分析； 微平面； 本构关系； 动态约束 中图分类号： TU313 文献标记码： A Comparison between microplane model and classical concrete theoretical model JU Jinxian, LU Hua, WANG Shuiping (Naval Representative Office, Shanghai Power Station Auxiliary Equipment Works Co., Ltd., Shanghai 200090, China) Abstract： As to the pr

3、oblems that the softening stage is not considered or considered unsuitable in nonlinear analysis of concrete by macroscopic constitutive model, and macroscopic constitutive model can not perform a good simulation on concrete stress states under multi-axial force-bearing situation, the nonlinear cons

4、titutive relationship program is fulfilled in MATLAB based on microplane theory. Then the theory and the computing results are compared with the classical constitutive model based on stress tensor and its invariant. The comparative results show that, there is no significant difference between the te

5、st results and the calculation results of classical theory model and the microplane model in elastic deformation stage; the difference is larger in non-elastic deformation stage, and the latter consists with the test results better. The constitutive relation of concrete material can be described mor

6、e accurately in three direction force-bearing and softening stage, and the constitutive relation is more simple and clearer. Key words： concrete; nonlinear analysis; microplane; constitutive relation; dynamic constraint 0 引 言 混凝土是土木工程中应用最多、也是最重要的建筑材料之一.近年来，混凝土结构非线性分析变得日益重要，而非线性分析结果的正确与否则依靠于描述混凝土非线性性

7、能的本构关系和破坏准则.迄今为止，已有的本构模型基本上分为两类1-3：(1)宏观现象学模型，包括非线性弹性模型、弹塑性模型、损伤力学模型和内蕴时间模型等，其本构关系用宏观连续介质中应力、应变张量及其不变量表示；(2)微观力学模型，包括微平面模型和粒子模型等，其本构关系用微观上的应力、应变关系表示.前者是经典方法，在混凝土本构模型的发展中发挥过巨大作用，但其发展进入瓶颈阶段，随着工程要求的提高，经典模型很难满意工程精度的要求，因此很多学者将探讨方向转向微观力学模型. 混凝土应力-应变全曲线反映混凝土基本力学特性，是探讨钢筋混凝土结构强度和变形的主要依据之一，特殊是应力-应变曲线的下降段(或称为软

8、化段)，对构件的弹塑性全过程分析、极限状态下的截面应力分布、抗震结构的延性和复原力特性等有较大影响.1目前，在混凝土结构非线性分析宏观本构模型中未考虑混凝土的软化阶段或考虑不当，同时，在多轴受力状况下也不能很好地模拟混凝土的受力状态.BAANT等3-8用微平面理论模拟混凝土的应力-应变全曲线，可以克服上述问题，并且该理论建立的微平面模型采纳统一的公式描述应力-应变全曲线，具有概念简洁、可变材料参数少等特点. 1 微平面模型简介 微平面指在宏观连续介质中任一点的无穷小领域内，垂直于任一方向的平面.微平面理论的实质是：任一点的宏观应力-应变关系是该点全部微平面上的应力-应变关系的叠加，描述微观结构

9、不同方向平面上微裂缝的相互作用4-5.混凝土微观结构见图1，二十面体微平面模型见图2. 微平面理论的最初思想来自于Taylor，他建议多晶体金属的本构关系可以通过材料中不同平面上的应力-应变关系确定，在静约束或动约束下，宏观应力或应变张量等于全部平面上应力或应变矢量之和.11013年，BAANT等4首次将Taylor的思想用于混凝土等具有软化特征的脆性材料，并提出微平面模型. 由于混凝土等准脆性材料的非弹性性能是由微裂缝而不是由塑性滑移确定的，对混凝土引入动态约束(不同平面上的应变为宏观应变张量的投影重量)构造本构关系4，使得微平面模型可以比较简单地模拟软化问题.经过不断改进和发展，混凝土微平

10、面本构模型演化出5代模型. 微平面理论的5个假设3基础为： (1)任何微平面上的应变等于宏观应变张量ij的投影重量，表示动态约束； (2)微平面既承受正应变N，又承受剪切应变Ti，且Ti与Ti方向相同； (3)每个微平面上的响应都与平均侧向应变有关； (4)在不出现卸载的状况下，每个微平面的应力-应变关系与路径无关； (5)在每个微平面上，体积应力与体积应变的关系VV，偏应力与偏应变的关系DD以及剪切应力与剪切应变的关系TT互不耦合. 2 微平面模型与经典理论模型的理论比较 2.1 微平面模型本构关系与经典模型本构关系的异同依据假设(1)(5)，微平面上的应力-应变关系3为： V=CV(V)V

11、 D=CD(D)D T=CT(T,C)T 式中：CV，CD和CT为弹性常数；CV(V)表示CV是V的函数；C为宏观应力不变量，对剪切刚度有较大影响，其值的大小与其次、第三主应力值有关.由于微平面上的V与V，D与D和T与T互不耦合，而经典理论模型在非弹性阶段各应力-应变理论上相互耦合2，因此表达式相对困难. 2.1.1 弹性阶段应力-应变关系 BAANT等3,6-7依据微平面上增量型应力-应变关系，在体积重量与偏量重量分别的状况下，应用虚功原理，令微观和宏观上的应力在各自对应的虚应变上所做的功相等，得到增量型体积重量同偏量重量分别的宏观应力-应变与微平面上应力-应变之间的联系，经过积分并与胡克定

12、律相比较，最终得到各个微平面上的弹性常数 C0V=E(12) C0D=C0T=E(1+) 式中：C0V，C0D和C0T表示弹性常数；E为弹性模量；G为剪切弹性模量；为泊松比. 上述理论推导说明，在弹性阶段微平面上的应力-应变关系与宏观理论模型一样，都是用弹性模量和泊松比表征应力-应变的关系，并且宏观表征的应力-应变关系相同. 2.1.2 非弹性阶段应力-应变关系 11016年，BAANT等6首先提出应力-应变边界的概念；文献7给出M4模型微平面上应力-应变边界，包括水平边界、正应力边界、体应力边界、偏应力边界和剪切应力边界等. 图 3 应力-应变边界的作用方式 Fig.3 Action mod

13、e of stress-strain boundary应力-应变边界的作用方式见图3.在边界内，将每个加载增量步看作是弹性的，但规定应力不能超出边界；只有当前应变增量与应力同为正值或同为负值时，加载才沿边界流淌；否则，发生卸载.当应力达到边界时斜率急剧下降，但由于不同的微平面在不同时刻进入加载或卸载阶段，因此宏观响应连续. 宏观经典模型一般是在应力达到某一值时材料屈服，即为应力确定的屈服极限；而应力-应变边界的作用事实上可看作是应变确定的屈服极限，这种方法在经典的宏观塑性表述中很难实现.边界的表达式中牵涉到宏观应力-应变的很多重量，而在微平面上只会出现很少的几个应力、应变重量,因此，应力-应变

14、边界简化为弹性阶段到软化破坏阶段的分析. 2.2 微平面模型和宏观经典模型的计算流程 微平面模型和经典模型的计算流程6见图4. 经典理论模型干脆给定宏观应力和宏观应变之间的关系以建立材料的本构关系，而微平面模型需经过动态约束微观应变微观应力虚功原理这个困难过程，目的是解决宏观经典模型不能精确描述应变软化阶段的混凝土应力-应变关系的问题.微平面可以看作是骨料间的接触面，事实上为“薄弱面”，在充分考虑混凝土中裂缝的基础上模拟应变软化阶段的应力-应变关系. 2.3 微平面模型的优缺点 与经典的基于张量不变量的本构模型相比，微平面模型的优势8为: (1)基于矢量而不是张量的本构关系，概念更清楚、简洁，

15、且在微平面上仅有少数的应力和应变重量，表达式相对于宏观本构模型简洁得多.另外，微平面上的应力-应变关系只针对一个平面，不会发生旋转，因此微平面模型不像宏观模型那样必需满意张量不变性条件，而是在对全部微平面方向进行积分时自动满意. (2)虽然每个微平面的应力-应变关系与路径无关，但是宏观应力-应变关系是全部微平面应力-应变关系的叠加.不同的宏观加载路径将导致微平面应力-应变关系的不同叠加，所以宏观应力-应变关系表现为与路径有关，且不同微平面上的加载过程和卸载过程进行组合，会产生困难的路径相关性和多样性，该过程比宏观模型操作起来更便利. (3)微平面模型采纳应力-应变边界的概念，易于考虑当前的屈服

16、极限与应变重量的相关性，而不用像宏观模型那样采纳标量的硬化、软化参数. (4)与张量式塑性本构模型相比，塑性微平面模型等价于大量同时活动的屈服曲面，且每个屈服曲面都有明确的形式. 微平面理论的缺点在于，大量的微平面和边界条件增加计算的负担，但随着计算机性能的不断提高，该缺点将得到解决. 3 经典理论模型与微平面模型计算结果 3.1 有限元模型 采纳Marc分析混凝土在静水压力状态下的应力-应变关系.为便于分析计算，模型选取简洁的长条形试件，几何尺寸为1 m1 m1 m，采纳三维八节点实体单元，单元尺寸为0.25 m. 3.2 经典计算理论 (1)von Mises 理论9.破坏面在空间上表现为

17、一个圆柱体，且与静水压力轴平行. (2)Buyukozturk准则10.破坏面为 f=3J2+3J1+J212=0 式中：为固定值0.2；一般取3；为等效应力；J1和J2为偏应力常量. 3.3 静水压力状况计算结果 3.3.1 经典理论计算结果 11013年，Green和Swanson进行混凝土在静水压力作用下的试验.11试验中，混凝土的弹性模量E=35 163 MPa，强度f c=48.4 MPa.在Marc中取相同的弹性模量和强度值，得到计算结果与试验数据比较见图5. 3.3.2 微平面模型计算结果 用MATLAB编写二十面体微平面非线性本构关系程序，采纳单向受拉和单向受压受力状况验证程序

18、的正确性.对于单向拉伸状况，与Peterson的单向拉伸状况下混凝土应力-应变关系的试验数据11进行对比，见图6. 微平面模型中的参数k1是特别重要的自由参数，可通过变更其值来对应大多数不同强度的混凝土材料.以单向受拉状况为例，固定其他自由参数，通过变更k1的值分析随k1改变的应力-应变之间的改变规律，k1取不同值时对应的应力-应变曲线见图7.可知，k1值的变更基本不影响应力-应变曲线的形态，只是变更应力峰值的大小和应力峰值对应的应变值，对于依据不同的混凝土材料调整参数值进行计算特别有利. 对于单向受压状况，弹性模量E=32 173 MPa，与van Mier的单向压缩状况下混凝土的应力-应变

19、关系试验数据11进行对比，见图8. 通过对单向受拉和单向受压状况的验证，初步证明微平面理论的正确性，并且可以看出微平面的本构关系可较精确地模拟混凝土的软化阶段.进一步给出静水压力下的微平面应力-应变关系，给定增量步的应变，分别给定轴向3个方向的应变，利用显式算法模拟静水压力下的载荷工况.给定自由参数k2=160，k3=10，k4=150，计算与试验数据相同弹性模量(E=35 163 MPa)的应力-应变关系，静水压力试验数据与微平面模型计算结果比较见图9. 3.3.3 结果分析 (1)对于弹性变形阶段，依据胡克定律，体应变与体应力关系为 V=E12V (4) 可知，经典理论模型与微平面模型计算

20、结果和试验相同，无明显区分. (2)对于非弹性变形阶段，经典理论模型与微平面理论模型有较大差异. von Mises理论的屈服面沿静水压轴的方向是开口的，不存在帽盖；而静水压力加载是沿负静水压轴方向进行的，导致加载恒久不会达到屈服面.所以，采纳von Mises理论计算的结果表明，混凝土的反应始终处于线弹性阶段，应力-应变呈现线性关系；Buyukozturk准则的屈服面在子午面上的投影为 一个封闭椭圆，其与负静水压力轴的交点即为静水压力加载下的屈服载荷，因此当采纳Buyukozturk准则时，计算结果表明混凝土在很小的静水压力作用下很快屈服破坏.由试验可知，混凝土在静水压力下，初始表现为线性反

21、应，之后表现为非线性反应，因此上述两种状况都与试验数据不符，这也说明现有经典理论的问题. 微平面模型的计算结果显示，应力-应变的曲线明显分为两段：一段呈线性关系；另一段随着应变的增加呈非线性关系，且急剧上升，几乎呈现指数上升趋势，这个应力-应变的改变趋势与试验吻合得很好.从试验数据可知，混凝土在静水压力条件下可承受很大的作用力，计算结果也与该现象吻合得很好.另外，须要说明的是，应力-应变曲线的一些特别拐点与材料的性能亲密相关，对于不同的材料有不同的结果，需调整可变参数k1和弹性模量E，有时须要调整其他的可变参数来拟合实际状况. 总之，用二十面体模拟微平面模型得到的静水压力下应力-应变关系的改变

22、趋势与试验值吻合得较好，更加详细的数值模拟则完全可以通过调整可变参数和不同材料弹性模量的值来实现.虽然调整参数的过程相对繁琐，但是完全可以得出与试验数据相符的结果. 4 结束语 与传统的宏观张量本构模型不同，微平面本构模型属于细观本构模型.传统的宏观本构模型采纳宏观连续介质的应力、应变张量并用其不变量表示材料的本构关系，而细观本构模型通常采纳一些合理的简化模型描述应力-应变关系.本文通过理论与计算结果的比较可知，微平面本构模型不但能较好地描述混凝土材料应变软化阶段和三向受力状况下的力学行为，而且不须要考虑张量不变量的限制(积分时自动满意)，并且微平面模型用不同方向微平面上少量的应力、应变重量之

23、间简洁的关系定义材料的本构特性，使建模更便利、概念更清楚，为微平面理论的工程应用供应参考. 参考文献： 1 过镇海, 时旭东. 钢筋混凝土原理和分析M. 北京: 清华高校出版社, 2003: 5-37. 2 张盛东. 混凝土本构理论探讨进展与评述J. 南京建筑工程学院学报: 自然科学版, 2002(3): 44-53. ZHANG Shengdong. State of arts on constitutive equations of concreteJ. J Nanjing Architectural & Civil Eng Inst: Nat Sci, 2002(3): 44-53. 3

24、 BAANT Z P, PRAT P C. Microplane model for brittle plastic material: I. TheoryJ. J Eng Mech, 11018, 114(10): 1689-1732. 4 BAANT Z P, OH B H. Microplane model for fracture analysis of concrete structuresC/Proc Symp Interact NonNuclear Munitions Struct Colorado Springs, Colo, Air Force Academy, 11013:

25、 4953. 5 BAANT Z P, OH B H. Microplane model for progressive fracture of concrete and rockJ. J Eng Mech, 11015, 114(4): 559582. 6 BAANT Z P, XIANG Y, PRAT P C. Microplane model for concrete. I: stressstrain boundaries and finite strainJ. J Eng Mech, 11016, 122(3): 245254. 7 BAANT Z P, CANER F C, ADL

26、EY M A, et al. Microplane model M4 for concrete. I: formulation with workconjugate deviatoric stressJ. J Eng Mech, 2000, 126(9): 944953. 8 BAANT Z P, CANER F C. Microplane model M5 with kinematic and static constraints for concrete fracture and anelasticity. I: theoryJ. J Eng Mech, 2022, 131(1): 314

27、0. 9 陈惠发, 萨里普A F. 混凝土和土的本构方程M. 余天庆, 译. 北京: 中国建筑工业出版社, 2004: 2325. 10 BUYUKOZTURK O. Nonlinear analysis of reinforced concrete structuresJ. Computers & Structures, 11017(7): 149156. 11 CANER F C. Computational modeling of damage and fracture in concreteD. Evanston: Northwestern Univ, 2000. (编辑 陈锋杰) 第14页 共14页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页第 14 页 共 14 页