本发明属于模拟方法领域,具体是一种混凝土重力坝材料参数空间变异性的模拟方法。
背景技术:
混凝土重力坝是水利水电工程拦河大坝常用的一种坝型。混凝土是一种由砂石骨料和水泥浆搅拌混合而成的非均质材料。在修建混凝土重力坝的过程中,一般采用分段、分块方式浇筑混凝土。由于浇筑仓面大、浇筑施工质量不均匀,以及混凝土材料自身的非均质性,因此,混凝土材料的物理力学性质具有空间变异性,即对于同类混凝土,其在坝体不同位置处的物理力学参数是不一样的。混凝土材料参数的空间变异性对重力坝的应力变形和结构安全会产生影响,因此,有必要研究材料参数空间变异性对混凝土重力坝结构特性的影响。这就需要解决混凝土重力坝结构分析中材料参数空间变异性的模拟技术问题。
目前混凝土重力坝结构分析中材料参数空间变异性的模拟方法为:将坝体断面离散为有限个网格单元,通过给网格单元赋予不同的材料参数来模拟材料参数的空间变异性,并假定不同网格单元的材料参数服从某种概率分布(如weibull分布、正态分布等),这种方法反映了材料参数在空间中的随机分布,但没有考虑材料参数的空间相关性,即相邻位置(或相邻网格单元)的材料参数具有某种程度的正相关性,也就是相邻位置(或相邻网格单元)的材料参数不会有太大差别。目前的模拟方法没有考虑这种空间相关性,从而导致相邻网格单元的材料参数出现较大差异的不合理情况。
技术实现要素:
本发明的目的在于解决上述现有技术中存在的问题,提供一种混凝土重力坝材料参数空间变异性的模拟方法,所述模拟方法更加符合工程实际,能够提供更加准确的数据。
为实现上述发明目的,本发明所采用的技术方案是:
一种混凝土重力坝材料参数空间变异性的模拟方法,包括以下步骤:
(1)建立重力坝模型,对所述模型进行网格划分,获取每个网格单元的编号和形心坐标,组成矩阵
(2)利用matlab软件生成随机矩阵am×n,其中m为每个网格单元参数的随机抽样次数,n为网格单元的个数;
(3)利用自相关函数计算所述矩阵en×3的相关系数矩阵
(4)对所述矩阵ρn×n进行cholesky分解,即ρn×n=bn×nt·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)中所述自相关函数为负指数自相关函数,所述元素
优选的,步骤(6)中所述满足对数正态分布的参数μlnx和σlnx通过下式计算,
本发明的有益效果是:本发明提供的混凝土重力坝材料参数空间变异性的模拟方法基于随机场理论,采用自相关函数考虑混凝土材料参数的空间相关性,提出一种基于空间变异性且满足已知相关关系的相关随机场的模拟方法;本发明所述模拟方法更加符合工程实际,能够在后期结果中为重力坝的安全运行提供更加准确的数据支撑;该模拟方法能够用于重力坝有限元分析中材料参数空间变异性的模拟,为考虑材料参数空间变异性的重力坝结构分析提供技术支持。
附图说明
图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,根据公式
4.对所述矩阵ρ4957×4957进行cholesky分解,即
ρ4957×4957=bt4957×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,则抗拉强度随机场
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。