混凝土拱坝坝体优化方法技术

本发明专利技术公开了一种混凝土拱坝坝体优化方法。凝土拱坝坝体优化能够在保证工程安全的前提下，提出经济可靠的坝体设计方案。传统的拱坝坝体优化方法重复设计法的工作量大、过程繁琐和效率低，因此提出一种高效可靠的混凝土拱坝坝形优化方法是十分必要的。而在控制拱坝坝体优化的因素中，应力控制是最为重要的一个。本发明专利技术采用快速多级边界元法对优化的混凝土拱坝坝体进行应力控制分析，相较于采用有限元法进行应力控制分析的常规拱坝坝体优化方法，在同等自由度的前提下，本发明专利技术耗时更少，计算效率更高。此外，快速多级边界元法无需对坝体进行三维网格划分，只需要对坝体外表面划分二维单元，过程更为简单，网格划分质量要求较低，更为便捷。更为便捷。更为便捷。

【技术实现步骤摘要】
混凝土拱坝坝体优化方法


[0001]本专利技术属于坝体优化
，具体涉及一种混凝土拱坝坝体优化方法。

技术介绍

[0002]所谓拱坝体型优化设计，就是在已知的给定参数条件下，寻找满足工程的施工、运用、强度、稳定性等约束条件，并使目标函数取最小值的体型方案。因此，优化设计的首要一步是将工程实际问题表示成数学表达式，即建立数学模型，包括选定设计变量、选择目标函数、建立约束方程，然后采用适宜的优化方法，求出问题的解。
[0003]关于拱坝体型优化的研究，国外是从20世纪60年代开始的。1969年，R.Sharp提出了用数学规划法对拱坝进行体型优化设计，从而使拱坝体型优化设计作为一个专门课题，引起学术界和工程界的注意；1973年，加拿大学者W.Stensch提出了自由型拱坝设计的模型；1974年G.A.Harley和G.M.Mcneice形成了有限元分析用的网格自由划分系统，进一步改进了W.Stcnsch的工作；1975年R.E.Rhichetts和D.C.Liewicz发表了《混凝土坝的形状优化设计》，其中也研究了拱坝；1982年，在上海举行的国际有限元会议上，K.Wassermann发表了一篇拱坝体型优化的论文等。
[0004]我国拱坝体型优化设计的研究是在20世纪70年代末期由朱伯芳院士提出来的，尽管起步较晚，但发展很快。我国拱坝优化设计研究已经进行了30多年，在优化理论和实践上均取得了显著的成绩，目前我国在拱坝优化领域已处于国际领先地位。从基于材料力学优化发展到有限元优化，从静力优化发展到动力优化，从线性优化发展到非线性全过程分析，拱圈的线性从单心圆发展到二次曲线、混合线型等，逐步发展和完善了我国提出的高混凝土坝优化理论与方法。
[0005]目前，有限元法作为一种很受欢迎的数值分析方法，逐步成为坝型优化应力计算的主要方法，但一般而言，拱坝为超大规模水利工程，有限元分析需要划分大量的三维体网格，计算量巨大，耗时长，需要高额的计算成本。

技术实现思路

[0006]针对现有技术存在的不足，本专利技术提出了一种混凝土拱坝坝体优化方法，基于快速多级边界元法，模型重建过程更为简单，同等自由度下，精度更高。
[0007]为了解决上述技术问题，本专利技术提供的混凝土拱坝坝体优化方法，其特征在于：包括以下步骤：
[0008]步骤1：选择用以确定拱坝体型的基本参数，如拱冠梁处拱圈中心线的曲率半径、拱端中心角，以及确定坝体厚度的基本参数如拱冠梁厚度、拱端厚度等作为优化设计变量；
[0009]步骤2：选择适当的目标函数，一般以安全性为关键性指标，即以应力为目标函数；
[0010]步骤3：为保证工程安全和施工方便，确定坝体几何尺寸和工作性态两方面的约束，主要包括几何约束条件、应力约束条件和稳定约束条件；
[0011]步骤4：将拱坝的体型抽象成几何模型并用一组参数来完整地描述，基于快速多级
边界单元进行坝体网格自动化划分以及荷载的自动施加。
[0012]步骤5：基于快速多级边界元法自动化计算应力，利用等效应力法求出拱坝上下游面的三个应力分量，然后由该处微分体的平衡条件求得另外三个应力分量，进一步决定拱坝上下游面的主应力，进而确定出最优拱坝坝形。
[0013]作为优选，所述步骤3中输入应力约束条件为坝轴线位置移动范围的限制，坝顶最小厚度的限制，倒悬度约束，保凸性要求。
[0014]进一步地，所述步骤3中输入几何约束条件为坝体内部的最大主拉应力及最大主压应力均小于坝体混凝土的允许强度，即
[0015]|σ
t
|
max
≤[σ
t
],|σ
c
|
max
≤[σ
c
]ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0016]式中：[σ
t
]为容许拉应力，[σ
c
]为容许压应力。
[0017]进一步地，所述步骤4中输入拱坝几何模型参数包括拱冠梁和水平拱圈的参数化描述，拱冠梁参数主要包括冠梁中心线曲线方程Yc(z)和拱冠梁厚度Tc(z)，一般将Yc(z)和Tc(z)都设为z的3次方程，故只需有4个控制高程(z＝z1、z2、z3、z4)处对应的Yc和Tc便可按插值方法确定Yc(z)和Tc(z)方程中的系数。这样拱冠梁剖面便可用以下12个参数来描述：
[0018][0019]上述Yci和Tci分别为控制高程zi处的拱冠梁中心线y坐标及拱冠梁厚度。水平拱圈只要确定了拱圈中心线的方程及拱冠至拱端的厚度变化规律便可确定水平拱圈的形状。
[0020]进一步地，所述步骤5中的利用快速多级边界元法自动化计算应力，具体基本计算公式为：
[0021]c
ij
(x)u
i
(x)＝∫
Γ
[U
ij
(x,y)t
j
(y)
‑
T
ij
(x,y)u
j
(y)]dΓ(y)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0022]其中2π
·
c
ij
(x)为点x在边界内部的角，当点x所在的边界处平滑时，有c
ij
(x)＝0.5。u
i
(x)和t
i
(y)为位移和面力，以及
[0023][0024]其中i,j＝1,2,3，r为点x和y的距离，代表r关于x
i
的偏导数，v为泊松比，μ为剪切模量，以及
[0025][0026]进一步地，所述步骤5中的利用的快速多级加速边界元法，快速多级算法的具体计算步骤为：
[0027](5.1)对求解域进行计算域边界的离散；(5.2)尽可能采用最小的立方体包含求解域，将其命名为根立方体，也就是第0层立方体，设定每个叶子立方体中最多包含的节点个
数，构造八叉树结构；(5.3)指定多级展开系数，对于每个层数大于1的立方体，定义存储多极矩的矩阵。对于每个层数大于2的立方体，定义存储局部展开系数的矩阵。对于每一个叶子立方体，定义矩阵存储系数。定义矩阵存储每个近点直接计算的系数；(5.4)对于每个叶子立方体分别计算其相关系；(5.5)计算每个叶子立方体的多极矩；(5.6)对于处于第L层每一个非叶子立方体计算其多极矩，这个多极矩由其子立方体传递而来，这个过程可以称为向上传递。(5.7)对于第L层的每一个立方体，其中L从2到最终层，通过其相互作用列表中的立方体传递而来；(5.8)分别计算近点和远点的贡献系数；(5.9)更新迭代向量；(5.10)不断重读以上过程直到达到迭代算法的控制残差，进而求得应力和位移的未知量。
[0028]进一步地，所述步骤5中的等效应力法为：
[0029][0030]τ
tz
＝τ
zx
cosα+τ...

【技术保护点】

【技术特征摘要】
1.一种混凝土拱坝坝体优化方法，其特征在于：包括以下步骤：步骤1：选择用以确定拱坝体型的基本参数以及确定坝体厚度的基本参数，作为优化设计变量；所述拱坝体型的基本参数包括拱冠梁处拱圈中心线的曲率半径和拱端中心角；所述坝体厚度的基本参数包括拱冠梁厚度、拱端厚度；步骤2：选择适当的目标函数，以安全性为关键性指标，即以应力为目标函数；步骤3：为保证工程安全和施工方便，确定坝体几何尺寸和工作性态两方面的约束，包括几何约束条件、输入应力约束条件和稳定约束条件；步骤4：将拱坝的体型抽象成几何模型并用一组参数来完整地描述，基于快速多级边界单元进行坝体网格自动化划分以及荷载的自动施加；步骤5：基于快速多级边界元法自动化计算应力，利用等效应力法求出拱坝上下游面的三个应力分量，然后由该处微分体的平衡条件求得另外三个应力分量，进一步决定拱坝上下游面的主应力，进而确定出最优拱坝坝形。2.根据权利要求1所述的混凝土拱坝坝体优化方法，其特征在于：所述步骤3中，输入应力约束条件包括坝轴线位置移动范围的限制、坝顶最小厚度的限制、倒悬度约束和保凸性要求。3.根据权利要求1或2所述的混凝土拱坝坝体优化方法，其特征在于：所述步骤3中，输入几何约束条件为坝体内部的最大主拉应力及最大主压应力均小于坝体混凝土的允许强度，即|σ
t
|
max
≤[σ
t
],|σ
c
|
max
≤[σ
c
]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)；式中：[σ
t
]为容许拉应力，[σ
c
]为容许压应力。4.根据权利要求1或2所述的混凝土拱坝坝体优化方法，其特征在于：所述步骤4中，输入拱坝几何模型参数包括拱冠梁和水平拱圈的参数化描述；所述拱冠梁参数包括冠梁中心线曲线方程Yc(z)和拱冠梁厚度Tc(z)；将Yc(z)和Tc(z)都设为z的3次方程，故只需有4个控制高程(z＝z1、z2、z3、z4)处对应的Yc和Tc便可按插值方法确定Yc(z)和Tc(z)方程中的系数；拱冠梁剖面用以下12个参数来描述：上述Yci和Tci分别为控制高程zi处的拱冠梁中心线y坐标及拱冠梁厚度；水平拱圈只要确定了拱圈中心线的方程及拱冠至拱端的厚度变化规律便能确定水平拱圈的形状。5.根据权利要求3所述的混凝土拱坝坝体优化方法，其特征在于：所述步骤4中，输入拱坝几何模型参数包括拱冠梁和水平拱圈的参数化描述；所述拱冠梁参数包括冠梁中心线曲线方程Yc(z)和拱冠梁厚度Tc(z)；将Yc(z)和Tc(z)都设为z的3次方程，故只需有4个控制高程(z＝z1、z2、z3、z4)处对应的Yc和Tc便可按插值方法确定Yc(z)和Tc(z)方程中的系数；拱冠梁剖面用以下12个参数来描述：
上述Yci和Tci分别为控制高程zi处的拱冠梁中心线y坐标及拱冠梁厚度；水平拱圈只要确定了拱圈中心线的方程及拱冠至拱端的厚度变化规律便能确定水平拱圈的形状。6.根据权利要求1或2或5所述的混凝土拱坝坝体优化方法，其特征在于：所述步骤5中，利用快速多级边界元法自动化计算应力，具体公式为：c
ij
(x)u
i
(x)＝∫
Γ
[U
ij
(x,y)t
j
(y)
‑
T
ij
(x,y)u
j
(y)]dΓ(y)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)；其中，x,y为研究域内的点，i,j＝1,2,3，x,2π
·
c
ij
(x)为点x在边界内部的角，当点x所在的边界处平滑时，有c<...

【专利技术属性】
技术研发人员：高宇，刘彪，李通盛，王桥，周伟，岳强，黄诚斌，常晓林，
申请(专利权)人：大唐宣威水电开发有限公司，
类型：发明
国别省市：