土的抗剪强度指标的随机场统计方法

土的抗剪强度指标的随机场统计方法 //.paper.edu 土的抗剪强度指标的随机场统计方法1 闫澍旺，朱红霞，刘润 (天津大学建工学院，天津 300072) 摘 要:在利用随机场理论进行土工可靠度计算时，关键问题之一在于如何考虑土的抗剪强 度指标的自相关性及互相关性。本文以抗剪强度指标的传统统计方法为基础，结合随机场理 论，提出抗剪强度指标的随机场统计方法，在抗剪强度指标统计时既考虑单个指标的自相关 性，又考虑指标之间的互相关性;并通过一工程算例，进一步说明此统计方法可使土工可靠 度计算更加合理、简便。 关键词:抗剪强度指标，随机场，可靠度，相关性 1(前言 随着概率统计理论在岩土工程中的应用，如何选择合理的土性指标统计方法 成为地基可 靠度分析面临的主要问题之一。这是因为，与人工 材料 关于××同志的政审材料调查表环保先进个人材料国家普通话测试材料农民专业合作社注销四查四问剖析材料 不同，土的性质既有天然的变异性又 有一定的相关性，因此岩土工程中的可靠度分析必须要考虑土性指标的相关性。这种相关性 包括两方面，即同一土性指标的自相关性及不同土性指标间的互相关特性。 在关于土的自相关性研究中，Vanmarcke(19><>976)首先将随机场理论引入岩土工程的可靠 度分析之中，建立了土性剖面的随机场模型，即将土性剖面看作在空间上分布的随机场。这 是符合土的性状的，这样处理可使问题大为简化。 对于土的抗剪强度指标的相关性问题，不同的学者有不同的考虑方法。有些学者建议在 指标统计时考虑自相关性;而对互相关性，则在计算过程中加以考虑[1][2]。抗剪强度指标统 计时可采用传统法或高大钊提出的综合回归方法[3]。传统法即分组求出内聚力 和内摩擦角c ??tg 后再计算均值和子样方差，然后根据随机场理论求得均值方差;综合回归方法是在各试 验样本独立条件下，将所有试验数据作为一个样本，一次求得两个指标的均值与均值方差。 无论是采用传统法还是综合回归方法，在实际应用中都有一定的困难。因为传统法在由子样 方差求均值方差时，必然要研究 c、 ??tg 的自相关函数及互相关函数，这无疑增加了许多工 作量;若采用综合回归方法，只有各抗剪强度试验子样相互独立时，回归得到的才是内聚力 和内摩擦角的均值方差;且研究 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 明，此法适用于求τ 、 P 的回归系数的估计值 ??tg 、 的 方差情况，不能得到反映土的 c ??tg 、 变异性的方差c [4]。 1 本课题得到国家自然科学基金(批准号:5037<>9034)及高等学校博士学科点专项科研基金(批准号: 20020056012)资助。 - 1 - //.paper.edu 1<>98<>9 年黄传志、孙万禾根据 ??tg 、 c相关的概念考虑相关因素，提出了简化相关法。此 法与现今港口工程地基 规范 编程规范下载gsp规范下载钢格栅规范下载警徽规范下载建设厅规范下载 采用的τ 平均法是一致的[5]。这种方法在假定 ??tg 、 独立的情 况下，将两个土性指标的互相关性的影响通过回归自动反映在子样方差中。因此，在进行土 工可靠度计算时，就不用考虑指标间的互相关性。但是，这种方法没有考虑抗剪强度指标的 自相关性。 c 本文首先以抗剪强度指标的传统统计方法为基础，引入随机场理论，提出抗剪强度指标 的随机场统计方法，使之不必分别研究 、c ??tg 的自相关函数及互相关函数，而是通过回归 直接得到均值方差，使计算大为简化;同时，将随机场理论引入简化相关法，在抗剪强度指 标统计时既考虑单个指标的自相关性，又考虑指标之间的互相关性，从而使土工可靠度计算 更加合理、简便。 2(抗剪强度指标的随机场统计方法 2.1 随机场理论的引入 土是自然历史的产物，地质成因条件的影响使土的性质与位置有关。为了使模拟土的性 质的概率模型更符合实际，有必要引入随机场理论，把土性指标看成是随空间位置而变化的 一族随机变量，即在空间上分布的随机场。适用这种需要而在近年由 Vanmarcke(1<>976)发展 的土性剖面的随机场模型，为解决这类问题的概率模拟，提供了一个有效的工具。 工程实践中应用的随机场可看作齐次随机场，它的均值 是一个与坐标无关的常 数 )]([ PYE ?? ， 为随机函数。)(PY P 表示空间的一个点，相关函数 是一个仅依赖于向量 )])()()([(),( <>9;<>9; ???? ????= PYPYEPPR <>9;PP 的函数，而与 P 的位置无关，即有 )(),( <>9;<>9; PPRPPR = 。为简单起见， 只讨论一维各向同性的齐次随机场，这是符合土的实际情况的。这时， 可 用 来表 示，而且总是假定 。 )(PY )(zY 0)]([ =zYE 一维齐次随机场在 上的随机积分为 ],[ hzz + ? += hzzh dzzYhzY )(1)( 该随机积分的均值: ? ?+ + === hzz hzzh dzzYEhdzzYhEzYE 0)]([1])(1[)]([ (1) 该随机积分的方差: - 2 - //.paper.edu ])()1(2[ )()1(2)()1(2 )(1 )]()([1 ])(1[])(1[)]([ 0 2 0 2 0 0 0 21212 21212 2 ? ?? ? ? ? ? ? ? ??= ??=??= ??= = == + + + + h hh h h hz z hz z hz z hz z h d hh d hh dR hh zdzzzR h dzdzzYzYE h dzzY h EdzzY h VarzYVar ττρτσ ττρσττττ (2) 令 [ ] ? ??==Γ hh dhhzYVarh 022 )()1(2)()( ττρτσ (3) 式(2)建立了局部空间平均的方差与子样的点方差 之间的关系。2σ )(τρ 为相关函数， 是[0,1]之间的一个系数，表明空间平均方差小于点方差，因此可称为方差 折减函数。 方差减小的程度取决于 和相关函数 )(2 hΓ h )(τρ ，当 足够大时， ，h δ=Γ?? )(2 hh δ 称为相关距离， 在该距离内，土性强烈相关，大于该相关距离可认为基本不相关。在实际计 算中，一般采用 如下简化 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 [6]: ???? ???? ?? ? ? =Γ )( )(1 )(2 u u u h h h h δδ δ (4) 2.2 抗剪强度指标按随机场统计的传统法 将随机场理论引入抗剪强度指标的传统统计方法，在指标统计时考虑自相关性的影响。 首先取 n(i=1~n)组抗剪强度数据，每组数据对应 k(j=1~k)级荷载。即对某些特定的压力 Pj，j=1，2，„„，k，给出相应的若干组抗剪强度τ 的试验值: ijτ ，i=1，2，„„，n。按 照摩尔强度包线: θτ Pc += ，这里为方便起见，记 ??θ tg= 。 2.2.1 抗剪强度指标的均值统计 由式(1)可知，随机场的局部平均 的均值和原来的随机场 的均值一样。所 以， 抗剪强度指标的均值按传统法、综合回归方法或 )(zYh )(zY τ 平均法统计均可，研究表明，这三种统计 方法的均值相同[4]。下面给出 ??tg 、 均值的计算公式: c ? = = n i itg tgn 1 1 ???? ?? ???? tg —— ??tg 的平均值; (5) - 3 - //.paper.edu ? = = n i ic cn 1 1?? c?? —— c的平均值; (6) 式中: ? ? = = ?? ?? = k j Pj k j ijPj i P P tg 1 2 1 )( )( ?? τ?? ?? itg?? ——每组试验的 ??tg 的回归值; iPi tgc i ??????τ ??= ——每组试验的 c的回归值; ic ? = = k j jP Pk 1 1?? P?? ——每组试验各级垂直压力 Pj(j=1~k)的 均值; ? = = k j ijki 1 1 τ??τ iτ?? ——每一组试验(i=1~n)各级压力 (j=1~k)下抗 剪强度 ijτ 的平均值。 2.2.2 抗剪强度指标的方差统计 设 ??tg 、 是平稳随机过程，且是联合平稳的，则在c ],[ hzz + 上相应 的随机积分分别为: ? += hzzh dzzhz )(1)( ττ ? += hzzh dzzchzc )(1)( ? += hzzh dzzhz )(1)( θθ 对任一个确定的 P 值，有τ 也是平稳过程，且: )()()( zPzcz hhh θτ += 取方差运算，得: )()]()([)()( 222222222 hPhhPhh cccccc θθθθθθττ σσσσσ Γ+Γ+Γ+Γ=Γ (7) 式中: ? ??=Γ h dhhh 02 )()1(2)( ααρα ττ ? ??=Γ h cc dhhh 02 )()1(2)( ααρα - 4 - //.paper.edu ? ??=Γ h dhhh 02 )()1(2)( ααρα θθ ? ??=Γ h cc dhhh 02 )()1(2)( ααρα θθ ? ??=Γ h cc dhhh 02 )()1(2)( ααρα θθ 2τσ 、 、 分别为2cσ 2θσ τ 、 、c θ 的“点方差”; cc θθ σσ = 为 与c θ 的“点协方差”; τρ 、 cρ 、 θρ 分别为τ 、 、c θ 的相关函数; θρ c 、 cθρ 分别为 和c θ 、θ 和 的互相关函数。 c 通常按照传统法，对任一个 i，由 ),( ijjP τ ， j=1，2，„„，k，可 确定出一个 和ic iθ 。 因此从理论上说，由 和ic iθ ，i=1，2，„„，n 可确定 、 、2cσ 2θσ θσ c 、 cθσ 及 )(αρ c 、 )(αρθ 、 )(αρ θc 、 )(αρθc ，进一步可得到方差函数 、 及协方差函数 、 。这样做到了在抗剪强度指标统计时考虑自相关性，但要想在可靠度计算中考虑 互相关性，就不可避免地要去研究互相关函数。这无疑会给计算增加许多工作量。 Γ )(2 hcc θθσ Γ )(22 hcc Γσ )(22 hθθσ )(2 hcc θθσ Γ 一般来说，实际问题中的 z、h 总是确定的，因此 、)(zch )(zhθ 就退化为随机变量。因 此，从实用的角度考虑，可不必分别获得 、 、2cσ 2θσ θσ c 、 cθσ 及 )(αρ c 、 )(αρθ 、 )(αρ θc 、 )(αρθc ，而只需根据式(7)，直接回归得出 、 及 就可 用于可靠度的计算。 )(22 hcc Γσ )(22 hθθσ Γ )()( 22 hh cccc θθθθ σσ Γ+Γ 实际计算中，可先按照考虑自相关性的方法求出 k 个 ，然后利用数据 ，j=1，2，„„，k，按式(7)进行二项式的回归计算，这用计算机程序很 容 易实现。 )(22 h jj ττσ Γ ))(,( 22 hPj Γ jj ττσ 值得注意的是，在统计 时，所用到的2 jτσ ijτ 应该是修正过的值，即用每组回归出的 、ic iθ (i=1~n)，把每组数据 ijτ (j=1~k)都修正到 ijiij Pc θτ += 直线 上去。认为每组数据都严格服 从库仑定律，此时，式(7)才成立。 下面讨论如何求 。根据随机场理论， 可按式(4)求得。那么理论上首先 要根据试验所得的特定的压力 Pj 下对应的抗剪强度 )(2 h jτΓ )(2 hjτΓ ijτ ，分别求出对应的相关距离 juδ 。但 是考虑到这 n(i=1~n)组试验样本应来自于同一土层母体，因此用不同压力 Pj(j=1，2，„„， - 5 - //.paper.edu k)下对应的抗剪强度 ijτ 求得的相关距离值应该是相等的，即只有一个 相关距离 uδ 。进而利 用式(4)求得的 也应该是相等的，即 。实际计算中，相关距离)(2 h jτΓ )()( 22 hhj ττ ΓΓ , uδ 可近 似地取为平均值，即 ? = = k j juu k 1 1 δδ 。 这种土性指标的随机场统计方法考虑了指标的自相关性，而互相关性可根据回归得到的 进一步在可靠度计算中考虑。 )()( 22 hh cccc θθθθ σσ Γ+Γ 2.3 抗剪强度指标按随机场统计的简化相关法 将随机场理论引入简化相关法(τ 平均法)，在考虑抗剪强度指标互相关性的同时，考 虑自相关性的影响。 抗剪强度指标均值的统计与 2.2.1 节的统计方法相同，下面讨论方差的统计。在用 JC 法计算可靠度时，考虑将 、)(zch )(zhθ 作为独立的随机变量，因此式(7)可化为: )()()( 2222222 hPhh cc θθττ σσσ Γ+Γ=Γ (8) 比较式(7)及式(8)，可看出两个土性指标的互相关性的影响可通过式(8) 的回归自 动反映在均值方差中。 先按照 2.2.2 节考虑自相关性的方法求出 k 个 ，然后利用数据 ， j=1，2，„„，k，按式(8)用最小二乘法进行回归计算，得到 和 如下: )(22 h jj ττσ Γ ))(,( 222 hPj jj ττσ Γ )(22 hcc Γσ )(22 hθθσ Γ ?? ? == = Γ??Γ??=Γ k j k j k j jj hPhPkh jjjj 1 22 1 1 222222 )]()([1)( ττττθθ σσσ )()(1)(1)( 22 1 1 22222 hP k h k h k j k j jcc jj θθττ σσσ Γ??Γ=Γ ? ? = = 式中， ? ? = = ??=?? k j k j jj PPk 1 2 1 24 )( 这样得到的 和 可作为相互独立的随机变量直接应用于可靠度的计算， 而不必在可靠度计算中再考虑土性指标间互相关性的影响。 )(22 hcc Γσ )(22 hθθσ Γ 3(可靠度算例 以天津新港北港池地区的土质统计资料为基础，用随机场方法统计其抗剪强 度指标，并 - 6 - //.paper.edu 对其上的一长方形基础的地基承载力作可靠度分析，以进一步说明抗剪强度 指标的随机场统 计方法如何应用于可靠度计算。为便于比较，本算例也用其他方法对土性指标进行统计，并 根据不同的相关性考虑方式分别对地基承载力作可靠度分析。 3.1 计算条件 长方形基础长度 L =10m，宽度 =4m，埋置深度 =1m。基础上作用的有效垂直荷载 =120KPa，变异系数 0.2;水平荷载 =15KPa，变异系数 0.125。持力层最大深度 Zmax=6.0m。 B D vp hp 3.2 土层强度统计特性 长方形基础以下 0,6m为粘土层，此土层的容重变异性很小，为简化起见暂不考虑，取 其均值为 17.74 kN/m3。采用固结快剪试验结果，试样间距 0.5m，用相关函数法求得抗剪强 度τ 的相关距离为 τδ =0.152m。 滑动面最大切割深度约为 4m，在此深度范围内对土性指标进行空间平均特 性统计，各 种统计方法得出的抗剪强度指标分别如表 1、表 2 所示: 表 1 传统法 未考虑自相关 考虑自相关 土性指标 平均值 “点” 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 差 “点”变异系数 均值标准差 均值变异系数 c (kPa) 18.65 6.6<>93 0.35<>9 3.864 0.207 υ(。) 14.728 2.41 0.164 1.37><92 0.0<>94 c 、υ协方差 -0.188 -0.063 表 2 简化相关法 原简化相关法 按随机场统计的简化相关法 土性指标 平均值 “点”标准差 “点”变异系数 均值标准差 均值变异系数 c (kPa) 18.65 6.333 0.34 3.656 0.1<>96 υ(。) 14.755 1.33 0.0<>9 0.76<>9 0.052 从表 1 和表 2 的抗剪强度指标统计结果可看到: (1)各种统计方法得到的抗剪强度指标的均值近似相等。 (2)若不考虑指标的自相关性，则原简化相关法的“点”标准差小于传统法统计的“点” 标准差，这是因为简化相关法考虑了 、c ??的相关性。 (3)抗剪强度指标按随机场统计的简化相关法得到的标准差小于原简化相关法统计的 - 7 - //.paper.edu “点”标准差，这是因为考虑了指标的自相关性，统计得到的是均值标准差。 (4)抗剪强度指标按随机场统计的简化相关法小于传统法统计的均值标准差，这也是 因为考虑了指标的互相关性。 3.3 地基承载力的可靠度计算 本例采用 Hansen 公式计算极限承载力，考虑的随机变量有水平荷载 ，竖直荷载 ， 内聚力 和内摩擦角 hp vp c ?? 。各随机变量都符合正态分布。极限状态方程为: vvh pctgciSdNctgcDiBSNPPcgZ ????????++== ????γγ?? γγγ )(2 1),,,( 式中，B 为基础宽度，D为基础埋深，γ 为地基土的容重，这三个参数都作为常量考虑; 为承载力系数，是内摩擦角γN qN ??的函数; 、 为与基础形状有关的、 系数; 、 为与基础底面作用合力的倾斜率有关的倾斜系数; 为与基础埋深有关的深度系数。 γS qS γi qi qd 采用国际“结构安全度联合委员会”(JCSS)推荐的 JC 法计算可靠度指标 β 。JC 法考 虑了极限状态方程中各基本变量的实际概率分布，将功能函数在验算点处用泰勒级数展开并 使之线性化，最后求解可靠度指标，是一种精度较高、较实用的近似概率分析方法。 JC法假定随机变量之间是相互独立的，本文提出的抗剪强度指标的随机场统计方法在土 性指标统计时就考虑了指标之间的相关性，所以可直接应用JC法计算可靠度指标。若抗剪强 度指标按传统法统计，若要在可靠度计算中考虑指标间的互相关性，可采用改进的JC法[7]。 根据上节中不同方法统计的抗剪强度指标，按照土性指标相关性的不同考虑方式，分六 种情形对本例中的地基承载力进行可靠度分析，结果见表 3。 表 3 计算结果汇总表 计算结果 序 号 抗剪强度指标 统计方法 考虑 自相关 考虑 互相关 β fP (%) K 1 1.<>9<>95 2.302% 2 ? 3.266 0.055% 3 ? 2.172 1.4<>93% 4 传统法 ? ? 3.354 0.040% 5 简化相关法 ? 2.173 1.48<>9% 6 按随机场统计的简化相关法 ? ? 3.537 0.020% 2.37 注:表中， β 为可靠度指标; 为失效概率;fP K 为按定值法求得的 安全系数。 - 8 - //.paper.edu 3.4 可靠度计算结果分析 根据表 3 中可靠度分析结果，可得: (1)当采用传统的抗剪强度指标统计方法时，若既不考虑指标的自相关性，又不考虑 指标的互相关性，则所得的可靠度指标偏小，失效概率偏大，达到 2.302%。与传统的安全 系数 2.37 相比，相差甚远，这在工程设计中是难以为人们接受的。 (2)当采用传统法时，若分别考虑指标的自相关性和互相关性，则所得的可靠度指标 都有一定的提高。相比较下，考虑抗剪强度指标的自相关性对可靠度指标的影响较之考虑互 相关性要大一些。 (3)采用按随机场统计的传统法，在土性指标统计时考虑自相关性，在可靠度计算中 考虑互相关性，得到的失效概率为 0.040%，与传统的安全系数较为匹配。 (4)原简化相关法只考虑了抗剪强度指标的互相关性，使计算所得的可靠度指标有一 定提高，但仍偏小;按随机场统计的简化相关法既考虑了自相关性又考虑了互相关性，所得 的失效概率最小，为 0.020%，与传统的安全系数相比，是较为合理的。 (5)上表中第 4 种情形和第 6 种情形(即本文提出的两种抗剪强度指标的随机场统计 方法)计算所得的可靠度指标很接近，说明虽然这两种方法虽然对相关性考虑方式不同，但 所得结果较为一致，在工程应用中都是可行的。 4(结语 (1)在应用随机场理论进行可靠度计算时，土性指标的自相关性和指标间的互相关性 不容忽视。 (2)本文提出的抗剪强度指标按随机场统计的传统法和按随机场统计的简化相关法， 分别从不同角度对相关性加以考虑，较之传统的统计方法，更为全面和合理。 (3)两种统计方法所得的可靠度指标接近，且与传统的安全系数相匹配，证明在土工 可靠度计算时，两种方法都可作为合理的抗剪强度指标统计方法。 参考文献 [1] 包承纲,黄卫峰,张庆华.随机场理论在重力式码头地基承载力计算中的应用.见:地基工程可靠度分析 方法研究[M].武汉:武汉测绘科技大学出版社,1<>9<>96:110-117. [2] 熊启东,高大钊.上海地区地基承载力的可靠度分析[J].岩土力学,1<>9<>97,18(1):73-77. [3] 高大钊.土的抗剪强度指标的统计与应用.见:地基工程可靠度分析方法研究[M].武汉:武汉测绘科技大 学出版社,1<>9<>96:<>98-104. [4] 孙万禾,黄传志,叶国良等.土的抗剪强度指标统计方法的分析. [J].港口工程,1<>9<>96,(3):5-16. [5] 港口工程地基规范(JTJ250-<>98),1<>9<>98,4. [6] Vanmarcke Erik H. Probabilistic Modeling of Soil Profiles[J]. Journal of the Geotechnical Engineering Division, ASCE,1<>977,103(GT11):1227-1246. - <>9 - //.paper.edu [7] 张庆华,包承纲.改进 JC 法及其在土工可靠度计算中的应用.见:地基工程可靠度分析方法研究[M].武 汉:武汉测绘科技大学出版社,1<>9<>96:105-10<>9. [8] 高大钊.土力学可靠性原理[M].北京:中国建筑工业出版社,1<>98<>9. [<>9] 闫澍旺,贾小黎,郭怀志.土性剖面随机场模型的平稳性和各态历经性验证[J].岩土工程学 报,1<>9<>95,17(3):1-<>9. [10] 刘润 ,闫澍旺 ,吕涛 .信息管理技术在海洋工程地基可靠度分析中的应用 [J]. 海洋技 术,2004,23(4):<>90-<>95. [11]刘润,闫澍旺.海工结构物浅基础稳定性的可靠度分析[J].海洋工程,2001,1<>9(2):23-28. The stochastic statistical methods for characteristics of shear strength YAN ShuWang ZHU HongXia LIU Run (School of Constructional Engineering, Tianjin University, Tianjin 300072, China) Abstract When applying random field theory to reliability analysis, one of the most important problems is how to consider the auto correlation and cross correlation of the shear strength of soils. In this paper, the stochastic statistical methods for characteristics of shear strength are presented based on the traditional statistical methods. Both the auto correlation of the individual index and the cross correlation between the indexes can be considered when using the stochastic statistical methods. It is proved that the methods can make the reliability analysis more reasonable and convenient by calculating an example of engineering. Keywords: shear strength, random field, reliability, correlation 作者简介: 闫澍旺:男，1<>950 年生，教授，主要从事岩土工程专业的研究和教学工 作。 - 10 -