还剩12页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
abaqus混凝土粘性系数取值范围混凝土非粘性塑性损伤模型的基本要点介绍如下应变率分解对率无关的模型附加假定应变率是可以如下分解的E=护+加,亡是总应变率,〃,是应变率的弹性部分,金加是应变率的塑性部分应力应变关系应力应变关系为下列弹性标量损伤关系a=(1-(£-Ep,)=D:(£-M).其中是材料的初始(无损)刚度,D=(l-,/)D3是有损刚度,”是刚度退化变量其值在0(无损)至!J1(完全失效)之间变化,与失效机制(开裂和压碎)相关的损伤导致了弹性刚度的退化在标量损伤理论框架内,刚度退化是各向同性的,它可由单个标量d来描述按照传统连续介质力学观点,有效应力可定义如下a=:(£一叫).Cauchy应力通过标量退化变量(d)转化为有效应力O=(1—d)a.对如任何一个给定的材料截面,因子(1―⑴代表承力的有效面积占总截面积的比重(总截面积剪除受损面积)在无损时d=0,有效应力等于cauchy应力然而,当损伤发生后,有效应力比cauchy应力更能代表实际情况,因为损伤后截面承力的是有效无损的面积因此,可以很方便的用有效应力来建立塑性相关公式正如后面将要谈论的那样,退化变量的演化是由一组硬化参数和有效应力控制的即d=献)硬化变量]cF{ff.£pl=G-3»+髀行max-7-^max一万式即<,
4.
5.2-11式中八和7是无量纲材料参数一
1.T〃=一弓I•J时有效静水压力,是Mises等效应力,S=〃I+G是有效应力张量的偏量部分,而沆心乂是行的代数最大主值,函数8洌=-Q-l+.仆©形式如下历4式中访和式分别为有效拉压内聚力在双轴受压时,万mm*=方程
4.
5.2T1就退化为Drucker-Prage屈服条件,材料系数可由单轴受压强度和双轴受压强度比值给出060—0cO・Q=-----------------------2T/»n—7/0一般材性试验给出的单双受压强度比值在
1.10-
1.16之间,那么C取值在
0.08-
0.12之间Lubliner etal,,1989系数:只在三维受压时才出现在公式中,它可以通过比较沿拉压子午线的强度比值得到根据定义拉子午线是满足主应力空间中抗U”=司的二-3的轨迹线,而压子午线是满足抗=益〉方2=方3的轨迹线其中包,i和是应力主值显然易求得,沿拉压子午线其表达式为抗uaxTM=初一力,1—Pl./^inax\\1=jo当抗…时,响应的屈服准则为:I+1G―6+3P=1-呢,TMQ}+1C-h+3a=1-a
7.CMcK=7+3令人[,〃为静水压力,那么就有一27+3;事实上大多数=GTM/7CZ,31-A=c试验也并没有证明八7是变化的,因此就可求出‘-2KLi:对于混凝土来说一般取Ac=小,那么=
3.当抗““X时-,沿拉压子午线的屈服函数就简化为|+1—3+3ap=1-ad.TMc+1G—/3+3ap=1-aa.CMc=3+3K同理令A7=r.\iA7wi,那么23+3在偏片面上典型的屈服面见图
4.
5.2-4,图
4.
5.2-5是平面应力时的屈服面图
4.
5.2-4对应于不同的A「值在片平面内的屈服面C.M.图
4.
5.2-5平面应力时的屈服面Auniaxial tension1人与何-3a「+而2=DAuniaxial compressionbiaxialtension1_一z(N q-3ap7T(q-3aP)=aoobiaxial compression流动法则本模型取的是非关联流动法则:塑性势G取为Drucker-Prager双曲函数的形式G=\J”o tan日产+1—ptan,
4.
5.2-12式中是p-q面内高围压时的膨胀角,内是单轴抗拉强度,是势函数偏心率,它描F述势函数向其渐近线逼近的速度(当偏心率趋于零时,流动势函数趋于直线)流动势函数的连续光滑性保证了流动方向的唯一性当围压很高时流动势函数渐近于线性Drucker-Prager势函数,且与静水轴的交角是90度在“Models forgranular orpolymerbehavior,“Section
4.
4.2,中对这个势函数有详细的讨论因为采用了非关联流动法则,刚度矩阵将会出现非对称粘塑性归一化在隐式分析程序里,当材料模型出现软化或刚度退化是往往难收敛有些收敛困难可以通过对模型的粘塑性归一化来解决本模型可用粘塑性归一化,因而就允许有效应力超出屈服面根据Duvaut-Lions归一化粘塑性应变定义为自=一非,).
(45213)式中“是粘性参数表征粘塑性系统随时间的松弛,和“是非粘性backbone model的计算塑性应变同理,在粘塑性体系里,粘性刚度退化变量为人:1(l=一(〃一%).(
4.
5.214)vd是非粘性backbone model的刚度退化变量,那么粘塑性模型的应力-应变关系就为a=(1-八)D:(£-非).(452-15)当〃〃T X粘塑性系统的解就趋近于非粘性情况其中代表时间当用较小的粘性系t数时就可以对其进行粘塑性归一化,通常能够改善模型在软化段内的收敛速度,而不会影响计算结果模型数值计算模型采用后退欧拉法进行计算,该方法在ABAQUS里对塑性计算用的很多在平衡迭代中采用了与积分算子相容的雅克比矩阵受拉和受压的损伤状态由两个独立的硬化变量滑和“描述,他们分别代表受拉和受压时的等效塑性应变硬化参数的演化由下式给出(下文将进一步讨论)ept=归力;浮=1】(访髀).加,混凝土的微裂纹和压碎由不断增大的硬化变量来描述这些硬化变量控制着屈服面和弹性刚度退化他们也与产生新裂纹面所要消耗的断裂能有密切的关系屈服函数屈服函数口无故)在有效应力空间内代表一个空间曲面,它决定了失效或损伤的状态屈服函数,至于本粘性无关的塑性损伤模型其屈服函数的具体形式稍后详细介绍F(万,加)<
0.流动法则c)a根据流动法则,塑性流动由塑性势G来确定,形式为do,式中为非负的流动因子,塑性势也是定义在有效应力空间里的其具体形式稍后介绍由于使用的是非相关联流动法则,所以刚度矩阵将会是非对称的b=D;/:(£-£叫G{团敢)<()},—二■%到”)•/.(
4.
5.21)小结总之,塑性损伤本构模型的混凝土弹塑性损伤是在有效应力空间和硬化变亍=D『(£一£叫G—,/()},—二h0洌)•洌.(
4.
5.21)加=叩一量来描述的」一加,式中、和F满足Kuhn-Tucker条件AF=0:A0:F
0.Cauchy是由刚度退化变量力.占)和有效应力按下式a=(1一炯(L
5.2-2)计算得到的从等式可以看出,弹塑性关系与刚度退化是非耦合的式
4.
5.2-2的优点在于他能方便计算机数值计算此处总结的非粘性塑性损伤模型可以很轻易地进行拓展就能考虑粘塑性影响了,只要允许有效应力超出屈服面然后对其归一化就可以了损伤和刚度退化硬化变量包”的演化规律可以很方便的先通过考虑单轴情况在推广到多轴情况来确定(但实际上从单轴到多轴的推广往往并不容易的,译者认为)4=4工,仇力,0f/1,心=心耳’・
①£.0f5d
1.e单轴情况演化:g=(曰君.仇力),452-3首先假定单轴应力-应变关系可以通过下式转化成应力-塑性应变关系:式中下表t c分别代表拉压J和7是拉压时的等效塑性应变率,j=Jd和=Jo弓’山是拉压等型塑性应变,0是温度,fi,(i=L
2....)是其它预定义常变量在单in uniaxialtension and452-4轴拉压情况下有效塑性应变率为in uniaxialcompn^ssion.这一节里面我们约定建是正数,它代表的是单压时的应力值,即正如在图
4.
5.2-1中显示的那样,当从应力-应变曲线的应变软化段卸载时,可以发现卸载的响应是退化了的,也就是说材料的弹性模量看起来变小了损伤了弹性刚度的损伤在拉压试验中表现是大不相同的但在拉压两种情况中,随着塑性变形的增加损伤效果都是越来越明显的混凝土的损伤响应由两个独立的单轴损伤变量,和心,控制,他们是塑性应变、温度和其它行变量的函数4=力工,仇加,0^1,、.
4.
5.2-0心=”工,仇£,0d
1.c图混凝土单轴拉和压应力-应变曲线Oto单轴刚度退化变量是等效塑性应变的非减函数,他们的取值范围在0(无损伤)到1(完全损伤)之间如果及表示材料的初始弹性刚度,那么在单轴拉压下的内=(1一4)瓦仁一工),应力-应变关系分别为屋=(1一/「)石(一一羽)・在单轴加载条件下,裂纹是沿着与应力垂直方向发展的裂纹的成核和扩展就造成了界面有效承载面积的减小,因此就导致了有效应力的增加在单轴压是这种承载面积减小的效果还要稍好一点,因为开始是裂纹基本上是平行于应力方向扩展的,但是当压碎发展到比较厉害时有效承载面积也将显著地减小那么有效单轴内聚力环和外形式如下%=K巴才、=Eo(々一工),(1-df)鼠==Eoin-年’).(1一人)有效单轴内聚力决定了屈服(破坏)面的大小单轴循环加载在单轴循环加载条件下,刚度退化机制比较复杂,它设计到预先存在裂纹的开闭问题和裂纹间的相互作用问题试验观察发现,但循环加载的应力符号变号是反向加载的刚度有所恢复这种刚度恢复也称之为“单边效应”它是混凝土循环加载的一个显著特点特别是当应力有拉变为压是,效应很明显,这时压应力是的受拉形成的裂纹闭合从而是受压刚度得到恢复混凝土塑性损伤模型假定弹性模量按标量减小变量1退化f=(1-
(1)E,QE是材料的初始(无损)模量这个关系式在拉压曲线中都是成立的,刚度减小变量d是应力状态和单轴损伤变量刈和心的函数,在单轴循环条件下ABAQUS假定下式成立(1—(I)=(1—,“心)(1—sdf),0St,s1,(
4.
5.2G)c c式中“和以应力状态的函数,引入他们是为了反应由于反向加载是刚度恢复效应,他们定义为=1—万11;0W wf1,s=1—w1—r*7n;0w1,c ec所if万ri0{/%1=H i=if Tu其中,权系数〃)和〃[这里假定为材料参数,他们分别控制应力反向是的刚度恢复能力举例来说,考虑图
4.
5.2-2荷载有拉变成压的情况假定材料没有初始预损伤,也就是羽=()及(=(),那么此时有1-d=1-sd=1-l-w l—「*力.c tc拉应力(万U>0)时1=1,正如预计的那样=,反之压应力(<711<)时/*=,=(1一
7..如果7=1那么〃=0,材料恢复到受压无损状态E=En,反之,若“「=时,〃=/,材料没有刚度恢复当〃7在0-1之间取值时表示刚度只能部分恢复图
4.
5.2-2受压刚度恢复参数〃7效应的示意图£r器=『*羽,452-7第,单轴循环加载时的等效塑性演化方程也可以进行推广如下:它在单拉或单压就退化为方程
4.
5.2-4的形式多轴情况有必要把硬化变量的演化规律推广到多轴情况下,在Lee andFenves1998的工作基础上,假定有效塑性应变率可由下式计算得到
15.2-8式中皿和分别是塑性应变率张量加的最大和最小主值,.方0rd1£3同是拉压应力权重系数,若有效应力张量三个主值全是正时为1,反之为0Macauley运算定义为《「〉=5Lrl+-r单轴加载情况下方程
4.
5.2-8op[[A退化为单轴定义式
4.
5.2-4和
4.
5.2-7,因为此时单拉时天=Si,单压时nJ rA《21=—若果对塑性应变率张量的主值进行排序如^pl二二二^pc一=^1为昂=4,in,那么多轴普通应力条件下等效塑性盈利率演化可以写成一下矩阵形式,「切].,el=7p/=h^,£p/-,E00h伉邛=0-1-「句弹性刚度退化混凝土塑性损伤模型认为混凝土的弹性刚度退化时各向同性的,且可以用一个单标量写成如下形式D〃=l_dD『;0rf
1.
4.
5.2-9式中的刚度退化标量变量d必须与单轴单调加载时的响应一致,同时还要能够反应在循环加载退化机制带来的复杂性对普通多轴加载情况ABAQUS假定,1—d=1—Sf/l—sdf.0Sq,s1,45210r cc形式上与单轴相同,只是现在通过应力权重系数将它推广到多轴情况了:St=1—w«ra;0Wt1,s=1—w l—ra;0wLc cc显然,很容易验证方程
4.
5.2-10的标量退化式与单轴加载时是一致的很多准脆性材料混凝土的试验表明,当拉应力换到压应力时由于裂纹闭合受压刚度将会恢复但是另一方面,当受压是的微裂纹压碎时,由受压换到受拉时的受拉刚度将不会恢复鉴于此,ABAQUS默认条件下,假定由=0及〃1即只有受压刚度恢复而没有受拉刚度恢复图
4.
5.2-3就是默认条件下的一个应力循环的曲线图图3默认条件下小,=0,〃7=
1.单轴应力循环曲线图拉-压-拉屈服条件本模型的屈服条件基于Lubliner等人1989建议的屈服函数,它综合了Lee andFenves您的修正以考虑拉压不同时强度的不同演化规律用有效应力表达时的屈服函数为。