拱坝计算书

1计算书目录：1、设计参数及控制指标2、拱坝体形3、应力计算4、拱肩稳定计算5、消能计算6、坝体细部及放空、取水孔设计1、设计参数及控制指标1.1设计参数坝体材料：200#砼，容重2.4t/m3,弹模1.7E6(坝体弹模考虑徐变的影响，取为瞬时弹性模量的0.6--0.7)，泊松比0.167，线胀系数1×10-5/℃，导温系数3m2/月。坝基：灰岩，容重3t/m3,弹模2E6，泊松比0.27，线胀系数1.4×10-5/℃，导温系数3m2/月。淤沙浮容重按1t/m3,内摩擦角14°。水文及地基f、c等有关各专业的基础资料请见附件1。温度荷载按规范(SD145-85)附录公式由程序动计算，封拱灌浆温度取8-12℃。1.2控制指标大坝拱肩稳定及应力控制指标均严格按照《混凝土拱坝设计规范》（SD145-85）执行，见表1-1、1-2。表1-1抗滑稳定安全系数表类别荷载组合情况安全系数抗剪断强度K’基本组合3特殊组合2.5表1-2大坝允许应力表类别荷载组合情况允许应力允许拉应力(kg/cm2)基本组合12特殊组合15允许压应力(kg/cm2)基本组合50特殊组合552、拱坝体形拱坝体形为双曲拱坝，拱圈平面曲线采用圆弧。因两岸地形不完全对称而采用两岸不同半径的双曲拱坝。2.1坝顶高程的拟定2.1.1已知:校核洪水位(p=0.2%)：848.87m2设计洪水位(p=2%)：848.35m正常蓄水位：848m2.1.2坝顶高程根据各种运行情况的水库静水位加上相应超高后的最大值确定。顶超高值Δh按下式计算(请见SD145-85《混凝土拱坝设计规范》第八章拱坝构造）Δh=2hl+h0+hc式中：Δh………坝顶距水库静水位高度(m)2hl………浪高(m)h0………波浪中心线至水库静水位的高度(m)hc………安全超高(m)：正常运用情况取0.4m，非常运用情况取0.3m。2.1.3波浪要素按“官厅——鹤地”公式计算：2hl=0.0166Vf5/4Df1/32L1=10.4(2hl)0.8h0=4πhl2/(2L1)式中：2L1………波长（m）；Df………吹程，由坝前沿水面至对岸的最大直线距离(km),取1Km。Vf………计算风速(m/s)，正常和设计情况取15m/s，校核情况取9.6m/s。2.1.4坝顶高程计算成果见下表：计算情况设计情况校核情况相应水位（m）848.35848.87Vf（m）159.6Df（km）11hc（m）0.30.42hl（m）0.4900.281h0（m）0.1280.066Δh（m）0.9180.747坝顶高程（m）849.268849.617结论：综合考虑后，取坝顶高程为849.7m河床底部高程为782.3m，按地质提供的基本资料，挖深按20.5m计，则坝高为87.4m，属于高坝，在坝顶高程849.7m开挖后的河谷水平宽度为193.63m，则河谷宽高比L/H=193.63/87.4=2.215，属“V”型河谷。2.2、拱坝体型设计2.2.1坝顶厚度Tc（参见《砌石坝设计》P137）3Tc=0.4+0.01(L+3H)式中Tc——坝顶厚度（m）H——最大坝高（m）=87.4L——顶拱弦长（m）=193.63经计算得Tc=0.4+0.01(193.63+3×87.4)=4.9583m，考虑交通等要求，取Tc=5m。2.2.2拱圈厚度（1）坝底厚Tb①参见《砌石坝设计》P138经验公式Tb=(0.132(L/H)0.269+(2H/1000))×H式中Tb——坝底厚度（m）H——最大坝高（m）=87.4L——顶拱弦长（m）=193.63经计算得Tb=(0.132(193.63/87.4)0.269+(2×87.4/1000))×87.4=29.57m。②由美国肯务局的经验公式式中Tb——坝底厚度（m）H——最大坝高（m）=87.4L1——顶拱弦长（m）=193.63L2——顶底以上0.15H处弦长（m）=42.2经计算得Tb=8.77m。综合考虑后，取坝底厚度为Tb=22m。水平拱圈厚度按T=Tc+(Tb–Tc)Y/H计算，即为T=5+(22–5)Y/85.9=5+0.1979Y2.2.3上游面曲线初拟参考其他工程经验及规范要求进行设计：设β1=Hd/Hβ2=Aud/HK=Aub/Tb一般情况下有：β1=0.3~0.4β2=0.135~0.175K=0.3~0.8因本工程河谷宽高比L/H=193.63/87.4=2.2154，综合考虑后取参数计算如下：122/213)122/(0012.0HHLHLTb4Hd=0.38H=0.37×85.9=31.783mAud=0.175H=0.16×85.9=15.0325mAub=0.7×20=14mm=Aud-Aub最终的成果取值为：(请见右图)Hd=31.8mAud=14.944mAub=13.5mm=14.944–13.5=1.444m由此可知，A、D、B三点的坐标为：点号XYA00D-14.94454.1B-13.585.9将地形图上可利用的基岩线确定后，再确定每一层拱圈对应的水平弦长，定出半径后（用半中心角控制，如顶拱的半中心角在80°～110°之间，底拱的半中心角在40°～80°之间，半径的连线为一光滑的曲线）经过应力控制的优化后得到拱坝体型。计算采用的程序为ADASO。设计过程中严格按照有关规范的要求来进行大坝体形的设计调整，最后确定大坝体形如下(请见下页)：101.392.3393.11010849.7843拱冠梁782.375.6557.7873.6745.8756.873440.122.0625.04827.2811.4795.6779.8763.853、应力计算3.1计算方法及荷载组合拱坝应力按多拱多梁混合法计算，利用北京水科院ADASO程序在PentiumⅡ微机上计算。计算初拟采用7层拱圈，其高程分别为：849.7、843.0、827.2、811.4、795.6、779.8、763.8m(拟合层为砼垫层762.3m)。荷载组合情况分基本组合及特殊组合两类。基本组合为：⑴正常蓄水位+相应尾水位+设计正常温降+自重+扬压力+泥沙压力+浪压力。特殊组合为：⑵校核洪水位+相应尾水位+设计正常温升+自重+扬压力+泥沙压力+浪压力。地区地震烈度小于Ⅵ度，故不作动力计算。3.2温度参数说明多年年平均气温14.6℃多年平均最低月(1月)气温3.8℃多年平均最高月(7月)气温24.7℃日照影响按2℃计温降变幅=年平均气温-最低月平均气温=14.6-3.8=10.8℃考虑日照影响后，设计正常温降=10.8+2=12.8℃，计算时间取1.5温升变幅=最高月平均气温-年平均气温=24.7-14.6=10.1℃考虑日照影响后，设计正常温升=10.1+2=12.1℃，计算时间取7.5库底水温因无实测资料，按规范(SD145-85)附录规定，Tkd可近似按最低3个月平均气温计，Tkd=(3.8+5.0+5.9)/3=4.9℃水表面年平均水温=年平均气温+日照影响=14.6+2=16.6℃表面水温年变幅=气温年变幅的一半=(24.7-3.8)/2=10.45℃温度计算的初相位取值为6.5（月）。3.3计算成果经计算，拱坝最大拉应力为1.42MPa，发生在▽795.6m高程右拱端的上游面，计算工况为校核洪水位+温升；最大压应力为4.03MPa，发生在▽795.6m高程右拱端的下游面，计算工况为校核洪水位+温升，均小于大纲中规定的允许应力。应力分布规律符合一般规律，较合理，计算原始数据和结果文件见附。64拱肩稳定计算4.1计算方法拱肩抗滑稳定按刚体极限平衡法计算，并根据《混凝土拱坝设计规范》规定(SD145-85)，拱肩稳定按平面分层抗剪断公式核算，若分层安全系数不满足要求时则按空间问题计算整体抗滑稳定安全系数，若均满足要求则不再按空间问题计算整体抗滑稳定安全系数。计算步骤：⑴绘制各层拱肩抗滑岩体图。⑵根据地质专业提供的资料，列出几组不同侧滑面（裂隙或岩体），在抗滑岩体图中分别量出侧滑面与拱端径向的夹角、底滑面面积、侧滑面长度。⑶根据前面量出的参数，各层拱端高程，计算工况的水位，各侧滑面、底滑面分别对应的抗剪断参数，程序计算得到的各层拱端作用力，代入抗剪断公式，分别计算各层、各侧滑面的K值，取其中每层最小K值为该层安全系数K。4.2已知条件4.2.1由拱坝应力计算得拱端力系计算结果如下表：拱端力系表（单位：t）工况拱端高程849.7843.0827.2811.4795.6779.8763.8正常蓄水位+温降左拱端推力-428.7383.11290.11746.72153.12082.4934.8拱端剪力11.6-25.9184.3480.1860.41494.92117.0右拱端推力-138.0531.11119.31642.82157.32231.61084.2拱端剪力-6.4293.8907.31645.22368.82414.71243.5校核洪水位+温升左拱端推力813.31137.81304.41699.22116.42007.1876.2拱端剪力-33.7-395.5-106.5351.2751.91347.51915.2右拱端推力1653.21591.71105.71582.32131.62038.5835.2拱端剪力1581.51480.21158.01595.92270.12167.8959.84.2.2水文：校核洪水位848.87m；4.2.3地质：①河谷呈“V”型，出露地层为T1m2薄至中厚层灰岩。②坝址岩层产状平缓，层面（岩层产状）为：左岸N5°E/∠6°；右岸N20°W/∠10°。主要发育裂隙有：E0°W/∠90°；S0°N/∠90°；断层主要有：f101：E0°W/N∠20~30°、破碎带宽0.3~0.5m；f102：N50°W/SW∠20~35°、破碎带宽0.3~0.5m；4.3计算假定根据地质报告，坝肩两岸岩层产状倾角较缓，陡倾角裂隙较发育，无其它大的结构面发育，拱肩稳定按平面分层抗滑计算，为简化计算，在计算中作如下假定：⑴地基为灰岩，因其产状倾角较缓（10°～20°），故可近似假定岩体层面为水平底滑动面。7⑵几组裂隙均为陡倾角（60°～80°），故可假定裂隙面即为可能的侧滑面。⑶若在拱端处岩体无裂隙面或裂隙面对抗滑无不利影响时，则可根据地形条件假定几组侧滑面，按岩体非结构面参数取用。⑷底滑面扬压力按单位高度计，取1m。⑸侧滑面渗透压力按实际水头计，即侧滑面上游端水头为实际H，下游面无水可视为0，侧滑面下游端水头为0。但考虑到侧滑面较长时，下游岸坡裂隙发育，渗水一般在坝后很快渗出地表，即零压力点不在侧滑面下游端与岸坡交点处。近似按拱下游面后一定距离为零计，如5倍拱厚、1倍岸坡厚等，为方便计算，计算中不直接确定零压力点而采用取综合折减系数的方法，综合折减系数取0.5计算渗透压力。认为U3=0。⑹为方便计算，认为裂隙面与拱端面为同一平面。⑺缓倾面认定为层面，作用力为1/3面积。⑻当剪力向上游时，取其值为0。⑼为偏于安全，未计梁底压力。4.4计算模型及公式两坝肩地基特点是层面倾角较缓、裂隙发育且倾角较陡，计算分层稳定时视层面为底滑面，裂隙面为侧滑面。计算参数请见附图。计算时按抗剪断K值公式计算，如下：K=[(N′-U1)f2+(Ws+W-U2)f1+C1Lab]/Q′式中各符号意义如下：N′:拱端作用力在侧裂面上的法向分力，Q′:拱端作用力在侧裂面上的切向分力。Ws:梁自重。W:抗滑岩体自重。U1:侧裂面渗透压力。U2:底滑面扬压力。f2:侧滑面抗剪系数。f1、C1:底滑面抗剪断系数。Lab:底滑面有效面积。为提高计算速度，编制了gbwd.bas程序，按上述公式计算各分层K值。该程序需输入的参数有拱端作用力N、Q，底滑面面积S，侧裂面与拱端径向夹角α及侧滑面长度L，以及底滑面、侧裂面的抗剪、抗剪断参数，各层拱高程及计算工况上游水位。为偏于安全计算，近似底滑面有效面积Lab等于底滑面面积S的1/3。U1按前假定规定折减。N,Q由ADASO程序计算结果而得。