Dynamic responses of hollow steel pipes directly buried in high-saturated clay to blast waves
-
摘要: 设计和实施了系列爆炸波作用下钢管的动态响应实验，获得了管道应变、振动速度和加速度及地表振动速度的时间历程。分析实验数据可知：在爆炸波近中场，峰值应变大小同管土相对刚度因数负相关，并随比例距离呈幂函数形式衰减，且在爆炸波不同分区衰减指数不同；地表峰值粒子速度、各管的峰值振动速度和测点峰值应变之间具有良好的线性关系。对各测试量做FFT (Fast Fourier Transformation)频谱分析得：振动信号主要集中在低频段，质心频率在10～60 Hz范围内，频谱同天然地震波谱有明显差异，动应变谱质心频率随药量增加亦呈幂函数形式衰减；考虑了爆腔因素以后，管道和地表振速频谱的质心频率和比例距离取对数后具有线性衰减关系。实验所得数据可应用于相似条件下管道的抗震计算，一些结论可作为深入研究爆炸波作用下埋地管道冲击振动机理的依据。Abstract: A series of experiments were designed and implemented to explore dynamic responses of steel pipes to blast waves The time histories of dynamic strains, vibration velocities and accelerations of the steel pipes were obtained, and the vibration velocity-time curves were gained. It is known by analyzing the experimental data that in the near and middle fields of the blast wave, the peak dynamic strains are negatively correlated with the relative stiffness coefficient of the pipe and the soil, and they follow the attenuation law of the power function with scaled distance; that the attenuation indexes are different for the blast wave propagating in the different sections of the field. The peak particle vibration velocities of the ground and the pipes have good linear correlations with the peak strains of the measuring points at the pipes. On the basis of the spectrum analysis by the fast Fourier transform on each test quantity, it is found that the spectrum energy of each test quantity is mainly concentrated in the low-frequency band, the centroid frequency is in 10−60 Hz, but there are obvious differences compared to the spectra of natural seismic waves. The centroid frequency of the dynamic strain spectrum is decayed in the exponent form of a power function with the increase of the explosive charge. By taking logarithm, there is a linear attenuation relationship between the scaled distance and the centroid frequencies of velocity spectra that the blasting cavity factor is considered. The test data can be directly applied to the seismic calculation under similar conditions, and some conclusions can be used as the theoretical basis for further study of the impact damage mechanism of buried pipelines.
-
Key words:
- blast wave /
- spectrum analysis /
- buried pipe /
- dynamic strain /
- vibration velocity /
- vibration acceleration
-
玻璃钢全称玻璃纤维增强树脂，是复合材料的典型代表，在实际工业与生活中有多种用途[1]。玻璃钢具有隔热、耐疲劳、阻燃、耐腐蚀、可设计性强等优点[2-6]，从20世纪60年代开始，英国、日本、德国就将玻璃钢应用于列车的结构与非结构件上[7]。玻璃钢由于结构的非均匀性，在力学性能上具有各向异性和层间剪切模量较低等复杂特性[8-9]。
针对玻璃钢复合材料力学行为方面的研究工作，国内外学者已经做了大量的工作。Reis等[10]针对玻璃钢材料进行了多种温度下的拉伸，并通过实验结果得到玻璃钢的冲击拉伸强度会受到温度影响的结论。Moradpour等[11]针对不同树脂基的玻璃钢进行了拉伸实验，发现不同的树脂基体会影响到玻璃钢的强度，通过对比实验结果确定了PMDI作为基体能获得最好的结果。刘华等[12]对玻璃钢进行了拉伸实验，并通过SEM观测破坏断面图发现玻璃钢的失效形式为纤维和树脂基体的组合破坏，可以通过改善树脂与玻璃纤维的界面粘结性能，获得性能更为优异的玻璃钢复合材料。张磊等[13]对玻璃钢进行了顺纤维方向和垂直纤维方向的拉伸实验，通过实验结果说明实际使用中应将玻璃钢材料沿纤维方向制成拉杆或管状形式，以便充分利用其优异的抗拉性能，将玻璃钢材料和输电杆塔中的钢材进行对比，发现了同等载荷下，玻璃钢材料制成的杆塔结构发生的变形更大，宜在使用时关注结构本身的变形控制。陈鲁等[14]对轻型结构做了大量实验，并与钢材进行对比，实验结果表明材料的强度和模量比钢材低很多。实际上，玻璃钢的强度确实比一般钢材低，但玻璃钢密度比钢材小，因此玻璃钢的比强度要比钢材更强。
数字图像相关技术(digital image correlation，DIC)是一种通过匹配物体表面不同状态下的数字化散斑图像上的几何点，跟踪点的运动获得物体表面变形信息的方法。通过DIC技术可以方便地获得实验试件的位移场。张燕南等在进行全场位移分析的时候，就用到了DIC数字图像相关技术，并得到了位移场呈倾斜的梯度分布的结论[15]。然而利用DIC技术对单向增强玻璃钢材料进行静/动态拉伸时的应变场以及精确的工程失效应变、强度和动态失效应变进行分析的有关研究还比较少。
本文针对单向增强玻璃钢复合材料，进行不同方向、不同应变率的静/动态拉伸实验，利用一整套静/动态拉伸装置以及高速摄影与DIC相结合的技术，采集材料拉伸过程中的全场应变信息，采用力传感器获得了材料拉伸过程的应力时域信号；结合应力及应变的时域信号结果，绘制材料不同方向、不同应变率的静/动态应力-应变曲线，根据DIC分析结果获得了材料不同方向、不同应变率破坏时刻的动态失效应变，以期更加精准地描述单向增强玻璃钢复合材料的动态失效行为。
1. 实验方案
1.1 实验设备
本实验针对单向增强玻璃钢材料进行了3种不同应变率的拉伸实验，分别是准静态（10−3 s−1）、10 s−1、100 s−1。其中准静态（10−3 s−1）拉伸实验采用万能试验机，另外两种动态拉伸实验采用中高速应变率材料试验机，如图1所示。中高速应变率材料试验机可以实现1～500 s−1应变率范围的材料拉伸试验。针对拉伸试件设计了专门的测力传感器，如图2所示。所设计的测力传感器不仅可以在拉伸过程中对力进行测量，还可以防止试件在拉伸过程中滑出。
为了便于DIC处理处应变和时间的信息，配备了高速相机来拍摄试验过程。高速相机型号为Photron SA1.1，在分辨率为1 024×1 024像素下，采样率可以达到5 400 s−1，通过降低分辨率，采样率最高可达675 000 s−1。实验选用了80 000 s−1的采样率，分辨率为896×120像素，既能准确地记录试验中每个试件不同时刻的散斑变形图片，还保证了后期处理软件对试件变形的处理精度。整套高速摄影装备如图3所示。
1.2 样品准备
实验材料采用青岛海铁船舰有限公司提供的玻璃纤维增强树脂。实验用的玻璃纤维增强树脂铺层方向为单向纤维增强方式，拉伸实验试件厚度为2 mm，切割方向选择0°、45°、90°三个方向，具体如图4所示。其中纤维增强方向，记为方向1，与纤维增强方向成45°的方向，记为方向2，与纤维增强方向垂直的方向，记为方向3。
图 4 纤维铺层方向及试件切割方向Figure 4. The layer directions of fibers and the cutting directions of specimens针对准静态（10−3 s−1）与动态拉伸试验不同的加载特点，设计了两种尺寸的拉伸试件。准静态（10−3 s−1）拉伸实验试件尺寸如图5(a)所示。试件厚度为2 mm，动态拉伸试件尺寸如图5(b)所示，试件厚度为2 mm。在拉伸试件的中间区域喷涂了人工散斑，用于玻璃钢试件全场变形量的测量。
1.3 结果采集
在试件表面喷涂散斑，拉伸过程如图6所示。利用高速摄影对试件表面的形貌进行拍摄得到散斑在整个动态拉伸过程中的形貌变化。采用DIC数字图像处理技术处理高速摄影拍摄的照片得到该工况下试件标距段位移的时空分布，对位移-时间曲线进行处理得到应变的时空分布。利用力传感器得到拉伸过程试件的载荷-时间曲线，对载荷-时间曲线进行处理得到应力-时间曲线。将应力-时间曲线和应变的时空分布进行联合处理得到试件的应力-应变曲线和动态失效应变。
相较于传统的引伸计测量应变方法，采用DIC图像处理技术的一个明显的优势为可以方便地获得应变的时空分布。在全场应变不均匀的情况下，可以针对性分析断口附近的应变-时间曲线，而不是用标距段平均应变来代替。
2. 结果与讨论
2.1 结果分散性讨论
试验针对特定方向和应变率的每一种工况都进行了三次重复试验，对得到的三组原始应力-应变曲线进行分散性分析，选取其中更接近平均值的应力-应变曲线作为代表来研究材料的力学性质。具体实验情况如表 1所示。
表 1 不同工况重复实验次数Table 1. Times of repeated experiments under different conditions方向 实验重复次数 准静态 10 s−1 100 s−1 方向 1 3 3 3 方向 2 3 3 3 方向 3 3 3 3 针对典型的方向1下，100 s−1应变率工况进行分散性讨论。三次重复试验得到的应力应变曲线和试件破坏图分别如图 7和图8所示。
图 7 典型工况下重复试验的应力-应变曲线Figure 7. Stress-strain curves of repeated tests under typical working conditions在三次重复试验下，试件最终的断裂形貌均发生在标距段，属于有效实验。针对数据分散性的讨论，选取了三个参考指标：弹性模量Einitial、工程失效应变和强度。分析方法为将指标项的三个结果分别和平均值作比较。经过比较发现，弹性模量的误差为8%，工程失效应变为10%，强度的误差为1%。在后续的分析中，将每种工况下的应力应变曲线均进行了相似的比较，所有的工况下最大误差均控制在10%的范围内。后续分析中全部应力-应变曲线的选取原则确定为各项参考指标数值居中的曲线。
2.2 典型工况的应变时空分布
对于动态拉伸过程中应变场的分析，选取了方向1应变率为100 s−1的情况下进行DIC图像处理。针对整个拉伸过程进行应变场的分析，得到玻璃钢试件不同时刻试件标距段形貌以及标距段全场应变如图9和图10所示。
图 10 方向1动态拉伸全场应变分布Figure 10. Strain field in direction 1 under dynamic tension condition通过整个拉伸过程的试件形貌图，可以发现在拉伸初期图10(a)～(b)由于内部玻璃纤维的单向增强分布导致材料的各向异性，材料表面应变分布并不均匀。拉伸中期图10(c)～(d)，材料标距段整体应变分布不均匀性进一步增强。断裂前期图10(e)，试件形貌方面部分区域颜色变白明显，由于方向1起增强作用的玻璃纤维是主要的受力部分，夹持端出现了两条白色带状区域。通过图10(e)也可以看出，试件局部应变显著增大。在整个断裂前期，拉伸方向位移场呈不均匀分布。断裂失效后图10(f)～(g)，材料发生破坏，观察试件形貌能明显地看到标距段左端出现破坏，标距段未断裂部分应变逐渐恢复为0。
2.3 典型工况的应力-应变曲线
方向1应变率为100 s−1的拉伸工况下获取的应力-应变曲线如图11所示。曲线中的时间标识与图9和图10所指代的相同。
图 11 方向1在100 s−1应变率下应力-应变曲线Figure 11. Stress-strain curve of the direction 1 at 100 s−1 strain rate观察整个动态拉伸应力-应变曲线，加载过程为一条近似通过原点的斜直线，加载中前期曲线的线性比较好，等到了后期，在即将断裂失效的时候，即图11中的刚度变化点N，直线的斜率变小。在N之前，主要是玻璃纤维在维系着复合材料的弹性变形；当达到N点的应力之后，玻璃纤维部分断裂，剩余纤维以及树脂基体共同承担载荷，树脂基的贡献使该复合材料增加了塑性应变的成分，纤维的断裂以及塑性应变成分的增加导致材料的应力-应变曲线产生斜率上的变化。随着纤维断裂增多，剩余纤维以及树脂基体不能再承受断裂纤维传递的载荷，试件破坏。变化前的刚度通过弹性模量Einitial来衡量，刚度变化后的模量若记为Echanged，则该材料从35.1 GPa的弹性模量Einitial降低为7.5 GPa的Echanged。联系拉伸实验过程中陆续听到几次“啪”的纤维断裂响声。由此推测此时材料内部部分纤维开始断裂。随着纤维断裂增多，邻近纤维不能再承受断裂纤维传递的载荷，试件开始破坏。后期材料内部的纤维的损伤导致材料的刚度发生了变化。
图12为准静态（10−3 s−1）拉伸情况下材料在方向1、2、3上的应力-应变曲线。如图12所示，准静态情形下三个方向的强度分别为608、75、90 MPa，工程失效应变分别为0.025、0.021、0.022，图中曲线的选取原则以最接近平均结果为准。在弹性阶段可获得3个方向的弹性模量Einitial分别为27.7、3.6、4.1 GPa。同时方向1曲线工程失效应变与强度是最大的。内部的玻璃纤维束设计的增强方向主要在方向1，但实际生产过程中需要在方向3也增加少量的玻璃纤维对方向1的玻璃纤维束进行固定，因此方向3的强度比方向2要高20%。此材料内部的玻璃纤维束不仅可以显著改善材料的强度，还使得材料的工程失效应变有小幅度的增加。方向1的工程失效应变较方向2增加了19%。在即将断裂的时候，方向1曲线在图中刚度产生了变化。类似于图11中典型工况应力-应变曲线的刚度变化点N。在N点之后，玻璃纤维部分断裂，剩余纤维以及树脂基体共同承担载荷，树脂基的贡献使该复合材料增加了塑性应变的成分，变化前的刚度通过弹性模量Einitial来衡量，刚度变化后的模量若记为Echanged，则该材料从27.7 GPa的弹性模量Einitial减小为18.7 GPa的Echanged。
图 12 准静态拉伸不同方向应力-应变曲线Figure 12. Stress-strain curves in different directions under quasi-static condition图13～15为材料在3种不同应变率下不同方向的应力-应变曲线。
图 13 方向1不同应变率拉伸应力-应变曲线Figure 13. Tensile stress-strain curves of the direction 1 under different strain rates
图 15 方向3不同应变率下拉伸应力-应变曲线Figure 15. Tensile stress-strain curves of the direction 3 under different strain rates
图 14 方向2不同应变率下拉伸应力-应变曲线Figure 14. Tensile stress-strain curves of the direction 2 under different strain rates如图13所示三种应变率（准静态（10−3 s−1）、10 s−1、100 s−1）下的强度分别为608、967、1 123 MPa，工程失效应变为0.025、0.066、0.040。结果表明，随着应变率的提高，材料的强度也得到了显著增强，具备应变率效应。而工程失效应变虽然在两种动态拉伸情况下都比准静态（10−3 s−1）拉伸情况下的工程失效应变要大，但是应变率为10 s−1的情况下的工程失效应变却比应变率为100 s−1情况下的要大，因此，与强度不同，工程失效应变并没有表现出随着应变率的增大而增大的特性。
在材料的拉伸过程后期，会听到纤维断裂的“啪”声，纤维的断裂导致材料的应力-应变曲线产生斜率上的变化。随着纤维断裂增多，邻近纤维不能再承受断裂纤维传递的载荷，试件开始破坏。纤维的陆续断裂宏观上表现为材料刚度发生变化，准静态（10−3 s−1）应变率拉伸下由27.7 GPa的弹性模量Einitial减小为18.7 GPa的Echanged，10 s−1应变率拉伸下由35.1 GPa的弹性模量Einitial减小为13.7 GPa的Echanged，100 s−1应变率拉伸下由35.1 GPa的弹性模量Einitial减小为7.5 GPa的Echanged。三种应变率情况下的初始弹性模量都在30 GPa左右。材料的刚度在刚度变化点N之后产生了下降，改变后的Echanged按照应变率从低到高分别只有初始弹性模量Einitial的67.5%、39.0%、21.4%。可见，不同应变率拉伸情况下，尤其是两种动态工况下，树脂基的贡献导致的材料刚度变化更加明显。且变化率随应变率的增加而增加。
如图14所示，方向2在准静态（10−3 s−1）、10 s−1、100 s−1三种应变率下的强度分别为75、67、58 MPa，工程失效应变为0.020、0.011、0.013，强度上没有表现出正的应变率效应。方向2与纤维增强方向成45°，拉伸时候主要是基体承受载荷。且整个拉伸过程仅在最后断裂时刻听到一次剧烈响声，没有纤维陆续断裂的过程。整个拉伸过程，刚度不断变化，不能找到明显的刚度变化点N。
如图15所示，方向3在准静态（10−3 s−1）、10 s−1、100 s−1三种应变率下的强度分别为90、151、221 MPa，工程失效应变为0.022、0.043、0.038，强度上表现出应变率效应，工程失效应变类似于方向1，呈现出先增大后减小的特点。方向3与纤维增强方向成90°, 内部的玻璃纤维束设计的增强方向主要在1方向，但实际生产过程中需要在方向3也增加少量的玻璃纤维对方向1的玻璃纤维束进行固定，因此方向3的在三个方向上的强度比方向2都大，且强度具备应变率效应。方向3拉伸时候主要是基体和少量玻璃纤维承受载荷。整个拉伸过程，不能听到明显的纤维断裂声音，仅在最后断裂时刻听到一次剧烈响声。整个拉伸过程没有找到明显的刚度变化点N。
三个方向在不同应变率下的强度及工程失效应变对比，以及1方向刚度变化点N前后不同阶段的刚度大小，如表2所示。
表 2 不同应变率、不同方向的强度、工程失效应变及方向1刚度Table 2. Tensile strength and engineering failure strain of different strain rates in different directions and the stiffness in 1 direction应变率 强度/MPa 方向 1 Einitial/GPa 方向 1 Echanged/GPa 失效应变 方向 1 方向 2 方向 3 方向 1 方向 2 方向 3 准静态 608 75 90 27.7 18.7 0.025 0.020 0.022 10 s−1 967 67 151 35.1 13.7 0.066 0.011 0.043 100 s−1 1 123 58 221 35.1 7.5 0.040 0.013 0.038 在不同方向上的强度方面，结合此前张硕等[16]对单向纤维增强的玻璃钢进行的不同方向的准静态拉伸实验：沿纤维方向的强度值达到了垂直纤维方向强度值的10倍左右。实验进一步在动态拉伸实验下得到了相似的结论，沿纤维方向的强度值比其他方向不在同一个数量级，方向1的拉伸强度比方向3的拉伸强度高一个数量级，比方向2的拉伸强度高两个数量级。可见在玻璃钢的实际应用中，拉伸载荷主要是由玻璃纤维来承担的，而基体树脂更多的起到了一个固定的作用。
玻璃钢材料由于内部的纤维铺层排列，在宏观上表现出了力学的各向异性。方向1和方向3，强度上都表现出了应变率效应。因为方向1是纤维的增强方向，而为了固定1方向上的纤维束，需要在方向3也添加少量的玻璃纤维，因此在方向3也表现出了应变率效应。在方向2上，载荷主要由树脂基体来承受，强度上不具备应变率效应。在纤维增强的1方向上，材料的弹性模量Einitial在应力-应变曲线上具有一个刚度变化点N。材料的刚度在点N之后产生了下降，改变后的Echanged在准静态（10−3 s−1）、10 s−1、100 s−1三种应变率下的值只有初始弹性模量Einitial的67.5%、39.0%、21.4%。不同应变率拉伸情况下，尤其是两种动态工况下，材料的刚度的变化更加明显。且刚度的变化率随应变率的增加而增加。
2.4 动态失效参数分析
利用高速摄影与DIC相结合的方法，获得了100 s−1应变率下，不同铺层方向破坏时刻的动态失效参数。x和y方向分别表示拉伸方向和面内与拉伸方向垂直向，最终方向1断裂瞬间的试件形貌、应变场分布、x和y方向的应变随动态拉伸过程变化如图16所示，其中图(a)中选框部分的应力云图见图(b)，图(b)椭圆区域处的应力应变曲线见图(c)。
图 16 方向1在100 s−1应变率下动态失效区域Figure 16. Dynamic failure region in the direction 1 at 100 s−1 strain rate由图可以获得材料在100 s−1应变率动态拉伸下，最终断裂处的εx=0.387，εy=−0.013 24，鉴于二维试样，假设此时满足εy=εz。应力、应变都是二阶张量，当忽略了球张量部分，只考虑偏张量部分，就可以引入等效应力和等效应变的概念。断裂局部的等效应变即为动态失效应变
εeff 。εeff=√23√(εx−εy)2+(εy−εz)2+(εz−εx)2 将相应数值代入上式，得：εeff=0.267。
方向2和3的断裂瞬间的试件形貌、应变场分布、x和y方向的应变随动态拉伸过程变化如图17～18所示，其中图(a)中选框部分的应力云图见图(b)，图(b)椭圆区域处的应力应变曲线见图(c)。
图 17 方向2在100 s−1应变率下动态失效区域Figure 17. Dynamic failure region in the direction 2 at 100 s−1 strain rate
图 18 方向3在100 s−1应变率下动态失效区域Figure 18. Dynamic failure region in the direction 3 at 100 s−1 strain rate同理可得，在100 s−1应变率动态拉伸下，方向2和3的动态失效应变为0.078和0.099。结果表明，方向2、3的动态失效参数较为接近，与方向1相差较大。
综上，表3给出了最终得到此种单向铺层玻璃钢材料在应变率为100 s−1情况下动态拉伸参数。
表 3 不同方向的动态失效应变及其它相关参数(应变率：100 s−1)Table 3. Dynamic failure strain and other relevant parameters in different directions (strain rete: 100 s−1)方向 Einitial/GPa σf/MPa εf εeff 1 27.7 1 123 0.040 0.267 2 3.6 58 0.013 0.078 3 4.1 221 0.038 0.099 由于材料在动态加载过程中，标距段应变场分布不均匀，材料断口处的动态失效应变与常规的工程失效应变存在很大差别，工程失效应变在方向1、2、3上仅有动态失效应变的15.0%、16.7%、38.4%。而确定材料的动态失效应变能更为精确的了解材料的动态性能，在实际结构分析中可以更好地模拟结构在动态加载情况下的破坏行为。
3. 结 论
本次实验针对单向增强玻璃钢复合材料的动态力学性能的研究，提出了静/动态拉伸实验结合DIC的方法：设计一系列针对单向增强玻璃钢复合材料的静/动态实验，利用一整套静动态拉伸装置以及高速摄影与DIC相结合技术，获得了材料不同方向、不同应变率的应力-应变曲线以及材料在不同方向上的动态失效应变，精确的描述了材料的动态失效行为。通过实验结果还可以得到以下结论。
（1） 方向1在准静态（10−3 s−1）、10 s−1、100 s−1三种应变率下的拉伸应力-应变曲线均存在一个刚度变化点N。在N之前，主要是玻璃纤维在维系着复合材料的弹性变形；当达到N点的应力之后，玻璃纤维部分断裂，剩余纤维以及树脂基体共同承担载荷，树脂基的贡献使该复合材料增加了塑性应变的成分，纤维的断裂以及塑性应变成分的增加导致材料的应力-应变曲线产生斜率上的变化。随着纤维断裂增多，剩余纤维以及树脂基体不能再承受断裂纤维传递的载荷，试件破坏。这一过程宏观上表现为材料的刚度发生变化。三种应变率情况下初始弹性模量Einitial为：27.7、35.1、35.1 GPa，N点之后的模量Echanged为：18.7、13.7、7.5 GPa。三者的初始弹性模量较为接近，均为30 GPa左右，刚度变化点之后的Echanged在三种应变率下的值只有初始弹性模量Einitial的67.5%、39.0%、21.4%。不同应变率拉伸情况下，尤其是两种动态工况下，材料的刚度的变化更加明显。且刚度的变化率随应变率的增加而增加。由于只有方向1是纤维增强方向，方向2、3没有一个纤维陆续断裂的过程，故没有找到明显的刚度变化点N。
（2）此材料在不同应变率（10−3、10、102 s−1）拉伸情况下，纤维增强的方向1上强度最高（分别为608、967、1 123 MPa），方向2强度最低（分别为75、67、58 MPa），方向3强度较低（分别为90、151、221 MPa）。材料具备明显的各向异性，沿纤维方向的强度值比其他方向高一个数量级。在玻璃钢的实际应用中，拉伸载荷主要是由玻璃纤维来承担的，而基体树脂更多的起到了一个固定的作用。同时材料在方向1和3上也表现出了应变率效应，材料在102 s−1拉伸应变率下强度值都能达到准静态拉伸（10−3 s−1）下的2倍左右。方向1和3的工程失效应变还表现出来随着应变率增加先增大后减小的特点。
（3）利用高速摄影和DIC技术，对单向增强玻璃钢材料的三个不同方向进行了静/动态拉伸实验，得到了三个方向在100 s−1应变率的动态失效应变分别为0.267、0.078、0.099、100 s−1应变率下材料断口处的动态失效应变与常规的工程失效应变存在很大差别，工程失效应变在方向1、2、3上仅有动态失效应变的15.0%、16.7%、38.4%。这种方法获得的材料动态失效应变更加准确，可以更好的描述此单向增强玻璃钢复合材料的动态失效行为。
-
图 4 不同测点的动应变波形对比
Figure 4. Comparison of dynamic strain waveforms at different measured points
图 5 S333应变衰减曲线(R=2.7 m，H=1.5 m)
Figure 5. Strain decay curves of S333 (R=2.7 m, H=1.5 m)
图 6 S233应变衰减曲线(R=2.7 m，H=1.0 m)
Figure 6. Decay curves of S332 strain (R=2.7 m, H=1.0 m)
图 7 S331Z最大拉和压应变(H=1.5 m)
Figure 7. The maximum tension and compression strains of S331Z (H=1.5 m)
图 11 动应变质心变频率衰减曲线(R=2.7 m, H=1.5 m)
Figure 11. The decay curves of dynamic strain centroid frequencies (R=2.7 m, H=1.5 m)
图 12 S3管道振动的峰值粒子速度和地表振动的峰值粒子速度随比例爆距的衰减
Figure 12. Attenuation of peak particle velocity with scaled explosion distance for vibrations of S3 pipe and ground
图 13 S3管道振动的峰值粒子速度(PPVPV)和地表振动的峰值粒子速度(PPVGV)的关系
Figure 13. Relations between peak particle velocities of S3 pipe vibrations (PPVPV) and ones of ground vibrations (PPVGV)
图 16 S3管振速质心频率衰减
Figure 16. Centroid frequency attenuation of pipe S3 vibration velocity
图 17 地面振速质心频率衰减
Figure 17. Centroid frequency attenuation of ground vibration velocity
图 18 S3管道的应变和地表振动峰值粒子速度（PPVGV）的关系
Figure 18. Relationships between strain of pipe S3 and peak particle velocity of ground vibration (PPVGV)
表 1 管道几何参数
Table 1. Geometrical parameters of pipes
管道编号 D/mm δ/mm L/m S1 110 4.24 4.6 S2 160 4.70 4.6 S3 300 4.40 4.6 
下载: 导出CSV
表 2 实验方案
Table 2. Experimental schemes
H/m R/m Q/g 2.0 2.2 50 75 100 125 150 175 200 2.7 50 75 100 125 150 175 200 1.5 2.2 50 75 100 125 150 175 200 2.7 50 75 100 125 150 175 200 3.2 50 75 100 125 150 175 200 1.0 2.7 50 75 100 125 150 175 200
下载: 导出CSV
表 3 实验参数
Table 3. Experimental parameters
H/m R/m Q/g ˉR/(m·kg−1/3) 管道编号 1.0 2.7 75 6.402 S3 1.0 3.2 125 6.4 S3
下载: 导出CSV
表 4 校验结果
Table 4. Results of check
测点编号 相关系数 可信概率 测点编号 相关系数 可信概率 S323H 0.836 0.872 S332Z 0.855 0.893 S311H 0.786 0.834 S341Z 0.897 0.910
下载: 导出CSV
表 5 S3管不同测点的最大峰值应变(R=2.2 m)
Table 5. The maximum peak strains at different measured points of pipe S3 (R=2.2 m)
Q/g 最大峰值应变/10−6 S321H S341Z S323H S323Z S331Z S332H S332Z S333H S333Z S334Z 50 −31.1 35.1 −50.7 −31.7 104.3 53.8 −66.3 48.3 44.5 −58.5 75 −32.1 39.4 −51.8 −35.8 103.9 62.7 −75.4 43.8 50.4 −65.9 100 −30.7 52.9 −63.7 −35.9 97.3 61.1 −77.8 45.2 52.9 −61.6 125 −26.4 44.3 −45.2 −33.3 122.5 74.5 −81.2 60.3 54.3 −58.2 150 −42.7 60.7 −98.1 56.5 174.9 90.6 −117.7 69.6 76.4 −78.1 175 −78.5 72.2 −71.7 −64.3 187.8 93.6 −118.6 79.9 79.2 −79.3 200 −53.9 86.3 −89.2 −84.1 228.9 103.4 −131.3 100.5 90.4 −87.4 注：压应变按绝对值最大给出，下同。
下载: 导出CSV
表 6 S3管不同测点的最大峰值应变(R=2.7 m)
Table 6. The maximum peak strains at different measured points of pipe S3 (R=2.7 m)
Q/g 最大峰值应变/10−6 S321H S341Z S323H S323Z S331Z S332H S332Z S333H S333Z S334Z 50 −16.5 7.2 −33.7 −14.4 −26.1 44.5 14.6 21.0 7.6 −29.8 75 −22.4 10.4 −52.0 −20.2 −34.9 61.1 19.4 28.1 11.2 −28.2 100 −28.2 14.9 −67.2 −28.3 −49.4 80.7 28.8 36.2 17.9 −31.3 125 −30.5 17.6 −69.6 −30.9 −58.3 88.7 33.6 39.1 23.9 −42.0 150 −36.5 21.0 −91.0 −37.1 −71.3 102.9 42.9 48.2 30.1 −48.4 175 −38.2 33.2 −103.2 −48.5 −94.1 121.3 52.3 60.1 34.5 −57.0 200 −42.1 28.9 −103.7 −44.6 −85.7 118.9 49.3 54.5 34.3 −48.5
下载: 导出CSV
表 7 S3管不同测点的最大峰值应变(R=3.2 m)
Table 7. The maximum peak strains at different measured points of pipe S3 (R=3.2 m)
Q/g 最大峰值应变/10−6 S321H S341Z S323H S323Z S331Z S332H S332Z S333H S333Z S334Z 50 −30.2 18.6 −36.9 −22.1 −49.1 39.3 −26.9 27.6 −15.5 34.4 75 −16.2 17.4 −32.0 −22.3 −45.9 31.8 −25.5 26.2 −16.8 32.8 100 −18.6 22.7 −36.6 −26.8 −52.8 35.4 −32.9 26.5 −19.6 33.3 125 −19.1 22.5 −38.9 −32.7 −34.1 45.3 −30.7 31.1 −20.9 37.9 150 −19.3 24.7 −41.5 −37.5 −61.1 52.3 −35.1 33.0 −22.1 38.2 175 −23.0 28.0 −43.0 −36.6 −71.8 53.0 −41.3 34.8 −28.6 40.9 200 −25.9 29.5 −47.7 −38.7 −90.8 58.1 −43.0 38.7 −38.1 42.5
下载: 导出CSV
表 8 S3管道振动的峰值粒子速度和地表振动的峰值粒子速度(R=2.7 m, H=2.0 m)
Table 8. Peak particle velocities of S3 pipe and ground vibrations (R=2.7 m, H=2.0 m)
Q/g 峰值粒子速度/(cm∙s−1) 比值/% Q/g 峰值粒子速度/(cm∙s−1) 比值% S3管道振动 地表振动 S3管道振动 地表振动 50 4.3 13.8 31 150 9.4 25.1 37 75 3.7 11.4 32 175 12.1 30.6 39 100 4.9 15.0 33 200 14.5 35.2 41 125 5.2 13.7 38
下载: 导出CSV
表 9 管道S3和地面振动速度质心频率(R=2.7 m, H=2.0 m)
Table 9. Centroid frequencies of pipe S3 and ground vibration velocities (R=2.7 m, H=2.0 m)
Q/g fc/Hz Q/g fc/Hz 管道 地表 管道 地表 X方向 Y方向 Z方向 X方向 Y方向 Z方向 X方向 Y方向 Z方向 X方向 Y方向 Z方向 50 24.3 18.0 27.5 18.3 22.3 18.3 150 22.6 26.9 28.8 29.7 16.0 29.8 75 25.1 19.7 25.8 14.9 19.6 16.4 175 20.4 23.0 22.7 32.9 16.5 39.9 100 25.4 22.4 28.8 19.1 18.0 15.5 200 33.2 24.5 26.9 36.3 15.7 45.8 125 26.2 26.2 36.8 28.6 17.0 15.2
下载: 导出CSV
表 10 S2管道振动的峰值加速度(H=1.5 m, R=2.7 m)
Table 10. Peak particle vibration acceleration of pipe S2 (H=1.5 m, R=2.7 m)
Q/g 加速度/(m·s−2) X方向 Y方向 Z方向 合成 Max Min Max Min Max Min 50 4.7 −11.1 3.7 0.3 6.7 −3.7 13.5 75 6.7 −12.9 4.0 0.0 8.1 −3.5 15.8 100 4.1 −10.4 2.2 −1.2 7.7 −2.5 13.2 125 4.7 −16.1 1.5 −2.9 8.4 −7.0 17.8 150 3.1 −15.5 0.9 −3.1 6.7 −6.5 17.1 175 1.7 −14.7 1.0 −3.3 6.8 −7.1 16.6 200 3.0 −18.2 0.3 −4.0 6.9 −7.4 20.1
下载: 导出CSV
表 11 加速度质心频率(H=1.5 m, R=2.7 m)
Table 11. Centroid frequencies of acceleration (H=1.5 m, R=2.7 m)
Q/g fc/Hz X方向 Y方向 Z方向 50 55.9 40.1 32.9 100 41.3 41.9 33.0 200 41.1 39.1 36.8
下载: 导出CSV
表 12 S3管加速度、速度和位移峰值(H=2.0 m, R=2.2 m)
Table 12. Peak acceleration, peak velocity and peak displacement of S3 (H=2.0 m, R=2.2 m)
Q/g dmax/mm vmax/(cm·s−1) amax/(m·s−2) (vmax⋅d−1max)/Hz (amax⋅v−1max)/Hz 50 0.99 7.45 12.73 75.60 170.85 75 0.82 5.89 11.23 71.69 190.78 100 1.02 7.21 15.17 70.50 210.32 125 1.33 9.33 23.59 70.40 252.82 150 1.57 10.19 20.53 64.92 201.37 175 1.88 11.85 25.31 62.91 213.55 200 2.54 16.68 41.92 65.58 251.26
下载: 导出CSV
-
[1] 李强, 陈德利, 屈洋. 爆破对输气管道本体影响的监测 [J]. 油气储运, 2015, 34(2): 190–194. DOI: 10.6047/j.issn.1000-8241.2015.02.017.LI Q, CHEN D L, QU Y. Monitoring on effects of blasting on gas pipeline body distortion [J]. Oil & Gas Storage and Transportation, 2015, 34(2): 190–194. DOI: 10.6047/j.issn.1000-8241.2015.02.017. [2] 程围峰, 梁旭, 王振宇. 隧道爆破施工对邻近输油管道的影响评价 [J]. 石油工程建设, 2011, 37(4): 44–46; 9. DOI: 10.3969/j.issn.1001-2206.2011.04.013.CHENG W F, LIANG X, WANG Z Y. Assessment of blasting impact in tunnel excavation on nearby oil pipeline [J]. Petroleum Engineering Construction, 2011, 37(4): 44–46; 9. DOI: 10.3969/j.issn.1001-2206.2011.04.013. [3] 张紫剑, 赵昌龙, 张黎明, 等. 埋地管道爆破振动安全允许判据试验探究 [J]. 爆破, 2016, 33(2): 12–16. DOI: 10.3963/j.issn.1001-487X.2016.02.003.ZHANG Z J, ZHAO C L, ZHANG L M, et al. Experimental investigation of blasting vibration safety criterion on buried pipeline [J]. Blasting, 2016, 33(2): 12–16. DOI: 10.3963/j.issn.1001-487X.2016.02.003. [4] 戴联双, 张海珊, 孟国忠, 等. 在役油气管道周边爆破作业风险分析 [J]. 油气储运, 2012, 31(11): 801–803; 887. DOI: 10.6047/j.issn.1000-8241.2012.11.001.DAI L S, ZHANG H S, MENG G Z, et al. Risk analysis of blasting operations around in-service oil and gas pipelines [J]. Oil & Gas Storage and Transportation, 2012, 31(11): 801–803; 887. DOI: 10.6047/j.issn.1000-8241.2012.11.001. [5] 王振洪, 侯雄飞, 边明, 等. 爆破对天然气长输管道振动影响的安全判据 [J]. 油气储运, 2016, 35(8): 813–818. DOI: 10.60 47/j.issn.1000-8241.2016.08.003.WANG Z H, HOU X F, BIAN M, et al. Safety criterion for impacts of bursting on long-distance gas pipeline [J]. Oil & Gas Storage and Transportation, 2016, 35(8): 813–818. DOI: 10.60 47/j.issn.1000-8241.2016.08.003. [6] SISKIND D E, STAGG M S, WIEGAND J E, et al. Surface mine blasting near pressurized transmission pipelines [R]. Minneapolis, MN: US Department of the Interior Bureau of Mines, 1994: 2−3; 22−23. [7] DOWDING C H. Blast vibration monitoring and control [M]. New Jersey: Prentice-Hall, 1985: 167−171. [8] ESPARZA E D, WESTINE P S, WENZEL A B. Pipeline response to buried explosive detonations: Vol. 2 [R]. Arlington, USA: American Gas Association, 1981: 146−152. [9] KOURETZIS G P, BOUCKOVALAS G D, GANTES C J. Analytical calculation of blast-induced strains to buried pipelines [J]. International Journal of Impact Engineering, 2007, 34(10): 1683–1704. DOI: 10.1016/j.ijimpeng.2006.08.008. [10] ABEDI A S, HATAF N, GHAHRAMANI A. Analytical solution of the dynamic response of buried pipelines under blast wave [J]. International Journal of Rock Mechanics and Mining Sciences, 2016, 88: 301–306. DOI: 10.1016/j.ijrmms.2016.07.014. [11] 钟冬望, 龚相超, 涂圣武, 等. 高饱和黏土中爆炸波作用下直埋聚乙烯管的动力响应 [J]. 爆炸与冲击, 2019, 39(3): 033102. DOI: 10.11883/bzycj-2017-0334.ZHONG D W, GONG X C, TU S W, et al. Dynamic responses of PE pipes directly buried in high saturated clay to blast wave [J]. Explosion and Shock Waves, 2019, 39(3): 033102. DOI: 10.11883/bzycj-2017-0334. [12] 冯双, 马郧, 徐光黎. 武汉地区基坑一般黏性土的力学特性与其物理指标间的相关性研究 [J]. 安全与环境工程, 2016, 23(5): 149–154. DOI: 10.13578/j.cnki.issn.1671-1556.2016.05.025.FENG S, MA Y, XU G L. Correlation between mechanical properties and physical indicators of general cohesive soil in foundation pits in Wuhan area [J]. Safety and Environmental Engineering, 2016, 23(5): 149–154. DOI: 10.13578/j.cnki.issn.1671-1556.2016.05.025. [13] 库特乌佐夫. 爆破工程师手册 [M]. 刘清泉, 译. 北京: 煤炭工业出版社, 1992: 137–138 [14] 卢文波, 张乐, 周俊汝, 等. 爆破振动频率衰减机制和衰减规律的理论分析 [J]. 爆破, 2013, 30(2): 1–6; 11. DOI: 10.3963/j.issn.1001-487X.2013.02.001.LU W B, ZHANG L, ZHOU J R, et al. Theoretical analysis on decay mechanism and law of blasting vibration frequency [J]. Blasting, 2013, 30(2): 1–6; 11. DOI: 10.3963/j.issn.1001-487X.2013.02.001. [15] 陈勇, 杨庆华. 埋地管道在地震波作用下的变形分析 [J]. 世界地震工程, 2015, 31(4): 275–279.CHEN Y, YANG Q H. Deformation analysis of buried pipelines under seismic waves [J]. World Earthquake Engineering, 2015, 31(4): 275–279. 期刊类型引用(6)
1. 李翔云，万军，冯俐，杨小乐，郝铭. 海洋平台玻璃钢挡风墙应用研究. 应用科技. 2023(01): 51-55 . 
百度学术2. 乔井彦，李金柱，张羲黄，姚志彦，申仕良. 环氧树脂玻璃钢的动静态拉伸力学特性. 高压物理学报. 2023(03): 29-38 . 百度学术3. 王昊康. 玻璃纤维增强复合材料力学性能研究进展. 合成纤维. 2022(11): 40-43+68 . 百度学术4. 肖志斌，武丽丽，裘雄伟，柯贤朝，蔡亮. 数字图像相关法在复合材料研究中的应用进展. 理化检验(物理分册). 2021(05): 39-45+49 . 百度学术5. 宋肖肖，李柯，周律，宋先捷，陈亚军. VIC-3D技术在材料力学单轴拉伸实验教学中的应用研究. 实验室科学. 2021(05): 45-47+52 . 百度学术6. 王付胜，欧阳俊杰，孔繁淇，邸则坤，陈亚军. 基于DIC研究应力比对T700碳纤维复合材料疲劳性能的影响. 复合材料科学与工程. 2020(11): 54-60 . 百度学术其他类型引用(5)
-
下载:
百度学术
计量
- 文章访问数: 5042
- HTML全文浏览量: 1586
- PDF下载量: 82
- 被引次数: 11