《ANSYS概率设计PDS讲义》由会员分享,可在线阅读,更多相关《ANSYS概率设计PDS讲义(47页珍藏版)》请在金锄头文库上搜索。
1、概率设计概率设计1概率设计概率设计 前言:可靠度基本理论前言:可靠度基本理论第一节:基于有限元的概率设计技术第一节:基于有限元的概率设计技术 1.1 基于有限元的概率设计基于有限元的概率设计(PDS)简介简介 1.2 PDS的基本概率与过程数据流的基本概率与过程数据流 1.3 PDS中的参数分布函数及其选用中的参数分布函数及其选用 1.4 Monte Carlo法法 1.5 响应面法响应面法第二节:基于有限元的概率设计基本过程第二节:基于有限元的概率设计基本过程 2.1 创建分析文件创建分析文件 2.2 初始化概率设计分析及参数初始化概率设计分析及参数 2.3 进入进入PDS并指定分析文件并指
2、定分析文件 2.4 定义概率设计模型定义概率设计模型 2.5 选择概率设计方法或工具选择概率设计方法或工具 2.6 执行概率设计分析执行概率设计分析 2.7 拟合和使用响应面拟合和使用响应面 2.8 概率设计结果后处理概率设计结果后处理第三节:概率设计分析的实例第三节:概率设计分析的实例 承受横向集中力板的承受横向集中力板的LHS抽样抽样MCS概率设计实例概率设计实例 三根杆桁架系统的直接抽样三根杆桁架系统的直接抽样MCS概率分析实例概率分析实例2l结构的极限状态:整个结构的一部分超过某一特定状态就不能满足设计规定结构的极限状态:整个结构的一部分超过某一特定状态就不能满足设计规定的某一功能要求
3、。结构的极限状态实质上是结构工作状态的一个阀值,如果的某一功能要求。结构的极限状态实质上是结构工作状态的一个阀值,如果工作状态超过这一阀值,则结构处于不安全、不耐久或不适用的状态;若工工作状态超过这一阀值,则结构处于不安全、不耐久或不适用的状态;若工作状态没超过这一阀值,则结构处于安全、耐久、适用的状态作状态没超过这一阀值,则结构处于安全、耐久、适用的状态l用用 表示结构的工作状态,称作结构的功能函数。则结构的工作状态可表示结构的工作状态,称作结构的功能函数。则结构的工作状态可表示为:表示为: 结构的工作状态结构的工作状态可靠度基本理论可靠度基本理论31. 用定义计算结构可靠度2.用统计分析计
4、算可靠度式中,试验的总次数;实验中的次数。 结构的工作状态结构的工作状态可靠度基本理论可靠度基本理论41.1 基于有限元的概率设计基于有限元的概率设计(PDS)简介简介利用概率设计方法可以帮助用户确定利用概率设计方法可以帮助用户确定“失效失效”情况发生的可能性,这情况发生的可能性,这样就使得用户可以改进设计直到满足用户可以接受的样就使得用户可以改进设计直到满足用户可以接受的“极限极限”即可。即可。概率设计技术是用来评估输入参数的不确定性对于系统响应的影响概率设计技术是用来评估输入参数的不确定性对于系统响应的影响行为及其特性。行为及其特性。输入参数包括几何尺寸、加工误差、材料、载荷等不确定因素。
5、输入参数包括几何尺寸、加工误差、材料、载荷等不确定因素。响应参数包括温度、应力、位移等。响应参数包括温度、应力、位移等。有限元分析技术与概率设计技术相结合,就是基于有限元的概率设有限元分析技术与概率设计技术相结合,就是基于有限元的概率设计,即计,即ANSYS程序提供的程序提供的PDS技术(技术(Probabilistic Design System).51.1 基于有限元的概率设计基于有限元的概率设计(PDS)简介简介当有限元模型的输入参数不确定时,有限元结当有限元模型的输入参数不确定时,有限元结果的不确定程度有多大?响应参数的置信度有果的不确定程度有多大?响应参数的置信度有多高?多高?输入参
6、数的不确定性决定响应参数的不确定性,输入参数的不确定性决定响应参数的不确定性,目标产品满足设计要求的概率有多大?工作失目标产品满足设计要求的概率有多大?工作失效概率有多大?效概率有多大?在所有不确定的输入参数中哪个参数的不确定在所有不确定的输入参数中哪个参数的不确定性对于响应参数的影响程度最大,或者说对于性对于响应参数的影响程度最大,或者说对于目标产品最容易引起其工作失效?响应参数对目标产品最容易引起其工作失效?响应参数对输入参数变化的灵敏度多大?输入参数变化的灵敏度多大?ANSYS提供的基于有限元的概率设计系统(提供的基于有限元的概率设计系统(PDS)的主要应用方向:)的主要应用方向:6 1
7、.2 PDS的基本概率与过程数据流的基本概率与过程数据流随机输入参数随机输入参数(RVsrandom input variables ) 又称设计驱动参数,直接影响分析结果,需指定分布类型以特征参数又称设计驱动参数,直接影响分析结果,需指定分布类型以特征参数相关性相关性(Correlation) 指两个(或多个)随机输入参数之间存在统计上的关联性指两个(或多个)随机输入参数之间存在统计上的关联性随机输出变量随机输出变量(RPsrandom output parameters) 指有限元分析结果指有限元分析结果 RP是是RV的函数的函数概率设计参数概率设计参数 (probabilistic de
8、sign variables) RV和和RP统称为概率设计参数,在定义时必需指定统称为概率设计参数,在定义时必需指定样本(样本(Sample) 一个样本就是一序列确定的随机输入参数值一个样本就是一序列确定的随机输入参数值仿真(仿真(Simulation)分析文件(分析文件(Analysis file) 是一个是一个ANSYS输入文件,包含一个完整的分析过程,如前处理、求解和后处理等输入文件,包含一个完整的分析过程,如前处理、求解和后处理等 必须包含参数化自动建模的过程,所有输入和输出项,将可能被定义成随机输入参数和随机输出必须包含参数化自动建模的过程,所有输入和输出项,将可能被定义成随机输入参
9、数和随机输出参数参数基本概念:基本概念:基本概念:基本概念:7 1.2 PDS的基本概率与过程数据流的基本概率与过程数据流循环文件(循环文件(Loop file) *.loop文件,由文件,由ANSYS自动根据分析文件生成。自动根据分析文件生成。 利用该文件进行概率设计循环利用该文件进行概率设计循环概率设计模型(概率设计模型(Probabilistic model) 以分析文件形式存在,包括所有定义和设置:以分析文件形式存在,包括所有定义和设置:RVs、相关性、相关性、RPs、概率设计方法和相关、概率设计方法和相关 参数等参数等概率设计数据库(概率设计数据库(PDS database) 包括当
10、前设计的环境,包括包括当前设计的环境,包括RVs、相关性、相关性、RPs、概率设计方法、被执行的概率分析及存储其结果、概率设计方法、被执行的概率分析及存储其结果的各种文件、使用哪个概率设计分析中的哪个输出参数来拟合响应表面、拟合中所使用的回归模的各种文件、使用哪个概率设计分析中的哪个输出参数来拟合响应表面、拟合中所使用的回归模型、拟合结果等。型、拟合结果等。 可以被存储到,并且可重新读入。结果不存储在这个数据库中。拟合响应表面的样本即存储在数可以被存储到,并且可重新读入。结果不存储在这个数据库中。拟合响应表面的样本即存储在数据库中。据库中。均值(均值(Mean value)、中间值()、中间值
11、(Median value)、标准方差()、标准方差(Standard deviation)基本概念:基本概念:基本概念:基本概念:81.2 PDS的基本概率与过程数据流的基本概率与过程数据流ANSYS数据库文件SAVERESUM有限元模型数据库可靠性分析数据库PDSAVE可靠性分析数据库文件PDRESUM循环文件PDEXEPDEXE/EXIST分析文件数据流程数据流程数据流程数据流程91.3 PDS中的参数分布函数及其选用中的参数分布函数及其选用高斯分布高斯分布高斯分布高斯分布(GAUS)(GAUS)特征参数:特征参数:特征参数:特征参数: 均值均值均值均值 标准方差标准方差标准方差标准方差
12、 截断高斯分布截断高斯分布截断高斯分布截断高斯分布(TGAU)(TGAU)特征参数:特征参数:特征参数:特征参数: 均值均值均值均值 标准方差标准方差标准方差标准方差 截断下限截断下限截断下限截断下限XminXmin 截断上限截断上限截断上限截断上限XmaxXmax对数正态分布(对数正态分布(对数正态分布(对数正态分布(LOGLOG)特征参数:特征参数:特征参数:特征参数: 均值均值均值均值 标准方差标准方差标准方差标准方差 101.3 PDS中的参数分布函数及其选用中的参数分布函数及其选用三角分布三角分布三角分布三角分布(TRIA)(TRIA)特征参数:特征参数:特征参数:特征参数: 最小值
13、最小值最小值最小值XminXmin 可能值可能值可能值可能值XmivXmiv 可能值可能值可能值可能值XmaxXmax均匀分布均匀分布均匀分布均匀分布(UNIF)(UNIF)特征参数:特征参数:特征参数:特征参数: 截断下限截断下限截断下限截断下限XminXmin 截断上限截断上限截断上限截断上限XmaxXmax指数分布指数分布指数分布指数分布特征参数:特征参数:特征参数:特征参数: 衰减系数衰减系数衰减系数衰减系数 下限下限下限下限XminXmin111.3 PDS中的参数分布函数及其选用中的参数分布函数及其选用BETABETA分布分布分布分布(BETA)(BETA)特征参数:特征参数:特征
14、参数:特征参数: 形状参数形状参数形状参数形状参数r r 形状参数形状参数形状参数形状参数t t 下限下限下限下限XminXmin 上限上限上限上限XmaxXmax伽马分布伽马分布伽马分布伽马分布(GAMA)(GAMA)特征参数:特征参数:特征参数:特征参数: 衰减系数衰减系数衰减系数衰减系数 幂指数幂指数幂指数幂指数k k威布尔分布(威布尔分布(威布尔分布(威布尔分布(WEIBWEIB)特征参数:特征参数:特征参数:特征参数: 威布尔特征值威布尔特征值威布尔特征值威布尔特征值XchrXchr 威布尔指数威布尔指数威布尔指数威布尔指数mm 最小值最小值最小值最小值XminXmin12 1.4
15、Monte Carlo法法Monte Carlo法(Monte Carlo Method)又名随机模拟法或统计试验法 将随机变量将随机变量代入功能函数代入功能函数,得出一个函数值。若,得出一个函数值。若,则在计算机程序中记入一次失效的实现;若,则在计算机程序中记入一次失效的实现;若,则不记入。,则不记入。这样就完成了一次计算,再产生下一个随机数,重复上面的计算,直至完这样就完成了一次计算,再产生下一个随机数,重复上面的计算,直至完成预定的实验次数为止。此时,失效概率成预定的实验次数为止。此时,失效概率为 Monte Carlo法可选择法可选择直接抽样法直接抽样法、超拉丁方抽样超拉丁方抽样和用户
16、抽样处理。和用户抽样处理。13 1.4 Monte Carlo法法Monte Carlo法特点:法特点:Monte Carlo方法及其程序结构简单,较容易实现;方法及其程序结构简单,较容易实现;收敛的概率和收敛的速度与问题的维数无关;收敛的概率和收敛的速度与问题的维数无关;用模拟的方法计算结构系统的失效概率,不需考虑失效模式的相关用模拟的方法计算结构系统的失效概率,不需考虑失效模式的相关性;性;只要抽样次数足够多,该方法计算所得的结构可靠度的精度满足要只要抽样次数足够多,该方法计算所得的结构可靠度的精度满足要求,所以一般用来检验其他方法的计算结果。求,所以一般用来检验其他方法的计算结果。141
17、.4 Monte Carlo法法直接抽样直接抽样Monte Carlo模拟技术中最常用的基本方法,可直接模拟技术中最常用的基本方法,可直接用于模拟各种工程真实过程。可模拟零件在现实中任用于模拟各种工程真实过程。可模拟零件在现实中任何行为。何行为。效率不高,需做大量仿真循环。效率不高,需做大量仿真循环。对抽样过程没有对抽样过程没有“记忆记忆”功能,会出现重复抽样。功能,会出现重复抽样。需要指定随机输入参数的样本种值、仿真循环次数和需要指定随机输入参数的样本种值、仿真循环次数和循环终止准则(均值和标准方差精度等)循环终止准则(均值和标准方差精度等)151.4 Monte Carlo法法拉丁超立方抽
18、样(拉丁超立方抽样(LHS)比直接抽样法更先进、更有效。比直接抽样法更先进、更有效。对抽样过程有对抽样过程有“记忆记忆”功能,可避免直接抽样法数据点功能,可避免直接抽样法数据点集中而导致的仿真循环重复问题。集中而导致的仿真循环重复问题。强制抽样过程中抽样点必须离散分布于整个抽样空间。强制抽样过程中抽样点必须离散分布于整个抽样空间。LHS抽样法比直接抽样法要少抽样法比直接抽样法要少20%40%的仿真循环资的仿真循环资料。料。需要指定仿真循环次数、重复次数、样本分布位置、循需要指定仿真循环次数、重复次数、样本分布位置、循环终止准则(均值和标准方差精度等)和随机输入参数环终止准则(均值和标准方差精度
19、等)和随机输入参数样本种子值。样本种子值。161.5 响应面法响应面法比比Monte Carlo模拟需要的循环次数少。模拟需要的循环次数少。可进行非常低概率问题的分析。可进行非常低概率问题的分析。拟合系数表示近似函数与响应数值的近似程度拟合系数表示近似函数与响应数值的近似程度,可通过用户自己定义。可通过用户自己定义。单个循环之间相互独立,非常适用于并行计算。单个循环之间相互独立,非常适用于并行计算。循环次数取决于随机输入变量个数,变量数不能太多。循环次数取决于随机输入变量个数,变量数不能太多。要求输入变量与输入变量的函数平滑。如接触分析不可以用该方法。要求输入变量与输入变量的函数平滑。如接触分
20、析不可以用该方法。171.5 响应面法响应面法响应面法可选择三种方法:中心合成设计、响应面法可选择三种方法:中心合成设计、Box-Bchnken矩阵法和用户指定法。矩阵法和用户指定法。用数学函数(二次函数)表达随机输入变量用数学函数(二次函数)表达随机输入变量和随机输出变量之间的关系。使用回归分析和随机输出变量之间的关系。使用回归分析技术(通常是用最小二乘法)确定函数的各技术(通常是用最小二乘法)确定函数的各项系数。项系数。响应面法两个步骤:响应面法两个步骤:1,进行仿真循环计算,进行仿真循环计算对应随机输入变量空间样本点的随机输出变对应随机输入变量空间样本点的随机输出变量的数据;量的数据;2
21、,进行回归分析确定近似函数。,进行回归分析确定近似函数。181.5 响应面法响应面法中心合成设计抽样(中心合成设计抽样(中心合成设计抽样(中心合成设计抽样(CCDCCD) 中心合成设计抽样包括一个中心点、中心合成设计抽样包括一个中心点、N个轴线点和位于个轴线点和位于2N-f阶乘个阶乘个N维超立方体的顶点。维超立方体的顶点。式中,式中,N随机输入变量数目随机输入变量数目 f中心合成设计阶乘因子表达式中的一个参数。中心合成设计阶乘因子表达式中的一个参数。Box-BehnkenBox-Behnken矩阵抽样矩阵抽样矩阵抽样矩阵抽样 Box-Behnken矩阵抽样包括一个中心点矩阵抽样包括一个中心点N
22、维超立方体每边中维超立方体每边中心点。心点。19PDS基本过程基本过程ANSYS基于有限元的概率设计分析具体实现过程如下基于有限元的概率设计分析具体实现过程如下:2.1 创建分析文件创建分析文件.文件应该包括完整的分析过程,如:文件应该包括完整的分析过程,如:参数化有限元模型参数化有限元模型(PREP7);求解求解(SOLUTION);获取数据,用做随机输入参数和随机输出参数获取数据,用做随机输入参数和随机输出参数(POST1/POST26);2.2 建立概率有限元分析数据库和所有参数建立概率有限元分析数据库和所有参数2.3 进入进入PDS模块并指定分析文件模块并指定分析文件(PDS)。2.4
23、 定义随机输入变量和随机输出变量定义随机输入变量和随机输出变量(PDS)。2.5 选择概率设计工具或方法选择概率设计工具或方法(PDS)。2.6 执行概率设计分析所需要的循环执行概率设计分析所需要的循环(PDS)。2.7 拟合响应面(拟合响应面(PDS)2.8 观察概率设计结果(观察概率设计结果(PDS)。)。202.1 创建分析文件创建分析文件分析文件就是基于分析文件就是基于APDL的参数化有限元分析过程的参数化有限元分析过程创建分析文件方法:创建分析文件方法: (1)GUI方式,再经方式,再经LOG文件整理出来;文件整理出来; (2)在文本编辑器中直接编写。)在文本编辑器中直接编写。21具
24、体过程包括具体过程包括: (1)定义并初始化参数()定义并初始化参数(RV) *CREATE,FILENAME,MAC !创建宏文件!创建宏文件 如:如:A1=10 或用或用*SET,A1,10 (2)参数化创建有限元模型)参数化创建有限元模型 /PREP7 !进入前处理器!进入前处理器 ET,1,LINK1 !杆单元!杆单元 R,1,A1 !以!以RV为参数的实常数为参数的实常数 MP,EX,1,2.1E5 !定义材料定义材料 MP,PRXY,1,0.3 N,1,0,0,0 !创建节点创建节点N,2,10,0,0E,1,2 !创建单元创建单元FINISH !退出前处理器!退出前处理器(3)约
25、束、加载、求解)约束、加载、求解/SOLU !进入求解器!进入求解器D,1,ALL !约束!约束F,4,FY,-1000 !加载加载SOLVE !求解!求解FINISH !退出求解器!退出求解器2.1 创建分析文件创建分析文件22(续前页续前页)(4)提取结果数据并存储到参数中)提取结果数据并存储到参数中/POST1 !进入后处理器进入后处理器SET,FIRST !读入第一个结果序列读入第一个结果序列ETABLE,axst,LS,1 !将单元应力放入表将单元应力放入表AXST中中*GET,sig1,ELEM,1,ETAB,AXST !sig1=单元单元1的轴向应力的轴向应力*GET,sig2,
26、ELEM,2,ETAB,AXST *GET,sig3,ELEM,3,ETAB,AXST SSUM !将单元表格内数据求和将单元表格内数据求和*GET,TVOL,SSUM, ,ITEM,VOLU !提取结构总体积提取结构总体积FINISH !退出后处理器退出后处理器*END !完成宏定义!完成宏定义2.1 创建分析文件创建分析文件232. 初始化概率设计分析及参数初始化概率设计分析及参数GUI方式方式:Utility MenuFileRead Input from APDL方式:方式:/INP,FILENAME,MAC3.进入进入PDS并指定分析文件并指定分析文件GUI方式方式:Main Men
27、uProb DesignAnasis FileAssign APDL方式:方式: /PDS PDANL,FILENAME,MAC24(1)定义随机输入参数)定义随机输入参数GUI方式方式:Main MenuProb DesignProb Definitns Random Input APDL方式:方式: PDVAR, Name, Type, PAR1, PAR2, PAR3, PAR42.4 定义概率设计模型定义概率设计模型设置变量分布类型;设置变量分布类型;设置变量分布参数。设置变量分布参数。25绘制输入参数分布图绘制输入参数分布图GUI方式方式: Main MenuProb DesignP
28、rob DefinitnsPlot APDL方式:方式: PDPLOT, Name, PLOW, PUP2.4 定义概率设计模型定义概率设计模型26定义输入参数相关性定义输入参数相关性GUI方式方式: Main MenuProb DesignProb DefinitnsCorrelationAPDL方式:方式: PDCORR, Name1, Name2, CORR2.4 定义概率设计模型定义概率设计模型27(2)指定随机输出参数)指定随机输出参数GUI方式方式: Main MenuProb DesignRandom OutputAPDL方式:方式: PDVAR, Name, Resp, PAR
29、1, PAR2, PAR3, PAR4 2.4 定义概率设计模型定义概率设计模型选择输出变量。选择输出变量。 28(1)使用)使用PDS向导自动选择最合适的概率设计方法向导自动选择最合适的概率设计方法GUI方式方式: Main MenuProb DesignRandom MethodMethod Wizard2.5 选择概率设计方法或工具选择概率设计方法或工具 可以使用方法向导让可以使用方法向导让ANSYS程序帮助选择适合当前程序帮助选择适合当前问题的最有效随机设计方法。问题的最有效随机设计方法。 主要考虑因素包括:分析文件执行一次的时间、主要考虑因素包括:分析文件执行一次的时间、CPU数目、
30、随机输出参数响应特性、输出参数需要评估的项数目、随机输出参数响应特性、输出参数需要评估的项等。等。 29(1)使用)使用Monte Carlo法进行概率设计法进行概率设计GUI方式方式: Main MenuProb DesignProb MethodMonte Carlo Sims APDL方式:方式:PDMETH, MCS, DIR/LHS2.5 选择概率设计方法或工具选择概率设计方法或工具DIR直接抽样法直接抽样法设置仿真次数设置仿真次数设置终止条件设置终止条件设置随机输入参数的样本种子值设置随机输入参数的样本种子值LHS拉丁超立方抽样拉丁超立方抽样设置仿真次数设置仿真次数设置仿真重复次数
31、设置仿真重复次数指定区间内的样本分布位置指定区间内的样本分布位置设置终止条件设置终止条件设置随机输入参数的样本种子值设置随机输入参数的样本种子值30(2)使用响应面法进行概率设计)使用响应面法进行概率设计GUI方式方式:Main MenuProb DesignProb MethodResponse Surface APDL方式:方式:PDMETH,RSM, CCD/BBM2.5 选择概率设计方法或工具选择概率设计方法或工具设置概率下限值设置概率下限值设置概率上限值设置概率上限值31运行序列求解运行序列求解GUI方式方式:Main MenuProb DesignRunExec SerialRun
32、 Serial APDL方式:方式:PDEXE, Slab, SER, NFAIL, FOPT, Fname2.6 执行概率设计分析执行概率设计分析输入当前输入当前PDS求解结果的标识字符串求解结果的标识字符串可以选择保存样本文件可以选择保存样本文件32运行运行PDS并行分析并行分析GUI方式方式:Main MenuProb DesignRunExec Parallel APDL方式:方式:PDEXE, Slab, PAR, NFAIL, FOPT, Fname2.6 执行概率设计分析执行概率设计分析定义用来并行处理的远程主机定义用来并行处理的远程主机指定要复制到远程主机上的文件指定要复制到远
33、程主机上的文件开始并行分析开始并行分析332.7 拟合和使用响应面拟合和使用响应面运行求解完概率设计计算后,使用存储在结果文件中运行求解完概率设计计算后,使用存储在结果文件中的的PDS数据可以拟合响应表面。数据可以拟合响应表面。如果不能确定响应面法需要多少样本点,那么样本点如果不能确定响应面法需要多少样本点,那么样本点数至少是响应面系数的两倍。数至少是响应面系数的两倍。Monte Carlo法和响应面法模拟后都可以拟合响应面。法和响应面法模拟后都可以拟合响应面。342.7 拟合和使用响应面拟合和使用响应面拟合响应表面拟合响应表面GUI方式:方式:Main MenuProb DesignResp
34、onse SurfFit Resp SurfAPDL方式:方式:RSFIT, RSlab, Slab, Name, Rmod, Ytrans, Yval, Xfilt, CONF指定存储名称字符串指定存储名称字符串用于拟合响应面的用于拟合响应面的PDS计算文件计算文件选择要拟合响应面的随机输出参数对象选择要拟合响应面的随机输出参数对象选择拟合响应面的回归模型类别,包括:选择拟合响应面的回归模型类别,包括:线性、没有交叉项的二次函数和带交叉项线性、没有交叉项的二次函数和带交叉项的二次函数的二次函数响应面输出结果数据的数学变换处理。响应面输出结果数据的数学变换处理。确定在进行回归计算时是否过滤不相
35、关项确定在进行回归计算时是否过滤不相关项352.7 拟合和使用响应面拟合和使用响应面绘制响应表面绘制响应表面GUI方式:方式:Main MenuProb DesignResponse SurfPlt Resp SurfAPDL方式:方式:RSPLOT, RSlab, YName, X1Name, X2Name, Type, NPTS, PLOW, PUP 选择要拟合的随机输出参数对象选择要拟合的随机输出参数对象选择拟合响应面的第一个随机输入参数选择拟合响应面的第一个随机输入参数选择拟合响应面的第二个随机输入参数选择拟合响应面的第二个随机输入参数响应面显示方式:响应面显示方式:2D,3D各坐标轴
36、上的刻度点数各坐标轴上的刻度点数概率上下限概率上下限362.7 拟合和使用响应面拟合和使用响应面输出响应表面输出响应表面GUI方式:方式:Main MenuProb DesignResponse SurfPrn Resp SurfAPDL方式:方式:RSPRNT, RSlab, YName, Xout输出文件包括信息有:输出文件包括信息有: 使用的变换函数使用的变换函数 回归项回归项 回归系数回归系数 拟合逼近程度拟合逼近程度设置内容包括:设置内容包括:选择拟合得到的响应面函数选择拟合得到的响应面函数选择要拟合响应面的输出对象选择要拟合响应面的输出对象选择输出方式选择输出方式372.7 拟合和
37、使用响应面拟合和使用响应面根据响应表面生成根据响应表面生成Monte Carlo样本样本GUI方式:方式:Main MenuProb DesignResponse SurfRS SimulationAPDL方式:方式:RSSIMS, RSlab, NSIM, Seed用拟合得到的响应面函数代替有限元模型进行用拟合得到的响应面函数代替有限元模型进行Monte Carlo模拟,模拟效模拟,模拟效率高。率高。模拟次数及计算精度。模拟次数及计算精度。382.8 概率设计结果后处理概率设计结果后处理在完成在完成PDS分析之后,对分析之后,对PDS结果进行适当的后处理操作。结果进行适当的后处理操作。1.统
38、计分析统计分析 (statistics) 对单个变量的评估,允许对随机输出变量或随机输对单个变量的评估,允许对随机输出变量或随机输入变量做统计分析。入变量做统计分析。样本历史样本历史柱状图柱状图累积分布函数累积分布函数概率概率逆概率逆概率2.趋势分析趋势分析(trends) 对多个变量的评估。对多个变量的评估。灵敏度灵敏度散点图散点图相关矩阵相关矩阵3.自动报告自动报告(report)设置报告选项设置报告选项输出输出HTML报告报告392.8 概率设计结果后处理概率设计结果后处理样本历史样本历史查看循环结果和循环次数的函数关系,查看变量的分布参数。查看循环结果和循环次数的函数关系,查看变量的分
39、布参数。GUI方式:方式:Main MenuProb DesignProb ResultsStatisticsSampl HistoryAPDL方式方式:PDSHIS, Rlab, Name, Type, CONF选择抽样数据文件选择抽样数据文件选择要分析的变量选择要分析的变量选择图形显示方式(均值、标准方差、最值)选择图形显示方式(均值、标准方差、最值)指定置信度水平指定置信度水平402.8 概率设计结果后处理概率设计结果后处理直方图直方图查看变量的离散情况。查看变量的离散情况。如图抽样次数足够多,直方图比较光滑,没有大的间隙。如图抽样次数足够多,直方图比较光滑,没有大的间隙。GUI方式:方
40、式:Main MenuProb DesignProb ResultsStatisticsHistogramAPDL方式方式:PDHIST, Rlab, Name, NCL, Type选择抽样数据文件选择抽样数据文件选择要分析的变量选择要分析的变量指定绘制直方图数据点数指定绘制直方图数据点数指定柱状图类型。指定柱状图类型。412.8 概率设计结果后处理概率设计结果后处理累积分布函数累积分布函数查看可靠性或失效概率的工具。函数一点的数值等于数据出现在该点下的概率。查看可靠性或失效概率的工具。函数一点的数值等于数据出现在该点下的概率。GUI方式:方式:Main MenuProb DesignProb
41、 ResultsStatisticsCumulativeDFAPDL方式方式:PDCDF, Rlab, Name, Type, CONF, NMAX 选择抽样数据文件选择抽样数据文件选择要分析的变量选择要分析的变量指定柱状图类型指定柱状图类型设置置信度水平设置置信度水平绘制曲线数据点总数绘制曲线数据点总数422.8 概率设计结果后处理概率设计结果后处理输出概率输出概率查看在给定数据点的失效概率。查看在给定数据点的失效概率。GUI方式:方式:Main MenuProb DesignProb ResultsStatisticsProbabilitiesAPDL方式方式:PDPROB, Rlab,
42、Name, Relation, LIMIT, -, CONF选择抽样数据文件选择抽样数据文件选择要分析的变量选择要分析的变量逻辑关系逻辑关系比较数值比较数值设置置信度水平设置置信度水平432.8 概率设计结果后处理概率设计结果后处理灵敏度灵敏度输出量和设计变量的相关程度输出量和设计变量的相关程度;如何提高零件质量如何提高零件质量;如何在质量保证条件下节约资金。如何在质量保证条件下节约资金。GUI方式:方式:Main MenuProb DesignProb ResultsTrendsSensitivitiesAPDL方式方式:PDSENS, Rlab, Name, Chart, Type, SL
43、EVEL选择抽样数据文件选择抽样数据文件选择要分析的变量选择要分析的变量绘制灵敏度图的类型绘制灵敏度图的类型计算灵敏度的相关性系数类型计算灵敏度的相关性系数类型显著性水平显著性水平442.8 概率设计结果后处理概率设计结果后处理报告设置和生成报告报告设置和生成报告包括上述后处理所有结果包括上述后处理所有结果GUI方式:方式:Main MenuProb DesignProb ResultsReportReport Options Main MenuProb DesignProb ResultsReportGenerate ReportAPDL方式方式:PDROPT, RVAR, CORR, ST
44、AT, SHIS, HIST, CDF, SENS, PDWRITE, File, Fnam, Lnam45 如图所示,两边固定方板承受集中力载荷模型。其尺寸和材料属性均是不确定的输入参数。随机条件如下:方板边长100mm,板厚1mm,服从均匀分布;材料弹性模量,服从高斯分布。标准方差是均值的倍;密度均值8000kg/mm3,集中载荷只能是正值,服从LOG1分布,标准方差为均值的10%;,在上述条件下,板的最大变形和固定边界的最大等效应力的输出为随机行为,具体研究内容如下:检查统计结果,确定PDS是否执行了足够多的仿真循环计算数目;确定最大变形低于指定值的概率;计算随机响应结果相对于随机输入参
45、数的灵敏度值;生成输出参数相对于最重要输入参数的离散图; 承受横向集中力板的承受横向集中力板的LHS抽样抽样MCS概率设计实例概率设计实例46如图所示,三根铰接杆承受集中力载荷模型。其尺寸和材料属性均是不确定的输入参数。随机条如图所示,三根铰接杆承受集中力载荷模型。其尺寸和材料属性均是不确定的输入参数。随机条如图所示,三根铰接杆承受集中力载荷模型。其尺寸和材料属性均是不确定的输入参数。随机条如图所示,三根铰接杆承受集中力载荷模型。其尺寸和材料属性均是不确定的输入参数。随机条件如下:件如下:件如下:件如下:截面积截面积截面积截面积 A1A1均值为均值为均值为均值为10mm210mm2,标准差为,
46、服从高斯分布,标准差为,服从高斯分布,标准差为,服从高斯分布,标准差为,服从高斯分布截面积截面积截面积截面积A2A2最小值为最小值为最小值为最小值为10,10,最可能的值为最可能的值为最可能的值为最可能的值为11,11,最大值为最大值为最大值为最大值为1212,服从三角分布,服从三角分布,服从三角分布,服从三角分布截面积截面积截面积截面积A3A3最小值为最小值为最小值为最小值为9,9,最大值为最大值为最大值为最大值为1111,服从均匀分布,服从均匀分布,服从均匀分布,服从均匀分布定义输入变量定义输入变量定义输入变量定义输入变量A1A1与与与与A3A3之间的关系之间的关系之间的关系之间的关系, ,相关系数为相关系数为相关系数为相关系数为在上述条件下,杆件的最大轴向应力的输出在上述条件下,杆件的最大轴向应力的输出SIG1、SIG2、SIG3为随机行为,具体研究内容如下:为随机行为,具体研究内容如下:观察变量的抽样过程,确定观察变量的抽样过程,确定PDS是否执行了足够是否执行了足够多的仿真循环计算数目;多的仿真循环计算数目;绘制绘制SIG1响应历史曲线;响应历史曲线;绘制绘制SIG2的分布柱状图;的分布柱状图;对对VTOT进行灵敏度分析;进行灵敏度分析; 三根杆桁架系统的直接抽样三根杆桁架系统的直接抽样MCS概率分析实例概率分析实例47