中国学术杂志网

层状土中支护结构稳定方法

 论文栏目:建筑勘察论文     更新时间:2012-9-3 8:44:18   

土钉支护结构整体稳定性分析是判定土钉支护结构是否正常工作的一个重要依据,在土钉支护结构设计计算和研究中具有极其重要的地位。整体稳定性分析即可以确定土钉设计所选各个参数的是否合理,是否可行,同时也可以为其安全性和适用性提供保证。这样,如何确定土钉支护结构的最小安全系数和其对应的最危险滑动面位置,就成为土钉支护设计工作中的首要任务。对于土钉支护结构整体稳定性分析,学者们基于大量的试验研究,提出了以极限平衡方法为基础的分析方法。采用极限平衡理论研究土钉支护结构的整体稳定性,通常首先假定一个潜在的滑动面,对于假设滑动面形状的不同,稳定性的分析方法可以分为楔体稳定分析法和圆弧稳定分析法。其中,最危险滑动面的形状和位置的确定成为土钉支护设计的关键。文献[1]假定土钉墙潜在破裂面为圆弧形考虑土钉抗拔力与土条自身的阻滑力的耦合作用利用混沌优化方法搜索最危险滑动面。文献[2]则利用复合形法寻求滑面几何控制参数,采用分步黄金分割法寻找最危险滑裂面圆弧,由此计算出与之对应的最小安全系数值。文献[3]采用在分析土钉墙整体稳定性时,潜在滑裂面为对数螺旋曲线,并借助复合形法优化理论确定最危险滑裂面的位置。但对于目标函数是多峰的、或者搜索空间很不规则的优化问题,在搜索最优解时很容易陷于局部最优解对于深基坑开挖边坡土层分布均匀的土钉墙,其最危险滑裂面搜索可归结为一个目标函数为凸函数的优化问题,用上述方法可以得到全局最优解。但如果开挖边坡土层分布不均匀时,得到的目标函数往往是具有多峰的非凸函数,那么采用上述方法得到的就很可能是局部最优解。而遗传进化算法的优点恰好擅长搜索全局最优解,即使在所定义的适应函数(或目标函数)是不连续的、非规则的情况下,它也能以很大的概率找到整体最优解[4]。文献[5]成功地将遗传算法应用到以圆弧滑动简单条分法为基础的土钉支护结构整体稳定性分析中。但遗传算法同样存在过早的收敛到局部最优解和在最优值附近收敛速度慢等缺点,为此,文献[6]将单纯形法引入经验遗传,提出了一种针对连续设计变量的经验遗传-单纯形算法。文献[7]将人工神经网络与遗传算法相结合,利用神经网络算法建立起优化设计变量和安全系数最小值之间的非线性映射关系,从而提高了运算速度和搜索效率。本文根据极限平衡原理和遗传算法理论[8],提出了采用动态自适应技术和非标准的遗传操作算子改进遗传算法的新算法,并将其引入到土钉支护结构整体稳定性分析中去,建立了一种能同时确定土钉支护最危险滑动面和最小安全系数的动态自适应遗传算法模型,为土钉支护结构的优化设计提供新方法。

1层状土中土钉支护整体稳定性分析模型简单条分法是土钉支护整体稳定性分析的最常用的一种方法。在分析过程中,假定滑动面为圆弧,土体为刚塑形材料,不考虑土条之间的相互作用力和力矩。本文研究了层状非均质土中土钉支护的整体稳定性的计算方法。

1.1基本假设1)如图1[9],最危险滑动面为坡脚圆弧,通过土钉墙地面角点A。2)最危险滑动面圆心位于土钉墙上方,即位于BC上方。

1.2模型的建立1)土钉稳定性分析的坐标系。采用圆弧滑动简单条分法进行土钉墙整体稳定性验算,首先建立直角坐标系,以基坑底坡脚A点位坐标原点,基坑深度方向为y轴,与基坑深度垂直的方向为x轴。2)滑动面圆弧半径。如图1所示,圆弧滑动面的圆心坐标为O(xo,yo),圆弧的半径为R,基坑深度为H,那么滑动面的方程为[10-11]:对于非均质土,特别是各层土的性质相差较大时,不能简单的进行加权平均后代入上面的公式计算土钉其中:γk为整体滑动分项系数,可取1.3;γ0为基坑侧壁重要性系数;cik为第i分条滑动面处土体粘聚力;φik为第i分条滑动面处土体内摩擦角;Sh为土钉的水平间距。

1.3滑动面圆心范围在滑动面搜索的过程中,圆心的范围取值应尽可能的大,从而使得最危险滑动面的圆心落在范围内,但取值范围过大会造成搜索效率降低和计算时间增长,同时还有可能产生一些不合理的滑动面。文献[3]中,圆心的横坐标向坡脚A左右分别取4倍坡高的范围,纵坐标自坡肩向上4倍坡高,坡肩向下到坡底的范围。而文献[5]中则认为,圆心取定的范围应该是x∈[0,4H],y∈[H,5H]。根据转动理论滑动体中每一点的速度将垂直圆心和坡脚的连线,当圆心位于A点的右侧,如图5所示的O2,坡脚A点速度vA2是倾斜向上的,说明整个滑动体将发生斜向上的滑动,这与土钉支护的边坡失稳时的实际情况不符。而对于圆心位于O3,滑动面与坡面交于D的情况同样与实际情况相违背,这是由于vD2的方向是斜后方[12]。图中O1才是滑动面圆心正确的位置,无论vA1还是vB1的方向都是斜前方,因此,它的取值范围应该为:

2基于DAGA法土钉结构最危险滑

动面搜索土钉支护结构的整体稳定性分析的关键在于寻求最危险滑动面以及与之对应的最小安全系数。而寻求最危险滑动面的过程实际上就是一个动态搜索的过程,它在数学上是一个动态优化的过程。因此,提出的动态自适应遗传算法搜索滑动面方法,能够克服收敛到局部最优解和在最优值附近收敛速度慢的缺点,提高运算速度和搜索效率,从而得到全局优化解。

2.1目标函数对于最危险滑动面可有圆点O和坡脚A两点坐标确定,因此O点的纵、横坐标则成为染色体的基因表现型[13]。根据上面推导的非均质土的安全系数表达式,在给定圆点O的搜索范围的情况下,安全系数的目标函数可表示为:.2编码方式确定和初始化种群遗传算法染色体编码通常可采用二进制编码、格雷编码和浮点数编码等。对于土钉整体稳定性分析问题,为了提高遗传算法的精度,同时改善计算的复杂性和运算效率,采用浮点向量法编码更加方便。浮点数编码是指将个体的每个基因值用某一范围内的一个浮点数来表示,个体编码长度等于其决策变量的位数[14-15]。在土钉整体稳定性分析模型中,滑动面所对应的圆心坐标(xo,yo)设计变量,采用一个浮点向量V=[x0y0]表示。当进化代数计数器t←0时,随机产生n个染色体构建成初始种群,其中第i个染色体为Vi=[(x0)i(y0)i],其中(x0)i和(y0)i为Vi的基因[16]。
2.3适应度函数的设计适应度是用来衡量种群中各个个体在优化计算中能达到、接近于或有助于得到最优解的好坏程度。而度量个体的适应度的函数称之为适应度函数。适应度函数是区分种群中个体优良程度的标准,同时也是算法演化的动力。由于滑动面圆心的位置存在一定的范围,因此土钉整体稳定性分析是一个有约束的优化问题,在构建适应度函数时,采用将约束以动态方式合并到其中,即形成一个具有变化的惩罚项的适应度函数。罚函数的作用是将在解空间中无对应可行解的个体计算适应度时,处以一个惩罚,从而降低其个体的适应度,使该个体被遗传到下一代群体中的概率减小,便于指导搜索,从而提高了找到全局最优解的概率。由于土钉结构整体稳定性分析中安全系数计算是一个有约束的极小值的问题,所以必须将目标函数进行转换才能建立适应度函数。因此,转换后的适应度函数为:通过适应度函数将土钉整体稳定性分析的约束最优化问题转化为无约束最优化问题。式(18)表明,当点X不满足约束条件时,说明它在可行域之外,需对目标函数的值加以惩罚;或者当点X位于约束边界附近时,λΦ(X)将趋于无穷大,迫使迭代点只能在可行域内移动。

2.4竞争、选择操作根据Vi所对应的适应度的大小采用最佳保留策略来进行优胜劣汰操作,将迄今为止最佳个体直接传递到下一代或至少等同于前一代,确保最优个体不会被交叉、变异等遗传运算所破坏,从而保证遗传算法的收敛性。

2.5交叉、变异操作交叉操作是遗传算法中最重要的遗传操作,是区别于其他进化运算的重要特征。算术交叉运算是将种群内两个相配对的染色体按照线性的方式相互交换它们之间的部分基因,从而产生新的染色体。变异操作是以一个很小的概率随机地改变染色体上的某一或某些基因而产生新的个体,即以变异概率随机改变某一染色体Vi=[(x0)i(y0)i]中一个或者几个基因,从而产生新的染色体。为了降低进化过程中出现局部最优解的可能性,分别对交叉算子和变异算子进行改进。其方法是,对于适应度较小的个体,将交叉算子pc和变异算子pm取较大值,这样可以提高搜索速度;而对于适应度较大的个体,交叉算子pc和变异算子pm则取较小值,这样可以降低收敛于局部最优解的可能性[17-18]。其中:pc1、pc2为交叉概率的上、下限;pm1、pm2为变异概率的上、下限;f′是交叉的2个个体中适应度较大值;favg、fmax分别表示当前种群的平均适应度和个体的最大适应度[19]。

2.6迭代终止循环执行2.3—2.5步的操作,直到目标函数f(X,λ)达到满意值或达到预先设定的代数时,终止计算。这时f(X,λ)所对应的设计变量组合为最危险滑动面。

3工程应用

一个尺寸为110m×78m的高层建筑基坑,深度为7.8m。其土层分布最上层为1.9m填土,第2层为0.8m的细砂层,第3层为2.7m的粉质粘土层,第4层为1.8m的粉土,以下为强分化砂岩,其厚度大于12m。各土层相关物理力学参数见表1。基坑采用土钉支护结构,土钉墙面与水平面的夹角为80°。土钉层数为6层,其水平间距和垂直间距均为1.2m,土钉钻孔直径为φ110mm,土钉入射角为10°,钢筋采用Ⅱ级φ25螺纹钢。土钉墙面板采用厚度为100mm的喷射混凝土。混凝土设计强度等级为C20,钢筋网采用φ8@250mm×250mm。图6为该基坑的土层分布情况以及土钉支护形式。采用本文方法求取的最小安全系数为1.2628所对应的圆心坐标为(-6.8,16.8)。在搜索最危险滑动面时,一些文献有时为了计算的简便将层状土层转化为均质土考虑,即将土层参数均值化。均值化方法分为2种,即直接平均和加权平均,前者是将所有土层简单数学平均,后者则根据土层厚度的贡献平均[20]。本文也采用加权平均法对该实例进行了计算。将非均质的层状土转化为均质土后,即对土层重度、粘聚力、摩擦角、土的摩阻力分别进行加权平均后得到的数据进行计算。结果表明,当圆心位于(-5.8,13)时,安全系数达到最小值1.3221,与分层计算时所得最小安全系数为1.2628存在一定的差别。从而验证了土性参数的取值对工程计算的影响很大,这说明将层状土转化为均质土的计算结果是欠准确的。

4结论

1)运用简单条分法对层状非均质土中土钉支护结构进行了稳定性分析,得出了最危险滑动面的搜索模型和最小安全系数的计算公式。通过实例对比分析了非均质层状土和等效均质土计算方法,分析结果表明采用等效均质土的计算明显提高了土钉支护结构的最小安全系数,对于其整体稳定性是不利。2)为了解决遗传算法中过早的收敛到局部最优解和在最优值附近收敛速度慢等缺点,采用动态自适应技术改进遗传算子,不仅大大提高空间搜索能力,而且提高了种群中表现优良个体的交义率和变异率,从而提高获得优化的全局最优解可能性。3)将DAGA法应用到土钉支护结构的整体稳定性分析领域,建立了一种能有效搜索土钉支护结构的最危险滑动面极其对应的安全系数的全局优化算法。算例分析表明,与传统的遗传算法相比,本文提出的DAGA法在运算速度和搜索效率更具优势。


学术网收录7500余种,种类遍及
时政、文学、生活、娱乐、教育、学术等
诸多门类等进行了详细的介绍。


建筑勘察论文
@2008-2012 学术网
主机备案:200812150017
值班电话
0825-6697555
0825-6698000

夜间值班
400-888-7501

投诉中心
13378216660
咨询电话
唐老师:13982502101
涂老师:18782589406
文老师:15882538696
孙老师:15982560046
何老师:15828985996
江老师:15228695391
易老师:15228695316
其它老师...
咨询QQ
89937509
89937310
89903980
89937302
89937305
89937307
89937308
业务
综合介绍
论文投稿
支付方式
常见问题
会员评价
关于我们
网站简介
版权声明
友情链接
人员招聘
联系我们