Adao

拱坝分析与优化软件系统ADAO （ADCAS ＆ ADOPT） 使 用 说 明 书 浙江大学水工结构研究所 刘国华 主编 2008年1月 目 录 4§1 ADAO拱坝分析与优化软件系统概要 6§2 ADAO 输入数据与软件使用的总体说明 14§3　应力分析输入数据文件.CAS填写说明 14§3.1　标 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 行 14§3.2　主控行 18§3.3 分析计算中所用拱梁网格的拱圈高程行。 19§3.4 虚结点位置描述行 19§3.5 基岩材料性质描述主参数行 20§3.6 基岩材料性质描述附加参数行 21§3.7 坝体材料性质描述主参数行 22§3.8 坝体材料性质描述附加参数行 25§3.9 响应谱法地震应力分析所需的场地特性与地震强度参数 25§3.10 用于描述体形的拱圈层数行 26§3.11 拱圈参数描述行 26§3.11.1 用上下游面圆弧描述的单心等厚圆拱圈参数行(IAD=-1) 26§3.11.2 用上下游面圆弧描述的五心变厚圆拱圈参数行(IAD=-2) 27§3.11.3 用中心线圆弧描述的单心等厚圆拱圈参数行(IAD=1) 28§3.11.4 用中心线圆弧描述的五心变厚拱圈参数行(IAD=2) 29§3.11.5 用中心线描述的抛物线型拱圈参数行(IAD=3) 30§3.11.6 用中心线描述的对数螺旋线拱圈参数行(IAD=4) 31§3.11.7 用中心线描述的悬链线型拱圈参数行(IAD=5) 32§3.11.8 用中心线描述的椭圆线型拱圈参数行(IAD=6) 33§3.11.9 用中心线描述的混合线型拱圈参数行(IAD=7) 35§3.11.10 非直立拱冠梁剖面的拱冠位置与方向角的描述数据 36§3.12 坝体扬压力径向分布折点位置描述行 36§3.13 气温、水温资料描述行 38§3.14 封拱次数与封拱温度控制行 39§3.15　封拱温度数据行 40§3.16 封拱条件行 41§3.17 总工况数及地震烈度 41§3.18 各工况荷载数据行 42§3.19 裂缝分布描述控制行（针对有已知裂缝的已建拱坝） 42§3.19.1 坝体裂缝总体情况描述行 43§3.19.2 上游坝面的特殊裂缝描述 44§3.19.3 下游坝面的特殊裂缝描述 46§3.19.4 特殊类型裂缝的拱区间的描述 47§3.19.5 特殊类型裂缝的梁区间的描述 47§3.20 超载计算条件描述行 49§4 拱坝优化输入数据文件.OPT填写说明 49§4.1 主控行 51§4.2 约束指标行之一（应力指标） 51§4.3 约束指标行之二（倒悬度、拱厚比、曲率半径变化率指标） 52§4.4 约束指标行之三（拱冠梁厚度下限） 52§4.5 约束指标行之四（拱冠梁厚度上限） 53§4.6 约束指标行之五（拱端厚度下限） 53§4.7 约束指标行之六（拱端厚度上限） 54§4.8 约束指标行之七（左半中心角上限） 54§4.9 约束指标行之八（右半中心角上限） 54§4.10 约束指标行之九（与左岸基岩等高线夹角之下限） 55§4.11 约束指标行之十（与右岸基岩等高线夹角之下限） 55§4.12 约束指标行之十一（坝厚指数限值，坝体方量期望值） 56§4.13 约束指标行之十二（左岸拱端顺河向期望位置） 56§4.14 约束指标行之十三（右岸拱端顺河向期望位置） 56§4.15 约束指标行之十四（左岸拱端顺河向位置允许偏移量） 57§4.16 约束指标行之十五（右岸拱端顺河向位置允许偏移量） 57§4.17 左岸坝肩可利用基岩等高线描述行 57§4.17.1 棱形河谷左岸坝肩可利用基岩等高线描述行 58§4.17.2 非棱形河谷左岸坝肩可利用基岩等高线描述行 59§4.18 右岸坝肩可利用基岩等高线描述行 59§4.18.1 棱形河谷右岸坝肩可利用基岩等高线描述行 59§4.18.2 非棱形河谷右岸坝肩可利用基岩等高线描述行 60§4.19 约束指标行之十六（封拱温度下限） 61§4.20 约束指标行之十七（封拱温度上限） 61§4.21 约束指标行之十八（拱冠处拱圈曲率半径之下限） 62§5 坝肩稳定分析输入数据文件.STB填写说明 62§5.1 左坝肩稳定分析主控行 63§5.2 左岸坝肩各层拱圈的拱端下游侧端点有效嵌深 64§5.3 左岸坝肩有效岩体等高线高程描述行 64§5.4 用于表示左岸坝肩有效岩体等高线的折线的折点数 65§5.5 用于表示左岸坝肩有效岩体等高线的折线的折点数 65§5.6 右坝肩稳定分析输入数据 §1 ADAO拱坝分析与优化软件系统概要 　　ADAO软件系统是拱坝应力分析ADCAS和拱坝优化ADOPT的集成系统，既可用于拱坝应力分析，也可用于拱坝优化 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 ，以及计算机辅助下的手工体形调整，已应用于许多拱坝工程的设计中，有助于提高设计进度和设计质量，具有显著的社会、经济效益。95年被列入国家级科技成果重点推广 计划 项目进度计划表范例计划下载计划下载计划下载课程教学计划下载 ，98年荣获国家科技进步奖。 　　软件系统还提供了配套的前处理模块，用于交互地输入拱坝设计和优化的基本参数；以及一个后处理模块，用于观看分析结果。 　　ADCAS（Arch Dam Complete Adjustment Stress Analysis System ）应用多拱梁分载法的基本原理编制而成，它集三向调整法、全调整法等多种方法于一体，可选择下列五种不同方法进行计算： 　　三向调整法 进行拱梁径向、切向线变位与水平扭转（绕竖向）变位的调整，仅考虑与上述三个变位分量相对应的三个分载分量的影响。 　　三向调整＋2M法 在三向调整法的基础上，将梁水平扭转分载的２倍近似地作为梁的竖向弯曲（绕切向）分载施加于梁上。 　　四向调整法 在三向调整法的基础上，利用剪应力互等关系推求梁的竖向弯曲分载并施加于梁上。 　　五向调整法 在三向调整法的基础上，利用剪应力互等关系推求梁的竖向分载和竖向弯曲分载并施加于梁上，同时计入坝体竖向变位沿切向的变化率对梁切向变位计算的影响。相当于美国垦务局所采用的考虑竖向变位影响的全调整法。 　　全调整法 在上述五向调整法的基础上，以径向纤维直线假定取代传统的平截面假定来进行拱梁变形分析，进一步考虑绕径向弯曲分载与绕径向扭曲变位对梁变位计算的影响，能满足所有的拱梁变位协调条件、径向条形微体的平衡条件与坝体的边界条件，是传统试载法的提高与完善，研制者向用户推荐使用本方法。 　　ADOPT（Arch Dam Optimization System）采用优秀的约束变尺度法进行拱坝优化设计，在一维搜索中运用监控（Watch dog ）技术等智能方法选择新迭代点，用临时扩展可行域以及其它方法处理优化迭代中可能出现的失效情况，提高了优化算法的有效性和计算效率。应用强有力的应力分析系统ADCAS生成全面、可靠的应力约束条件；在几何约束方面，除一般考虑的坝厚、中心角、倒悬度等约束之外，还制定了一整套严格的常规体形要求的约束条件，如坝面凸性约束，拱圈厚度变化约束，和梁厚非颈缩约束等。 ADAO 概要 ──试载法模式下的各种分析方法齐全 　　·三向调整法 　　·三向调整＋2M法 　　·四向调整法 　　·五向调整法 ★全调整法 ──静动力、非线性分析功能齐备 　　·常规静力分析 　　·常规静力分析 + 拟静力法地震应力分析 　　★响应谱法地震应力分析 　　★分期浇筑、分期封拱、分期蓄水条件下的全过程应力分析 ★静力非线性分析（含超载条件下的非线性全过程分析） ──拱型曲线库完备 　　·单心、双心、三心等厚、变厚圆拱圈 　　·抛物线拱圈 　　·对数螺旋线拱圈 　　·悬链线拱圈 ·椭圆拱圈 ★新型的混合线型拱圈（圆锥曲线族、对数螺旋线和悬链线的混合形式） ──优化方法先进有效 ·高效的约束变尺度法 ·一维搜索的监控技术 ·约束条件全面合理 ★用可性域扩展技术和智能选点方法处理失效的迭代过程 ──体形优化可用于棱形河谷，也可适用于非棱形河谷 ──坝顶开口条件下的应力分析 ──最大网格密度：十拱二十三梁(根据用户需要可以增加到十五拱以上) ──可考虑拱梁之间泊桑比效应的影响 ──可考虑坝体扬压力的影响 ──能够按照混凝土拱坝设计 规范 编程规范下载gsp规范下载钢格栅规范下载警徽规范下载建设厅规范下载 （1985规范、2003规范）所推荐的公式计算温度荷载 ──可以对坝体和基础的不同区域指定不同的材料性质 ──输入数据量少，坝基的Vogt系数等均由软件自动计算 ──便利的后处理功能 ·分析结果图形化、汉化、表格化 ·可显示并绘制结构图、主应力图、等应力线图等 ★可输出详尽的施工放样点坐标 ──适用机型及操作系统：奔腾(pentium)等微机，Windows操作系统 显示器宜设置为像素800x600，16位增强彩色。 §2 ADAO 输入数据与软件使用的总体说明 　　ADAO 在运行时由对话框的选项或输入数据，控制、选择软件的执行功能，既可用于拱坝应力分析，也可用于拱坝优化设计，以及计算机辅助下的手工体形调整。 　　ADAO推出了两种不同形式的版本，其一是一次性转让版本，用户需以相对较高的价格才能获得该版本的使用权，而今后在使用过程中则不再受到限制；其二是付费使用版本，用户可用相当低的价格获得该版本的使用权，而相应地，在今后的使用过程中则会受到一定的限制，即所应用的每一个拱坝工程均需向软件研制者通报，索取相应的使用口令，并交纳软件使用费。 　　在运行软件的拱坝应力分析功能时，需要读入一个数据文件，用于描述拱坝体形、材料特性、荷载条件和计算要求等。该文件的文件名由用户任取，而文件的扩展名则由软件系统规定为.CAS。 　　在运行软件的拱坝优化功能时，也需要以同样的格式读入.CAS文件，其中所描述的体形为优化的初始体形；另外还需读入一个文件名与.CAS相同、规定扩展名为.OPT的文件，用于选定优化的体形模式、约束指标和坝肩地形条件等。 系统提供了EXAMPLE.CAS和EXAMPLE.OPT两个数据文件，供用户填写.CAS和.OPT文件时参考。 　　程序系统的应力分析结果存放于.STS文件中，应力等图形的AutoCAD 接口文件存放于.DXF文件中，优化 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 新体形的设计变量存放于.VAR文件中，与优化体形方案相应的.CAS文件以及优化体形的轴线位置存放于.NEW文件中，优化的迭代过程存放于.ITR文件中，输入数据的复核检查文件存放于.CHK文件中。这些输出文件的文件名均与输入数据文件名相同，扩展名由软件系统取定。对于付费使用的版本，系统会自动地将用户所使用过的口令记录于文件ADAO.PWD中，以确保用户对每个应用工程仅需回答一次使用口令。 在.CAS和.OPT数据文件中，程序系统规定凡是第一列为＃的行均被视为注释行，在数据输入时被略去。 用户在使用程序系统之前，必须先根据工程条件准备好.CAS文件（用于优化时还须.OPT文件）。 ADCAS ＆ ADOPT各种数据文件的含义总结如下： .CAS──应力分析的输入文件 .STB──坝肩稳定分析的输入文件 .OPT──优化的附加输入文件 .STS──应力分析结果文件 #分层小块体法坝肩稳定分析_左岸.TXT──分层小块体法左坝肩稳定分析成果 #分层小块体法坝肩稳定分析_左岸.TXT──分层小块体法右坝肩稳定分析成果 #放样坐标.XYZ──坝体（横缝）施工放样坐标 .DXF──应力等图形的AutoCAD接口文件 .VAR──优化体形设计变量文件 #优化.CAS──相应于优化体形方案的.CAS文件 .ITR──优化迭代过程文件 .ECH──输入数据的复核检查文件 （以上文件均为文本文件） ADAO.PWD──记录已被系统接受的用户使用口令（仅限于付费使用的用户） 程序系统规定的坐标系如下：整体坐标系的X轴沿水平面指向右岸，Y轴铅直向下，Z轴沿水平面指向下游。应力分析的输入／输出数据所取用的整体坐标系以拱冠梁剖面作为其OYZ坐标平面。因OXYZ坐标系随坝轴线的变动而变动，故在优化中另取坐标系O'X'Y'Z'，O'X'Y'Z'各坐标轴的基本朝向同上，此坐标系用于描述坝肩地形和确定坝体位置，应力分析的输出结果数据大多由局部坐标系oxyz所描述，x 轴沿坝体局部点上的拱圈中心线切向指向右岸，y轴铅直向下，z轴沿坝体局部点上的拱圈中心线法向指向下游。 软件打开运行时，以十分醒目的方式，郑重提醒用户： 软件并不能替代设计人员的经验， 软件只能起到辅助的作用， 使用者必须要有自己的设计经验和理性判断： 包括——对设计规范的全面理解； ——对设计参数的合理取值； ——对坝基嵌深的恰当选择； ——对合理体型的正确把握，等等…… 软件开始运行时，请先在“文件”菜单中选择工作目录及文件名，然后在“功能选择”菜单中选择所需的功能： “功能选择”菜单的子菜单： “设计方案”、“优化方案”分析子菜单的对话框： “体形优化设计”子菜单的对话框： 软件数据输入的对话框如下，对话框具有动态帮助的功能。对话框中的“输入项”代表一组输入变量。在输入时动态地显示输入项及输入变量的有关说明及建议值等。对话框中还有“本页输入：详尽→简化”或“本页输入：简化→详尽”的按钮，可以在“简化输入”与“详尽输入”之间进行切换。在“简化 输入”情况下，软件自动地给出某些建议值，同时相应的输入框被置为灰色。 软件在运行应力分析功能时，将首先弹出对话框提问输出量的级别（建议值为1）。输出量的级别的含义如下： =1（希望输出简略时的选择）。此时屏幕上仅有少量输出信息，输出文件.STS的内容比较简洁，仅给出坝体集合信息及主应力成果； =2（希望输出较详细时的选择）。此时屏幕上仅有少量输出信息，输出文件.STS中的内容齐全，包括坝体几何信息、坝面主应力、坝面拱梁应力、结点温度荷载、拱梁内力、位移、拱梁分载等； =3（希望输出详尽图形或作输入数据检查时的选择）： 输出文件.STS中的内容齐全，屏幕上可显示坝体剖面图、应力图和位移图等，并可形成AutoCAD的图形接口文件.DXF， 系统会有若干提问以控制或拷贝图形； =4（输入数据检查时的选择）： 有时输入数据出错导致软件运行中断而未能发现错在何处，此选择可以显示系统的一些计算过程，可帮助用户找到错误之所在，以及供研制者调试程序时使用，系统将会有若干次暂停，以便使用者有足够的时间观看屏幕上的内容。 对于体形优化和计算机辅助下的手工体形调整，系统不提问输出量的级别而自动取为1。 在完成了上述的交互输入之后，对于付费使用的版本，系统将从输入数据中抽取工程的特征参数，将它与文件ADAO.PWD中已有的口令所对应的特征参数进行匹配。若不能匹配，则系统向屏幕输出该工程的特征参数，在这种情况下，用户需要向系统研制者通报这些特征参数，并索取相应的使用口令以回答系统的提问，否则系统将无法正常运行。 作应力分析时，坝面上最大最小主应力结果的摘要将显示于屏幕上，其中封拱阶段0代表第一次封拱时的初始状态；方括号中的U代表最大主应力出现在上游坝面，D代表最大主应力出现在下游坝面，(R代表最大主应力出现在第(层拱圈，拱圈序号以顶拱为1向下编排，(C代表最大主应力出现在第(条梁，梁的序号以拱冠梁为0，左岸侧为负、右岸侧为正依次编排。 作体形优化时，在开始时将出现如下的三个对话框，给出优化初始方案的设计变量值，并允许对这些设计变量值进行变更，其中选择框“√”是否选上，代表相应的项是作为设计变量还是作为设计参数（设计参数在优化过程中取为定值）。这些设计变量或设计参数决定了沿高程作线性、二次或三次变化的体型参数。必要时，可改选B样条来拟合各高程上的设计变量或设计参数，此时，沿高程所选取的设计变量或设计参数的个数超过4个。 下图为优化迭代时屏幕显示的界面，其中“当前迭代超限应力约束概要”中给出超限应力的位置及大小：第一个数代表工况号（0代表自重应力）；“U”、“D”和“I”分别代表“上游坝面”、“下游坝面”和“内结点”（即非基础结点），“@”所在行给出独立坝段抗倾覆的安全系数；第二个数代表梁号；第三个数代表拱号；第四个数代表应力值。 屏幕右上方给出迭代序列的概要，其中最紧约束值代表在当前迭代点上违反最严重的约束值，且这些约束值按约束的最大最小允许值作归一化。在“窗口”菜单中，提供了“约束向量”和“分析优化模块”两种窗口可供切换，下图所示的是“分析优化模块”窗口。在“约束向量”窗口中给出了当前排序靠前的一些约束值及相应的约束名称与位置（约束名称后面的数字，2位数时表示拱号，3位数时表示结点号，结点号从左岸起自上而下、自左而右依次编排）。约束向量的次序除了第一个固定为“坝体方量约束”外，其余是根据初始方案的约束违反程度排列的，迭代过程中依据一定的准则对约束向量重新排序。 §3　应力分析输入数据文件.CAS填写说明 §3.1　标题行 　　本行的第1至第80列被读入作为该计算项目的标题，中西文均可。 分析结果的报告文件中将包含此标题，以利于识别。 §3.2　主控行 主控行中的各主控参数将控制着程序的运行功能和后续输入数据项的取舍。 【数据格式】 主控行中包含14个整数： NARCH，NCANTI，MODEL，IAD，ITH，NVIRT，ITEM，IPSN，IWP，IDFBL，IUNIT，INON，IDYN，IOPT 格式说明：自由格式。即规定本输入项可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 1、 NARCH： ·分析计算中所选用的拱梁网格的拱圈数目（不含拱冠梁底部 高程上的拱圈），≤15。 2、NCANTI： ·可以填1、3、5或7，整数。建议值：3。 · =1、3、5或7时，各边梁是从各条拱的拱端向上剖切而来的，其拱梁网格的悬 臂梁数目为2*(NARCH-1)+NCANTI。 · =1 时规定在河床上只设置一条梁，即拱冠梁。 · =3、5或7时，规定在河床上设置3、5或7条梁，中间一条为拱冠梁。 3、MODEL： ·可填1～5，用以选择分析方法。 =1 三向调整法。 =2 三向调整＋2M法。 =3 四向调整法。 =4 五向调整法。 =5 全调整法。 · 建议值为5。 4、IAD： ·可填-2～7，用以选择拱圈线型。 =－1 用上下游面圆弧描述的单心等厚圆拱圈。 =－2 用上下游面圆弧描述的五心变厚圆拱圈。 =1 用中心线圆弧描述的单心等厚圆拱圈。 =2 用中心线圆弧描述的五心变厚圆拱圈。 =3 用中心线描述的抛物线型拱圈。 =4 用中心线描述的对数螺旋线型拱圈。 =5 用中心线描述的悬链线型拱圈。 =6 用中心线描述的椭圆线型拱圈。 =7 用中心线描述的混合线型拱圈。 5、ITH： ·可填1～4，11～14，用以选择拱圈厚度变化规律。仅当IAD≥3时有效。            （ITH=11～14选项可使岸坡不光滑时的拱厚变化更为光滑平顺） =1 用公式Ti=Tc+(Ta-Tc)*(1-cosφi)/(1-cosφa)计算半中心角φi处的拱厚，其中Tc为拱冠梁处的拱厚，Ta为拱端的拱厚，φa为拱端的半中心角。 =2 用公式Ti=Tc+(Ta-Tc)(Si/Sa)1.75计算中心线弧长为Si（从拱冠梁起算）处的拱厚，其中Sa为拱端处的中心线弧长（从拱冠梁起算），其余同上。 =3 用公式Ti=Tc+(Ta-Tc)(φi/φa)Tindex计算拱厚，指数Tindex在拱圈参数描述行中提供。 =4 用公式Ti=Tc+(Ta-Tc)(Si/Sa)Tindex计算拱厚，指数Tindex在拱圈参数描述行中提供。         =11用公式Ti=Tc+(Te-Tc)*(1-cosφi)/(1-cosφe)计算半中心角φi处的拱厚，其中Tc为拱冠梁处的拱厚，Te为“似拱端”的拱厚，φe为“似拱端”的半中心角，其中“似拱端”非真正之拱端也， “似拱端”的位置由“似拱端”半中心角φe确定，φe由顶拱和底拱的拱端半中心角按高程线性内插而得；         =12用公式Ti=Tc+(Te-Tc)(Si/Se)**1.75计算中心线弧长为Si（从拱冠梁起算）处的拱厚，其中Se为“似拱端”处的中心线弧长（从拱冠梁起算），其余同上；         =13 用公式Ti=Tc+(Te-Tc)(φi/φe)**Tindex计算拱厚，指数Tindex在拱圈参数描述行中提供；         =14 用公式Ti=Tc+(Te-Tc)(Si/Se)**Tindex计算拱厚，指数Tindex在拱圈参数描述行中提供； 6、NVIRT： ·虚结点数，用于描述坝顶开口。 本程序系统用虚拱单元和虚梁单元（即刚度为零的拱、 梁单元）来模拟坝顶开口。 设置虚拱单元和虚梁单元实质上使得拱与梁在同虚拱与虚梁单元相连结的地方被切断。虚梁单元的长度应同开口的高度相同，虚拱单元的长度也宜尽量同开口的切向长度相同（如图1所示）。虚拱与虚梁单元相交叉的点称之为虚结点，如果用拱号（以顶拱为1号拱，向下顺序编号）和梁号（以拱冠梁为0号梁，以向左岸方向递减、向右岸方向递增的顺序编号）两个数来表示一个结点，则图1中有3 个虚结点，可分别表示为（1，-1）、（1，0）、（1，1）。从图中可以看出， 每个与虚结点相连接的拱单元和梁单元就是虚拱和虚梁单元。 本程序系统通过指定虚结点便可以方便地描述坝顶开口。在此主控行中读入虚结点数NVIRT，虚结点的拱梁号在后续的行片中读入。 ·建议值为0，表示没有或不考虑坝体开口。 　梁号：-8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 拱 号 1 ─┬─┬─┬─┬─┬─┐┈┈┈┈┈┈┈┌─┬─┬─┬─┬─┬─ ： │ │ │ │ │ │ ┊ ┊ ┊ │　│ │ │ │ │ 2 └─┼─┼─┼─┼─┼─┬─┬─┬─┼─┼─┼─┼─┼─┘ │ │ │ │ │ │ │ │ │ │ │ │ │ 3 └─┼─┼─┼─┼─┼─┼─┼─┼─┼─┼─┼─┘ │　│ │ │ │ │ │ │ │ │ │ 4 └─┼─┼─┼─┼─┼─┼─┼─┼─┼─┘ │ │ │ │ │ │ │ │ │ 5 （左岸） └─┼─┼─┼─┼─┼─┼─┼─┘　　 　（右岸） │ │ │ │ │ │ │ 6 └─┼─┼─┼─┼─┼─┘ │ │ │ │ │ 7 └─┼─┼─┼─┘ │ │ │ 8 └─┴─┘ 图1 坝顶与坝身开口示意图 7、ITEM： ·选择温度荷载的计算方法。 =1 按混凝土拱坝设计规范（1985年版）中的公式计算，有关参数在后续的数据项中读入。 =2 按混凝土拱坝设计规范（2003年版）中的公式计算，有关参数在后续的数据项中读入，其中坝体上下游水温按规范建议的公式由气温资料自动推求； =-2 按混凝土拱坝设计规范（2003年版）中的公式计算，有关参数在后续的数据项中读入，其中坝体上下游表面水温由用户输入，水温沿水深的变化规律按规范建议的公式自动推求； =3 按经验公式±47/(T+3.39)计算，其中各期封拱的工况一、二取负值，工况三、四取正值，Ｔ为坝体厚度。 =4 按经验公式±57.57/(T+2.44)计算，其中各期封拱的工况一、二 取负值，工况三、四取正值。 =5 不考虑温度荷载。 =6 输入各拱圈高程上的温度荷载Tm和Td，在§1.15里输入。其中TmCLOS作为各拱圈高程的平均温升荷载，TdCLOS作为各拱圈高程的线性温差荷载。各期封拱条件下的第三、第四工况的温度荷载取-TmCLOS和-TdCLOS。 · 建议值为2或1。 8、IPSN： ·选择是否考虑拱梁之间泊桑比效应的影响。 =1 考虑。 =2 不考虑。 · 建议值为2。 9、IWP： ·选择是否考虑坝体扬压力影响。 =1 考虑。 =2 不考虑。 · 建议值为2。 10、IDFBL： ·选择基础结点分载的取用方式。 =1 顶拱基础结点分载取用沿拱向的相邻结点分载，其余结点分载取用沿梁向的 相邻结点分载。 =2 顶拱基础结点分载取用沿拱向的相邻结点分载，其余结点分载取用沿梁向的 两相邻结点分载之线性外延值。 =3 基础结点分载取用拱向两相邻结点分载之线性外延值与梁向两相邻结点分载 之线性外延值的平均值。 · 建议值为1。 11、IUNIT: ·选择输入数据的单位制。 =1 选用国际惯用单位制，如力(重量)、长度、角度、应力(弹模、变模)、 温度、时间 分别以千牛顿、米、度、MPa、摄氏度、秒为单位，其中导温系数以米2／月 为单位。 =2 选用工程单位制，如力(重量)、长度、角度、应力(弹模)、 温度、时间分别以吨、米、度、吨/米2、摄氏度、秒为单位，其中导温系数以米2／月为单位。其中应力输出单位用MPa。 · 建议值为2。 12、INON ·选择线性或非线性分析。 =0 线性分析。 =1 非线性分析（无初始缺陷）。 =2 非线性分析（针对坝体有已知裂缝的已建拱坝）。 · 建议值为0。 13、IDYN ·选择动力或静力分析。 =0 静力分析（包括拟静力法地震应力分析）。 =1 响应谱法动力分析。 =2 时间历程法动力分析（时间历程法动力分析功能暂缺）。 ·建议值为0。 14、IOPT ·选择应力分析或体形优化。 =0 应力分析。 =1 体形优化。 =2 计算机辅助下的手工体形调整。 =3 演示优化过程。 ·建议值为0。 15、ISHAPE ·特殊体型描述。 =0 一般体型。 =1或3 第一至第二层拱圈之间上游面直立。 =2或3 第NARCH至第NARCH+1层拱圈作为基础垫层，垫层的上下游面直立，垫层体形尺寸由最下层拱圈决定。 =100 非直立拱冠剖面的拱坝体型（目前不推荐！）。 (100 非直立拱冠剖面的拱坝体型（目前不推荐！）。101、102、103分别对应上述的1、2、3。 ·建议值为0。 §3.3 分析计算中所用拱梁网格的拱圈高程行。 　　本输入行根据控制参数NARCH， 依次读入从顶拱开始直到拱冠梁底部高程的共（NARCH+1）个拱圈的高程值，单位m。在拱坝分析中以这些指定高程划分拱梁网格。这些高程可以与“拱圈参数描述行”所取的拱圈高程不同。参见§1.10中的说明。 【数据格式】 　　输入实数：（HIGHT(I),I=1～NARCH+1）。 格式说明：自由格式。即规定本输入项读入NARCH+1个实数， 可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 §3.4 虚结点位置描述行 当控制参数NVIRT=0时免去本批参数行的输入。 NVIRT＞0时，从本批数据读入NVIRT个虚结点的描述数据。 每个虚结点用该点所在的拱号和梁号两个数来表示。虚结点描述的顺序不限。 【数据格式】 　　输入整数：（NRVIRT(I),NCVIRT(I),I=1～NVIRT），即输入顺序为：虚结点 1拱号，虚结点1梁号，虚结点2拱号，虚结点2梁号，……。 格式说明：自由格式。即规定本输入项读入2*NVIRT个整数， 可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 参见图1及其文字说明。 §3.5 基岩材料性质描述主参数行 本输入行用于提供河床及左右岸一定范围内的基岩材料性质。 【数据格式】 　　输入两个实数、两个整数：RE，RMU，NLFT，NRGT。 格式说明：自由格式。即规定本输入项读入2个实数和2个整数，可用一行或一行以上的数据输入，同一行的各数据之间可用逗号或空格作为其分隔符。 【参数说明】 RE 河床及左右岸一定范围内的基岩变形模量，单位MPa（IUNIT=1）或T／m2（IUNIT=2）。 RMU 河床及左右岸一定范围内的基岩泊桑比。 NLFT 基岩材料性质与河床处材料性质（RE，RMU）不同的左岸坝基分区数，建议值为0，NLFT=0表示整个左岸坝基材性与河床处相同。 NRGT 基岩材料性质与河床处材料性质（RE，RMU）不同的右岸坝基分区数，建议值为0，NRGT=0表示整个右岸坝基材性与河床处相同。 §3.6 基岩材料性质描述附加参数行 当(NLFT+NRGT)=0时免去本批参数行的输入。 本批数据给出非均一的坝基材料特性的描述。 【数据格式】 本批数据共(NLFT+NRGT)组，每组输入三个实数：HABUT(I) ， ABUTRE( I) ， ABUTRMU(I)； 格式说明：自由格式。即规定本输入项分为(NLFT+NRGT)组，每组读入3个实数，每组数据需另起一行并可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 HABUT 表示高程。 ABUTRE、ABUTRMU 表示相应的变形模量和泊桑比。 其中：前NLFT行按高程从高向低排列，表示各相应高程以上部分的左岸基岩材性；后NRGT行也按高程从高向低排列，表示各相应高程以上部分的右岸基岩材性。 【例】 　　如果基岩材性主参数行是： # RE RMU NLFT NRGT 2000000，0.2， 0， 2 基岩材性副参数行是： # HABUT ABUTRE ABUTRMU 200, 1500000, 0.3 150, 1800000, 0.2 则它们表示： 左岸、河床及右岸150高程以下的坝基变形模量为2000000，泊桑比为0.2； 右岸150～200高程之间的坝基变形模量为1800000，泊桑比为0.2； 右岸200高程以上的坝基变形模量为1500000，泊桑比为0.3。 §3.7 坝体材料性质描述主参数行 本批输入项用于提供坝体主体部分的材料性质，包括线性特性描述行和非线性特性描述行。 系统假定坝体弹模等材性沿水平径向是相同的。当坝体沿水平径向材性不同时可按等效弹模等提供给程序。系统允许沿坝面分区以指定不同的材性。 坝体材料性质描述主参数行──常规特性部分 【数据格式】 　　线性特性描述行由六个实数和三个整数组成：CE，CMU，RO，ALFAT， a，ALFAWP，NDIFC，NDIFCAE，NDIFCBE； 格式说明：自由格式。即规定本输入项读入6个实数和3个整数，可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 CE 坝体弹模，单位MPa（IUNIT=1）或T／m2（IUNIT=2）。 CMU 坝体泊桑比系数。 RO 坝体容重，单位kN／m3（IUNIT=1）或T／m3（IUNIT=2）。 ALFAT 坝体热线胀系数，单位1/℃。 a 坝体导温系数，单位m2／月。参考值3。 ALFAWP 坝体扬压力折减系数（坝体径向扬压力分布图上折点处的扬压力定义为WP=H2+ALFAWP*(H1-H2)，其中H1为上游水头，H2为下游水 头。扬压力折点位置由“坝体扬压力径向分布折点位置描述行”描述） NDIFC 坝体材料分区控制参数。 ＝0 时表示不分区。 ≠0时表示除了本数据行所描述的坝体主体部分的材性之外，其它坝体分区还具有NDIFC种特殊的材性。 NDIFCAE 属于特殊材性分区的拱单元的描述行数，参见“坝体材料性质描述附加参数行”的说明。 NDIFCBE 属于特殊材性分区的梁单元的描述行数，参见“坝体材料性质描述附加参数行”的说明。 ·建议值为0。 坝体材料性质描述主参数行──非线性特性部分 当控制参数INON=0时免去非线性特性描述行的输入。 【数据格式】 INON≠0时，非线性特性描述行由十个实数组成：FT0，FT，CKI，FC，FU，EC，EU，F_CRACK，H_DAMAGE，V_DAMAGE。 格式说明：自由格式。即规定本输入项读入10个实数， 可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 FT0 坝面混凝土产生初裂的主拉应力值，单位MPa（IUNIT =1）或T／m2（IUNIT=2）。 FT 混凝土单轴受力状态下的极限抗拉强度，单位同上。 CKI 混凝土断裂强度因子，单位MPa·m1/2（IUNIT=1）或T／m2·m1/2（IUNIT=2）。 FC 混凝土单轴受力状态下的极限抗压强度，单位MPa（IUNIT=1）或 T／m2（IUNIT=2）。 FU 混凝土单轴受力状态下压碎时的强度，单位同上。 EC 混凝土单轴受力状态下，对应于极限抗压强度时的应变值。 EU 混凝土单轴受力状态下压碎时的应变值。 F_CRACK 开裂面的抗滑摩擦系数。 H_DAMAGE 水平施工层面的抗拉强度折减系数。 V_DAMAGE 径向竖直施工层面（横缝）的抗拉强度折减系数。 【特别说明】 本输入项的应力强度和应变指标值不论拉、压均以正值输入。 §3.8 坝体材料性质描述附加参数行 当NDIFC=0或(NDIFCAE+NDIFCBE)=0时免去本批参数行的输入。 本批数据用于描述近基岩处等各坝体特殊分区的材性以及分区的范围，由常规材性、非线性材性（当INON≠0即非线性分析时）和分区范围描述三部分组成。 前两部分需NDIFC组，用于描述NDIFC组特殊材料的材性，依输入的顺序编为材料特性第二组、第三组……（“坝体材料性质描述主参数行”中输入的材性编为第一组），每组一行：“坝体材料性质描述附加参数行──常规特性部分”（当INON=0时）或两行：“坝体材料性质描述附加参数行──常规特性部分”和“坝体材料性质描述附加参数行──非线性特性部分”（当INON ≠ 0 时）； 后一部分需 ( NDIFCAE+NDIFCBE)行，用于描述各材料分区的范围。 坝体材料性质描述附加参数行──常规特性部分 【数据格式】 　　本行输入2个实数：CE，CMU。 格式说明：自由格式。即规定本输入项读入2个实数， 可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 参数意义参见“坝体材料性质描述主参数行”说明。 坝体材料性质描述附加参数行──非线性特性部分 当INON=0时免去本参数行的输入。 【数据格式】 　　当INON≠0时，本输入项输入七个实数：FT0，FT，CKI，FC，FU，EC，EU。 　　格式说明：自由格式。即规定本输入项读入7个实数，可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 参数意义参见“坝体材料性质描述主参数行”说明。 坝体材料性质描述附加参数行──特殊材料分区的拱梁单元描述 【数据格式】 本批数据共有(NDIFCAE+NDIFCBE)行，前NDIFCAE行用于指定属于特殊材料分区的拱单元，每行上有四个整数IAEB，IAEE，ICAE1，ICAE2；后NDIFCBE 行用于指定属于特殊材料分区的梁单元，每行上有四个整数IBEB，IBEE，ICBE1，ICBE2。 　　格式说明：自由格式。即规定本输入项分(NDIFCAE+NDIFCBE)组读入，每组4个整数，每组数据需另起一行并可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 IAEB 某分区的起始拱单元号。 IAEE 某分区的终止拱单元号。 ICAE1 该分区的材性组号之一。 ICAE2 该分区的材性组号之二。 IBEB 某分区的起始梁单元号。 IBEE 某分区的终止梁单元号。 ICBE1 该分区的材性组号之一。 ICBE2 该分区的材性组号之二。 【特别说明】 拱单元的编号规则：拱单元从1#编起，在拱梁网格中按从上到下、从左岸到右岸的顺序依次编号（拱冠梁底部高程上的拱圈不作为拱单元，不参加编号）； 梁单元的编号规则：梁单元从1#编起，在拱梁网格中按从上到下、从左岸到右岸的顺序依次编号； 材性组号的编号规则：“坝体材料性质描述主参数行”中输入的材性编为第一组，本节中输入的材性依输入的顺序编为材料特性第二组、第三组……。 前NDIFCAE行中的每一行表示：拱单元号为IAEB～IAEE的拱单元， 其材料特性取分区材性组号ICAE1和ICAE2两种材性的平均值（当该单元介于两个分区之间时可照此取用），如果ICAE2=0，则仅取分区材性组号ICAE1的材性。系统允许将IAEB填为负值，并赋予特殊的含义：即规定其绝对值代表起始的拱单元号，此时该描述行所指定的是拱单元号│IAEB│～│IAEE│之间的与基础相连接的那部分拱单元。 类似地，后NDIFCBE行中的每一行表示：梁单元号为IBEB～IBEE的梁单元， 其材料特性取分区材性组号ICBE1和ICBE2两种材性的平均值（当该单元介于两个分区之间时可照此取用），如果ICBE2=0，则仅取分区材性组号ICBE1的材性。系统允许将IBEB填为负值，并赋予特殊的含义：即规定其绝对值代表起始的梁单元号，此时该描述行所指定的是梁单元号│IBEB│～│IBEE│之间的与基础相连接的那部分梁单元。 其它未指定分区材性组号的拱梁单元，均取用“坝体材料性质描述主参数行”中所描述的材性，即第一组材性。 §3.9 响应谱法地震应力分析所需的场地特性与地震强度参数 当控制参数IDYN = 0时免去本参数项的输入。 IDYN = 1时，本输入项读入的5个数为响应谱法地震动力分析的有关参数。 【数据格式】 　　本行输入5个实数：BETAM，DDTG，GKHz，GKHx，GKv。 格式说明：自由格式。即规定本输入项有3个实数，可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 BETAM 设计地震响应谱曲线的放大系数峰值βmax，一般取值2.25~2.5。取较大值时，相应所求得的动应力较大。 DDTG 设计地震响应谱曲线的地基卓越周期Tg。一般取值0.2。 GKHz 顺河向水平地震加速度系数（重力加速度的倍数）。 GKHx 横河向水平地震加速度系数（重力加速度的倍数）。 GKv 竖向地震加速度系数（重力加速度的倍数）, 含遇合系数。 §3.10 用于描述体形的拱圈层数行 坝体体形是通过一层一层拱圈的参数来描述的。为了便于灵活地布置拱梁网格，本程序系统允许分析计算中所用拱梁网格的拱圈高程与描述坝体体形的拱圈高程不同，也允许两者所取用的拱圈层数不同。在两者拱圈高程、层数取用不同的情况下，前者的拱圈参数通过插值而得。 本行输入用于描述体形的拱圈层数NARCH0。 当NARCH0＝3时，分析计算中所用的拱圈参数通过二次插值而得； 当NARCH0＝4时，分析计算中所用的拱圈参数通过三次插值而得； 当NARCH0＞4时，分析计算中所用的拱圈参数通过三次B样条插值而得。 【数据格式】 　　本行输入一个整数：NARCH0。 格式说明：自由格式。即规定本输入项有1个整数，占用一个输入行。 【参数说明】 NARCH0 用于描述体形的拱圈层数。 §3.11 拱圈参数描述行 本批数据输入拱圈的参数，用于描述坝体的体形。每一层拱圈的参数为一组，共NARCH0组，各组的前后顺序按高程从高到低排列。各组数据的内容根据主控行所选定的拱圈线型按下述各相应线型的要求填写。 §3.11.1 用上下游面圆弧描述的单心等厚圆拱圈参数行(IAD=-1) 【数据格式】 　　本输入项分NARCH0组，每组输入六个实数：HIGHT0，ZCROWNU，THICK， Rd ，PHIl，PHIr； 格式说明：自由格式。即规定本输入项分NARCH0组读入，每组6个实数， 每组数据需另起一行并可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 HIGHT0 拱圈高程，单位m。 ZCROWNU 拱冠梁上游面的Ｚ坐标（坐标正向指向下游），单位m。 THICK 圆拱圈厚度，单位m。 Rd 圆拱圈下游弧半径，单位m。 PHIl 圆拱圈左岸拱端半中心角，单位°。 PHIr 圆拱圈右岸拱端半中心角，单位°。 其中，半中心角以拱冠梁剖面为零值，不论左右岸均以正值输入。 §3.11.2 用上下游面圆弧描述的五心变厚圆拱圈参数行(IAD=-2) 【数据格式】 　　本输入项分NARCH0组，每组输入十二个实数：HIGHT0，ZCROWNU，Rmd，Rmu ，Rld，Rlu，Rrd，Rru，PHIl，PHIml，PHImr，PHIr。 格式说明：自由格式。即规定本输入项分NARCH0组读入，每组12个实数，每组数据需另起一行并可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 HIGHT0 拱圈高程，单位m。 ZCROWNU 拱冠梁上游面的Ｚ坐标（坐标正向指向下游），单位m。 Rmd 中间圆拱圈下游弧半径，单位m。 Rmu 中间圆拱圈上游弧半径，单位m。 Rld 左岸圆拱圈下游弧半径，单位m。 Rlu 左岸圆拱圈上游弧半径，单位m。 Rrd 右岸圆拱圈下游弧半径，单位m。 Rru 右岸圆拱圈上游弧半径，单位m。 PHIl 左岸圆拱圈上游弧拱端半中心角，单位°。 PHIml 中间圆拱圈左岸侧半中心角，单位°。 PHImr 中间圆拱圈右岸侧半中心角，单位°。 PHIr 右岸圆拱圈上游弧拱端半中心角，单位°。 其中，半中心角以拱冠梁剖面为零值，不论左右岸均以正值输入。 §3.11.3 用中心线圆弧描述的单心等厚圆拱圈参数行(IAD=1) 【数据格式】 　　本输入项分NARCH0组，每组输入六个实数：HIGHT0，ZCROWNU，THICK， R，PHIl，PHIr。 格式说明：自由格式。即规定本输入项分NARCH0组读入，每组6 个实数，每组数据需另起一行并可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 HIGHT0 拱圈高程，单位m。 ZCROWNU 拱冠梁上游面的Ｚ坐标（坐标正向指向下游），单位m。 THICK 圆拱圈厚度，单位m。 R 圆拱圈中心线半径，单位m。 PHIl 圆拱圈左岸拱端半中心角，单位°。 PHIr 圆拱圈右岸拱端半中心角，单位°。 其中，半中心角以拱冠梁剖面为零值，不论左右岸均以正值输入。 §3.11.4 用中心线圆弧描述的五心变厚拱圈参数行(IAD=2) 【数据格式】 　　本输入项分NARCH0组，每组输入十二个实数：HIGHT0，ZCROWNU， THICKm，THICKl，THICKr，Rm，Rl，Rr，PHIl，PHIml，PHImr，PHIr。 格式说明：自由格式。即规定本输入项分NARCH0组读入，每组12个实数，每组数据需另起一行并可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 HIGHT0 拱圈高程，单位m。 ZCROWNU 拱冠梁上游面的Ｚ坐标（坐标正向指向下游），单位m。 THICKm 中间圆拱圈厚度，单位m。 THICKl 左拱端拱圈厚度，单位m。 THICKr 右拱端拱圈厚度，单位m。 Rm 中间圆拱圈中心线半径，单位m。 Rl 左岸圆拱圈中心线（近似圆弧）半径，单位m。 Rr 右岸圆拱圈中心线（近似圆弧）半径，单位m。 PHIl 左岸圆拱圈中心线拱端半中心角，单位°。 PHIml 中间圆拱圈左岸侧半中心角，单位°。 PHImr 中间圆拱圈右岸侧半中心角，单位°。 PHIr 右岸圆拱圈中心线拱端半中心角，单位°。 其中，半中心角以拱冠梁剖面为零值，不论左右岸均以正值输入。 §3.11.5 用中心线描述的抛物线型拱圈参数行(IAD=3) 【数据格式】 　　本输入项分NARCH0组，每组输入十个实数：HIGHT0， ZCROWNU ， THICKc，THICKl，THICKr，Rcl，Rcr，PHIl，PHIr，Tindex。 格式说明：自由格式。即规定本输入项分NARCH0组读入，每组10个实数，每组数据需另起一行并可用一行或一行以上的数据输入，同一行各数据之间可用逗号或空格作为其分隔符。 【参数说明】 HIGHT0 拱圈高程，单位m。 ZCROWNU 拱冠梁上游面的Ｚ坐标（坐标正向指向下游），单位m。 THICKc 拱冠梁处的拱圈厚度，单位m。 THICKl 左拱端拱圈厚度，单位m。 THICKr 右拱端拱圈厚度，单位m。 Rcl 拱冠梁处左岸侧的拱圈中心线曲率