ansys的可靠度分析

1、ansys分析可靠度2007-11-11 10:29:41|分类：Ansys特辑|举报|字号订阅关于ansys分析可靠度的问题,他有两种方法：monte-carlo和响应面法。在现在的可靠度分析中monte-carlo法有中心点抽样法、直接重要抽样法、更新重要抽样法、渐进重要抽样法、方向抽样法，这里的中心点抽样法是最古老、效率最低的一种，但ansys里只有这一种方法，只是在抽样选点时有不同的两种选择； 并且，monte-carlo在工程计算中只用于校合，不能用于工程实践； 中心点抽样法在计算中一般要进行计算次数的讨论：当可靠指标为1.0时，失效概率1.5866E-01；当可靠指标为2.0时，失

2、效概率2.275E-02；当可靠指标为3.0时，失效概率1.3499E-03；当可靠指标为4.0时，失效概率3.1671E-05；一般结构的可靠指标为2-4，假设计算结构的可靠指标为3.0，此时的最少有限元计算次数为1/1.3499E-03（由于在计算过程中的多维变量随机选点不理想等原因，实际的计算次数远大于此），这对于写论文还可以，对于实际复杂的体系可靠度而言，是没法完成的； 下面我们来讨论一下ansys响应面法以及构件可靠度和体系可靠度： 响应面法计算可靠度不需要monte-carlo那么多次的有限元计算，对于构件可靠度他是现在一个很热门的研究方法，但是，对于体系可靠度，他没有考虑体系可靠

3、度的失效模式； 现在对于体系可靠度有两种认识：一种认为体系可靠度是由构件可靠度构成的，只有先知道构件可靠度，才能知道体系可靠度，要知道体系失效，先知道构件失效及其失效路径，在这方面大连理工大学的许林博士和张小庆博士开发了一套程序（程序思想是以上面的体系可靠度的认识为理论基础），程序的流程如下： 利用经过二次开发生成的新的ANSYS，进行可靠度计算的具体运算过程为：1) 利用APDL建立结构分析文件和优化文件；2) 运行ANSYS的批处理方式，利用分析文件建立模型、进行结构分析与敏度分析；3) 进入用户优化模块完成可靠度分析的一次迭代过程；4) 重新利用分析文件建立模型、进行结构分析与敏度分析；

4、5) 根据结构分析函数值和敏度值，以及前一点的结构分析函数值，用前面介绍的近似曲面构造法寻求拟合误差最小的近似极限状态函数；6) 对上一步得到的近似函数进行可靠度分析；7) 比较两次计算结果收敛与否，是则结束迭代，否则转到第4步，进行下一轮迭代。下面我要说的是另外一种体系可靠度思想，就是认为体系可靠度是由构件组成，一个（无论那一个构件）失效，体系就失效，不考虑失效路径，ansys是可以分析这种体系的可靠度的。现阶段已经开发的可靠度计算软件及相关工作如下：4.1结构可靠度与有限元有限元法作为一种数值分析方法，其萌芽可追溯到18世纪的Euler，然而，只是在电子计算机开始普及使用后，它才真正受到人

5、们的重视，并以令人惊讶的速度迅速发展，至今已成为计算力学在上世纪最重要和最辉煌发展的一个方面。在实践上，从50年代中期飞机薄壁结构分析的初期应用开始，至今已扩展到几乎所有的科学技术领域，对推动当代科学技术和工程实践的发展，已起到并将继续发挥难以估量的作用。计算力学极大地增强了经典力学解决自然科学和工程问题的计算能力，实现了大量复杂力学问题的数值求解，扩展了力学研究的领域，逐渐成为与试验、理论并列的力学研究三大支柱之一；其次，它的发展极大地提高了力学为工程服务的能力，改变了整个工程设计的面貌。它不仅使许多过去无法实现的复杂工程分析成为现实，而且可以采用优化设计的方法能动地优选设计方案，提高设计水

6、平和产品性能，缩短设计周期，将力学与工程更紧密地联系在一起。这里谈及的力学中行之有效的一些数值计算方法均是以确定性方法进行分析的，也就是没有考虑实际问题中的不确定因素。然而不计及不确定因素的分析方法是无法正确地判断结构安全可靠度的。因此把考虑不确定因素的方法引入到有效的结构数值分析方法中自然就具有非常大的意义。把概率分析部分整合到有限元分析中，求解结构响应的概率分布特性及可靠度的方法有两种。一种是随机有限元法（Stochastic Finite Element Method）。它在随机场概念的基础上，把材料特性的变异性融入到有限元分析过程中，分析这些随机性对结果的影响情况。早期在有限元中处理随

7、机变异性的方法是摄动法（Hasselmen & Hart,19721; Vanmarcke & Grigoriu,19832;陈虬，刘先斌，19933；刘宁，20014。Vanmarcke et al.(1986)5综述了随机场相关的概念，以及摄动法建立随机有限元的基本列式，回顾了一系列相关的问题。张圣坤等人结合船舶结构的特点分析了用随机有限元方法计算可靠度的基本过程。采用Taylor技术展开只能得到小随机问题响应的前两阶摄动解，而且数值解不太稳定，Yamazaki et al.(1988)6利用Neumann级数展开推导出计算响应变异的有限元格式，提高了结果的精度和适应性。这种展开方法的缺点

8、是它需要结合Monte Carlo抽样计算响应的前两阶统计矩，而且更高阶矩的计算仍然非常困难（Spanos & Ghanem，19897）。为了克服Monte Carlo法需要大量抽样的缺陷，Spanos & Ghanem(1989)8提出了基于Galekin形式的改进Neumann级数展开,并且在实际工程中加以应用。采用这一方法的软件目前比较缺少。另外一种是把概率分析部分和现成的有限元程序相对独立地连接起来。这种做法有一个明显的优点，它接利用经过实践考验的已经较为完善的有限元软件，从而计算各种结构的可靠度问题。但其缺点是计算花费比较大，而且通常不能处理随机场问题。因为有限元分析部分一般作为基

9、本独立的模块，仅为概率计算提供功能函数值，功能函数对随机变量的梯度需要采用有限差分的方法得到。Maymon（1994）9在商用软件ANSYS平台上利用它的APDL语言成功的实现了利用确定性分析软件计算验算点位置。Borri & Speranzini(1997)10同样在ANSYS平台上利用它内含的优化模块计算验算点位置，继而分析结构的可靠度。Bucher et al.11，12利用ANSYS 的结构分析功能，结合了响应面分析、Monte Carlo抽样和方向抽样方法进行结构可靠度计算，开发了SLang（CAL-Reliability）在有限元程序FEAP上集成了一次可靠度算法、二次可靠度算法、

10、Monte Carlo模拟法。他既可以计算元件可靠度，也能够分析系统可靠度。STRUREL包含了STATREL、COMREL、SYSREL、COSTREL、PERMAS-RA。STATREL专门对数据进行统计分析、处理时间序列问题；COMREL可以使用多种方法对元件进行时变/时不变可靠度分析；SYSREL是系统可靠度分析模块；COSTREL是一个集成了可靠度计算和优化分析的模块，采用的算法是单层次优化方法，可以分析时变/时不变可靠度约束的优化问题；PERMAS-RA是一个有限元分析程序，它可以结合响应面分析。COSSAN（Computational Stochastic Structural

11、Analysis）提供了响应面分析和各种Monte Carlo抽样方法（重要抽样、自适应抽样），在分析方面它借助外界的有限元分析程序。ISPUD（Importance Sampling Procedure Using Design Points）是一个专门进行抽样模拟的程序。PROBAN（Probabilistic Analysis）是Veritas Sesam Systems公司推出的海洋工程领域著名概率分析软件，它的功能强大，但是价格也比较昂贵。它可以计算元件和系统的可靠度及其灵敏度、确定响应的概率分布，主要采用的分析方法有一次/二次可靠度算法、Monte Carlo模拟法、拉丁超立方抽样

12、和其他一些抽样方法。NESSUS可以计算累积概率分布函数、失效概率、结构可靠度、系统可靠度、故障树分析。COMPASS（Computer Methods for Probabilistic Analysis of Structures and System）是Martec公司开发的一个随机系统可靠度和风险分析软件，还可以计算疲劳累计损伤、随机断裂力学，建立复合失效准则。国内的结构可靠度基础研究工作也取得了可喜的进步，但在软件开发和使用方面仍然停留在研究阶段。国内西北工业大学曾经开发了一个结构静强度可靠度计算程序，其中包含了杆单元、四边形/三角形板单元、空间梁单元和三角形壳单元。元件失效定义为单

13、元形心处的应力超过许用应力，体系失效定义为刚度阵奇异、结构变形超过限值。浙江大学的金伟良等基于海洋平台结构分析软件SACS开发了具有针对性的可靠度计算程序，可以进行FORM、SORM和重要抽样法计算可靠度。 众所周知，在体系可靠度中，关于失效模式的可靠度计算存在两种观点：第一种观点：体系可靠度的含义是从概率意义上度量结构作为一个整体抵抗某一种外部作用的能力，它的功能函数可以写成一个隐函数形式，这种情况下任一种函数近似方法都可以使用，响应面法也不例外；第二种观点：体系可靠度强调组成一个系统的各元件以及它们之间的相关关系，认为只有首先得出各元件可靠度和它们间的相关系数才可以计算系统的可靠度。对于第

14、二种方法，现阶段的工作主要集中在关于失效模式寻找这一理论探索上，这方面的大型结构可靠度计算还处于初步阶段（没看这方面的文献，此处可能有冒昧之处），主要的问题也是失效模式的判断很难确定。 我计算的是桥梁可靠度，是以第一种体系可靠度的理论为基础的（利用ANSYS的单元生死法模拟各施工过程，本文采用极限状态设计法大桥进行了静力计算，为实现大型结构的体系可靠度计算，在后处理中运用ANSYS的参数化语言做二次开发，使各荷载步的计算结果实现了累加、排序等处理，得到桥梁的最后计算结果，为可靠度计算提供了静力方面的数据，并为可靠度主要失效模式的判断提供了依据。可靠度思想概括如下：运用ANSYS参数化语言编写可

15、靠度程序计算大型复杂结构，依据静力计算和经验判定结构的主要失效模式，利用其中的响应面分析法使结构不同失效模式的隐式极限状态方程显式化，结合一次二阶矩法得到结构体系可靠度，这种分析方法称为混合分析法。计算表明，它是一种快捷有效的方法，在工程实践中具有很强的实用价值）。我说的“采用几何法作可靠度分析”，并非利用ansys的可靠性模块构造响应面函数，进而求解可靠度。而是用梯度优化法（在结构可靠性分析及随机有限元法（武清玺）一书中有详细的推导）求解可靠度，采用APDL实现这一过程。其本质与一次二阶矩法相同，采用线形近似的方法计算可靠性指标，只是迭代计算的格式不同而已。至于变量的分析类型问题，我们更多简

16、单地处理为正态分布，当然，如果不是正态分步的话，需要在迭代计算之前对验算点作当量正态化处理，这个我想应该不是什么难点。至于系统可靠度，我主要做的是机械传动系统的分析，目前机械系统方面的可靠性分析不如结构体系可靠度作得那么深入，做了很多简化。其失效模式的寻找是个头痛的问题，是个真正难点。我现在做的主要是对系统作故障树分析，很多东西交给软件（relex）去做。你的方法行不通:1.ANSYS的有限元是确定性有限元;2.ANSYS的MONTE-CARLO法分析桥梁,对于一个大型桥梁,单元达到2000个,一个奔四3.0的电脑,一次循环至少要1-2分钟,2万次有限元计算就是要花费300小时,这样的计算到了

17、一定的时候计算机C盘的空间不够(大量的有限元数据需要处理),电脑自动死机,其实你可以用APDL语言做二次开发,用MONTE-CARLO一般重要抽样法、MONTE-CARLO一般重要方向抽样法、响应面法都不需要过多的循环就能达到相当好的结果基于ANSYS的结构可靠性分析-叶勇2011-01-14 23:45:29|分类：默认分类|举报|字号订阅 1。ansys的概率分析功能ansys是一款强大的有限元分析软件，其提供的概率分析功能可以解决以下问题，根据模型中，输入参数的布确定性计算模型中结果变量的不确定性，确定由于输入参数的不确定性导致结构失效概率的值，已知结构的失效概率确定结构容许的变形范围；

18、判断对输出结果和失效概率影响较大的参数，计算输出结果相对于输入参数的灵敏度；确定输入变量和输出变量之间的灵敏度。1,1 概率分析方法在ANSYS中进行结构可靠度分析时，通常有以下几步：生成分析文件、可靠度分析阶段和结构的后处理，首先建立结构的循环分析文件，它包括预处理模块、求解模块和结构提出模块等，在可靠性分析阶段主要包括，指定可靠性分析文件，选择和定义输入和输出变量之间的相关函数，确定各输入变量服从的分布类型和函数，选择分析工具和方法（响应面法和Montre-carle法），而后处理阶段包括抽样过程显示，绘制设计变量取值分布函数，绘制失效概率分布函数，确定结构可靠度分析中的相关系数的矩阵，假

19、定已知结构的失效概率寻求对应的输入变量、灵敏度分析、生成分析报告等。2,。具体实例（来自ansys使用帮助）桁架结构如图所示：各节点的坐标如下：1（0，0，0），2（10，0，0），3（20，0，0），4（10，-10，0）（单位：cm）。1、2、3节点完全固定，节点4受集中力作用，FX20000N，FY20000N。杨氏模量EX2.1e5N/cm2。单元1的横截面积A1N（10，0.5），单元2的横截面积A2服从三角分布，Xmin=10，Xmlv=11，Xmax=12；单元3的横截面积A3服从均匀分布，Xmin=9，Xmax=11；A1、A3的相关系数为0.2。试分析杆件横截面积的不确定性对

20、结构体积、杆单元轴向应力的影响。解析：本题采用link1单元，直接构件有限元模型求解，提取结果体积和个单元轴向力作为可靠度分析输出结果变量，选择各杆横截面面积作为输入变量，采用monter-carle直接抽样法进行分析，命令流如下：*creat, pds3bar*set,A1,10*set,A2,10*set,A3,10/PREP7ET, 1, 1 ！定义单元类型，et,itype,ename(单元名或标号)，kop1,kop2,inoprMP,EX, 1, 2.1E5 !定义材料属性，MP,lab(ex,nuxy,damp等)，mat(材料参考号)，C0，C1，C2,R, 1, A1 !定义

21、实常数，R，NSET(实常数组号),R1,（各类单元有不同的实常数值）R, 2, A2 R, 3, A3 N, 1 N, 2, 10N, 3, 20N, 4, 10, -10!生成有限元模型REAL, 1E, 1, 4 ! E,I,J,K,L,M,O,P（通过节点定义单元）REAL, 2E, 2, 4REAL, 3E, 3, 4FINISH！荷载施加/SOLUD, 1, ALL, 3F, 4, FX, 20000F, 4, FY, -20000SOLVEFINISH/POST1SET, 1ETABLE, VOLU, VOLU !将单元的体积放入表VOLU中 ！ETABLE,LAB(用户定义的单

22、元名)，iTEM,COMPETABle, AXST, LS, 1 !将单元轴向应力放入表AXST中*GET, SIG1, ELEM, 1, ETAB, AXST !SIG1单元1的轴向应力*GET, SIG2, ELEM, 2, ETAB, AXST ！*GET,Par（结果数据的参数名）,Entity（NODE, ELEM, KP, LINE, AREA, VOLU, PDS）,ENTNUM（单元的实体标签）,Item1,IT1NUM,Item2,IT2NUM*GET, SIG3, ELEM, 3, ETAB, AXSTSSUM !将单元表格内数据求和*GET,TVOL, SSUM, ITE

23、M, VOLU!提取结构的总体积FINISH*end !结束宏定义*use, pds3bar !执行宏操作/PDS !进入可靠性分析模块pdanl, pds3bar !指定分析文件!定义输入变量PDVAR, A1, GAUS, 10,0.5 !PDVAR,Name,Type(参数的概率分布类型，beta、expo、log1、weib、unif、gama、gaus),PAR1,PAR2,PAR3,PAR4（对应于不同的分布类型，各参数的意义不同)PDVAR, A2, tria, 10,11,12 ！tria三角分布，PAR1= Minimum value. Defaults to 0.0.PAR

24、2= Most Likely ！Value (MLV). Defaults to 0.5.PAR3= Maximum value. Defaults to 1.0.PAR4is ignored.PDVAR, A3, unif, 9, 11PDCOR, A1, A3, 0.2 ！pdcor,Name1,Name2,CORR ,定义name1和那么2的相关系数。!定义输出变量PDVAR, SIG1, respPDVAR, SIG2, respPDVAR, SIG3, respPDVAR, TVOL, respPDMETH, MCS, DIR !指定分析方法，mcs-monter-carle法，rs

25、m-响应面法，dir-直接或者粗糙的montercalrePDDMCS, 100, NONE, ALL, !设定分析循环次数！PDDMCS,NSIM（分析循环的次数）,-,Astop（auto、all）,ACCMEAN（当为all时缺省，一般缺省为0.01，）,ACCSTDEV（缺省时为0.02）,CHECK（一般要大于1，缺省时为10）,Seed（推荐使用）PDEXE, mcs3bar !执行分析,Slab（执行分析的文件名）,MRUN（Keyword for the processing of the simulation loopsSERExecutes the simulation l

26、oops of the probabilistic analysis using a single machine; serial processing (default).PARExecutes the simulation loops of the probabilistic analysis in parallel using the ANSYS parallel-processing tool.）,NFAIL（This parameter sets the number of loops that are allowed to fail before the execution of the probabilistic analysis is terminated. For example, if geometric parameters are used as random variables then for some simulations (loops) the loop will fail because meshing was not successful. This might happen infrequently which might be acceptable.