CN106407620B - 一种基于abaqus的工程结构响应面随机有限元分析处理方法 - Google Patents
一种基于abaqus的工程结构响应面随机有限元分析处理方法 Download PDFInfo
- Publication number
- CN106407620B CN106407620B CN201611013444.4A CN201611013444A CN106407620B CN 106407620 B CN106407620 B CN 106407620B CN 201611013444 A CN201611013444 A CN 201611013444A CN 106407620 B CN106407620 B CN 106407620B
- Authority
- CN
- China
- Prior art keywords
- parameter
- finite element
- engineering structure
- space
- standard normal
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
一种基于ABAQUS的工程结构响应面随机有限元分析处理方法,包括建立参数化有限元模型,确定随机分布特性,采用序列响应面法拟合极限状态方程,将随机分布参数当量正太化,转化标准正态空间等步骤,能够减小误差,提高工程结构可靠度和灵敏度。
Description
技术领域
本发明涉及工程数据处理应用技术领域,更为具体地,涉及一种基于ABAQUS的工程结构响应面随机有限元分析处理方法。
背景技术
计算机辅助工程(Computer Aided Engineering,简称为CAE)分析是以有限元法为基础发展起来的一种十分有效的计算机数值仿真与优化设计技术。有限元法是一种高效能、常用的数值计算方法,其基于确定性的参数进行计算,其中确定性的参数是指结构的材料性能、载荷及边界条件等都是确定的。但在工程实际中,结构的材料性能、载荷及边界条件等都存在着很大的统计随机性,而这些结构参数的统计随机性,会对结构的临界性能和可靠性有较大的影响,给分析结果带来误差。
发明内容
本发明的目的在于克服现有技术的不足,提供一种能够减小误差,提高工程结构可靠度和灵敏度的基于ABAQUS的工程结构响应面随机有限元分析处理方法。
本发明提供了一种基于ABAQUS的工程结构响应面随机有限元分析处理方法,具体包括如下步骤:
步骤S1:建立工程结构参数化有限元模型;
步骤S2:根据工程结构参数化有限元模型,通过数据采集单元采集工程结构中存在随机分布特性的参数,形成参数集合X=(x1,x2,...,xn),并根据工程标准、试验数据和实际经验,确定所述参数集合中每个参数的随机分布特性;
步骤S3:采用序列响应面法拟合极限状态方程g(X):
其中,a、bi、ci分别为极限状态方程的二次项系数,xi为原始空间的随机分布参数;
步骤S4:取参数集合在原始空间中的均值点为验算点的初值并通过Rackwits-Fiessler算法将随机分布参数xi当量正态化;其中,在验算点处,当量正态随机变量的累积分布函数与原随机变量的累积分布函数相等,以及当量正态随机变量的概率密度函数与原随机变量的概率密度函数相等;
式(1)中,为标准正态空间的随机分布参数;
式(2)中,为标准正态概率密度函数;
式(3)中,和分别为随机分布参数xi对应的近似正态分布函数的均值和方差;
根据式(1)~(3),求出近似正态分布函数的均值和方差
步骤S5:将极限状态响应面方程函数转化到标准正态空间变为:
其中,g'(U)为转化到标准正态空间的极限状态方程函数,α0、ri、λi分别为极限状态方程函数的二次项系数,ui为标准正态空间的随机分布参数;
步骤S6:取标准正态空间中拟合点的参数δ=0.1~0.5,形成标准正态空间中的2n+1个拟合点,2n+1个拟合点分别为:
步骤S7:将标准正态空间中的2n+1个拟合点换算成原始空间中的拟合点,原始空间的拟合点为:
步骤S8:将原始空间的拟合点代入参数化有限元模型获得极限状态响应值g1、g2、…、g2n、g2n+1;
步骤S9:将响应值代入标准正态空间响应面拟合方程获得联立方程式,求出联立方程式的系数a0、r1、λ1、…、rn、λn;其中,
步骤S10:将联立方程式的系数a0、r1、λ1、…、rn、λn分别带入标准正态空间响应面拟合方程,获得响应面方程式:
步骤S11:对响应面方程式采用JC算法求得工程结构的可靠度指标β、结构可靠度对随机分布参数的灵敏度系数αi和新验算点的初值其中,
JC算法的计算过程是先将极限状态方程展开为泰勒级数,在分别计算,由于JC算法是现有技术,故计算过程的细节在此不过多说明;
步骤S12:根据可靠度指标β和新验算点的初值重复步骤S7~S11直至可靠度指标β收敛,其中可靠度指标β收敛是指上次求出的β至上次求出的β≤允许误差ε,而允许误差ε根据工程结构的要求确定。
本发明的基于ABAQUS的工程结构响应面随机有限元分析处理方法,将可靠性计算方法与有限元法相结合建立的,能够求解工程复杂结构可靠度的问题,该方法不受随机参数变异性大小的限制,计算较为准确,具有高效、方便的特点。
附图说明
图1为基于ABAQUS的工程结构响应面随机有限元分析处理方法的流程示意图;
图2为悬臂梁的力学模型示意图;
图3为轴承座有限元网格模型;
具体实施方式
下面详细说明本发明的具体实施,有必要在此指出的是,以下实施只是用于本发明的进一步说明,不能理解为对本发明保护范围的限制,这些方面指示的仅仅是可使用本发明的原理的各种方式中的一些方式,本发明旨在包括所有这些方面以及它们的等同物,并且该领域技术熟练人员根据上述本发明内容对本发明做出的一些非本质的改进和调整,仍然属于本发明的保护范围。
本发明提供了一种基于ABAQUS的工程结构响应面随机有限元分析处理方法,首先,ABAQUS(有限元分析)是工程模拟的有限元软件,其解决问题的范围从相对简单的线性分析到复杂的非线性问题。本发明通过结合ABAQUS的二次开发技术,基于Python语言编程,实现工程结构的参数化建模,通过调用有限元分析结果建立工程结构的极限状态响应面方程,结构可靠度计算方法(优选为JC算法)来计算结构的可靠度和分析参数的敏感度,从而建立一种能有效地解决实际工程复杂结构可靠性分析问题的响应面随机有限元分析处理方法。
图1示出了根据本发明的基于ABAQUS的工程结构响应面随机有限元分析处理方法的流程示意图,具体包括如下步骤:
步骤S1:建立工程结构参数化有限元模型;
由于ABAQUS软件程序的主体框架均由Python语言构成,利用Python才可以对ABAQUS软件进行二次开发,所以需要基于Python语言建立参数化有限元模式;
步骤S2:根据工程结构参数化有限元模型,通过数据采集单元采集工程结构中存在随机分布特性的参数,形成参数集合X=(x1,x2,...,xn),并根据工程标准、试验数据和实际经验,确定所述参数集合中每个参数的随机分布特性;
采集工程结构中存在随机分布特性的参数,为采集工程结构的尺寸、材料参数、载荷及边界条件等存在着随机分布特性的参数,参数的随机分布特性是根据工程实际经验或试验测试统计获得的。确定随机分布特性是为了掌握参数的随机分布情况,在此基础上才能进行随机有限元分析,计算可靠性。
步骤S3:采用序列响应面法拟合极限状态方程g(X):
其中,a、bi、ci分别为极限状态方程的二次项系数,xi为原始空间的随机分布参数;
步骤S4:取参数集合在原始空间中的均值点为验算点的初值并通过Rackwits-Fiessler算法将随机分布参数xi当量正态化;其中,在验算点处,当量正态随机变量的累积分布函数与原随机变量的累积分布函数相等,以及当量正态随机变量的概率密度函数与原随机变量的概率密度函数相等;
式(1)中,为标准正态空间的随机分布参数;
式(2)中,为标准正态概率密度函数;
式(3)中,和分别为随机分布参数xi对应的近似正态分布函数的均值和方差;
根据式(1)~(3),求出近似正态分布函数的均值和方差
步骤S5:将极限状态响应面方程函数转化到标准正态空间变为:
其中,g'(U)为转化到标准正态空间的极限状态方程函数,α0、ri、λi分别为极限状态方程函数的二次项系数,ui为标准正态空间的随机分布参数;
步骤S6:取标准正态空间中拟合点的参数δ=0.1~0.5,形成标准正态空间中的2n+1个拟合点:
步骤S7:将标准正态空间中的2n+1个拟合点换算成原始空间中的拟合点,原始空间的拟合点为:
步骤S8:将原始空间的拟合点代入参数化有限元模型获得极限状态响应值g1、g2、…、g2n、g2n+1;
步骤S9:将响应值代入标准正态空间响应面拟合方程获得联立方程式,求出联立方程式的系数a0、r1、λ1、…、rn、λn;其中,
步骤S10:将联立方程式的系数a0、r1、λ1、…、rn、λn分别带入标准正态空间响应面拟合方程,获得响应面方程式:
步骤S11:对响应面方程式采用JC算法求得工程结构的可靠度指标β、结构可靠度对随机分布参数的灵敏度系数αi和新验算点的初值其中,
步骤S12:根据可靠度指标β和新验算点的初值重复步骤S7~S11直至可靠度指标β收敛,其中可靠度指标β收敛是指上次求出的β至上次求出的β≤允许误差ε,而允许误差ε根据工程结构的要求确定。
上述内容详细说明了本发明提供的基于ABAQUS的工程结构响应面随机有限元分析处理方法的流程,下面将以两个实例对上述方法可取得的有益效果进行说明。
实例一
以某悬臂梁为例,如图2所示,悬臂梁长为L,梁的截面形状为矩形,尺寸为a×b。悬臂梁上表面受均布载荷,载荷大小为F,材料屈服应力为σs。分析中考虑结构尺寸a、b、L及载荷F、屈服强度σs五个变量的随机性,具体分布如表1所示:
表1随机变量的分布
根据材料力学知识,该悬臂梁的极限状态方程为:
采用序列响应面拟合标准正态空间的极限状态方程为:
g=a0+r1A+λ1A2+r2B+λ2B2+r3l+λ3l2+r4f+λ4f2+r5s+r5s2;
其中,A、B、l、f、s均为标准正态随机向量。
将最后求得的β值与Monte-Carlo抽样法或其它解析方法计算的β值进行比较,以Monte-Carlo抽样法的β值为精确解,计算允许误差,比较结果的分析如表2所示:
表2可靠度指标β的结果分析
通过表2可以看出,由表中的数据可知,通过本发明提供的方法计算得到的β与精确解相比,允许误差只有1.997%,虽然比二次二阶矩法的精度低,但比一次二阶矩法的均值法要精确,而1.997%的相对误差在大多数工程实际中是可以接受的。同时,该方法只需根据变异性设定随机变量参数即可,不受随机变量变异性大小的限制,适用性广。
实例二
某轴承座通过4个安装孔进行固定,轴承孔的下半部分承受由轴传来的径向压力载荷P,轴承孔圆周上承受推力载荷P/5。轴承座材料为钢,弹性模量为E,泊松比μ=0.3,屈服强度为σs。分析中考虑径向压力载荷P1、推力载荷P2、弹性模量E、屈服强度σs四个变量的随机性。
根据轴承座的几何模型建立轴承座的有限元模型,建立的轴承座的有限元模型如图3所示。
将最后求得的β值与Monte-Carlo抽样法的β值进行比较,以MCSFEM法的β值为精确解,计算允许误差,比较结果的分析如表3所示:
表3可靠度指标β的结果分析
由表中的数据可知,本发明提供的方法与MCSFEM法(10000次抽样实验)计算得到的精确解相比误差只有1.93%,而计算时间只有MCSFEM法的1.5%,证明本发明通过的方法计算较为准确,计算效率很高,可用于工程实际。
尽管为了说明的目的,已描述了本发明的示例性实施方式,但是本领域的技术人员将理解,不脱离所附权利要求中公开的发明的范围和精神的情况下,可以在形式和细节上进行各种修改、添加和替换等的改变,而所有这些改变都应属于本发明所附权利要求的保护范围,并且本发明要求保护的产品各个部门和方法中的各个步骤,可以以任意组合的形式组合在一起。因此,对本发明中所公开的实施方式的描述并非为了限制本发明的范围,而是用于描述本发明。相应地,本发明的范围不受以上实施方式的限制,而是由权利要求或其等同物进行限定。
Claims (3)
1.一种基于ABAQUS的工程结构响应面随机有限元分析处理方法,其特征在于,具体包括如下步骤:
步骤S1:建立工程结构参数化有限元模型;
步骤S2:根据工程结构参数化有限元模型,通过数据采集单元采集工程结构中存在随机分布特性的参数,形成参数集合X=(x1,x2,...,xn),并根据工程标准、试验数据和实际经验,确定所述参数集合中每个参数的随机分布特性;
步骤S3:采用序列响应面法拟合极限状态方程g(X):
其中,a、bi、ci分别为极限状态方程的二次项系数,xi为原始空间的随机分布参数;
步骤S4:取参数集合在原始空间中的均值点为验算点的初值并通过Rackwits-Fiessler算法将随机分布参数xi当量正态化;其中,在验算点处,当量正态随机变量的累积分布函数与原随机变量的累积分布函数相等,以及当量正态随机变量的概率密度函数与原随机变量的概率密度函数相等;
式(1)中,为标准正态空间的随机分布参数;
式(2)中,为标准正态概率密度函数;
式(3)中,和分别为随机分布参数xi对应的近似正态分布函数的均值和方差;
根据式(1)~(3),求出近似正态分布函数的均值和方差
步骤S5:将极限状态响应面方程函数转化到标准正态空间变为:
其中,g'(U)为转化到标准正态空间的极限状态方程函数,α0、ri、λi分别为极限状态方程函数的二次项系数,ui为标准正态空间的随机分布参数;
步骤S6:取标准正态空间中拟合点的参数δ=0.1~0.5,形成标准正态空间中的2n+1个拟合点:
步骤S7:将标准正态空间中的2n+1个拟合点换算成原始空间中的拟合点,原始空间的拟合点为:
步骤S8:将原始空间的拟合点代入参数化有限元模型获得极限状态响应值g1、g2、…、g2n、g2n+1;
步骤S9:将响应值代入标准正态空间响应面拟合方程获得联立方程式,求出联立方程式的系数a0、r1、λ1、…、rn、λn;其中,
步骤S10:将联立方程式的系数a0、r1、λ1、…、rn、λn分别带入标准正态空间响应面拟合方程,获得响应面方程式:
步骤S11:对响应面方程式采用JC算法求得工程结构的可靠度指标β、结构可靠度对随机分布参数的灵敏度系数αi和新验算点的初值其中,
步骤S12:根据可靠度指标β和新验算点的初值重复步骤S7~S11直至可靠度指标β收敛,其中可靠度指标β收敛是指上次求出的β至上次求出的β≤允许误差ε,而允许误差ε根据工程结构的要求确定。
2.如权利要求1所述的方法,其特征在于:步骤S2中根据工程结构参数化有限元模型,通过数据采集单元采集工程结构中存在随机分布特性的参数是通过工程实际经验或试验测试统计获得。
3.如权利要求1或2所述的方法,其特征在于:通过数据采集单元采集工程结构中存在随机分布特性的参数为工程结构的尺寸、材料参数和/或载荷及边界条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611013444.4A CN106407620B (zh) | 2016-11-17 | 2016-11-17 | 一种基于abaqus的工程结构响应面随机有限元分析处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611013444.4A CN106407620B (zh) | 2016-11-17 | 2016-11-17 | 一种基于abaqus的工程结构响应面随机有限元分析处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106407620A CN106407620A (zh) | 2017-02-15 |
CN106407620B true CN106407620B (zh) | 2019-08-06 |
Family
ID=58068884
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611013444.4A Active CN106407620B (zh) | 2016-11-17 | 2016-11-17 | 一种基于abaqus的工程结构响应面随机有限元分析处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106407620B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107403037B (zh) * | 2017-07-04 | 2021-05-07 | 清华大学 | 一种开源有限元求解及优化分析方法 |
CN111159948B (zh) * | 2019-12-30 | 2022-08-16 | 西北工业大学 | 考虑随机不确定性的关节轴承的可靠性分析方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6301970B1 (en) * | 1998-08-06 | 2001-10-16 | The United States Of America The Secretary Of The Navy | Cumulative damage model for structural analysis of filed polymeric materials |
CN102332046A (zh) * | 2011-09-30 | 2012-01-25 | 北京工业大学 | 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 |
CN102663176A (zh) * | 2012-03-28 | 2012-09-12 | 北京航空航天大学 | 针对高可靠机械产品的主动可靠性分析评价方法 |
CN103136428A (zh) * | 2013-03-12 | 2013-06-05 | 上海交通大学 | 基于两种不确定的轿车车身结构稳健设计方法 |
CN103902785A (zh) * | 2014-04-14 | 2014-07-02 | 北京航空航天大学 | 一种基于多元不确定性的结构有限元模型修正方法 |
CN105787149A (zh) * | 2015-12-27 | 2016-07-20 | 北京航空航天大学 | 一种由弧齿锥齿轮传动系统轴上功率谱向齿面应力谱精确转换的方法 |
CN106021186A (zh) * | 2016-05-13 | 2016-10-12 | 北京航空航天大学 | 一种高效求解大规模非线性随机结构系统状态的多尺度迭代方法 |
-
2016
- 2016-11-17 CN CN201611013444.4A patent/CN106407620B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6301970B1 (en) * | 1998-08-06 | 2001-10-16 | The United States Of America The Secretary Of The Navy | Cumulative damage model for structural analysis of filed polymeric materials |
CN102332046A (zh) * | 2011-09-30 | 2012-01-25 | 北京工业大学 | 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 |
CN102663176A (zh) * | 2012-03-28 | 2012-09-12 | 北京航空航天大学 | 针对高可靠机械产品的主动可靠性分析评价方法 |
CN103136428A (zh) * | 2013-03-12 | 2013-06-05 | 上海交通大学 | 基于两种不确定的轿车车身结构稳健设计方法 |
CN103902785A (zh) * | 2014-04-14 | 2014-07-02 | 北京航空航天大学 | 一种基于多元不确定性的结构有限元模型修正方法 |
CN105787149A (zh) * | 2015-12-27 | 2016-07-20 | 北京航空航天大学 | 一种由弧齿锥齿轮传动系统轴上功率谱向齿面应力谱精确转换的方法 |
CN106021186A (zh) * | 2016-05-13 | 2016-10-12 | 北京航空航天大学 | 一种高效求解大规模非线性随机结构系统状态的多尺度迭代方法 |
Non-Patent Citations (1)
Title |
---|
基于可靠性的复合材料定向管优化设计;蔡德咏 等;《工程力学》;20120630;第29卷(第6期);第259-264页 |
Also Published As
Publication number | Publication date |
---|---|
CN106407620A (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | A computational method for automated detection of engineering structures with cyclic symmetries | |
CN108920786B (zh) | 一种基于切比雪夫多项式拟合的区间不确定性分析方法 | |
Wang et al. | A new special element for stress concentration analysis of a plate with elliptical holes | |
CN106597574A (zh) | 一种基于时变云模型的天气温度预测方法及装置 | |
US20190121838A1 (en) | A dynamically non-gaussian anomaly identification method for structural monitoring data | |
CN106407620B (zh) | 一种基于abaqus的工程结构响应面随机有限元分析处理方法 | |
Chen et al. | Improved response surface method for anti-slide reliability analysis of gravity dam based on weighted regression | |
Interian et al. | Tumor growth modelling by cellular automata | |
Wang et al. | Voronoi polygonal hybrid finite elements with boundary integrals for plane isotropic elastic problems | |
Kalita et al. | Unsteady separation leading to secondary and tertiary vortex dynamics: the sub--phenomena | |
CN104123471A (zh) | 重型施工机械重载结构安全评估方法 | |
Rocco et al. | Global sensitivity analysis in a multi-state physics model of component degradation based on a hybrid state-space enrichment and polynomial chaos expansion approach | |
Liu et al. | A smoothed finite element method (S-FEM) for large-deformation elastoplastic analysis | |
CN107727350B (zh) | 微纳卫星矢量振动试验方法 | |
CN113343423B (zh) | 基于强度空间变异性的随机裂隙网络生成方法 | |
CN112149311B (zh) | 一种基于数量规约的非线性多元统计回归测井曲线预测方法 | |
CN110334459B (zh) | 一种输电塔塔线体系快速精细化建模系统及方法 | |
Tvergaard | Discrete modelling of ductile crack growth by void growth to coalescence | |
CN106934729A (zh) | 建筑检测鉴定方法及装置 | |
CN106021186A (zh) | 一种高效求解大规模非线性随机结构系统状态的多尺度迭代方法 | |
Hernández-Rivera et al. | Anisotropic grain growth modeling under the spparks framework | |
Zhao et al. | Uncertainty static analysis of structures with hybrid spatial random and interval properties | |
CN110276087A (zh) | 一种非平稳随机扭转振动疲劳寿命预测方法 | |
Hrnjica et al. | Numerical calculation of J-integral using finite elements method | |
Nemade et al. | The mesh quality significance in finite element analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |