CN103577641A - 基于多相材料布局优化的坝料分区方法 - Google Patents
基于多相材料布局优化的坝料分区方法 Download PDFInfo
- Publication number
- CN103577641A CN103577641A CN201310540326.9A CN201310540326A CN103577641A CN 103577641 A CN103577641 A CN 103577641A CN 201310540326 A CN201310540326 A CN 201310540326A CN 103577641 A CN103577641 A CN 103577641A
- Authority
- CN
- China
- Prior art keywords
- optimization
- design
- energy density
- strain energy
- dam
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于多相材料布局优化的坝料分区方法。是一种通过设计多相材料布局进行过程优化布置的新方法。以实际工程作为输入条件;根据邓肯E-B模型确定材料初始弹性模量;确定垫层、主次堆石材料的弹性模量;然后通过有限元软件获取结果数据,并基于拓扑优化原理进行设计变量更新。当所有材料均达到各自的体积约束临界值时得到布局方案。与原设计相比,优化后的面板在水压力作用下变形小,分布均匀,有利于阻止面板发生断裂。
Description
技术领域
本发明涉及一种坝料分区的设计方法,更具体地说,涉及一种基于多相材料拓扑优化理论的面板堆石坝的坝料分区设计方法。
背景技术
水坝是水利建筑物的重要组成部分,是重要的挡水建筑物。其中,面板堆石坝因具有经济、可靠、广泛适应性和环保等显著优点而受到坝工界的欢迎。
长期以来,面板因坝体变形不均匀而产生裂缝极为常见。而控制坝体变形可从坝高、河谷形状、坝料物理特性及坝料分区、施工程序等角度入手。其中河谷形状和坝料物理特性为客观量,而坝高、料分区和施工程序为主观可变量。因此,在坝高给定的情况下,坝料分区和施工程序成为影响坝体变形的主要因素。传统的坝料分区设计方法仅对基于经验的原设计做适当的调整,而原设计存在的缺陷无法在优化后得到有效解决,结构的性能并不能得到明显改善。因此,在尺寸优化前进行拓扑优化和形状优化,有利于对设计从源头得到更好的优化。
发明内容
有鉴于此,本发明提供了一种基于拓扑优化理论的面板堆石坝的坝料分区设计方法。在满足工程要求的情况下,能够得到最优材料布局方案,使结构在现有工况下表现出最大刚度,保证开挖料得到充分利用。
为实现上述目的,本发明提供如下技术方案:
基于多相材料拓扑优化理论的面板堆石坝的坝料分区设计方法,包括以下步骤:
a) 建立结构的有限元分析模型,初始化算法各参数,令;
b) 分析材料集合为{Me,j}(j属于m-1,m,…,Nm)的结构,获得各单元的应变能密度;
c) 升序排列单元的应变能密度,将序列前端一定数量的低应变能密度单元的材料由Me,m改变为Me,m-1,并累计含Me,m材料单元的体积;
d) 若含Me,m材料单元的体积没有达到其约束临界值,则返回到第 b)步重分析,否则到第 e)步;
e) 若所有材料均达到各自的体积约束临界值,则到第 f)步;否则,令m=m-1,返回到第 b)步重分析;
f) 停止迭代;
在上述优化设计方法中,所涉及到的建模,优化软件均为ANSYS软件。
在上述步骤(f)中为确保算法稳定,本方法迭代过程中每次仅更新不超过 5%的总单元数的材料号。
本发明的有益效果是,可以在满足实际工程条件的情况下,优化后的面板在水压力作用下变形小,分布均匀,有利于阻止面板发生断裂。从而证明该方法可为此类工程实践提供具有理论基础的参考设计。
附图说明
下面结合附图和实施例对本发明专利进一步说明。
图1 :优化分区前的坝体剖面有限元模型。进行多相材料分区前的坝剖面有限元网格划分并作为优化设计的设计域。
图 2 :原始坝料分区方法设计的布局方案。
图 3 :使用拓扑优化方法优化后得到的多相材料布局方案。基于拓扑优化方法对设计域优化后得到的多相材料布局,对比图2,图3可以看出,新型设计方法得到的设计方案与原始设计方法有明显区别。材料分布更具规律及合理性。
Claims (6)
1.一种基于多相材料拓扑优化理论的面板堆石坝的坝料分区设计方法,其特征在于,包括步骤:
以实际工程作为输入条件;
根据邓肯E-B模型确定材料初始弹性模量;
确定垫层、主次堆石材料的弹性模量;
a) 建立结构的有限元分析模型,初始化算法各参数,令;
b) 分析材料集合为{Me,j}(j∈m-1,m,…,Nm)的结构,获得各单元的应变能密度;
c) 升序排列单元的应变能密度,将序列前端一定数量的低应变能密度单元的材料由Me,m改变为Me,m-1,并累计含Me,m材料单元的体积;
d) 若含Me,m材料单元的体积没有达到其约束临界值,则返回到第 b)步重分析,否则到第 e)步;
e) 若所有材料均达到各自的体积约束临界值,则到第 f)步;否则,令m=m-1,返回到第 b)步重分析;
f) 停止迭代;
结构安全性校核,并按需要进行局部调整;完成最终设计方案。
2.如权利要求1所述的优化设计方法,其特征在于,所述建模和前后处理、求解软件为ANSYS。
3.如权利要求1所述的优化设计方法,其特征在于,所述拓扑优化方法应用应变能密度排序算法通过 APDL 编写实现(二次开发)。
4.如权利要求1所述的优化设计方法,其特征在于,步骤1.b)中的所述材料初始弹性模量的确定方法为Ei=K*Pa(σ3/pa)n。
5.如权利要求1所述的优化设计方法,其特征在于,步骤1.d)中所述的荷载及边界约束条件具体为结构总体实际使用时的荷载和自由度。
6.如权利要求1所述的优化设计方法,其特征在于,步骤1.f)中所述的材料更新方法为应变能密度排序法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310540326.9A CN103577641A (zh) | 2013-11-05 | 2013-11-05 | 基于多相材料布局优化的坝料分区方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310540326.9A CN103577641A (zh) | 2013-11-05 | 2013-11-05 | 基于多相材料布局优化的坝料分区方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103577641A true CN103577641A (zh) | 2014-02-12 |
Family
ID=50049409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310540326.9A Pending CN103577641A (zh) | 2013-11-05 | 2013-11-05 | 基于多相材料布局优化的坝料分区方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103577641A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599354A (zh) * | 2016-11-08 | 2017-04-26 | 西安理工大学 | 一种内序列多相材料拓扑优化方法 |
CN107273580A (zh) * | 2017-05-22 | 2017-10-20 | 西安理工大学 | 一种确定多相双模量材料布局问题体积约束的方法 |
CN107491599A (zh) * | 2017-08-03 | 2017-12-19 | 华中科技大学 | 一种应力约束下多相材料柔性机构拓扑优化方法 |
CN109492261A (zh) * | 2018-10-15 | 2019-03-19 | 国家电网有限公司 | 一种高寒地区坝体变形控制方法、装置及控制器 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976280A (zh) * | 2010-10-14 | 2011-02-16 | 西北工业大学 | 多相材料结构构型设计方法 |
-
2013
- 2013-11-05 CN CN201310540326.9A patent/CN103577641A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976280A (zh) * | 2010-10-14 | 2011-02-16 | 西北工业大学 | 多相材料结构构型设计方法 |
Non-Patent Citations (6)
Title |
---|
JIAO SHI等: "Effects of Moduli Differences on Layout Optimization of a Continuum with Multiple Materials", 《ADVANCED MATERIALS RESEARCH》, vol. 217218, 30 March 2011 (2011-03-30), pages 1414 - 1418 * |
KUN CAI等: "Volumes Constrained Layout Optimization of a Continuum with Multiple Materials", 《KEY ENGINEERING MATERIALS》, vol. 480481, 30 June 2011 (2011-06-30), pages 619 - 623 * |
南洪: "拓扑优化理论及其在面板堆石坝优化设计中的应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, no. 11, 15 November 2008 (2008-11-15), pages 037 - 39 * |
南洪: "邓肯E-B模型原理及在成屏面板坝中的应用", 《人民长江》, vol. 39, no. 11, 5 June 2008 (2008-06-05), pages 72 - 74 * |
王汝军: "白沙面板堆石坝剖面优化设计", 《湖北水力发电》, no. 5, 30 September 2009 (2009-09-30), pages 48 - 52 * |
高照良等: "拓扑优化方法在面板堆石坝分区设计中的应用", 《农业工程学报》, vol. 29, no. 16, 15 August 2013 (2013-08-15) * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599354A (zh) * | 2016-11-08 | 2017-04-26 | 西安理工大学 | 一种内序列多相材料拓扑优化方法 |
CN106599354B (zh) * | 2016-11-08 | 2020-02-18 | 西安理工大学 | 一种内序列多相材料拓扑优化方法 |
CN107273580A (zh) * | 2017-05-22 | 2017-10-20 | 西安理工大学 | 一种确定多相双模量材料布局问题体积约束的方法 |
CN107273580B (zh) * | 2017-05-22 | 2020-11-17 | 西安理工大学 | 一种确定多相双模量材料布局问题体积约束的方法 |
CN107491599A (zh) * | 2017-08-03 | 2017-12-19 | 华中科技大学 | 一种应力约束下多相材料柔性机构拓扑优化方法 |
CN109492261A (zh) * | 2018-10-15 | 2019-03-19 | 国家电网有限公司 | 一种高寒地区坝体变形控制方法、装置及控制器 |
CN109492261B (zh) * | 2018-10-15 | 2023-05-12 | 国家电网有限公司 | 一种高寒地区坝体变形控制方法、装置及控制器 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103577641A (zh) | 基于多相材料布局优化的坝料分区方法 | |
Yammani et al. | Optimal placement and sizing of distributed generations using shuffled bat algorithm with future load enhancement | |
CN102043883B (zh) | 基于材料破损约束的连续体结构拓扑设计建模及优化设计方法 | |
CN112910012B (zh) | 一种配电系统弹性提升评估方法、存储介质及计算设备 | |
CN105279615A (zh) | 一种基于双层规划的主动配电网网架规划方法 | |
CN103795068A (zh) | 一种高压配网动态无功补偿装置容量的优化配置方法 | |
CN109829217A (zh) | 压裂性裂缝油藏产能模拟方法及装置 | |
CN103632046A (zh) | 一种电网潮流计算方法 | |
Ersavas et al. | Optimum allocation of FACTS devices under load uncertainty based on penalty functions with genetic algorithm | |
Chen et al. | Physics-shielded multi-agent deep reinforcement learning for safe active voltage control with photovoltaic/battery energy storage systems | |
CN107808200B (zh) | 配电网短时检修计划优化方法和装置 | |
Kreishan et al. | Allocation of dump load in islanded microgrid using the mixed-integer distributed ant colony optimization | |
Zhao et al. | Stress‐constrained multiscale topology optimization with connectable graded microstructures using the worst‐case analysis | |
CN117175548B (zh) | 一种基于随机潮流的连锁故障紧急控制方法 | |
Kumar et al. | Comparative analysis of particle swarm optimization variants on distributed generation allocation for network loss minimization | |
CN111274624B (zh) | 一种基于rbf代理模型的多工况异形节点拓扑优化设计方法 | |
Bachir et al. | An efficient chaotic cuckoo search framework for solving non-convex optimal power flow problem | |
CN107919985A (zh) | 一种复杂网络级联失效容量负载架构的应用 | |
CN103984832A (zh) | 一种铝电解槽电场仿真分析方法 | |
CN105631066A (zh) | 一种基于模态试验结果的有限元三维模型模态调整方法 | |
Yao et al. | Topology optimization design of pre-stressed plane entity steel structure with the constrains of stress and displacement | |
CN104793107B (zh) | 一种基于改进opa模型的电网连锁故障确定方法 | |
Wang et al. | Stability constrained opf in microgrids: A chance constrained optimization framework with non-gaussian uncertainty | |
CN104268312B (zh) | 用于抵抗长期高温和内压共同作用的安全壳结构计算方法 | |
Yin et al. | Reduced-dimensional skip-inception feature-aggregated classified proportional-integral-derivative for suppression of mixed-mode oscillations in hydropower units |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140212 |