CN112395791A - 机械结构静强度全工况评定方法、系统、设备及存储介质 - Google Patents
机械结构静强度全工况评定方法、系统、设备及存储介质 Download PDFInfo
- Publication number
- CN112395791A CN112395791A CN202011253385.4A CN202011253385A CN112395791A CN 112395791 A CN112395791 A CN 112395791A CN 202011253385 A CN202011253385 A CN 202011253385A CN 112395791 A CN112395791 A CN 112395791A
- Authority
- CN
- China
- Prior art keywords
- condition
- working
- full
- stress
- maximum stress
- 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.)
- Granted
Links
- 230000003068 static effect Effects 0.000 title claims abstract description 58
- 238000011156 evaluation Methods 0.000 title claims abstract description 48
- 239000000463 material Substances 0.000 claims abstract description 66
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000000007 visual effect Effects 0.000 claims abstract description 11
- 238000012800 visualization Methods 0.000 claims abstract description 9
- 238000003466 welding Methods 0.000 claims description 14
- 230000006872 improvement Effects 0.000 claims description 5
- 238000005457 optimization Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 12
- 230000003137 locomotive effect Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 6
- RUDATBOHQWOJDD-UZVSRGJWSA-N ursodeoxycholic acid Chemical compound C([C@H]1C[C@@H]2O)[C@H](O)CC[C@]1(C)[C@@H]1[C@@H]2[C@@H]2CC[C@H]([C@@H](CCC(O)=O)C)[C@@]2(C)CC1 RUDATBOHQWOJDD-UZVSRGJWSA-N 0.000 description 5
- 238000012545 processing Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 231100000817 safety factor Toxicity 0.000 description 2
- 229910000838 Al alloy Inorganic materials 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
Images
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种机械结构静强度全工况评定方法、系统、设备及存储介质，整个有限元模型的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况进行可视化显示，材料利用度可视化图能全面显示机械结构的静强度承载性能，可查看任何位置的材料利用度；全工况最大应力可视化图能全面显示机械结构的静强度全工况最大应力，可查看任何位置的全工况最大应力；载荷工况可视化图能全面显示机械结构的工况特性，可查看任何位置全工况最大应力对应的载荷工况；该方法能够从结构静强度所有载荷工况即全工况数据进行评定，能全面地评估结构所有静强度载荷工况的计算结果，形成直观的云图和详细的数据。
Description
技术领域
本发明属于轨道交通技术领域，尤其涉及一种机械结构静强度全工况评定方法、系统、设备及存储介质。
背景技术
安全可靠是机械结构功能的基础能力，机械结构设计需要满足产品的运行要求，静强度计算是机械承载结构设计中重要的环节。例如通过机车车体静强度计算和分析可以有效地评定车体的安全性和材料的力学承载情况，可以方便地找出车体的最危险位置及其安全系数。在机械结构设计时会根据规范或标准规定产品计算分析的多个静强度载荷工况，机械结构应用在不同产品上对应有不同静强度载荷工况。静强度均以所使用材料的屈服强度或屈服应力作为评价指标，当结构任意点的材料利用度(计算应力/屈服应力)>1时，结构的静强度就不能满足要求。
传统的静强度计算分析通常以各个工况的计算应力结果进行评价，计算应力以云图方式和颜色显示，颜色显示越红表示应力越大，这种显示方式只能针对高应力区域和关键部位，再将这些部位的计算应力转换为安全系数或材料利用度。当结构中有低屈服强度或屈服应力材料时，容易造成评估遗漏，如采用焊缝连接的大型承载结构，焊缝的屈服强度会比其母材低，特别是铝合金材料，焊缝的屈服强度最多为母材的50％，而连接焊缝在结构大范围的存在，这就容易造成焊接连接部分的评估遗漏。
传统的静强度计算分析，各个载荷工况间的应力云图、安全系数或材料利用度结果是分开的，不能在同一个云图或表格中显示，结构的静强度计算报告要对每个工况进行评定。因而，要分析整个结构的承载能力，不仅步骤多、图表多，而且还分析不全面。需要查看多个载荷工况，每个载荷工况需要顶面和底面的应力结果，结果查看工况多，一般只会查看每工况的最大应力和关注部位的结果，应力不高的焊缝部位常被遗漏。
传统的静强度计算分析，其工况之间的应力数据无联系，仅凭分析师的记忆，无法了解结构的静强度承载能力，只有大概的判断，无总体的强度数据。无法了解工况对结构的影响能力，绝大结构部位对工况没有直观的分析，工况越多，越不能把握结构的强度特性。
发明内容
本发明的目的在于提供一种机械结构静强度全工况评定方法、系统、设备及存储介质，以解决多个载荷工况的静强度计算结果数据处理量大，易遗漏部分需要评估的区域，无法评估工况对大部分结构区域的影响能力，以及静强度数据分析不全面等问题。
本发明独立权利要求的技术方案解决了上述发明目的中的一个或多个。
本发明是通过如下的技术方案来解决上述技术问题的：一种机械结构静强度全工况评定方法，包括以下步骤：
建立待评定机械结构的有限元模型；
将所述有限元模型的评定对象分成多个对象集合，每个所述对象集合内每个对象的屈服应力相同，并输出多个对象集合文件，每个所述对象集合文件的内容为对应对象集合的对象编号；
将每个所述对象集合文件分解成对应对象集合内每个对象的屈服应力，形成对象编号-屈服应力对应的评定指标文件；
按照对象编号顺序以文件的形式输出每个所述对象在所有载荷工况下的TOP面应力和BOTTOM面应力，形成所有载荷工况下对象的应力数据文件；
根据每个对象在所有载荷工况下的TOP面应力和BOTTOM面应力得到每个对象的全工况最大应力，并记录每个对象在全工况最大应力时对应的载荷工况，形成对象编号-全工况最大应力-全工况最大应力时的载荷工况文件；
检索所有对象的屈服应力和全工况最大应力，计算每个对象的材料利用度；
以所述有限元模型的单元和节点组成数据，形成结果输入文件，所述结果输入文件包含各对象的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况；
读取所述结果输入文件，对整个有限元模型的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况进行可视化显示，并得到模型的可视化结果文件。
本发明中，材料利用度可视化图能全面显示机械结构的静强度承载性能，可查看任何位置的材料利用度；全工况最大应力可视化图能全面显示机械结构的静强度全工况最大应力，可查看任何位置的全工况最大应力；载荷工况可视化图能全面显示机械结构的工况特性，可查看任何位置全工况最大应力对应的载荷工况；该方法能够从结构静强度所有载荷工况即全工况数据进行评定，能全面地评估结构所有静强度载荷工况的计算结果，形成直观的云图和详细的数据，快速、简洁、高效地反映了结构的整体静强度性能，避免了数据处理量大，容易遗漏部分区域的缺点。
优选地，所述评定对象为有限元模型的所有节点或所有单元。利用有限元分析软件构建的有限元模型是由多个节点和多个单元构成的，在进行全工况评定时，可以以模型的所有节点为评定对象，也可以以模型的所有单元为评定对象。
进一步地，每个所述对象的全工况最大应力为：
σi,max＝max(σi,j,T,σi,j,B),j＝1,2,3,…,n
其中，σi,max为第i个对象的全工况最大应力，σi,j,T为第i个对象在第j种载荷工况下的TOP面应力，σi,j,B为第i个对象在第j种载荷工况下的BOTTOM面应力，n为载荷工况的数量。
进一步地，每个所述对象的材料利用度的计算公式为：
其中，μi为第i个对象的材料利用度，σi,max为第i个对象的全工况最大应力，σi,s为第i个对象的屈服应力。
进一步地，所述评定方法还包括：对所述材料利用度大于第一临界值的对象进行全方位筛选识别，判断材料利用度大于第一临界值的对象是否满足安全要求，如果不满足安全要求，则对该对象进行优化改进。
每个对象对应机械结构的一个局部位置，该判断为机械结构局部位置的优化改进提供了判断依据，提高了机械结构的整体安全性。
优选地，所述第一临界值为1。第一临界值为1表示对象的应力为其屈服强度，达到了塑性变形的临界值。
进一步地，所述评定方法还包括：对所述材料利用度小于第二临界值的位置进行轻量化处理，以降低成本和运营能耗。
优选地，所述第二临界值为0.5。第二临界值为0.5表示对象的利用水平低，可以局部变薄以减少重量，达到结构轻量化目的。
进一步地，所述评定方法还包括：检索某个对象在所有载荷工况下的应力数据，输出该对象的载荷工况-应力曲线，可全面、直观、高效地进行机械结构在所有载荷工况下的整体性能评定，能够分析各载荷工况对某个对象所对应位置的影响能力。
进一步地，所述评定方法还包括：输出焊缝工艺文件编制所需的焊缝最大应力和材料利用度数据。
本发明还提供一种机械结构静强度全工况评定系统，包括：
模型建立单元，用于建立待评定机械结构的有限元模型；
集合及文件输出单元，用于将所述有限元模型的评定对象分成多个对象集合，每个所述对象集合内每个对象的屈服应力相同，并输出多个对象集合文件，每个所述对象集合文件的内容为对应对象集合的对象编号；
指标文件生成单元，用于将每个所述对象集合文件分解成对应对象集合内每个对象的屈服应力，形成对象编号-屈服应力对应的评定指标文件；
应力数据文件生成单元，用于按照对象编号顺序以文件的形式输出每个所述对象在所有载荷工况下的TOP面应力和BOTTOM面应力，形成所有载荷工况下对象的应力数据文件；
全工况最大应力文件生成单元，用于根据每个对象在所有载荷工况下的TOP面应力和BOTTOM面应力得到每个对象的全工况最大应力，并记录每个对象在全工况最大应力时对应的载荷工况，形成对象编号-全工况最大应力-全工况最大应力时的载荷工况文件；
利用度计算单元，用于检索所有对象的屈服应力和全工况最大应力，计算每个对象的材料利用度；
结果输入文件生成单元，用于以所述有限元模型的单元和节点组成数据，形成结果输入文件，所述结果输入文件包含各对象的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况；
可视化单元，用于读取所述结果输入文件，对整个有限元模型的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况进行可视化显示，并得到模型的可视化结果文件。
本发明还提供一种设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上所述机械结构静强度全工况评定方法。
本发明还提供一种存储介质，其上存储有计算机程序，该程序被处理器执行时实现如上所述机械结构静强度全工况评定方法。
有益效果
与现有技术相比，本发明所提供的一种机械结构静强度全工况评定方法、系统、设备及存储介质，材料利用度可视化图能全面显示机械结构的静强度承载性能，可查看任何位置的材料利用度；全工况最大应力可视化图能全面显示机械结构的静强度全工况最大应力，可查看任何位置的全工况最大应力；载荷工况可视化图能全面显示机械结构的工况特性，可查看任何位置全工况最大应力对应的载荷工况；该方法能够从结构静强度所有载荷工况即全工况数据进行评定，能全面地评估结构所有静强度载荷工况的计算结果，形成直观的云图和详细的数据，快速、简洁、高效地反映了结构的整体静强度性能，避免了数据处理量大，容易遗漏部分区域的缺点。
附图说明
为了更清楚地说明本发明的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一个实施例，对于本领域普通技术人员来说，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例1中一种机械结构静强度全工况评定方法流程图；
图2是本发明实施例2中车体有限元模型的材料利用度云图；
图3是本发明实施例2中某个节点的载荷工况-应力曲线图。
具体实施方式
下面结合本发明实施例中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
如图1所示，本实施例所提供的一种机械结构静强度全工况评定方法，包括以下步骤：
1、建立待评定机械结构的有限元模型
结构有限元模型的建立为现有技术，可采用ANSYS软件进行结构有限元模型的建立。
2、将有限元模型的评定对象分成多个对象集合，每个对象集合内每个对象的屈服应力相同，并输出多个对象集合文件，每个对象集合文件的内容为对应对象集合的对象编号。
有限元模型是由多个节点和多个单元构成的，在进行全工况评定时，可以以模型的所有节点为评定对象，也可以以模型的所有单元为评定对象，即将有限元模型的所有节点分成多个节点集合，每个节点集合内每个节点的屈服应力相同，并输出多个节点集合文件，每个节点集合文件的内容为对应节点集合的节点编号。或者，将有限元模型的所有单元分成多个单元集合，每个单元集合内每个单元的屈服应力相同，并输出多个单元集合文件，每个单元集合文件的内容为对应单元集合的单元编号。
将有限元模型的所有N个节点分成M个节点集合Jk。
其中，Jk为第k个节点集合或单元集合，k＝1,2,3,…,M，M为节点集合或单元集合的数量，N为有限元模型的节点数量或单元数量，Ni为第i个节点或单元的编号。
每个节点集合内每个节点的屈服应力相同，或每个单元集合内每个单元的屈服应力相同，即屈服应力相等的节点或单元构成了一个节点集合或单元集合。
对象集合与对象集合文件是一一对应的，例如，某个对象集合Jk是由编号为371101、371102、371103、371104、371105、371106的节点构成的，这6个节点的屈服应力均为460MPa，则与该对象集合Jk对应的对象集合文件的内容为371101、371102、371103、371104、371105、371106，即为这6个节点编号。
3、将每个对象集合文件分解成对应对象集合内每个对象的屈服应力σi,s，形成对象编号-屈服应力对应的评定指标文件。
σi,s＝YieldStress(Jk)
其中，YieldStress(Jk)表示取对象集合Jk中每个对象的屈服应力，由于同一对象集合中每个对象的屈服应力相同，因此，相当于取该对象集合的屈服应力。评定指标文件可以按照对象编号的顺序显示对象编号-屈服应力。
4、按照对象编号顺序以文件的形式输出每个对象在所有载荷工况下的TOP面应力σi,j,T和BOTTOM面应力σi,j,B，形成所有载荷工况下对象的应力数据文件。
σi,j,T为第i个对象在第j种载荷工况下的TOP面应力，σi,j,B为第i个对象在第j种载荷工况下的BOTTOM面应力。应力数据文件的内容为对象编号-第j种载荷工况-TOP面应力-BOTTOM面应力，应力数据文件包含所有对象在所有载荷工况下的TOP面应力和BOTTOM面应力。载荷工况是根据机械结构具体应用的产品来决定的，全工况即指某个机械结构应用于某种产品时的所有载荷工况。
5、根据每个对象在所有载荷工况下的TOP面应力σi,j,T和BOTTOM面应力σi,j,B得到每个对象的全工况最大应力σi,max，并记录每个对象在全工况最大应力σi,max时对应的载荷工况LCi,max，形成对象编号-全工况最大应力-全工况最大应力时的载荷工况文件。
σi,max＝max(σi,j,T,σi,j,B|j＝1,2,3,…,n)
LCi,max＝LC(σi,max)
其中，σi,max为第i个对象的全工况最大应力，σi,max即为第i个对象在所有载荷工况下的最大应力，σi,j,T为第i个对象在第j种载荷工况下的TOP面应力，σi,j,B为第i个对象在第j种载荷工况下的BOTTOM面应力，n为载荷工况的数量，LCi,max为第i个对象在全工况最大应力时所对应的载荷工况，LC表示取载荷工况。
6、检索所有对象的屈服应力和全工况最大应力，计算每个对象的材料利用度。
每个对象的材料利用度的计算公式为：
其中，μi为第i个对象的材料利用度，σi,max为第i个对象的全工况最大应力，σi,s为第i个对象的屈服应力。
7、以有限元模型的单元和节点组成数据，形成有限元软件的结果输入文件，结果输入文件包含各对象的全工况最大应力σi,max、材料利用度μi以及全工况最大应力时的载荷工况LCi,max。
8、通过有限元软件读取结果输入文件，对整个有限元模型的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况进行可视化显示，得到模型的可视化结果文件，展示机械结构的应力水平、静强度材料利用度等强度整体性能。
可视化结果文件对全工况结果进行图形显示，显示内容包括材料利用度、全工况最大应力、全工况最大应力时的载荷工况等信息，根据图形显示判断是否满足要求。
9、对材料利用度大于第一临界值的对象进行全方位筛选识别，判断材料利用度大于第一临界值的对象是否满足安全要求，如果不满足安全要求，则对该对象进行优化改进。
每个对象对应机械结构的一个局部位置，该判断为机械结构局部位置的优化改进提供了判断依据，提高了机械结构的整体安全性。
本实施例中，第一临界值为1。第一临界值为1表示对象的应力为其屈服强度，达到了塑性变形的临界值。
10、对材料利用度小于第二临界值的位置进行轻量化处理，以降低成本和运营能耗。
本实施例中，第二临界值为0.5。第二临界值为0.5表示对象的利用水平低，可以局部变薄以减少重量，达到结构轻量化目的。第一临界值和第二临界值均是根据产品需求来设置的。
11、检索某个对象在所有载荷工况下的应力数据，输出该对象的载荷工况-应力曲线，可全面、直观、高效地进行机械结构在所有载荷工况下的整体性能评定，能够分析各载荷工况对某个对象所对应位置的影响能力。
12、输出焊缝工艺文件编制所需的焊缝最大应力和材料利用度数据。
实施例2
以机车车体静强度全工况评定为例，说明实施例1方法的具体实施过程和在实际工程项目中的应用。
机车车体的静强度全工况评定在ANSYS通用有限元软件中进行，车体结构使用了Q460ME、Q355ME，16MnDR等多种材料，车体结构所用材料的最高屈服应力为Q460ME材料的屈服强度(或屈服应力)，为460MPa，最低屈服应力为16MnDR材料的屈服强度(或屈服应力)，为295MPa。车体各部分结构采用焊接连接，按照标准，焊缝的屈服强度为母材的屈服强度除以1.1。该车体结构的静强度载荷工况包括纵向压缩工况、纵向拉伸工况、起吊救援及架车工况、3向冲击和端部排障器等19种载荷工况。
为了说明比较简单的说明实施例1方法的实施方式，本实施例只列出Q460ME、Q355ME，16MnDR三种材料的静强度全工况评定。
该车体结构静强度全工况评定的具体实施过程为：
1、在ANSYS软件中建立车体的有限元模型。
2、在ANSYS软件中对该有限元模型的所有节点，定义了3个节点集合，3个节点集合分别为Q460，Q355，16MnDR，用NLIST命令对每个节点集合输出包含节点编号和节点坐标(X,Y,Z)的节点集合文件，3个节点集合文件分别为Q460.LIS，Q355.LIS和16MnDR.LIS。这三个节点集合文件均可用普通文本编辑器打开，此节点集合文件对于静强度全工况评定方法的有效信息为文件的首列节点编号。文件格式如下：
3、以每行进行一个节点号集合文件设置屈服应力的方法，将3个节点集合文件分解成每个节点集合内每个节点的屈服应力，形成节点编号-屈服应力对应的评定指标文件，并输出该评定指标文件。
节点集合文件屈服应力设置示例：
Q355.lis,418
Q460.lis,460
16MnDR.lis,295
每行的意义为该行文件中含有的节点编号，该节点编号的屈服应力为“，”后的数值。如Q460.lis,460，表示Q460.lis中首列节点，其屈服应力为460MPa。
评定指标文件示例：
此评定指标文件的首列表示节点编号，第二列表示该节点对应的屈服应力，第三列表示该节点处于评定指标文件的第2行的文件集合中。一个模型仅有一个评定指标文件，多个节点集合文件。本实施例中，有3个节点集合文件。
4、在ANSYS软件中对该有限元模型的所有节点，输出每个节点在所有载荷工况下的节点TOP和BOTTOM面应力。本实施例共输出了19种载荷工况的全部载荷工况数据文件，19种载荷工况的TOP面节点应力数据文件：seqvt1.txt～seqvt19.txt。19种载荷工况的BOTTOM面节点应力数据文件：seqvb1.txt～seqvb19.txt。
对于38个节点应力数据文件(19种载荷工况的全部载荷工况数据文件+19种载荷工况的TOP面节点应力数据文件)，计算每个节点的全工况最大应力，并记录各节点全工况最大应力时的载荷工况，按照节点编号顺序输出节点编号-全工况最大应力-全工况最大应力时的载荷工况文件；以节点编号顺序输出节点编号-载荷工况1最大应力-载荷工况2最大应力...载荷工况n最大应力...的全工况应力数据文件，其中载荷工况n最大应力表示该节点在载荷工况n下的TOP面应力和BOTTOM面的最大应力。
全工况最大应力文件示例：
全工况最大应力文件首列表示节点编号，第二列表示节点全工况最大应力，第三列表示全工况最大应力时的载荷工况编号。
全工况应力数据文件示例：
全工况应力数据文件首列表示节点编号，第二列表示载荷工况1的TOP和BOTTOM面的最大应力，第三列表示载荷工况2的TOP和BOTTOM面的最大应力，……，后续第n+1列表示载荷工况n的TOP和BOTTOM面的最大应力。
5、检索所有节点的屈服应力、全工况最大应力，并计算节点的材料利用度，以结构有限元模型的单元和节点组成数据，形成有限元软件的全工况结果输入文件。
ANSYS软件全工况结果输入文件示例：
DESOL,314260,371101,S,X,0.076,31.700,418.200,13,0,2,
DESOL,314260,371103,S,X,0.085,35.663,418.200,13,0,2,
DESOL,314261,371103,S,X,0.085,35.663,418.200,13,0,2,
DESOL,314261,371104,S,X,0.085,35.462,418.200,13,1,2,
DESOL,314262,371104,S,X,0.085,35.462,418.200,13,1,2,
......
全工况结果输入文件首列表示软件结果输入命令(本例为ANSYS)，第2列表示单元编号，第3列表示单元的节点编号，第6列表示第3列节点的全工况材料利用度，第7列表示第3列节点的全工况最大应力，第8列表示第3列节点的屈服应力，第9列表示第3列节点的最大应力工况。
6、通过有限元软件读入全工况结果文件，形成结果输入文件，所述结果输入文件包含各对象的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况，展示车体的应力水平、材料利用度等强度整体性能。
图2是车体的全工况材料利用度及查询的结果示例，表示结构在静强度全工况下的材料利用程度，图形以颜色表示材料利用度的大小，最大的材料度会显示为红色(本例为0.802)，利用度小于1表示结构的静强度载荷工况的最大应力小于其屈服应力，利用度大于1则表示结构会在静强度载荷工况发生永久变形，本例中以节点应力分量SX(同第5步的S，X输入定义)显示。通过图1可全面、直观、高效地进行结构的整体性能评定。全工况最大应力、载荷工况均可进行云图显示和查询，图3是某个节点的载荷工况-应力曲线图，通过图3可全面、直观、高效地进行该节点的所有载荷工况的整体性能评定。
以上所揭露的仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或变型，都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种机械结构静强度全工况评定方法，其特征在于，包括以下步骤：
建立待评定机械结构的有限元模型；
将所述有限元模型的评定对象分成多个对象集合，每个所述对象集合内每个对象的屈服应力相同，并输出多个对象集合文件，每个所述对象集合文件的内容为对应对象集合的对象编号；
将每个所述对象集合文件分解成对应对象集合内每个对象的屈服应力，形成对象编号-屈服应力对应的评定指标文件；
按照对象编号顺序以文件的形式输出每个所述对象在所有载荷工况下的TOP面应力和BOTTOM面应力，形成所有载荷工况下对象的应力数据文件；
根据每个对象在所有载荷工况下的TOP面应力和BOTTOM面应力得到每个对象的全工况最大应力，并记录每个对象在全工况最大应力时对应的载荷工况，形成对象编号-全工况最大应力-全工况最大应力时的载荷工况文件；
检索所有对象的屈服应力和全工况最大应力，计算每个对象的材料利用度；
以所述有限元模型的单元和节点组成数据，形成结果输入文件，所述结果输入文件包含各对象的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况；
读取所述结果输入文件，对整个有限元模型的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况进行可视化显示，并得到模型的可视化结果文件；
优选地，所述评定对象为有限元模型的所有节点或所有单元。
2.如权利要求1所述的机械结构静强度全工况评定方法，其特征在于：每个所述对象的全工况最大应力为：
σi,max＝max(σi,j,T,σi,j,B),j＝1,2,3,…,n
其中，σi,max为第i个对象的全工况最大应力，σi,j,T为第i个对象在第j种载荷工况下的TOP面应力，σi,j,B为第i个对象在第j种载荷工况下的BOTTOM面应力，n为载荷工况的数量。
3.如权利要求1所述的机械结构静强度全工况评定方法，其特征在于：每个所述对象的材料利用度的计算公式为：
其中，μi为第i个对象的材料利用度，σi,max为第i个对象的全工况最大应力，σi,s为第i个对象的屈服应力。
4.如权利要求1～3中任一项所述的机械结构静强度全工况评定方法，其特征在于，所述方法还包括：对所述材料利用度大于第一临界值的对象进行全方位筛选识别，判断材料利用度大于第一临界值的对象是否满足安全要求，如果不满足安全要求，则对该对象进行优化改进；
优选地，所述第一临界值为1。
5.如权利要求1～3中任一项所述的机械结构静强度全工况评定方法，其特征在于，所述方法还包括：对所述材料利用度小于第二临界值的位置进行轻量化处理；
优选地，所述第二临界值为0.5。
6.如权利要求1～3中任一项所述的机械结构静强度全工况评定方法，其特征在于，所述方法还包括：检索某个对象在所有载荷工况下的应力数据，输出该对象的载荷工况-应力曲线。
7.如权利要求1～3中任一项所述的机械结构静强度全工况评定方法，其特征在于，所述方法还包括：输出焊缝工艺文件编制所需的焊缝最大应力和材料利用度数据。
8.一种机械结构静强度全工况评定系统，其特征在于，包括：
模型建立单元，用于建立待评定机械结构的有限元模型；
集合及文件输出单元，用于将所述有限元模型的评定对象分成多个对象集合，每个所述对象集合内每个对象的屈服应力相同，并输出多个对象集合文件，每个所述对象集合文件的内容为对应对象集合的对象编号；
指标文件生成单元，用于将每个所述对象集合文件分解成对应对象集合内每个对象的屈服应力，形成对象编号-屈服应力对应的评定指标文件；
应力数据文件生成单元，用于按照对象编号顺序以文件的形式输出每个所述对象在所有载荷工况下的TOP面应力和BOTTOM面应力，形成所有载荷工况下对象的应力数据文件；
全工况最大应力文件生成单元，用于根据每个对象在所有载荷工况下的TOP面应力和BOTTOM面应力得到每个对象的全工况最大应力，并记录每个对象在全工况最大应力时对应的载荷工况，形成对象编号-全工况最大应力-全工况最大应力时的载荷工况文件；
利用度计算单元，用于检索所有对象的屈服应力和全工况最大应力，计算每个对象的材料利用度；
结果输入文件生成单元，用于以所述有限元模型的单元和节点组成数据，形成结果输入文件，所述结果输入文件包含各对象的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况；
可视化单元，用于读取所述结果输入文件，对整个有限元模型的全工况最大应力、材料利用度以及全工况最大应力时的载荷工况进行可视化显示，并得到模型的可视化结果文件。
9.一种设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：所述处理器执行所述程序时实现权利要求1～7中任一项所述机械结构静强度全工况评定方法。
10.一种存储介质，其上存储有计算机程序，其特征在于：该程序被处理器执行时实现权利要求1～7中任一项所述机械结构静强度全工况评定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011253385.4A CN112395791B (zh) | 2020-11-11 | 2020-11-11 | 机械结构静强度全工况评定方法、系统、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011253385.4A CN112395791B (zh) | 2020-11-11 | 2020-11-11 | 机械结构静强度全工况评定方法、系统、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112395791A true CN112395791A (zh) | 2021-02-23 |
| CN112395791B CN112395791B (zh) | 2024-02-06 |
Family
ID=74600979
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011253385.4A Active CN112395791B (zh) | 2020-11-11 | 2020-11-11 | 机械结构静强度全工况评定方法、系统、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112395791B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116050227A (zh) * | 2023-03-29 | 2023-05-02 | 上海波客实业有限公司 | 一种有限元结构力图显示方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH05288635A (ja) * | 1992-04-10 | 1993-11-02 | Kobe Steel Ltd | Ａｓｍｅ規格に沿った構造物の応力評価装置 |
| CN106570204A (zh) * | 2016-09-23 | 2017-04-19 | 西安交通大学 | 一种基于cpu+gpu异构并行计算的透平机械叶片静强度特性分析方法 |
-
2020
- 2020-11-11 CN CN202011253385.4A patent/CN112395791B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH05288635A (ja) * | 1992-04-10 | 1993-11-02 | Kobe Steel Ltd | Ａｓｍｅ規格に沿った構造物の応力評価装置 |
| CN106570204A (zh) * | 2016-09-23 | 2017-04-19 | 西安交通大学 | 一种基于cpu+gpu异构并行计算的透平机械叶片静强度特性分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| 刘斌;张开林;王秋实;何媛媛;: "复荷载单边止挡轴箱的结构强度评估", 机械, no. 07 * |
| 王博;蔡园武;李传龙;吴昌华;: "9600kW六轴电力机车转向架构架的静强度分析", 计算力学学报, no. 1 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116050227A (zh) * | 2023-03-29 | 2023-05-02 | 上海波客实业有限公司 | 一种有限元结构力图显示方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112395791B (zh) | 2024-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112465748B (zh) | 基于神经网络的裂缝识别方法、装置、设备及存储介质 | |
| CN105893337A (zh) | 用于文本压缩和解压缩的方法和设备 | |
| CN112395791A (zh) | 机械结构静强度全工况评定方法、系统、设备及存储介质 | |
| CN109543744B (zh) | 一种基于龙芯派的多类别深度学习图像识别方法及其应用 | |
| CN110414806B (zh) | 员工风险预警方法及相关装置 | |
| CN110909674B (zh) | 一种交通标志识别方法、装置、设备和存储介质 | |
| CN116049989A (zh) | 转向架数字孪生模型构建方法、系统、电子设备及介质 | |
| CN113609195A (zh) | 报告生成方法、装置、电子设备及存储介质 | |
| CN113780673A (zh) | 离职预测模型的训练方法、装置及离职预测方法、装置 | |
| CN112862345A (zh) | 隐患质检方法、装置、电子设备及存储介质 | |
| CN114312930B (zh) | 基于日志数据的列车运行异常诊断方法和装置 | |
| CN113705201B (zh) | 基于文本的事件概率预测评估算法、电子设备及存储介质 | |
| US20210318949A1 (en) | Method for checking file data, computer device and readable storage medium | |
| CN114091699A (zh) | 一种电力通信设备故障诊断方法及系统 | |
| KR102380021B1 (ko) | 통계추론모델에 기반한 특허 성과지표 출력장치 및 방법 | |
| CN1828579A (zh) | 专利分析系统 | |
| Wiecek et al. | Multi-scenario multi-criteria optimization in engineering design | |
| CN109584235A (zh) | 塑料泡沫餐盒的质量检测方法、装置、服务器及存储介质 | |
| CN110362896A (zh) | 基于非线性度的飞机结构静强度试验数据筛选方法及设备 | |
| CN110322772A (zh) | 一种基于bim的实验室桥梁监测演示教学系统 | |
| CN117315368B (zh) | 一种大规模数据中心智能运维巡检方法 | |
| CN110826942B (zh) | 转运中心效率分析监测方法、装置、设备和存储介质 | |
| WO2024095715A1 (ja) | 環境負荷評価装置および方法 | |
| JP2002056036A (ja) | 構造物の座屈解析法 | |
| Yanagimoto et al. | Working Load in a Container Stack Nonlinearity With Twist Lock Separation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |