堆石体的摩擦约束对面板温度应力的影响.docxVIP

堆石体的摩擦约束对面板温度应力的影响 影响板堆堆质量和安全的主要因素是板的裂缝。国内外已建和建筑工地混凝土板堆石坝中板的裂缝的例子很多。虽然小裂缝没有对水库有很大影响，但穿过的裂缝会导致漏水，从而影响水库的正常使用和总体安全。国内外许多科学家对板面裂缝的成因进行了深入的研究。结果表明，温度应力和干缩应力是板面裂缝的主要原因。板面的温度（或干缩）应力来自两个方面，即内外高差和均匀温度降。内外温差的应力可以形成表面裂缝，但不会形成贯穿裂缝。均匀温降会导致整个截面的拉应力，这将导致整个截面的拉力学。由于均匀温度的变化和水库对板的限制，混凝土板被建造在由散粒料组成的垫层中。当板变形较大时，可以在板与垫层之间发生不连续位移，即板可以沿水库移动。因此，水库不是基于连续变形的弹性限制，而是通过板与地板之间的摩擦来限制板的变形。传统的计算方法采用板与水库的相互分离，采用连续元素法或弹性基质梁法，这无疑夸大了堆垛体对板的限制，并夸大了流经板的温度应力的振幅。 本文将不连续分析方法中所使用的接触摩擦的处理方式引入到有限单元法,并模拟分析了混凝土面板和坝体之间的摩擦约束.用某设计中的混凝土面板堆石坝的实际数据,分析不同的摩擦系数对面板温度应力及可能形成的贯穿裂缝的影响,讨论面板产生贯穿裂缝的可能情况. 1 贯入距离dn 本文计算温度场及应力场时,采用常规有限单元法,此处仅介绍接触摩擦的模拟方法.接触模型主要用以处理不连续问题,目前主要有两种方法,嵌接处罚法(Penalty)和拉格郎日法(Lagrange).虽然拉格郎日法物理意义明确,但需要增加未知量,求解不易稳定.因此本文采用嵌接处罚法.在嵌接处罚法中,接触是通过点和面的相互关系来模拟的(图1).P1和P2P3的接触有两种状态:一种是锁定接触状态.当点和面之间的应力不满足摩尔-库仑准则,即剪切力小于摩擦力和凝聚力(当存在凝聚力)时为锁定状态,认为P1和P2P3之间既不能产生贯入也不能有相对的滑动,此时,需要在法向和切向用刚性较大的弹簧锁定.另一种为摩擦接触状态.当点和面之间的剪切力大于摩擦力时,则P1和P2P3之间不能产生贯入但允许产生相对滑动,此时为摩擦接触,需要在法向加弹簧,在切向加摩擦力. 1.1 法向接触当单元i中的P1点与单元j的线P2P3接触时,需要在P1与P2P3之间加法向嵌接.设xi\,yi(i=1,2,3)分别为P1、P2、P3的坐标,ui\,νi(i=1,2,3)分别为三点的位移,则P1与P2P3之间的贯入距离dn可用下式求出: dn=S0l+{Ηi}{δi}+{Gj}{δj}(1)dn=S0l+{Hi}{δi}+{Gj}{δj}(1) 式中:l=√(x2-x3)2+(y2-y3)2l=(x2?x3)2+(y2?y3)2??????????????????√;S0=|1x1y11x2y21x3y3|S0=∣∣∣∣111x1x2x3y1y2y3∣∣∣∣;{Ηi}=1l[Νi(x1,y1)]Τ{y2-y3x3-x2}?{Gj}=1l[Νj(x2,y2)]Τ{y3-y1x1-x3}+1l[Νj(x3,y3)]Τ{y1-y2x2-x1}. 由法向接触引起的附加刚度和附加荷载可根据弹簧的应变能求出,为 {[Κcn]=[[Κcnii][Κcnij][Κcnji][Κcnjj]]?[Κcnii]=p{Ηi}{Ηi}Τ?[Κcnij]=p{Ηi}{Gj}Τ?[Κcnji]=p{Gj}{Ηi}Τ?[Κcnjj]=p{Gj}{Gj}Τ?[Fcni]=-p(S0l){Ηi},[Fcnj]=-p(S0l){Gj}(2) 式中:p为法向弹簧刚度;Ni\,Nj分别为单元i和j的形函数. 1.2 切向接触当图1中P1和P2P3为锁定接触时,需要在切向施加一剪切弹簧以限制P1和P2P3之间的相对滑动.设P0为P1和P2P3在变形之前的接触点,则剪切弹簧沿P2P3方向连接P0和P1.当P0点变形到P1时,沿P2P3方向的剪切位移为: ds=S0l+{Ηi}Τ{δi}+{Gj}Τ{δj}(3) 式中:S0=(x1-x0y1-y0){x3-x2y3-y2}?{Ηi}=1l[Τi(x1,y1)]Τ{x3-x2y3-y2}?{Gj}=1l[Τj(x0,y0)]Τ{x2-x3y2-y3}. 由剪切弹簧引起的附加刚度和附加荷载同样可以由式(2)求出. 1.3 摩擦力荷载当接触点之间产生相对滑动时,会在相对滑动的点和面之间存在摩擦力.摩擦力的大小与法向接触力pdn及摩擦系数tan(?)成正比,?为摩擦角.即:f=p·dn·sgn·tan(?),其中sgn为摩擦力的符号,取决于→Ρ0Ρ1与→Ρ2Ρ3之间的夹角.当夹角小于90°时为正,大于90°时为负.方向为→Ρ2Ρ3的方向.根据产生相对滑动时摩擦力所做的功,可以求出摩擦力在