CN106407620B - 一种基于abaqus的工程结构响应面随机有限元分析处理方法 - Google Patents

一种基于abaqus的工程结构响应面随机有限元分析处理方法 Download PDF

Info

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
engineering structure
equation
response surface
finite element
space
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.)
Expired - Fee Related
Application number
CN201611013444.4A
Other languages
English (en)
Other versions
CN106407620A (zh
Inventor
蔡德咏
刘福军
马大为
朱忠领
任杰
李士彦
胡建国
仲健林
陈福红
贾风雨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
72465 Troops Of Ji'nan Military Region Of Chinese People's Liberation Army
Original Assignee
72465 Troops Of Ji'nan Military Region Of Chinese People's Liberation Army
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 72465 Troops Of Ji'nan Military Region Of Chinese People's Liberation Army filed Critical 72465 Troops Of Ji'nan Military Region Of Chinese People's Liberation Army
Priority to CN201611013444.4A priority Critical patent/CN106407620B/zh
Publication of CN106407620A publication Critical patent/CN106407620A/zh
Application granted granted Critical
Publication of CN106407620B publication Critical patent/CN106407620B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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的工程结构响应面随机有限元分析处理方法
技术领域
本发明涉及工程数据处理应用技术领域，更为具体地，涉及一种基于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所述的方法，其特征在于：通过数据采集单元采集工程结构中存在随机分布特性的参数为工程结构的尺寸、材料参数和/或载荷及边界条件。
CN201611013444.4A 2016-11-17 2016-11-17 一种基于abaqus的工程结构响应面随机有限元分析处理方法 Expired - Fee Related CN106407620B (zh)

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 Expired - Fee Related CN106407620B (zh) 2016-11-17 2016-11-17 一种基于abaqus的工程结构响应面随机有限元分析处理方法

Country Status (1)

Country Link
CN (1) CN106407620B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107403037B (zh) * 2017-07-04 2021-05-07 清华大学 一种开源有限元求解及优化分析方法
CN111159948B (zh) * 2019-12-30 2022-08-16 西北工业大学 考虑随机不确定性的关节轴承的可靠性分析方法
CN118094821B (zh) * 2024-04-17 2024-07-02 中国人民解放军海军工程大学 机械系统可靠度优化方法、装置、设备以及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
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 北京航空航天大学 一种高效求解大规模非线性随机结构系统状态的多尺度迭代方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
基于可靠性的复合材料定向管优化设计;蔡德咏 等;《工程力学》;20120630;第29卷(第6期);第259-264页

Also Published As

Publication number Publication date
CN106407620A (zh) 2017-02-15

Similar Documents

Publication Publication Date Title
CN111832102A (zh) 一种高维随机场条件下的新型复合材料结构优化设计方法
CN106407620B (zh) 一种基于abaqus的工程结构响应面随机有限元分析处理方法
CN106570204A (zh) 一种基于cpu+gpu异构并行计算的透平机械叶片静强度特性分析方法
CN106650141A (zh) 一种预测周期性材料性能的不确定性分析方法
CN116305522A (zh) 一种基于数字孪生的飞机结构可靠性仿真试验模型标定方法、设备、介质
CN113722860B (zh) 基于降阶模型的瞬态热力状态在线评估方法、设备及介质
KR101420304B1 (ko) 구조물의 신뢰성 해석 방법
CN112926163B (zh) 基于Workbench和Matlab构建螺栓疲劳寿命预测模型的方法及装置
CN114792020B (zh) 一种基于机器学习的建筑抗震韧性快速评估方法及系统
CN106772645A (zh) 基于一般先验信息约束的核磁共振数据反演方法和装置
Broggi et al. Reliability assessment of axially compressed composite cylindrical shells with random imperfections
Chu et al. Reliability based optimization with metaheuristic algorithms and Latin hypercube sampling based surrogate models
CN113836760B (zh) 一种涡轮盘蠕变疲劳寿命可靠性评估方法
CN114528741A (zh) 一种基于有限元分析法的大跨度钢结构指标变化预测方法
CN114611768A (zh) 一种配电网业扩配套容量时序建设规模预测方法
CN113343559A (zh) 迭代重加权最小二乘法极限学习机响应面可靠性分析方法
CN112800540B (zh) 一种基于高斯过程回归的航空发动机载荷谱任务段建模方法
CN117725859A (zh) 飞机疲劳损伤关键部位的疲劳裂纹扩展预测方法及系统
CN117688316A (zh) 一种用于飞机载荷预测的数据降维处理方法
CN117077513A (zh) 基于depinn求解中子扩散特征值问题噪声分析方法
CN116227045B (zh) 一种结构试件的局部应力应变场构造方法及系统
CN109635452B (zh) 一种高效的多峰随机不确定性分析方法
CN111400948B (zh) 一种应力释放螺栓所受应力的有限元计算方法
CN113155384A (zh) 用于减小结构阻尼比识别的不确定性的传感器布置方法
CN112115540B (zh) 一种融合实测数据的海工支撑结构疲劳分析方法及系统

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190806