Mechanical properties of airbag frame ground beams for slope support
-
摘要: 针对地震、崩塌等自然灾害频发而导致的边坡失稳问题,提出了可快速应对抢险救灾工作的气囊式锚杆框架地梁新型结构,并阐述了该新型结构的工作机理。分析气囊框架地梁受荷时玻璃钢与气囊的受力特性,给出了气囊框架地梁承载力计算方法;结合气压传动和双层弹性地基梁理论,建立了气囊框架地梁-锚杆-土体协调变形的计算模型,给出了气囊式框架地梁在边坡支护施工阶段与工作阶段的力学效应计算方法。结果表明:①气囊式框架梁结构承载力主要由框架梁板厚、板材强度、气囊内压和高度决定,增加框架板厚、板材强度、气囊内压和高度可以提高结构承载力;②气囊式框架梁支护边坡存在施工和工作两个阶段,工作阶段位于坡脚处的框架梁剪力和弯矩出现较大的上升,设计时应作为最不利因素着重考虑;③支护设计时还应综合考虑支护效果、结构承载力和经济性给出设计方案,确保结构在实际抢险救灾中的可靠性、时效性和经济性。研究成果可在边坡快速支护中为新型结构的设计应用提供理论依据和指导。Abstract: In view of the slope instability caused by frequent natural disasters such as earthquakes and collapses, a new type of ground beam structure with airbag anchor frame which can quickly respond to rescue and relief work is proposed, and the working mechanism of the new structure is described. The stress characteristics of FRP and airbag are analyzed when the ground beam of airbag frame is loaded, and the method for the bearing capacity of the ground beam of airbag frame is given. Combined with the theory of pneumatic transmission and double-layer elastic foundation beam, the model for calculating the coordinated deformations of the airbag frame ground beams, bolts and soils is established, and the relevant method for the mechanical effects of the airbag frame ground beams at the construction stage and working stage of slope support is given. The results show that: (1) The bearing capacity of the airbag frame beam structure is mainly determined by the thickness of frame beam, the strength of plate, and the internal pressure and height of the airbag. Increasing the thickness of frame plate, the strength of plate, and the internal pressure and height of airbag can improve the bearing capacity of the structure. (2) The airbag frame beams for slope support have two stages: construction and working. The bending moments and shear forces at the working stage are larger than those at the construction stage, and the shear forces and bending moments of the frame beams at the foot of the slope rise greatly at the working stage, which should be considered as the most unfavorable factor in the design. (3) The support design should also comprehensively consider the support effects, structural bearing capacity and economy, and give the design scheme to ensure the reliability, timeliness and economy of the structure in the actual rescue and disaster relief. The research results may provide a theoretical basis and guidance for the design and application of the new structure in rapid slope support.
-
0. 引言
深筒形基坑因其平面尺寸小而开挖深度大,具有整体性强、空间效应明显等特点。这类基坑常因为规划需要而建设在河道堤岸旁边,使得基坑工程性状及设计与施工问题变得非常复杂,因而引起诸多工程及研究人员的关注。Faheem等[1]运用有限元法分析了三维矩形基坑的稳定性问题。龚维明等[2]、戴国亮等[3]对单片墙和单室井筒式地下连续墙基础的模型试验结果表明,这些结构在水平力作用下呈现整体倾斜破坏。王洪新[4]根据基坑宽度与插入深度的比值大小将基坑分为宽基坑、窄基坑和一般基坑,分别推导了3种基坑支护结构抗倾覆稳定安全系数的计算式。虽然对于类似深筒形基坑支护结构的整体稳定问题进行了一些探讨,但目前仍主要基于平面假定对平坦场地条件下的情况开展研究,而对临近岸坡等复杂场地下的基坑整体稳定问题缺乏研究。对这种临岸深筒形基坑支护结构整体倾覆破坏的规律性尚缺乏认识,还没有提出临岸深筒形基坑支护结构整体倾覆稳定的计算模式。为探讨这一问题的规律性,本文通过室内模型试验与数值模拟对临岸深筒形基坑支护结构整体倾覆破坏形式进行了研究,并提出了支护结构整体抗倾覆稳定安全系数的计算方法。
1. 室内模型试验
1.1 模型与试验条件
临岸基坑支护结构模型试验在长3 m,宽1 m,高1.2 m的试验槽中进行。土体采用不均匀系数为2.2的细砂土,利用有机玻璃板模拟基坑支护结构。支护结构长0.2 m,宽0.2 m,高0.8 m,厚2 mm,与实际工程支护结构的相似比例大约为1∶50。模型试验布置见图1,砂土及有机玻璃板的特性参数见表1。利用液压千斤顶加载,量测结构A、B、C侧水平位移及土压力。
表 1 试验材料参数Table 1. Parameters of materials试验材料 弹性模量E/MPa 泊松比 ν 重度 γ /(kN·m-3)内摩擦角φ/(°) 孔隙比e 细砂 60 0.30 14.5 32.61 0.78 有机玻璃板 3000 0.37 12.0 — — 试验采用先分层开挖后加载的方法,总开挖深度0.5 m,支护结构顶部与岸坡相距0.1 m,坡比为1∶2,坡高0.6 m。
1.2 试验步骤及测量
试验步骤为:①分层人工填土至结构底部;②安放支护结构模型,在结构上设土压力盒,安装位移测量装置;③继续填土至支护结构顶部;④模拟基坑分层开挖;⑤在坑外A侧施加超载。测量每一步完成后各侧结构水平位移及墙前后的土压力。
1.3 试验结果及现象
图2为通过定点拍照得到的支护结构顶部位移情况照片。由图2中参考线与支护结构相对位置可见,开挖至0.5 m深时,支护结构四侧只有微小向坑内位移;随着超载,基坑支护结构顶部整体向着岸坡发生很大的侧移,最终发生了很明显的整体倾斜。
2. 有限元数值模拟
2.1 建模与计算条件
根据室内模型试验条件,利用Z-SOIL软件建立如图3所示对称结构的一半模型。土体本构模型为弹性理想塑性模型,采用Mohr–Coulomb屈服及破坏准则,具体计算参数见表1。数值模拟通过设置接触单元考虑了有机玻璃板墙体与细砂土之间的摩擦接触特性;接触单元的法向及切向刚度系数按照接触面单元临近土层的弹性模量、接触面临近单元法向尺寸及经验调节因子(默认10-4)等因素由软件计算确定。
2.2 模拟结果及分析
数值模拟得到的基坑支护结构整体侧倾形态如图4所示。计算发现开挖阶段整个结构发生了正常量级的变形,而加载阶段结构整体深度范围内均向岸坡方向侧移,基本达到了倾覆程度。
结构A、B、C侧土压力分布情况如图5所示。开挖完成时,坑各侧内外土压力减小,不仅开挖面以上减到接近于零,开挖面以下也有明显的减小。随超载的增加土压力增大,深度0.2 m范围内土压力受超载的影响大。结构B侧土压力的分布情况在未加载之前与A侧相似。结构C侧土压力分布规律未加载之前与A侧相似,当施加超载时开挖深度0.5 m以下的土压力增大,开挖面以上土压力受影响不大。
图6为砂土中临岸深筒形基坑支护结构在不同嵌固深度下结构倾覆的位移矢量图。观察矢量图发现,当基坑支护结构发生倾覆破坏时,转动中心o随插入比即插入深度与支护结构高度之比h0/h不同而发生变化,但是大致都在距离支护结构底面为0.05h~0.10h,转动中心平面位置即距离支护结构底面中心偏心距e约为0.15b~0.30b。
根据数值模拟得出的土压力分布形式,土压力在支护结构底部有明显回收,转动中心高度以下土压力相对较小,对结构转动的影响也较小。假定深筒形基坑倾覆时,转动中心在支护结构底部平面内,与底面中心轴距离为e=Keb,Ke为转动中心位置系数。
3. 抗倾覆稳定安全系数
3.1 抗倾覆稳定分析模型
通过室内模型试验及有限元模拟发现,临岸深筒形基坑支护结构绕底部转动中心产生整体倾覆破坏,基于平面土压力理论的单片墙倾覆稳定分析方法不再适用,故建议支护结构整体抗倾覆稳定计算方法。
针对图7所示的支护结构倾覆破坏分析模式,本方法简化假定为:①转动中心在支护结构底面距中心轴上Keb处;②筒内土芯不随结构转动,仅对结构倾覆提供水平抗力;③考虑B、D两侧对结构提供摩阻力;④不考虑结构柔性变形对稳定性的影响;⑤不考虑地下水影响。
参考图7所示结构及荷载布置,定义深筒形基坑支护结构抗倾覆稳定安全系数为
Kc=∑MER∑MES=(MEai+MEco+MFao+MFai+MFco+MFci+2Mfbo+ 2Mfbi+MG)⋅(MEao+MEci)−1={Eailai+Ecolco+[Eao(b+e)+Eai(b−2t+e)+Eco(b−e)+ Eci(b−2t−e)]⋅tanδ2+fbi(d−2t)+fbod+Ge}⋅(Eaolao+Ecilci)−1 (1) 式中,
MEai ,MEao ,MEci ,MEco 分别为结构A侧、C侧内部、外部土压力所产生的弯矩;MFai ,MFao ,MFai ,MFco 分别为结构A侧、C侧内部、外部侧摩阻力所产生的弯矩;MFbi ,MFao 分别为结构B侧内部、外部侧摩阻力所产生的弯矩;MG为支护结构自重产生的弯矩。Eai,Eao,Eci,Eco分别为结构A侧、C侧内部、外部土压力;lai,lao,lci,lco分别为合力Eai,Eao,Eci,Eco作用点至支护结构底端的距离;fbi,fbo分别为结构B侧内部、外部侧摩阻力;δ为结构的外摩擦角;b,d,h分别为支护结构外缘的长、宽、高;h0为筒内土芯高度;t为挡土结构厚度;e为转动中心与结构底部中心的偏心距。结构A侧外部土压力Eao按照考虑空间效应和超载作用的无限土体主动土压力计算,土压力及作用点计算式见式(2),(3);结构B侧、D侧外部土压力Ebo,Edo按照考虑空间效应的无限土体静止土压力计算,在结构倾覆破坏时提供摩阻力fbo、fdo,见式(4),摩阻力方向水平向左;结构C侧外部土压力Eco按照考虑空间效应的有限土体被动土压力计算,如图8所示,对于一侧放坡的有限土体来说,滑裂面可能会与P1P2或P2P3或P3P4任一段相交,对应滑楔体的形状可能为三角形、四边形或者多边形,被动土压力大小为式(5)~(7)的最小值,作用点仍取为距离墙底1/3墙高处;结构A侧内部土压力Eai、结构C侧内部土压力Eci按照有限土体被动、主动土压力计算,假定土压力三角形分布[5-6],土压力作用点高度为h0/3;结构B侧、D侧内部土压力Ebi,Edi按照考虑空间效应的有限土体静止土压力计算,结构倾覆破坏时提供摩阻力fbi,fdi。具体计算表达式如下:
Eao=2∫d/20h2(γh+2q)KaKspodz, (2) lao=h36γKa+h22qKah2(γh+2q)Ka=γh2+3qh3(γh+2q), (3) fbo=fdo=Ebotanδ=2∫b/2012γh2K0Kspotanδdx, (4) E1=2∫d/20γh22tanθcos(θ−ϕ)sin(θ−ϕ−δ)Kspodz, (5) E2=2∫d/20γ2(h2tanθ−(htanθ−bt)2tanθ+1tanε)cos(θ−φ)sin(θ−φ−δ)Kspodz, (6) E3=2∫d/20γ2[(h−bb)2tanθ+bb(2bt+bbtanε)]⋅cos(θ−φ)sin(θ−φ−δ)Kspodz (7) Eco=min(E1,E2,E3), (8) Eai=2∫d/2−t0(1−12mp)mpγh20KpKspidz, (9) mp=bh0(45∘+φ2)(当mp>1 时,取mp=1), (10) Eci=2∫d/2−t0(1−12ma)maγh20KaKspidz, (11) ma=bh0cot(45∘−φ2)(当ma>1 时,取ma=1), (12) fbi=fdi=Ebitanδ=2∫b/2−t012γh02K0Kspitanδdx 。 (13) 式中
γ 为砂土重度;q为基坑外离岸侧超载值;Ka,Kp,K0分别为主动、被动、静止土压力系数;θ为滑裂面与墙背的夹角;φ为砂土内摩擦角;ε为坡角;Kspo,Kspi分别为坑外、坑内空间效应影响系数,计算方法如图8,9所示;其他符号意义同前。如图9所示,考虑支护结构A、B、D侧坑外土体无限长,C侧临岸,四侧坑外土压力空间系数按照同一方法计算。作用在结构上坑外的土压力随着距离坑角的距离的减小而减小,假定在角部效应影响范围内的土压力呈抛物线分布。其中角部效应影响范围为
b1=(h−h0)tan(45∘−φ2)。 (14) 坑外空间效应影响系数按照式(15)或者(16)计算。沿x方向:
Kspo={k1√xb1 (x≤b1) k2−k1b−2b1(x−b1)+k1 (b1<x≤b−b1)k2√b−xb1 (b−b1<x≤b) ; (15) 沿z方向:
Kspo={k1√zb1 (z≤b1) k2−k1d−2k1(z−b1)+k1 (b1<z≤d−b1)k2√d−zb1 (d−b1<z≤d) 。 (16) 式中,k1,k2为空间效应影响系数的放大系数,其值的大小与基坑所处位置和加载情况相关。
支护结构内部土芯属于两墙之间的有限土体,假定坑内角部效应影响范围内土压力呈线性分布,如图10所示,其中角部效应的影响范围为
b2=h0tan(45°+φ2)。 (17) 坑内空间效应影响系数为
Kspi={1 z(或x)≥ b2z(或x)b2 z(或x) <b2。 (18) 3.2 算例参数分析
为说明本计算方法的可行性,按本方法编制成深筒形基坑整体抗倾覆稳定计算程序,计算各工况下基坑的抗倾覆稳定安全系数,并进行若干参数分析。
以尺寸长×宽×高×厚为5 m×5 m×15 m×0.2 m的钢板桩挡土结构为例。其基坑内部土芯高度h0为6 m,基坑顶部距离岸坡bt为1 m,坡角
ε 为18°,岸坡深度hb为12 m;取转动中心位置系数ke=0.2;砂土内摩擦角φ为20°,重度γ 为16 kN·m-3,外摩擦角δ为内摩擦角的1/4,超载q为30 kPa/m。以这些数据为基本组,通过改变各因素参数观察其对抗倾覆安全系数和抗滑移安全系数的影响。保持其他参数不变,只改变砂土重度和内摩擦角,计算结果如图11。由图11可知,基坑整体抗倾覆稳定安全系数随土的内摩擦角和重度增大而增大。
图12,13分别为其他参数不变,改变岸坡相关参数得到的岸坡各因素对于结构安全系数的影响。抗倾覆安全系数随岸坡坡深的增大而减小,当岸坡深度大于基坑挡土结构高度时,岸坡深度的变化对两个安全系数不再造成影响。岸坡坡角
ε 越小,安全系数越大;挡土结构距岸坡越远,安全系数越大。超载对安全系数的影响如图14所示。超载越大,结构的抗倾覆稳定安全系数越小。
4. 结论
通过室内模型试验与有限元数值模拟来分析支护结构的整体倾覆破坏模式,并提出了砂土中支护结构整体抗倾覆稳定安全系数计算方法,并通过初步算例分析,得出以下3点结论。
(1)临岸深筒性形基坑支护结构整体性强,空间效应明显,结构在离岸一侧加载过大时倾向于发生整体倾覆破坏。
(2)支护结构整体抗倾覆稳定分析方法考虑了支护结构空间效应、坑内及坑外临岸坡侧有限土压力的影响。
(3)基坑支护结构整体抗倾覆稳定安全系数随砂土重度、内摩擦角的增大而增大;结构距离岸坡越远、坡角越小、坡深越浅,结构越安全;远离岸坡侧坑外超载值越大,结构越不安全。
-
表 1 土层参数
Table 1 Soil parameters
黏聚力/kPa 内摩擦角/ (°) 泊松比 重度/(kN⋅m-3) 弹性模量/MPa 19 23 0.3 18 35 表 2 框架梁、锚杆基本参数
Table 2 Basic parameters of frame beam and anchor bolt
立柱/横梁截面面积/m2 气囊初始气压/MPa 锚杆自由段长度/m 自由段弹性模量/MPa 锚杆锚固段长度/m 锚杆锚固段弹性模量/MPa 0.16 0.128 5 2×105 4 2.8×104 -
[1] 单斌, 熊熊, 郑勇, 等. 2013年芦山地震导致的周边断层应力变化及其与2008年汶川地震的关系[J]. 中国科学: 地球科学, 2013, 43(6): 1002-1009. https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201306008.htm SHAN Bin, XIONG Xiong, ZHENG Yong, et al. Changes of peripheral fault stress caused by Lushan earthquake in 2013 and its relationship with Wenchuan earthquake in 2008[J]. Scientia Sinica (Terrae), 2013, 43(6): 1002-1009. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201306008.htm
[2] 邓起东, 张培震, 冉勇康, 等. 中国活动构造基本特征[J]. 中国科学(D辑: 地球科学), 2002, 32(12): 1020-1030, 1057. https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200212006.htm DENG Qidong, ZHANG Peizhen, RAN Yongkang, et al. Basic characteristics of active structures in China[J]. Science in China, SerD, 2002, 32(12): 1020-1030, 1057. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200212006.htm
[3] 史培军, 张欢. 中国应对巨灾的机制: 汶川地震的经验[J]. 清华大学学报(哲学社会科学版), 2013, 28(3): 96-113, 160. https://www.cnki.com.cn/Article/CJFDTOTAL-QHDZ201303014.htm SHI Peijun, ZHANG Huan. China's response to the catastrophe: the experience of Wenchuan earthquake[J]. Journal of Tsinghua University (Philosophy and Social Sciences), 2013, 28(3): 96-113, 160. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-QHDZ201303014.htm
[4] 徐玖平, 卢毅. 地震救援·恢复·重建的组织技术及理论模式[J]. 系统工程理论与实践, 2011, 31(增刊1): 107-119. https://www.cnki.com.cn/Article/CJFDTOTAL-XTLL2011S1017.htm XU Jiuping, LU Yi. Organizational technology and theoretical pattern of earthquake rescue, recovery and reconstruction[J]. Systems Engineering-Theory & Practice, 2011, 31(S1): 107-119. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XTLL2011S1017.htm
[5] 张磊. 面向地震灾情时序变化的应急救援物资需求动态预测研究[J]. 灾害学, 2018, 33(3): 161-164. doi: 10.3969/j.issn.1000-811X.2018.03.031 ZHANG Lei. Research on dynamic demand prediction of emergency relief materials oriented to the temporal and spatial change of earthquake disaster losses[J]. Journal of Catastrophology, 2018, 33(3): 161-164. (in Chinese) doi: 10.3969/j.issn.1000-811X.2018.03.031
[6] 何琳, 赵应龙. 舰船用高内压气囊隔振器理论与设计[J]. 振动工程学报, 2013, 26(6): 886-894. doi: 10.3969/j.issn.1004-4523.2013.06.011 HE Lin, ZHAO Yinglong. Theory and design of high-pressure and heavy-duty air spring for naval vessels[J]. Journal of Vibration Engineering, 2013, 26(6): 886-894. (in Chinese) doi: 10.3969/j.issn.1004-4523.2013.06.011
[7] 陈静, 闫澍旺, 孙立强, 等. 隧道气囊在外压作用下的变形特性及试验验证[J]. 土木建筑与环境工程, 2018, 40(5): 16-26. https://www.cnki.com.cn/Article/CJFDTOTAL-JIAN201805003.htm CHEN Jing, YAN Shuwang, SUN Liqiang, et al. Model test on deformation characteristics of large diameter airbag in tunnel under external pressure[J]. Journal of Civil, Architectural & Environmental Engineering, 2018, 40(5): 16-26. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JIAN201805003.htm
[8] 邹炳燕. 液压与气压传动[M]. 北京: 中国铁道出版社, 2020. ZOU Bingyan. Hydraulic and Pneumatic Transmission[M]. Beijing: China Railway Publishing House, 2020. (in Chinese)
[9] LI Y W, NAHON M, SHARF I. Airship dynamics modeling: a literature review[J]. Progress in Aerospace Sciences, 2011, 47(3): 217-239. doi: 10.1016/j.paerosci.2010.10.001
[10] 张利国, 张嘉钟, 贾力萍, 等. 空气弹簧的现状及其发展[J]. 振动与冲击, 2007, 26(2): 146-151, 183. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ200702035.htm ZHANG Liguo, ZHANG Jiazhong, JIA Liping, et al. Future and development of air springs[J]. Journal of Vibration and Shock, 2007, 26(2): 146-151, 183. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ200702035.htm
[11] 廖航, 竺梅芳, 雷江利, 等. 大质量航天器气囊着陆缓冲过程研究[J]. 航天返回与遥感, 2020, 41(1): 28-38. https://www.cnki.com.cn/Article/CJFDTOTAL-HFYG202001005.htm LIAO Hang, ZHU Meifang, LEI Jiangli, et al. Airbag landing research of massive spacecraft[J]. Spacecraft Recovery & Remote Sensing, 2020, 41(1): 28-38. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HFYG202001005.htm
[12] 东南大学, 天津大学, 同济大学. 混凝土结构学习指导[M]. 3版. 北京: 中国建筑工业出版社, 2020. Southeast University, Zhejiang University, Hunan University. Learning Guidance of Concrete Structure[M]. 3rd ed. Beijing: China Architecture & Building Press, 2020. (in Chinese)
[13] 李怀志. 预应力锚索格构梁受力分析[D]. 贵阳: 贵州大学, 2009. LI Huaizhi. Stress Analysis of Prestressed Anchor Cable Lattice Beam[D]. Guiyang: Guizhou University, 2009. (in Chinese)
[14] 周勇, 朱彦鹏. 框架预应力锚杆柔性支护结构坡面水平位移影响因素[J]. 岩土工程学报, 2011, 33(3): 470-476. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract13964.shtml ZHOU Yong, ZHU Yanpeng. Influencing factors of horizontal displacement of wall facing of grillage flexible supporting structure with prestressed anchors[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(3): 470-476. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract13964.shtml
[15] SELVADURAI A P S. Elastic analysis of soil-foundation interaction[M]. New York: distributors for the United States and Canada, Elsevier-North Holland, 1979.
[16] 朱彦鹏, 董建华. 柔性支挡结构的静动力稳定性分析[M]. 北京: 科学出版社, 2015. ZHU Yanpeng, DONG Jianhua. Static and Dynamic Stability Analysis of Flexible Retaining Structure[M]. Beijing: Science Press, 2015. (in Chinese)
[17] 梁瑶, 周德培, 赵刚. 预应力锚索框架梁支护结构的设计[J]. 岩石力学与工程学报, 2006, 25(2): 318-322. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200602016.htm LIANG Yao, ZHOU Depei, ZHAO Gang. Design of support of frame beam and prestressed anchor[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(2): 318-322. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200602016.htm
[18] 王其昌. 铁路新型轨下基础应力计算[M]. 北京: 中国铁道出版社, 1987. WANG Qichang. Stress Calculation of Foundation under New Railway Track[M]. Beijing: China Railway Publishing House, 1987. (in Chinese)
[19] 建筑边坡工程技术规范: GB 50330—2013[S]. 北京: 中国建筑工业出版社, 2014. Technical Code for Building Slope Engineering: GB 50330—2013[S]. Beijing: China Architecture & Building Press, 2014. (in Chinese)
[20] 顾晓鲁. 地基与基础[M]. 3版. 北京: 中国建筑工业出版社, 2003. GU Xiaolu. Foundation and Foundation[M]. 3rd ed. Beijing: China Architecture & Building Press, 2003. (in Chinese)
[21] POLI︠A︡NIN A D, ZAĬT︠S︡EV V F. Handbook of Exact Solutions for Ordinary Differential Equations[M]. 2nd ed. Boca Raton: Chapman & Hall/CRC, 2002.
[22] 张兴龙. 纤维增强复合材料结构宏细观多尺度力学性能研究及应用[D]. 西安: 西安电子科技大学, 2015. ZHANG Xinglong. The Macro and Micro Mechanical Properties and Application of the Fiber Reinforced Composite[D]. Xi'an: Xidian University, 2015. (in Chinese)
-
其他相关附件