ABAQUS中钢‑聚丙烯混杂纤维混凝土塑性本构参数的取值方法与流程

本发明属于建筑材料
技术领域：
，具体涉及一种基于ABAQUS的钢-聚丙烯混杂纤维混凝土塑性本构模型参数的取值方法。
背景技术：
：本世纪以来，纤维混凝土得到了迅猛的发展，为了准确确定其真实的本构模型，国内外学者进行了大量的研究，提出了各种荷载条件下的应力-应变数学表达式，并将其运用到有限元仿真分析当中，但模拟结果却往往不能很好地反映纤维混凝土的非线性受力行为，如基体多轴强度、体积变形、滞回能耗等力学特性在纤维加入后的变化无法得到较好的体现。目前，大型通用有限元软件ABAQUS以其强大的仿真功能而被广泛应用，其中，其内置的混凝土损伤塑性模型(CDPM)在普通混凝土结构的力学性能分析中起着关键作用，但对于纤维混凝土，尤其是混杂纤维混凝土，该模型并不十分适用，需要进行进一步的修正，但相关工作迄今为止并未见诸任何报道，这在很大程度上限制了该领域课题的研究深度，不利于混杂纤维混凝土材料的进一步推广应用。技术实现要素：为了解决上述的技术问题，本发明提供了一种基于大型通用有限元软件ABAQUS适用于钢-聚丙烯混杂纤维混凝土(下文中如无做特殊说明，均简称为“混杂纤维混凝土”)塑性本构模型参数的取值方法。ABAQUS中钢-聚丙烯混杂纤维混凝土塑性本构参数的取值方法，包括以下步骤：步骤1：提出屈服准则修正的总体形式：所述钢-聚丙烯混杂纤维混凝土塑性本构模型的屈服准则为：F=11-αhf(q‾-3αhfp‾+βhf<σ‾max>-γhf<-σ‾max>)-σc0≤0---(1)]]>其中，有效静水压力Mises等效应力“:”表示张量的点积；有效偏应力张量为应力张量；I为单位矩阵；为最大有效应力值；〈·〉为取正符号，x表示一个数值；参数αhf，βhf，γhf的计算公式如下：αhf=σb0hf/σc0-12(σb0hf/σc0)-1,βhf=σc0σt0(1-αhf)-(1+αhf),γhf=3(1-Kchf)2Kchf-1---(2)]]>式中，为混杂纤维混凝土双轴抗压强度，σc0为普通混凝土单轴抗压强度，定义了屈服面在平面应力状态下的形状；σt0为普通混凝土单轴抗拉强度；规定了当静水压力时混杂纤维混凝土的拉、压子午线之比步骤2：修正硬化规律：所述钢-聚丙烯混杂纤维混凝土塑性本构模型的硬化规律定义如下：(1)钢-聚丙烯混杂纤维混凝土单轴拉应力-应变关系为：σhft＝(1-dt)Ehftεhft(10)dt=1-ρt[α1+(1.5-1.25α1)x+(0.25α1-0.5)x5]x≤11-ρtαt(x-1)1.7+xx>1---(11)]]>α1＝1.2(1+0.265λsf+0.277λpf)1.2≤α1≤2(12)αt=0.312(ft0)21+3.366λsf+3.858λpf,0≤αt≤1.5---(13)]]>fhft0=ft0(1+0.379λsf+0.02λsfλpf),ϵhft0=ϵt0(1+0.498λsf+0.697λpf)---(14)]]>x=ϵhft/ϵhft0,ρt=fhft0/(Ehftϵhft0)---(15)]]>式中，Ehft、σhft、εhft和dt分别表示混杂纤维混凝土受拉弹性模量、受拉应力、受拉应变和受拉损伤值；ft0分别为混杂纤维混凝土和普通混凝土的单轴抗拉强度；则分别对应上述两种材料在单轴受拉条件下的峰值应变，λsf，λpf分别表示钢、聚丙烯纤维的特征值；α1为混杂纤维混凝土单轴拉应力-应变关系上升段系数，αt为混杂纤维混凝土单轴拉应力-应变关系下降段系数；ρt为混杂纤维混凝土单轴抗拉强度与峰值应变和弹性模量的比值；(2)钢-聚丙烯混杂纤维混凝土单轴压应力-应变关系为：y=ax+(3-2a)x2+(a-2)x30≤x≤1y=xb(x-1)2+xx>1---(16)]]>a=28.2283-23.2771(fhfc0)0.0374+0.4772λsf-0.4917λpf---(17)]]>b=1+0.3688(fhfc0)-0.2846-λsf-λpf---(18)]]>fhfc0=fc0(1+0.206λsf+0.388λpf),ϵhfc0=263.3fhfc0×10-6---(19)]]>x=ϵhfc/ϵhfc0,y=σhfc/fhfc0---(20)]]>式中，σhfc和εhfc分别表示混杂纤维混凝土受压应力和受压应变；分别为混杂纤维混凝土和普通混凝土的单轴抗压强度；为混杂纤维混凝土在单轴受压条件下的峰值应变，λsf，λpf分别表示钢、聚丙烯纤维的特征值；a为混杂纤维混凝土单轴压应力-应变关系上升段系数，b为混杂纤维混凝土单轴压应力-应变关系下降段系数；需要指出的是，所述(1)、(2)步骤中分别给出了单轴拉、压情况下的应力-应变全曲线，其中的应变包括了弹性应变和塑性应变两部分，而在使用ABAQUS软件进行有限元分析时只需提供应力-塑性应变全曲线，因此，再输入应变值时需要将弹性应变先行减去；步骤3：修正流动法则：所述钢-聚丙烯混杂纤维混凝土塑性本构模型的剪胀角定义如下：ψhf=arctan(|(φ1hfq‾2-φ2hfq‾1||(φ1hf)2(φ2hf)2(q‾22-q‾12)|)---(27)]]>式中，ψhf为混杂纤维混凝土的剪胀角，φhf为混杂纤维混凝土塑性流动角，Mises等效应力各参数下标“1”、“2”分别对应两种不同应力状态下的值；根据公式(27)的定义，即可计算得到混杂纤维混凝土的剪胀角与纤维特征值之间的对应关系。而且，公式(2)中，所述混杂纤维混凝土拉压子午线之比的计算方法为：Kchf=q‾(TM)hfq‾(CM)hf=q‾(TM)·ktq‾(CM)·kc=Kc·ktkc---(3)]]>式中，Kc为普通混凝土拉子午线值与压子午线值之比，kt、kc分别为混杂纤维对普通混凝土拉、压子午线的影响系数，计算方法如下：kt＝1+0.08λsf+0.132λpf，kc＝1+0.056λsf(4)式中，λsf，λpf分别表示钢、聚丙烯纤维的特征值。而且，公式(2)中，所述混杂纤维混凝土双轴抗压强度与普通混凝土单轴抗压强度之比的计算方法如下：(1)根据弹塑性理论，等压双轴抗压强度点位于拉子午线上，令混杂纤维混凝土的双轴抗压强度为(压为负)，应力状态为σ1，σ2，σ3分别为x、y、z方向上主应力，该应力状态在Haigh-Westergaard坐标系下对应的静水压力值ξ与偏应力值分别为：ξ=I13=-2σb0hf3---(5)]]>ρthf=2J2=2·16·[(0+σb0hf)2+(0+σb0hf)2+(-σb0hf+σb0hf)2]=23σb0hf---(6)]]>式中，I1为主应力张量第一不变量，I1＝σ1+σ2+σ3；J2为偏应力张量第二不变量；(2)不同静水压力下拉子午线上混杂纤维混凝土的偏应力值建议如下：ρthf=-a1-a12-4a2(a0-ξσc0)2a2·σc0·kt---(7)]]>式中，a0、a1、a2分别为拟合系数，将式(5)带入式(7)，联解式(6)和(7)，可得：8a223kt2·(σb0hfσc0)2+(42a1a2+8a2kt3kt)·(σb0hfσc0)+4a2a0=0---(8)]]>式中各拟合系数a0、a1、a2的取值建议分别为a0＝0.1775，a1＝-1.4554，a2＝-0.1576，将他们代入式(8)中，整理得到σb0hfσc0=-(0.749/kt-0.728)+(0.749/kt-0.728)2+0.03/kt20.132/kt2---(9)]]>式中，kt为混杂纤维对普通混凝土拉子午线的影响系数，计算方法如下：kt＝1+0.08λsf+0.132λpf，式中，λsf，λpf分别表示钢、聚丙烯纤维的特征值。而且，步骤2中，对于单轴受拉情况，弹性极限点应变取混杂纤维混凝土达到峰值强度时对应的应变值，因此，其塑性应变计算公式如下：ϵ~tpl=0σt≤fhft0ϵ~tpl=ϵt-σt/Ehftσt>fhft0---(21)]]>对于单轴受压情况，弹性极限点应变取峰值强度的1/3处对应的应变值，因此，其塑性应变计算公式为：ϵ~cpl=0σc≤fhfc0/3ϵ~cpl=ϵc-σc/Ehfcσc>fhfc0/3---(22)]]>式中，为单轴受拉情况下受拉塑性应变，为单轴受压情况下受压塑性应变；εt为单轴受拉情况下受拉总应变，εc为单轴受压情况下受压总应变；σt为单轴受拉情况下拉应力，σc单轴受拉情况下压应力；Ehfc为混杂纤维混凝土受压弹性模量。而且，步骤3中，公式(27)的推导过程为：(1)流动法则采用非关联流动法则其中，塑性流动势函数ghf采用D-P双曲线方程ghf=(ehfσt0tanψhf)2+q‾2-p‾tanψhf=0---(23)]]>式中，ψhf为混杂纤维混凝土的剪胀角，描述了该材料的屈服面在子午平面内与静水压力轴的夹角；ehf为混杂纤维混凝土塑性势能函数的偏心率，定义了与线性Drucker-Prager塑性势能函数的渐进程度，本方法中ehf取0.1；ehfσt0则代表了该塑性势双曲线及其渐近线在静水压力轴截距之间的距离；(2)假定混杂纤维混凝土的塑性流动在偏平面任意点处的方向一致，且总塑性应变流动方向与塑性应变增量流动方向相同，即混杂纤维混凝土塑性流动角φhf与Lode角不相关；将塑性势函数式(23)两边分别对求偏导，得：φhftanψhf(ehfσt0tanψhf)2+q‾2=-q‾---(24)]]>取两种不同应力状态对塑性势参数ψhf进行标定:φ1hftanψhf(ehfσt0·tanψhf)2+q‾12=-q‾1---(25)]]>φ2hftanψhf(ehfσt0tanψhf)2+q‾22=-q‾2---(26)]]>式中，ψhf为混杂纤维混凝土的剪胀角，ehf为混杂纤维混凝土塑性势能函数的偏心率，φhf为混杂纤维混凝土塑性流动角，Mises等效应力各参数下标“1”、“2”分别对应两种不同应力状态下的值；联立方程组(25)和(26)，消除ehf并求解ψhf可得到：ψhf=arctan(|(φ1hfq‾2-φ2hfq‾1||(φ1hf)2(φ2hf)2(q‾22-q‾12)|)---(27)]]>式中，ψhf为混杂纤维混凝土的剪胀角，φhf为混杂纤维混凝土塑性流动角，Mises等效应力各参数下标“1”、“2”分别对应两种不同应力状态下的值。而且，公式(27)中φhf=12·[(ϵ1pl-ϵ2pl)2+(ϵ1pl-ϵ3pl)2+(ϵ2pl-ϵ3pl)2]/ϵ1pl+ϵ2pl+ϵ3pl3---(28)]]>式中，分别为x、y、z方向塑性应变；将两种不同应力状态下的值分别代入式(28)，即可得到对应的和值。式中的计算公式如下：ϵ1pl=ϵ1-[σ1-v(σ2+σ3)]/Eϵ2pl=ϵ2-[σ2-v(σ1+σ3)]/Eϵ3pl=ϵ3-[σ3-v(σ1+σ2)]/E---(29)]]>式中，ε1，ε2，ε3分别为x、y、z方向总应变；σ1，σ2，σ3分别为x、y、z方向上主应力；v为泊松比；E为弹性模量。而且，根据公式(24)～(29)中混杂纤维混凝土的剪胀角与纤维特征值之间的变化规律，建议一个简洁的计算公式ψhf＝ψ0(1-0.89λsf-0.196λpf)(30)其中，ψ0为普通混凝土的剪胀角，取值范围为30°～40°之间。而且，钢纤维的特征值λsf＝Vsf(lsf/dsf)，聚丙烯纤维的特征值λpf＝Vpf(lpf/dpf)；Vsf，Vpf分别为钢、聚丙烯纤维的体积参数；lsf，lpf分别为钢、聚丙烯纤维的长度；dsf，dpf分别为钢、聚丙烯纤维的直径。本发明的优势在于：本发明基于目前广泛使用的大型通用有限元软件ABAQUS中内置的混凝土损伤塑性模型(CDPM)，充分考虑了不同纤维类型、体积掺量及长径比等因素对混凝土力学性能的影响，对模型的相关参数进行了适当地修正，使其吻合钢-聚丙烯混杂纤维混凝土的屈服准则、硬化规律和流动法则，从而能够更为精确的预测混杂纤维混凝土材料的多轴强度、体积变形，以及相应构件的承载力、滞回耗能等力学性能，为纤维混凝土结构在复杂应力状态下的精细化非线性分析提供了一个很好的参考。附图说明图1：本发明实施例的方法总体流程图；图2：本发明实施例的材料试验压应力-非弹性应变关系对比图；图3：本发明实施例的材料试验拉应力-开裂应变关系对比图；图4：本发明实施例的材料试验单轴受压循环加卸载曲线对比图；图5：本发明实施例的构件试验柱子滞回曲线对比图对比图；图6：本发明实施例的构件试验柱子裂缝实际分布对比图；图7：本发明实施例的构件试验柱子应力模拟结果对比图。具体实施方式为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限定本发明。一种ABAQUS中钢-聚丙烯混杂纤维混凝土塑性本构参数的取值方法，包括以下步骤：步骤1钢-聚丙烯混杂纤维混凝土塑性本构模型的屈服准则为：F=11-αhf(q‾-3αhfp‾+βhf<σ‾max>-γhf<-σ‾max>)-σc0≤0---(1)]]>其中，有效静水压力Mises等效应力“:”表示张量的点积；有效偏应力张量为应力张量；I为单位矩阵；为最大有效应力值；<·>为取正符号，即x表示一个数值；参数αhf，βhf，γhf的计算公式如下：αhf=σb0hf/σc0-12(σb0hf/σc0)-1,βhf=σc0σt0(1-αhf)-(1+αhf),γhf=3(1-Kchf)2Kchf-1---(2)]]>式中，为混杂纤维混凝土双轴抗压强度，σc0为普通混凝土单轴抗压强度，定义了屈服面在平面应力状态下的形状；σt0为普通混凝土单轴抗拉强度；规定了当静水压力时混杂纤维混凝土的拉、压子午线之比步骤1.1混杂纤维混凝土拉压子午线之比的计算方法为：Kchf=q‾(TM)hfq‾(CM)hf=q‾(TM)·ktq‾(CM)·kc=Kc·ktkc---(3)]]>式中，Kc为普通混凝土拉子午线值与压子午线值之比，kt、kc分别为混杂纤维对普通混凝土拉、压子午线的影响系数，计算方法如下：kt＝1+0.08λsf+0.132λpf，kc＝1+0.056λsf(4)式中，λsf＝Vsf(lsf/dsf)，λpf＝Vpf(lpf/dpf)分别表示钢、聚丙烯纤维的特征值；Vsf，Vpf分别为钢、聚丙烯纤维的体积参数；lf，lf分别为钢、聚丙烯纤维的长度，df，df分别为钢、聚丙烯纤维的直径。步骤1.2混杂纤维混凝土双轴抗压强度与普通混凝土单轴抗压强度之比的计算方法如下：(1)根据弹塑性理论，等压双轴抗压强度点位于拉子午线上，令混杂纤维混凝土的双轴抗压强度为(压为负)，应力状态为σ1，σ2，σ3分别为x、y、z方向上主应力，该应力状态在Haigh-Westergaard坐标系下对应的静水压力值ξ与偏应力值分别为：ξ=I13=-2σb0hf3---(5)]]>ρthf=2J2=2·16·[(0+σb0hf)2+(0+σb0hf)2+(-σb0hf+σb0hf)2]=23σb0hf---(6)]]>式中，I1为主应力张量第一不变量，I1＝σ1+σ2+σ3；J2为偏应力张量第二不变量；(2)不同静水压力下拉子午线上混杂纤维混凝土的偏应力值建议如下：ρthf=-a1-a12-4a2(a0-ξσc0)2a2·σc0·kt---(7)]]>式中，a0、a1、a2分别为拟合系数；将式(5)带入式(7)，联解式(6)和(7)，可得：8a223kt2·(σb0hfσc0)2+(42a1a2+8a2kt3kt)·(σb0hfσc0)+4a2a0=0---(8)]]>式中各拟合系数a0、a1、a2的取值建议分别为a0＝0.1775，a1＝-1.4554，a2＝-0.1576，将他们代入式(8)中，整理得到σb0hfσc0=-(0.749/kt-0.728)+(0.749/kt-0.728)2+0.03/kt20.132/kt2---(9)]]>值得注意的是，当kt＝1时，即不考虑混杂纤维对屈服面在平面应力状态下形状的影响，该式则退化为计算结果为与普通混凝土的建议取值相同，表明本公式的适用性强，对普通混凝土也成立。步骤2钢-聚丙烯混杂纤维混凝土塑性本构模型的硬化规律定义如下：(1)钢-聚丙烯混杂纤维混凝土单轴拉应力-应变关系为：σhft＝(1-dt)Ehftεhft(10)dt=1-ρt[α1+(1.5-1.25α1)x+(0.25α1-0.5)x5]x≤11-ρtαt(x-1)1.7+xx>1---(11)]]>α1＝1.2(1+0.265λsf+0.277λpf)1.2≤α1≤2(12)αt=0.312(ft0)21+3.366λsf+3.858λpf,0≤αt≤1.5---(13)]]>fhft0=ft0(1+0.379λsf+0.02λsfλpf),ϵhft0=ϵt0(1+0.498λsf+0.697λpf)---(14)]]>x=ϵhft/ϵhft0,ρt=fhft0/(Ehftϵhft0)---(15)]]>式中，Ehft、σhft、εhft和dt分别表示混杂纤维混凝土受拉弹性模量、受拉应力、受拉应变和受拉损伤值；ft0分别为混杂纤维混凝土和普通混凝土的单轴抗拉强度；则分别对应上述两种材料在单轴受拉条件下的峰值应变，λsf，λpf分别表示钢、聚丙烯纤维的特征值；α1为混杂纤维混凝土单轴拉应力-应变关系上升段系数，αt为混杂纤维混凝土单轴拉应力-应变关系下降段系数；ρt为混杂纤维混凝土单轴抗拉强度与峰值应变和弹性模量的比值；其他参数的物理意义与前文相同。(2)钢-聚丙烯混杂纤维混凝土单轴压应力-应变关系为：y=ax+(3-2a)x2+(a-2)x30≤x≤1y=xb(x-1)2+xx>1---(16)]]>a=28.2283-23.2771(fhfc0)0.0374+0.4772λsf-0.4917λpf---(17)]]>b=1+0.3688(fhfc0)-0.2846-λsf-λpf---(18)]]>fhfc0=fc0(1+0.206λsf+0.388λpf),ϵhfc0=263.3fhfc0×10-6---(19)]]>x=ϵhfc/ϵhfc0,y=σhfc/fhfc0---(20)]]>式中，σhfc和εhfc分别表示混杂纤维混凝土受压应力和受压应变；分别为混杂纤维混凝土和普通混凝土的单轴抗压强度；为混杂纤维混凝土在单轴受压条件下的峰值应变；λsf，λpf分别表示钢、聚丙烯纤维的特征值；a为混杂纤维混凝土单轴压应力-应变关系上升段系数，b为混杂纤维混凝土单轴压应力-应变关系下降段系数；其他参数的物理意义与前文相同。需要指出的是，上述(1)、(2)步骤中分别给出了单轴拉、压情况下的应力-应变全曲线，其中的应变包括了弹性应变和塑性应变两部分，而在使用ABAQUS软件进行有限元分析时只需提供应力-塑性应变全曲线，因此，再输入应变值时需要将弹性应变先行减去。(3)对于单轴受拉情况，弹性极限点应变取混杂纤维混凝土达到峰值强度时对应的应变值，因此，其塑性应变计算公式如下：ϵ~tpl=0σt≤fhft0ϵ~tpl=ϵt-σt/Ehftσt>fhft0---(21)]]>对于单轴受压情况，弹性极限点应变取峰值强度的1/3处对应的应变值，因此，其塑性应变计算公式为：ϵ~cpl=0σc≤fhfc0/3ϵ~cpl=ϵc-σc/Ehfcσc>fhfc0/3---(22)]]>式中，为单轴受拉情况下受拉塑性应变，为单轴受压情况下受压塑性应变；εt为单轴受拉情况下受拉总应变，εc为单轴受压情况下受压总应变；σt为单轴受拉情况下拉应力，σc单轴受拉情况下压应力；Ehfc为混杂纤维混凝土受压弹性模量。步骤3钢-聚丙烯混杂纤维混凝土塑性本构模型的剪胀角定义如下：(1)该流动法则采用非关联流动法则其中，塑性流动势函数ghf采用D-P双曲线方程，即ghf=(ehfσt0tanψhf)2+q‾2-p‾tanψhf=0---(23)]]>式中，ψhf为混杂纤维混凝土的剪胀角，描述了该材料的屈服面在子午平面内与静水压力轴的夹角；ehf为混杂纤维混凝土塑性势能函数的偏心率，定义了与线性Drucker-Prager塑性势能函数的渐进程度，本方法中ehf取0.1；ehfσt0则代表了该塑性势双曲线及其渐近线在静水压力轴截距之间的距离。(2)本方法假定混杂纤维混凝土的塑性流动在偏平面任意点处的方向一致，且总塑性应变流动方向与塑性应变增量流动方向相同，即混杂纤维混凝土塑性流动角φhf与Lode角不相关。将塑性势函数式(23)两边分别对求偏导，得:φhftanψhf(ehfσt0tanψhf)2+q‾2=-q‾---(24)]]>取两种不同应力状态对塑性势参数ψhf进行标定:φ1hftanψhf(ehfσt0·tanψhf)2+q‾12=-q‾1---(25)]]>φ2hftanψhf(ehfσt0tanψhf)2+q‾22=-q‾2---(26)]]>式中，各参数下标“1”、“2”分别对应两种不同应力状态下的值。联立方程组(25)和(26)，消除ehf并求解ψhf可得到：ψhf=arctan(|(φ1hfq‾2-φ2hfq‾1||(φ1hf)2(φ2hf)2(q‾22-q‾12)|)---(27)]]>式中，ψhf为混杂纤维混凝土的剪胀角，φhf为混杂纤维混凝土塑性流动角，Mises等效应力各参数下标“1”、“2”分别对应两种不同应力状态下的值；其中φhf=12·[(ϵ1pl-ϵ2pl)2+(ϵ1pl-ϵ3pl)2+(ϵ2pl-ϵ3pl)2]/ϵ1pl+ϵ2pl+ϵ3pl3---(28)]]>式中，分别为x、y、z方向塑性应变；将两种不同应力状态下的值分别代入式(28)，即可得到对应的和值。式中的计算公式如下：ϵ1pl=ϵ1-[σ1-v(σ2+σ3)]/Eϵ2pl=ϵ2-[σ2-v(σ1+σ3)]/Eϵ3pl=ϵ3-[σ3-v(σ1+σ2)]/E---(29)]]>式中，ε1，ε2，ε3分别为x、y、z方向总应变；σ1，σ2，σ3分别为x、y、z方向上主应力；v为泊松比；E为弹性模量。(3)根据公式(24)～(29)的定义，即可计算得到混杂纤维混凝土的剪胀角与纤维特征值之间的对应关系。为进一步地方便本方法的应用，根据上述两者间的变化规律，本方法建议了一个简洁的计算公式，即ψhf＝ψ0(1-0.89λsf-0.196λpf)(30)其中，ψ0为普通混凝土的剪胀角，取值范围为30°～40°之间。材料试验结果验证：将通过上述修正方法得到的各参数输入ABAQUS软件中，计算混杂纤维混凝土在单轴单调受拉、受压、循环拉压荷载作用下材料的响应，并将其与试验结果对比，如附图2、3、4所示，对比结果表明数值结果与试验结果吻合良好。构件试验结果验证：将通过上述修正方法得到的各参数输入ABAQUS软件中，计算混杂纤维混凝土柱在低周反复荷载作用下的构件响应，并将其与试验结果对比，如附图5、6、7所示，对比结果表明数值结果与试验结果吻合良好。本发明基于目前广泛使用的大型通用有限元软件ABAQUS中内置的混凝土损伤塑性模型(CDPM)，充分考虑了不同纤维类型、体积掺量及长径比等因素对混凝土力学性能的影响，对模型的相关参数进行了适当地修正，使其吻合钢-聚丙烯混杂纤维混凝土的屈服准则、硬化规律和流动法则，从而能够更为精确的预测混杂纤维混凝土材料的多轴强度、体积变形，以及相应构件的承载力、滞回耗能等力学性能，为纤维混凝土结构在复杂应力状态下的精细化非线性分析提供了一个很好的参考。应当理解的是，本说明书未详细阐述的部分均属于现有技术。应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。当前第1页1 2 3