高压条件下石英砂中的静力触探试验及大变形模拟.pdf
《高压条件下石英砂中的静力触探试验及大变形模拟.pdf》由会员分享,可在线阅读,更多相关《高压条件下石英砂中的静力触探试验及大变形模拟.pdf(7页珍藏版)》请在咨信网上搜索。
1、第 45 卷第 4 期2023 年 8 月Vol.45 No.4Aug.2023土 木 与 环 境 工 程 学 报(中 英 文)Journal of Civil and Environmental Engineering高压条件下石英砂中的静力触探试验及大变形模拟孔壮壮,王栋,张民生,裴会敏(中国海洋大学 山东省海洋环境地质工程重点实验室,山东 青岛 266100)摘要:在深部土层中进行原位试验时,以往针对中低应力水平的数据解译公式不再适用。为探究砂土在高应力条件下的力学特性,利用自行研制的高压标定罐,改变围压水平,进行不同密实度石英砂中的静力触探试验,基于任意拉格朗日欧拉格式的大变形有限元方
2、法,模拟静力触探仪的完整贯入过程,有效避免了锥尖周围的网格扭曲。大变形模拟中采用修正摩尔库伦模型描述砂土的内摩擦角和剪胀角随累积塑性剪应变的变化,并通过修正弯曲元试验,确定合理的土体剪切模量。标定罐试验和大变形有限元方法获得的锥尖阻力结果能很好地相互验证;在此基础上,建立高应力条件下静力触探锥尖阻力与砂土相对密实度之间的关系,与已有的低应力水平试验结果进行对比,证明建立的经验公式能通过锥尖阻力很好地预测砂土的相对密实度。关键词:标定罐;静力触探;锥尖阻力;深部土层;力学特性;砂土中图分类号:TU411.93 文献标志码:A 文章编号:2096-6717(2023)04-0049-07High-
3、pressure cone penetration tests in silica sand using calibration chamber and large deformation analysisKONG Zhuangzhuang,WANG Dong,ZHANG Minsheng,PEI Huimin(Shandong Provincial Key Laboratory of Marine Environment and Geological Engineering,Ocean University of China,Qingdao 266100,Shandong,P.R.China
4、)Abstract:The previous data interpretation equation for medium and low stress is no longer applicable when in-situ tests are carried out in deeply seated soils.In order to explore the mechanical properties of sand under high stress conditions,in this paper,the cone penetration test in silica sand wi
5、th different relative densities and confining pressure levels was performed in a self-developed high-pressure calibration chamber,and the Arbitrary Lagrangian-Eulerian method,a large deformation finite element approach,was used to study cone penetration in silica sands.Frequent mesh generations were
6、 conducted to avoid the distortion of soil elements around the cone tip.A modified Mohr-Coulomb constitutive model was introduced to describe the internal friction and dilatancy angle varied with the plastic shear strain in silica sands.A method for determining shear DOI:10.11835/j.issn.2096-6717.20
7、22.061收稿日期:20220329基金项目:中央高校基本业务经费(841962013);国家自然科学基金(42025702、U1806230)作者简介:孔壮壮(1996-),男,主要从事海洋岩土工程研究,E-mail:。王栋(通信作者),男,教授,博士生导师,E-mail:。Received:20220329Foundation items:Fundamental Research Funds for the Central Universities(No.841962013);National Natural Science Foundation of China(No.42025702
8、,U1806230)Author brief:KONG Zhuangzhuang(1996-),main research interest:offshore geotechnical engineering,E-mail:.WANG Dong(corresponding author),professor,doctorial supervisor,E-mail:.开放科学(资源服务)标识码OSID:第 45 卷土 木 与 环 境 工 程 学 报(中 英 文)modulus by correcting the bending element tests was proposed.Numeric
9、al results of cone tip resistance agree reasonably well with calibration chamber tests.An empirical equation for cone tip resistance varied with different relative densities was established under high stress condition,and the comparison with the existing low stress test results show that the establi
10、shed equation can predict well the relative density of sand by the cone tip resistance.Keywords:calibration chamber;cone penetration test;cone tip resistance;deep soils;methanical properties;sand深基础在陆上和海洋工程中的应用日益增多,例如,近海导管架平台的桩基埋深可达 150 m 以上,桩基底端的上覆有效应力水平约在 12 MPa 之间1;中国南海第二轮天然气水合物试采的储层埋深为 208 m,水平开
11、采套管承受的上覆有效应力约为 1.5 MPa2。工程实践中常采用原位试验确定土层的物理和力学性质,但大多针对中等或低应力水平土层,已有的原位试验数据解译方法也很少考虑高应力水平的影响。静力触探(CPT)试验是海洋工程勘察中最常用的原位试验方法,可通过量测的锥尖阻力经验性地预测砂层的相对密实度和内摩擦角3-5。为建立预测公式,需要进行一定数量的室内标定罐试验:在标定罐中制备均匀砂样,对砂样施加预定的上覆压力或围压,量测施加压力后砂样的相对密实度;将静力触探仪贯入砂样中,记录贯入阻力随深度的变化;改变砂土的相对密实度和压力水平,进行多个试验。Lunne等6总结了不同地区石英砂的标定罐试验和现场测试
12、数据,建立起锥尖阻力 qc和砂土相对密实度Dr值及平均有效应力 p 之间的关系,但该依托试验的 p 值不超过 500 kPa。Ahmadi等7开展了不同 K0条件下的标定罐试验,考虑不同 K0状态对锥尖阻力的影响,也提出过类似公式,但其研究中的 p 值不超过 200 kPa。目前公开报道的高应力条件下砂土标定罐试验和现场试验的数据极少。标定罐试验砂样准备时间长、成本高,数值模拟可以作为有效辅助手段,以完善锥尖阻力与相对密实度之间的关系。然而,传统的小变形或有限变形框架内的有限元方法很难模拟 CPT,原因是探头在土中的贯入导致周围土体网格的严重扭曲,使得计算无法进行。近年来,多种大变形数值方法已
13、被用于追踪触探仪在砂土中的贯入过程,如任意拉格朗 日欧 拉 方 法(Arbitrary Lagrangian Eulerian method,简称 ALE)8、耦合欧拉拉格朗日方法9和物质点法10。ALE 方法结合了拉格朗日和欧拉步骤,允许网格独立于材料移动,保持剩余的网格拓扑不变。与传统的拉格朗日分析相同,满足平衡方程、边界条件、外部载荷和接触条件,然后重新定位节点,并将所有变量从旧网格映射到新网格。因此,当发生大变形时,可以保持高质量的网格。Kouretzis等11采用 ALE 法模拟砂土中静力触探试验,并与离心机试验结果相互验证,建立了相对密实度和锥尖阻力之间的关系,但其试验中p 值不超
14、过340 kPa。笔者采用高围压标定罐试验和 ALE 大变形有限元方法,进行砂土中静力触探贯入,建立高应力条件下锥尖阻力与石英砂相对密实度的关系。大变形模拟中采用修正摩尔库伦模型(MMC)本构模型描述砂土内摩擦角和剪胀角随累积塑性剪应变的变化。1标定罐试验采用某种青岛石英砂进行标定罐试验,砂的颗粒级配曲线如图 1 所示,平均粒径 D50=0.173 mm,不均匀系数 Cu=1.70,最大孔隙比 emax=0.949,最小孔隙比 emin=0.490。早期标定罐侧壁是刚性的,当罐体直径不够大时,容易出现明显的边界效应。目前,标定罐都采用柔性双壁式标定罐,在刚性外壁和柔性内壁之间施加水头形成围压,
15、从而有效消除刚性边界对贯入阻力的增强效应。笔者使用的高压标定罐(图 2)也属于柔性双壁式:试样直径为 600 mm、高度为 750 图 1青岛某石英砂颗粒级配曲线Fig.1Particle size distribution curve of a silicaQingdao sandmm,试样四周被柔性囊包裹,围压最大为 5 MPa。采用砂雨法制备试样,罐口上部安装水平往复运动的砂斗,运动速度自动控制。砂斗底部开口宽度及其与砂面的相对高度为手动调节,试样的相对密实度取决于撒砂高度、开口宽度和运行速度12。砂中静力触探仪的贯入基本处于排水条件,Kluger 等13进行了饱和砂和干砂中的 CPT
16、比较试验,认为当应力水平很低时,干砂样的锥尖阻力会略高于饱和砂样;但当平均有效应力 p 增加到 200 kPa时,干砂样和饱和砂样的锥尖阻力无明显差异。因此,采用干砂雨法制备不同相对密实度的试样:为了制备两种不同相对密实度的砂样,砂样下落高度分别采用 0.5、1.3 m,每次撒砂的砂斗开口保持 3 mm 不变,运行速度控制在 3 mm/s。在标定罐内部每隔 0.1 m 进行标记,撒砂厚度一旦达到标记线,即升高砂斗 0.1 m,以保证砂样的下落高度基本不变。施加有效围压 3分别为 0.5、1.0、2.0 MPa,采用等向固结,即静止土压力系数 K0=1。提前在罐底布设 4 个体积为 200 mm
17、3的带底环刀,用于试验后确定砂样的相对密实度。施加围压后砂样大致可分为两组,相对密实度分别约为 46%和 60%。共完成 6 个不同相对密实度和围压条件下的标定罐试验,对应的条件见表 1。采用的 CPT 探头直径 D=20 mm,锥角为 60,探头贯入速度为 20 mm/s。2有限元模拟2.1本构模型及参数标定砂土在剪切时表现出的应变硬化或软化行为取决于砂土的相对密实度与应力水平。传统摩尔库伦模型规定砂土内摩擦角和剪胀角为定值,不能描述砂土的软化和剪胀性。这里采用 Hu 等14提出的 MMC 模型,通过引进内摩擦角和剪胀角随累积塑性剪应变 的变化来描述砂土的软化和剪胀性,其变化规律如图 3 所
18、示,图中 i、p和 cv分别为初始、峰值和临界内摩擦角。因为摩尔库伦类模型无法描述剪缩,所以假定当 1时,剪胀角为 0。此后剪胀角迅速线性增大,在=2时达到峰值剪胀角 p;剪胀角保持峰值至=3,然后线性减小,直至在临界应变 4处减小到 0。当砂样表现为硬化和剪缩性时,MMC 模型实质上退化为传统的摩尔库伦模型:内摩擦角保持为cv,剪胀角始终为0。对于一般石英砂,按照 Hu等14和 Zheng等15的建议,取 i=cv、1=1%、2=1.2%、3=5%和 4=15%。MMC 模 型 包 含 的 参 数 有:泊 松 比,取 定 值0.35;黏聚力 c 为 0;剪切模量 G;峰值内摩擦角 p;临界内
19、摩擦角 cv和峰值剪胀角 p。通过弯曲元试验推算 G,利用三轴排水剪切试验获得 p、cv和 p。2.1.1弯曲元试验弯曲元试验用于确定土体在微小应变(应变约为 10-6量级)时的最大剪切模量Gmax。将弯曲元量测设备安装在三轴仪上,在三轴仪中分别制备初始孔隙比为 0.694 和 0.764 的砂样,逐级施加 1001 500 kPa 的有效围压。每一级的固结完成后,首先测定排水量,用于计算固结后的孔隙比 e,然后进行弯曲元试验,获得微小应变水平的 Gmax。Gmax是 e和 p 的函数16。图 4为弯曲元试验得到的归一化最大剪切模量Gmax/pa随 p/pa的变化,试验数据可拟合为表 1标定罐
20、中砂样相对密实度Table 1Relative density of sand sample in calibration chamber图 2中国海洋大学研制的高压标定罐Fig.2High-pressure calibration chamber developed in Ocean University of China图 3MMC模型中内摩擦角和剪胀角的变化Fig.3Variation of friction and dilation angles of modified Mohr-Coulomb model50第 4 期孔壮壮,等:高压条件下石英砂中的静力触探试验及大变形模拟mm,试样
21、四周被柔性囊包裹,围压最大为 5 MPa。采用砂雨法制备试样,罐口上部安装水平往复运动的砂斗,运动速度自动控制。砂斗底部开口宽度及其与砂面的相对高度为手动调节,试样的相对密实度取决于撒砂高度、开口宽度和运行速度12。砂中静力触探仪的贯入基本处于排水条件,Kluger 等13进行了饱和砂和干砂中的 CPT 比较试验,认为当应力水平很低时,干砂样的锥尖阻力会略高于饱和砂样;但当平均有效应力 p 增加到 200 kPa时,干砂样和饱和砂样的锥尖阻力无明显差异。因此,采用干砂雨法制备不同相对密实度的试样:为了制备两种不同相对密实度的砂样,砂样下落高度分别采用 0.5、1.3 m,每次撒砂的砂斗开口保持
22、 3 mm 不变,运行速度控制在 3 mm/s。在标定罐内部每隔 0.1 m 进行标记,撒砂厚度一旦达到标记线,即升高砂斗 0.1 m,以保证砂样的下落高度基本不变。施加有效围压 3分别为 0.5、1.0、2.0 MPa,采用等向固结,即静止土压力系数 K0=1。提前在罐底布设 4 个体积为 200 mm3的带底环刀,用于试验后确定砂样的相对密实度。施加围压后砂样大致可分为两组,相对密实度分别约为 46%和 60%。共完成 6 个不同相对密实度和围压条件下的标定罐试验,对应的条件见表 1。采用的 CPT 探头直径 D=20 mm,锥角为 60,探头贯入速度为 20 mm/s。2有限元模拟2.1
23、本构模型及参数标定砂土在剪切时表现出的应变硬化或软化行为取决于砂土的相对密实度与应力水平。传统摩尔库伦模型规定砂土内摩擦角和剪胀角为定值,不能描述砂土的软化和剪胀性。这里采用 Hu 等14提出的 MMC 模型,通过引进内摩擦角和剪胀角随累积塑性剪应变 的变化来描述砂土的软化和剪胀性,其变化规律如图 3 所示,图中 i、p和 cv分别为初始、峰值和临界内摩擦角。因为摩尔库伦类模型无法描述剪缩,所以假定当 10D 时,qc略微减小或基本保持不变。qc略微减小的原因可能是在试样上放置的顶盖使得浅层试样的密实度略高于深部。高应力条件下锥尖阻力仍满足随相对密实度和围压增大而增大的现象。对比 Dr约为 4
24、6%的砂样,围压 1.0 MPa 或 2.0 MPa 时的锥尖阻力大约为0.5 MPa 时的 1.76 倍和 2.66 倍;对比 Dr约为 60%的砂样,围压 1.0 MPa 或 2.0 MPa 时的锥尖阻力大约为 0.5 MPa 时的 1.58 倍和 2.66 倍。这表明对于不同相对密实度的砂样,高应力条件下的锥尖阻力仍依赖围压水平。如果围压相同,Dr约为 60%砂样的锥尖阻力是 Dr约为 46%砂样的 1.171.3倍。整体来说,ALE 方法能够较好地模拟高围压石英砂中的 CPT贯入。从图 6中可以看出:1)对于不同相对密实度的砂样,当贯入深度H10D 时,有限元计算得到的 qc基本不再变
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 高压 条件下 石英砂 中的 静力 试验 变形 模拟
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。