CN108629111B - 一种混凝土重力坝材料参数空间变异性的模拟方法 - Google Patents
一种混凝土重力坝材料参数空间变异性的模拟方法 Download PDFInfo
- Publication number
- CN108629111B CN108629111B CN201810410957.1A CN201810410957A CN108629111B CN 108629111 B CN108629111 B CN 108629111B CN 201810410957 A CN201810410957 A CN 201810410957A CN 108629111 B CN108629111 B CN 108629111B
- Authority
- CN
- China
- Prior art keywords
- matrix
- lnx
- gravity dam
- grid cell
- concrete
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种混凝土重力坝材料参数空间变异性的模拟方法,包括以下步骤:建立重力坝模型,对所述模型进行网格划分,获取每个网格单元的编号和形心坐标,组成矩阵En×3;利用MATLAB软件生成随机矩阵Am×n;利用自相关函数计算所述矩阵En×3的相关系数矩阵ρn×n;对所述矩阵ρn×n进行Cholesky分解,得到上三角矩阵Bn×n;利用所述矩阵Bn×n对所述矩阵Am×n进行线性变换,得到矩阵Dm×n;由所述矩阵Dm×n计算得到相关随机场Xm×n;本发明所述模拟方法更加符合工程实际,能够在后期结果中为重力坝的安全运行提供更加准确的数据支撑。
Description
技术领域
本发明属于模拟方法领域,具体是一种混凝土重力坝材料参数空间变异性 的模拟方法。
背景技术
混凝土重力坝是水利水电工程拦河大坝常用的一种坝型。混凝土是一种由 砂石骨料和水泥浆搅拌混合而成的非均质材料。在修建混凝土重力坝的过程中, 一般采用分段、分块方式浇筑混凝土。由于浇筑仓面大、浇筑施工质量不均匀, 以及混凝土材料自身的非均质性,因此,混凝土材料的物理力学性质具有空间 变异性,即对于同类混凝土,其在坝体不同位置处的物理力学参数是不一样的。 混凝土材料参数的空间变异性对重力坝的应力变形和结构安全会产生影响,因 此,有必要研究材料参数空间变异性对混凝土重力坝结构特性的影响。这就需 要解决混凝土重力坝结构分析中材料参数空间变异性的模拟技术问题。
目前混凝土重力坝结构分析中材料参数空间变异性的模拟方法为:将坝体 断面离散为有限个网格单元,通过给网格单元赋予不同的材料参数来模拟材料 参数的空间变异性,并假定不同网格单元的材料参数服从某种概率分布(如 Weibu l l分布、正态分布等),这种方法反映了材料参数在空间中的随机分布, 但没有考虑材料参数的空间相关性,即相邻位置(或相邻网格单元)的材料参 数具有某种程度的正相关性,也就是相邻位置(或相邻网格单元)的材料参数 不会有太大差别。目前的模拟方法没有考虑这种空间相关性,从而导致相邻网 格单元的材料参数出现较大差异的不合理情况。
发明内容
本发明的目的在于解决上述现有技术中存在的问题,提供一种混凝土重力 坝材料参数空间变异性的模拟方法,所述模拟方法更加符合工程实际,能够提 供更加准确的数据。
为实现上述发明目的,本发明所采用的技术方案是:
一种混凝土重力坝材料参数空间变异性的模拟方法,包括以下步骤:
(1)建立重力坝模型,对所述模型进行网格划分,获取每个网格单元的 编号和形心坐标,组成矩阵其中n为网格单元的个 数,xi,yi分别为第i个网格单元形心的横、纵坐标,1≤i≤n;
(2)利用MATLAB软件生成随机矩阵Am×n,其中m为每个网格单元参数 的随机抽样次数,n为网格单元的个数;
(3)利用自相关函数计算所述矩阵En×3的相关系数矩阵其中元素ρi,j为第i和j个网格单元中参数的相关 系数,1≤i≤n,1≤j≤n;
(4)对所述矩阵ρn×n进行Cholesky分解,即ρn×n=Bn×n T·Bn×n,得到上三角矩 阵Bn×n;
(5)利用所述矩阵Bn×n对所述矩阵Am×n进行线性变换,即Dm×n=Am×n·Bn×n, 得到矩阵Dm×n;
(6)相关随机场Xm×n的计算:通过公式lnXm×n=σlnx×Dm×n+μlnx计算得到 lnXm×n,则相关随机场Xm×n=eσlnx×Dm×n+μlnx,其中μlnx和σlnx为满足对数正态分布的 参数。
优选的,步骤(1)中所述划分网格的尺度为1~1.5m。
优选的,步骤(3)中所述自相关函数为负指数自相关函数,所述元素其中θx和θy分别为横、纵坐标方向上混凝土材料的 相关距离,xi,yi分别为第i个网格单元形心的横、纵坐标,xj,yj分别为第j个 网格单元形心的横、纵坐标,1≤i≤n,1≤j≤n。
优选的,步骤(6)中所述满足对数正态分布的参数μlnx和σlnx通 过下式计算,其中μx为混凝土材 料参数的均值,αx为混凝土材料参数的变异系数。
本发明的有益效果是:本发明提供的混凝土重力坝材料参数空间变异性的 模拟方法基于随机场理论,采用自相关函数考虑混凝土材料参数的空间相关性, 提出一种基于空间变异性且满足已知相关关系的相关随机场的模拟方法;本发 明所述模拟方法更加符合工程实际,能够在后期结果中为重力坝的安全运行提 供更加准确的数据支撑;该模拟方法能够用于重力坝有限元分析中材料参数空 间变异性的模拟,为考虑材料参数空间变异性的重力坝结构分析提供技术支持。
附图说明
图1.本发明所述模拟方法的流程图;
图2.重力坝模型的剖面尺寸图;
图3.重力坝模型网格划分图;
图4.抗拉强度随机场的云图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面本发明中的技术方 案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例, 而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有 作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,一种混凝土重力坝材料参数空间变异性的模拟方法,包括 以下步骤:
1.某一重力坝模型的剖面尺寸如图2所示,坝高150m,坝底宽130m, 上游坝面折坡坡率为1:0.2,下游坝坡坡率为1:0.8,已知其抗拉强度服从对 数正态分布,其均值为2.15MPa,变异系数为0.1;如图3所示,将所述重力 坝模型划分了4957个网格单元,获取该模型中每个网格单元的编号和形心坐 标,组成矩阵
2.利用MATLAB软件中的normrnd(0,1,m,n)函数随机产生独立标准正态 分布随机变量抽样序列A10×4957,即单元的个数为4957,每个单元参数的抽样 随机组为10组,生产的随机矩阵
3.利用负指数自相关函数计算所述矩阵E4957×3的相关系数矩阵ρ4957×4957,选 取横、纵坐标方向上混凝土材料的相关距离θx和θy分别为10m和5m,根据 公式得到相关系数矩阵ρ4957×4957;
4.对所述矩阵ρ4957×4957进行Cholesky分解,即
ρ4957×4957=BT 4957×4957·B4957×4957,
其中矩阵B4957×4957为上三角矩阵,
5.对所述矩阵A10×4957进行线性变换,即D10×4957=A10×4957·B4957×4957,得到矩阵
6.首先计算得到抗拉强度满足对数正态分布的参数μlnx和σlnx分别为 0.7605和0.09975,再通过公式lnX10×4957=σlnx×D10×4957+μlnx计算得到lnX10×4957,则 抗拉强度随机场图4是10组 抗拉强度随机场中的四组云图展示结果,从模拟数据和云图来看,本发明所 述模拟方法能够考虑抗拉强度的空间变异性。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限 制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员 应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其 中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的 本质脱离本发明各实施例技术方案的精神和范围。
Claims (2)
1.一种混凝土重力坝材料参数空间变异性的模拟方法,其特征在于,包括以下步骤:
(1)建立重力坝模型,对所述模型进行网格划分,获取每个网格单元的编号和形心坐标,组成矩阵其中n为网格单元的个数,xi,yi分别为第i个网格单元形心的横、纵坐标,1≤i≤n;
(2)利用MATLAB软件生成随机矩阵Am×n,其中m为每个网格单元参数的随机抽样次数,n为网格单元的个数;
(3)利用自相关函数计算所述矩阵En×3的相关系数矩阵其中元素ρi,j为第i和j个网格单元中参数的相关系数,1≤i≤n,1≤j≤n;所述自相关函数为负指数自相关函数,所述元素其中θx和θy分别为横、纵坐标方向上混凝土材料的相关距离,xi,yi分别为第i个网格单元形心的横、纵坐标,xj,yj分别为第j个网格单元形心的横、纵坐标,1≤i≤n,1≤j≤n;
(4)对所述矩阵ρn×n进行Cholesky分解,即ρn×n=Bn×n T·Bn×n,得到上三角矩阵Bn×n;
(5)利用所述矩阵Bn×n对所述矩阵Am×n进行线性变换,即Dm×n=Am×n·Bn×n,得到矩阵Dm×n;
(6)相关随机场Xm×n的计算:通过公式lnXm×n=σlnx×Dm×n+μlnx计算得到lnXm×n,相关随机场其中μlnx和σlnx为满足对数正态分布的参数,通过下式计算:
其中μx为混凝土材料参数的均值,αx为混凝土材料参数的变异系数。
2.根据权利要求1所述的一种混凝土重力坝材料参数空间变异性的模拟方法,其特征在于,步骤(1)中所述划分网格的尺度为1~1.5m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810410957.1A CN108629111B (zh) | 2018-05-02 | 2018-05-02 | 一种混凝土重力坝材料参数空间变异性的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810410957.1A CN108629111B (zh) | 2018-05-02 | 2018-05-02 | 一种混凝土重力坝材料参数空间变异性的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108629111A CN108629111A (zh) | 2018-10-09 |
CN108629111B true CN108629111B (zh) | 2019-09-10 |
Family
ID=63695394
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810410957.1A Active CN108629111B (zh) | 2018-05-02 | 2018-05-02 | 一种混凝土重力坝材料参数空间变异性的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108629111B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110046400B (zh) * | 2019-03-22 | 2022-11-01 | 昆明理工大学 | 一种基于h-p型有限元法对混凝土重力坝进行静力特性分析的高精度数值模拟计算方法 |
CN110110394B (zh) * | 2019-04-18 | 2022-11-18 | 西安理工大学 | 一种堆石坝材料参数的空间随机场模拟方法 |
CN111539097A (zh) * | 2020-04-17 | 2020-08-14 | 安徽理工大学 | 一种深部煤层底板岩体参数随机场建模方法 |
CN111914327B (zh) * | 2020-07-28 | 2023-07-04 | 同济大学 | 一种考虑地层参数随机场模型的长隧道纵向抗震设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021853A (zh) * | 2016-05-09 | 2016-10-12 | 青岛理工大学 | 一种空间变异特性下边坡稳定可靠度分析软件的开发方法 |
CN106442153A (zh) * | 2016-09-23 | 2017-02-22 | 河海大学 | 模拟混凝土重力坝坝踵水力劈裂的试验装置及试验方法 |
CN107862146A (zh) * | 2017-11-15 | 2018-03-30 | 中国地质大学(武汉) | 序贯高斯模拟岩土参数空间变异结构的有限元分析方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7853045B2 (en) * | 2007-10-31 | 2010-12-14 | Saudi Arabian Oil Company | Geostatistical analysis and classification of core data |
US8447237B2 (en) * | 2010-04-12 | 2013-05-21 | Telefonaktiebolaget L M Ericsson (Publ) | Interference avoidance in white space communication systems |
CN104778331B (zh) * | 2015-04-24 | 2017-12-05 | 浙江工业大学 | 一种大跨桥梁监测数据空间插值方法 |
CN107449827B (zh) * | 2017-08-22 | 2020-02-18 | 中国水利水电科学研究院 | 一种堆石混凝土结构层面质量检测方法及评价方法 |
-
2018
- 2018-05-02 CN CN201810410957.1A patent/CN108629111B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021853A (zh) * | 2016-05-09 | 2016-10-12 | 青岛理工大学 | 一种空间变异特性下边坡稳定可靠度分析软件的开发方法 |
CN106442153A (zh) * | 2016-09-23 | 2017-02-22 | 河海大学 | 模拟混凝土重力坝坝踵水力劈裂的试验装置及试验方法 |
CN107862146A (zh) * | 2017-11-15 | 2018-03-30 | 中国地质大学(武汉) | 序贯高斯模拟岩土参数空间变异结构的有限元分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108629111A (zh) | 2018-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108629111B (zh) | 一种混凝土重力坝材料参数空间变异性的模拟方法 | |
CN205958224U (zh) | 一种采用振动格栅实时调控湍流度的风洞试验装置 | |
CN110046400A (zh) | 一种基于h-p型有限元法对混凝土重力坝进行静力特性分析的高精度数值模拟计算方法 | |
CN105631101B (zh) | 一种地震作用下高耸升船机塔柱结构动态分布系数分析方法 | |
CN110765695B (zh) | 一种基于高阶有限元法获取混凝土重力坝裂纹扩展路径的模拟计算方法 | |
CN110929443B (zh) | 一种基于高精度地形概化的二维洪水模拟方法 | |
CN102799713A (zh) | 堆石坝心墙水力劈裂的数值模拟方法 | |
KR200467341Y1 (ko) | 다단으로 분리한 태양광발전장치의 풍압저감 지지구조 | |
CN107563064A (zh) | 一种海啸波越浪过程的二维数值模拟方法 | |
CN107237321B (zh) | 径向小高度加翼大直径钢管桩及其水平承载力计算方法 | |
CN108897940A (zh) | 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 | |
CN107992715B (zh) | 一种基于虚拟边界力法的弧板式防波堤受力计算方法 | |
CN108842712A (zh) | 利用滚水坝促进过饱和do释放的方法、实验装置及释放规律的研究方法 | |
CN103823948B (zh) | 一种机场水泥混凝土道面板防翘曲变形设计方法 | |
CN115906682A (zh) | 一种基于水动力模型的鱼礁设计优化方法 | |
Major et al. | Numerical static analysis of the curtain wall with light steel structure | |
CN107657127A (zh) | 基于频域等能量插值的一维多变量随机过程高效模拟方法 | |
CN103643652A (zh) | 一种河工动床模型试验加糙方法及加糙载体 | |
CN110188412A (zh) | 大坝施工截流钢筋石笼起动流速试验及流场模拟方法 | |
CN110750865A (zh) | 一种钢板筒仓群仓优化布局方法及*** | |
Farouk et al. | Effect of elastic soil structure interaction on modulus of subgrade reaction | |
CN114003995B (zh) | 一种m型人工鱼礁空间布局优化方法 | |
Wong et al. | Pollutant dispersion over two-dimensional idealized street canyons: a large-eddy simulation approach | |
Cunhong et al. | New development in the numerical simulation of the tidal bore | |
Zang et al. | Influence of the Scorer parameter profile on the wavelength of trapped lee waves |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |