专利名称：：堆石体密度测定的附加质量法的制作方法
技术领域：
：本发明涉及一种堆石体密度的测定方法，具体地说是一种利用附加质量法测定水工建筑物堆石坝的堆石体密度的方法。在水工建筑物堆石坝的施工过程中，堆石密度监测是确保坝体质量的关键性工作，一般来讲，堆石材料的粒径级配是自然形成的，是极不均匀的，粒径小者到砂子，毫米级，大到巨石，1米左右。在这种情况下测定堆石的密度是比较困难的。工程中常用的方法有如下几种，规范的方法是挖坑、称重、量体积，即在测点挖坑，并称出所挖出砂石的质量，同时用灌水(铺塑料布在坑中灌水)或灌砂的办法测量所挖坑的体积，从而求出相应测点的密度值。这种方法的优点是可以直接得到坑测的密度值，比较直观，缺点是遇到大块石开挖困难，存在避重就轻等人为因素，所测密度值与实际情况有偏差，有时可能在巨石处产生较大偏差，另一方面，由于挖坑破坏了原堆石结构，需要将所挖之坑重新填充压实，费工费时，不便布置足够数量的测点，密度数据的代表性受到一定限制。第二种方法为压实降观测法，即观测压实前后堆石的高差以判断其密度是否满足设计要求，第三种方法为振动碾装加速度计法。第四种方法为控制压实遍数。第五种方法为静弹模法。第六种方法为动弹模法。后五种方法都不能直接测出密度数据，是一种间接的、相对的、定性的、宏观的方法。八十年代日本推出了GR-810面波探查机之后，很多研究者利用面波技术测定散填体或土基的密度，利用面波测堆石或地基土密度的方法基本上可分两条思路一条是频散曲线拟合法。由于波速的大小与介质的物理参数，如密度、剪切模量、压缩模量、泊松比密切相关，因此可以通过对实测资料如频散曲线的反演拟合得到岩土的横波速度、纵波速度、密度等参数。但拟合本身是有多解性的，即不同解的任何一种组合都有可能得到与实测资料的最佳拟合，没有唯一解。因此，这种方法至今未在工程中运用；文献《堆石坝压实密度快速无损检测新技术》(孙继增等，《水利水电技术》，1996年第1期)提出的另一条思路是利用数理统计方法，建立面波速度VR与堆石密度ρ的相关关系，该文献认为VR与ρ是线性相关的，这样利用所测VR很快即可推算出ρ。但从波动理论研究中可知VR、ρ是非线性关系，因此从测定VR去推算ρ的精度是很难保证的。中国专利ZL93107958.6公开了一种附加质量法测定地基承载力的方法，该方法采用由附加质量—承压板—地基土体组成的振动体系，用重锤击振体系使其振动，借助仪器测出体系相应的自振频率f，利用自振频率f、体系振动园频率ω、圆周率π之间的关系，代入公式D＝1/(2·π·f)2利用地基动刚度与相关函数D、参振质量△m的关系作D—△m曲线，求出地基动刚度，再建立动静关系式，从而得出地基承载力。该方法很好地解决了地基土的参振质量mo，为研究堆石体密度的附加质量法测定打下了良好的基础。本发明的目的是用附加质量法对堆石体密度进行原位、快速地测定。本发明的技术方案是这样实现的一种堆石体密度测定的附加质量法，其测定步骤为a在测点由附加质量块—承压板—地基坝体组成振动体系，用附加质量法测出振动体系相应的自振频率f，将f代入公式D＝1/ω2(ω＝2πf)作D—△m曲线，(△m为附加质量块)，找出D—△m曲线与横坐标交点到原点之间的绝对值即为地基坝体的参振质量mo；b在测点处布置检波点(点距≤0.5m)，地震仪接收对体系锤击所激发的弹性波，作地震时距曲线得直达波段曲线反斜率Vp，将Vp代入公式λ＝Vp/fo(fo为mo相应的体系自振频率)，得介质由地面向深部传播的纵波波长λ；c在测点采用挖坑、取样、称重、量体积的方法取得地基坝体点位的密度值ρ，将mo、λ、ρ代入公式ρ＝mo/(Ao·ho)得系数K(Ao为承压板面积，ho＝λ/(2K))；在同一地基的不同点位利用附加质量法测出各点位的参振质量mo、压板下介质纵波波长λ，利用所率定的系数K回代公式ρ＝mo/(Ao.ho)求得地基坝体测点位的密度。率定系数K时，所挖坑数至少为10，以保证有足够的样本数量，为采用算术平均法或线性回归法率定系数K提供充分的数据。为保证测试精度所采用的承压板为钢板，面积至少为2m2，厚度至少为25mm。本发明的积极效果是1.简便、快速、定量，在测量前率定系数K时只需开挖有限的坑，因而无损坝体结构的完整性；2.该方法适用广泛，不仅适用于堆石坝在施工碾压后的密度检测，对于未碾压后的天然或人工堆石、砂基、砂砾石地基以及各类土壤地基的检测亦能适用。本发明是建立在无阻尼模型理论基础上的，根据有关资料(《机械振动》，李文美、方同、陈松洪，科学技术出版社，1995)介绍，阻尼对振幅影响甚大，对频率影响甚微，因此阻尼对系统频率的影响非常小，在作频率测量时可以略去不计。本发明在作动参数测量时，只用频率，不用振幅，为了简化测量工作及资料处理工作，拟选择无阻尼模型为对象还是可行的。本发明在求解体系参振质量mo时采用至少为两级的附加质量法，这是因为(1)地基刚度及参振质量mo的测定必须在半无限地基上圈定出一个范围，而后才能利用质弹模型去量测，这个范围在放置承压板之前是未知的，然而用附加质量法解mo至少要有两个不同的体系自振频率，相应两级有不同的△m，承压板也有一定的质量，也应当算一级附加质量，承压板以下所覆盖的范围即为所圈定的测量范围；(2)在作体系的频率测量时往往有干扰振动及测量误差，在这种情况下，用多级附加质量法的曲线(D—△m)可以在一定范围和一定程度上消除这些因素的影响；(3)在单自由度质弹体系模型中，体系自振频率f有随△m增加而减小的理论关系，利用这种关系可以在多频振动的干扰中识别系统固有频率，这一点是多级附加质量法所具有的优点。在多种干扰下去识别系统固有频率，在瞬态激振中往往是一个非常困难而没有把握的事情，多数情况下凭经验估计，而D—△m曲线排除了干扰，解决了这一问题。(4)利用D—△m曲线求解mo方法简便，不需要大量分析计算工作，物理概念、几何概念非常明确。密度解析式的推导如下参振质量mo的动能To用动能公式表示为To＝mo·(Vo)2/2这里Vo为质点振动速度的幅值，设质点振动位移为Z，其幅值为Zo，质点振动速度为V，系统振动园频率为ωo，Z可以表示为Z＝Zo·sin(ωo·t)由动力学原理知V＝dZ/dt＝Zo·ωo·cos(ωo·t)当Z与V取幅值时cos(ωo·t)＝1，此时V＝Vo＝Zo·ωo因此有To＝mo·(Zo·ωo)2/2设地表下深度Z处，厚度为dz介质薄片的振动动能为dTz，则dTz＝dmz·(Vz)2/2其中dmz＝Az·ρ·dz，Vz、Az分别为相应该深度的薄片振动速度及薄片面积，将薄片振动动能等效为质点的动能，设深度Z处质量为dmz，质点振动幅值为Zz，振动园频率为ωz，类似地有质点振动速度的幅值Vz＝Zz·ωz，根据《振动计算与隔振设计》(姜俊平等，中国建筑工业出版社，1976)，深度Z与Zz、Zo、纵波波长λ及振动沿竖向衰减系数K之间的关系为Zz＝Zo·e-K-z/λ由于研究深度范围不大(一般＜2m)，故Z深度处的振动园频率ωz可近似等于ωo，综上有Vz＝Zo·ωo·e-K·z/λ，从而有dTz＝(1/2)·Az·ρ·(Zo·ωo)2·e-2·K·z/λ·dz令压板下介质振动的总动能To′等于介质薄片动能的积分即To′=∫o∞dTz]]>积分结果为To′＝ρ·(Zo·ωo)2·Ao·λ/(4·K)由于研究深度较小，压板下介质振动的扩散角对介质薄片面积的影响随深度变化不大(根据小浪底密度测试资料该扩散角＜1°)，忽略扩散角的影响有以上积分结果。令压板下参振质量mo的动能To等于介质薄片动能的积分To′，两式相等有ρ＝mo/(Ao·ho)其中ho＝λ/(2·K)为等效深度。K的率定方法有二，一是算术平均法，二是线性回归法。算术平均法即求取各已知密度(坑测而得)点位的K值之算术平均值；线性回归法则是利用密度解析式的变形形式K·mo＝Ao·λ·ρ/2令等式右边为N，则N与mo成线性关系，以N为纵坐标，mo为横坐标，将各点位由N与mo值组成的坐标点标在坐标系中，用最小二乘法拟合该直线，直线的斜率即为线性回归所率定的K值。下面结合实施实例对本发明作进一步详细介绍图1附加质量法求取mo现场工作布置示意2多级附加质量法D-△m曲线示意31#点位附加质量法D-△m曲线图4地震直达波法现场工作布置示意51#点位地震时距曲线实施例在测试现场选测点1#，第一，平整、铺2-3cm厚的砂土；第二，如图1所示，将2m2的承压板3徐徐平放在铺平砂土后的测点；第三，将拾振器2用石膏粘合剂藕合在承压板3中央并与中科院声学所研制的SL-II型振动信号分析仪5联接，用锤4击震承压板边侧堆石，使承压板3—堆石体组成的振动体系振动，此时，从振动信号分析仪5上显示出体系自振频率f，然后于承压板3上加一质量为1×103kg的附加质量1，再用锤4击震承压板3边侧堆石，使附加质量1—承压板3—堆石体组成的振动体系振动，此时，从振动信号分析仪5上显示出相应体系的自振频率，这样连续四次附加，每附加一次体系振动一次，振动信号分析仪5上显示一次体系振动频率f，所测各次数据列于表1。表1由公式D＝1/ω2＝1/(2π·f)2计算出相应的D值，其相应数据列于表1，然后利用表1中的△m值为横坐标，以相应D值为纵坐标在坐标系中找出各坐标点，取坐标点群中心线为D—△m曲线如图3所示，在曲线上任选一点，找出该点的横坐标△m＝4.0×103kg时，D＝2.252×10-5S2，△m＝1.0×103kg时，D＝1.667×10-5S2，由公式k＝△m/△D得曲线之反斜率k＝(4.0-1.0)/(2.252-1.667)＝5.128×108(kg·s-2)由曲线与纵坐标的交点得Do＝1.473×10-5S2，按公式mo＝Do·k，得参振堆石体质量mo＝1.473×10-5×5.128×108＝7.554×103kg如图4所示，在上述测点布置地震排列，采用地面直达波法测试弹性纵波速度Vp，在测点附近约2m处选定激发点放置激发板1，然后以检波点距0.5m安好速度检波器系列3，并与美国R24浮点式工程地震仪4联接，以大锤2击震激发板1，产生弹性波，弹性波向周围传播，被速度检波器系列3接收并传至地震仪4，量出激发点至各检波点之间的距离1.0m、1.5m、……，地震仪4所接收到的弹性波初至时间t及激发接收距离x如表2所示。表2</tables>由表2数据，以x为横坐标，t为纵坐标，绘制地震时距曲线，找出各坐标点，取坐标点群中心线为直达波曲线，如图5所示，其反斜率即为Vp，Vp＝(6.0-2.0)/(11.75×10-3-3.83×10-3)＝505m/s，由公式Do＝1/(2π·f)2，当Do＝1.473×10-5s2时，得fo＝41.469Hz，代入公式λ＝Vp/fo得λ＝505/41.469＝12.178m。在测点开挖一坑，用磅称称出所挖质量m＝3.550×103kg，将不透水塑料布放于坑中灌水，灌水体积为V＝1.642m3，根据公式ρ＝m/V得该点位密度为ρ＝3.550/1.621＝2.162×103kg·m-3，再将开挖前用附加质量法所测该点位的参振质量mo、纵波波长λ及压板面积Ao代入公式ρ＝mo/(Ao.λ/2K)求得该点位系数K＝3.485。类似地，对于该坝体2#、3#、4#、5#、……、12#点位，同样可获得各点位的系数K，如表3所示。该表列出了各点位系数K计算过程中的各个参数。k为附加质量法所得D—△m曲线的反斜率；Do为该曲线与纵坐标的交点；mo为k与Do相乘的结果，即参振质量；利用公式Do＝1/(2π·f)2，由Do可得fo；Vp为地震直达波法所得纵波波速；利用λ＝Vp/fo可得波长λ值；每级附加质量△m为1×103kg；土工试验原位密度为坑测值；利用公式ρ＝mo/(Ao·ho)由坑测密度值ρ及mo、Ao可得等效深度ho；再利用公式ρ＝mo/(Ao·λ/2K)可求得各点位的系数K值。采用线性回归法(或算术平均法)，可求得该坝体的率定系数K值为3.60(算术平均法率定系数K为3.614)。然后再按照附加质量法依次测定坝体未知密度点位的参振质量mo、纵波波长λ，利用在坝体中所率定的系数K＝3.60回代公式ρ＝mo/(Ao·λ/2K)即得该点位密度值ρ。权利要求1.一种堆石体密度测定的附加质量法，其测定步骤为a在测点由附加质量块—承压板—地基坝体组成振动体系，用附加质量法测出振动体系相应的自振频率f，将f代入公式D＝1/ω2(ω＝2πf)作D—△m曲线，(△m为附加质量块)，找出D—△m曲线与横坐标交点到原点之间的绝对值即为地基坝体的参振质量mo；b在测点处布置检波点(点距≤0.5m)，地震仪接收对体系锤击所激发的弹性波，作地震时距曲线得直达波段曲线反斜率Vp，将Vp代入公式λ＝Vp/fo(fo为mo相应的体系自振频率)，得介质由地面向深部传播的纵波波长λ；c在测点采用挖坑、取样、称重、量体积的方法取得地基坝体点位的密度值ρ，将mo、λ、ρ代入公式ρ＝mo/(Ao·ho)得系数K(Ao为承压板面积，ho＝λ/(2K))；在同一地基的不同点位利用附加质量法测出各点位的参振质量mo、压板下介质纵波波长λ，利用所率定的系数K回代公式ρ＝mo/(Ao.ho)求得地基坝体测点位的密度。2.根据权利要求1所述的堆石体密度测定的附加质量法其特征在于率定系数K时，所挖坑数至少为10。3.根据权利要求2所述的堆石体密度测定的附加质量法其特征在于所采用的承压板为钢板，面积至少为2m2，厚度至少为25mm。全文摘要本发明公开了一种堆石体密度测定的附加质量法，步骤a、附加质量法测地基的参振质量m文档编号E02D1/00GK1165292SQ96117338公开日1997年11月19日申请日期1996年12月24日优先权日1996年12月24日发明者李丕武申请人:黄河水利委员会勘测规划设计研究院物探总队