工程可靠分析.docVIP

工程结构是由钢、木、砖石、混凝土及钢筋混凝土等建造的各种建筑物和构筑物。工程结构在相当长的使用期内需要承受设备、人群、车辆等使用荷载，经受风、雪、冰、雨、日照或波浪、水流、土压力、地震等环境荷载的作用，它们安全可靠与否，不仅影响工农业生产，而且还常常关系到人身安危。20世纪80年代以来，北美、欧洲和亚洲的许多国家和地区，相继发生了一些大型桥梁、高层建筑和海上石油钻井平台的突然性断裂和倾覆事件，这些灾难性事故不仅引起了公众舆论的高度关注，也使各国政府和科研机构意识到，对大型重要工程结构的可靠性状况进行分析和评估已刻不容缓。 大型结构的体系可靠度分析过程中，伴随着大量的构件可靠度分析，而一般情况下，实际工程中的结构构造非常复杂，既使对其进行确定性分析，都需要借助于有限元等数值分工具。在这种情况下，进行构件可靠度分析时，常不能给出功能函数的明确表达式，因此对大型结构的体系可靠度分析，必须首先解决隐式功能函数的可靠度分析问题。对于功能函数不能直接明确表达的结构可靠度问题，可以用蒙特卡罗法结合有限元方法加以解决，也可以用随机有限元或响应面法进行求解。 蒙特卡罗法(Monte-Carlo)是一种随机模拟方法，被广泛用来处理具有随机性的问题。此方法易于被人们理解、接受，并有概率论中的中心极限定理来保证其可靠性与准确性，因而很自然地被引入到以概率论为基础的结构可靠度理论中，用来计算结构的失效概率。 蒙特卡罗法的优点是:在结构可靠度的数值模拟中，模拟的收敛速度与基本变量的维数、极限状态函数的复杂程度以及模拟过程无关;并且不管随机变量服从正态分布还是其它分布，不管是构件可靠度问题还是体系可靠度问题，其求解过程同样简单。蒙特卡罗法由于具有这些优越特性，因而引起了人们的广泛关注。但是，蒙特卡罗法需要成千上万次的模拟，对于隐式功能函数问题，也就需要成千上万次的有限元分析计算，工作量大，很不经济，因此蒙特卡罗法常被用于检验和校核各种近似方法的计算精度，很少用于大型结构的可靠度分析。 由概率定义可知，某事件的发生概率可以用大量试验中该事件发生的频率来估算。因此在结构可靠度计算中，可以通过对随机变量进行大量的随机抽样，然后把这些抽样值一组一组地代入结构功能函数，根据计算得到的功能函数值，确定结构安全与否，最后根据事件发生的次数，计算结构的可靠度或失效概率。 用Monte Carlo法分析结构可靠度问题，关键是产生已知分布变量的随机数。首先在区间(0，1)内产生均匀分布随机数，然后再变换成给定分布变量的随机数，通常在区间(0，1)内产生随机数的方法一般有数学方法和物理方法两种。产生随机数有很多方法，例如取中法、加同余法、乘同余法、混合同余法和组合同余法等等。乘同余法以其统计性质好、周期长等优点而被人们广泛采用。乘同余法的算式为 : (1) 式中a、c和m为正整数。表示以m为模数的同余式，+1为以m除(+c)后得到的余数。具体计算时，可引入参数，即令 (2) 式中lnt为取整符号，式(1)成为 (3) 若将+1除以模数m，可得标准化的随机数+1 ，即 (4) 用Monte Carlo法计算结构的失效概率时，需要进行大批抽样，需要确定最低抽样数N。显然，抽样数N与计算结果的精度有关。有文献建议，以95％的置信度保证Monte Carlo法计算精度，其允许误差ε为： 由上式可见，结构计算次数N越大，误差ε越小。因此，要达到一定的精度，N必须取得足够大.为简便起见，有文献建议N必须满足： N ≥100 / (5) 式中Pf为预先估计的失效概率。在工程实际中，一般要求失效概率Pf很小，这就需要计算的次数很多。例如，工程结构的失效概率一般在0.1％以下，要求计算的次数须达十万次以上。这个要求使实际计算不是遇到困难，就是耗时过多，Monte Carlo法的使用也因此而受到限制。 随机有限元法(Stochastic Finite Element Method，简记为SFEM)是一种将有限元与概率论相结合的可靠度分析方法。在各类工程结构中，存在着很多不确定性因素的影响。这些因素可以采用随机变量、随机过程或随机场来描述。随机分析理论与传统的有限元法相结合而产生的随机有限元法，是分析结构可靠性的最有效的