CN103207926A - 使用显式有限元分析、采用子循环和质量缩放的组合技术数值模拟产品的结构性能 - Google Patents
使用显式有限元分析、采用子循环和质量缩放的组合技术数值模拟产品的结构性能 Download PDFInfo
- Publication number
- CN103207926A CN103207926A CN2012104229025A CN201210422902A CN103207926A CN 103207926 A CN103207926 A CN 103207926A CN 2012104229025 A CN2012104229025 A CN 2012104229025A CN 201210422902 A CN201210422902 A CN 201210422902A CN 103207926 A CN103207926 A CN 103207926A
- Authority
- CN
- China
- Prior art keywords
- time step
- finite element
- length
- fea
- finite
- 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
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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了使用显式FEA、采用子循环和质量缩放的组合技术数值模拟产品的结构性能的方法和系统。指定代表一产品的FEA模型、以及指定用于该产品的时间推进模拟的最小时间步长(Δtmin)。建立N个元素组,其中第一元素组与作为所需的最小时间步长的Δtmin相关联,而其他元素组根据公式Δti=2i-1Δtmin,n=1,2,...,N与整数倍数的Δtmin相关联。基于有限元各个新的Δtcr将有限元周期性再分类为N个元素组。为了确保在用户指定的最小时间步长下处理FEA模型中的所有有限元,将适当数量的质量缩放应用于在第一元素组中已经变得过小以至于不能维持稳定解的那些有限元。
Description
技术领域
本发明总地涉及计算机辅助工程分析,更具体地,本发明涉及连同子循环和质量缩放的组合技术(combined technique)使用显式有限元分析(finiteelement analysis,FEA)数值模拟产品(例如,汽车、飞机等)的结构性能的方法及系统。
背景技术
有限元分析(FEA)是广泛应用于工业、以建模和解决与复杂系统有关的工程问题的计算机化方法,例如三维非线性结构设计和分析。FEA的名称起源于为考虑中的物体指定几何形状的方式。随着现代数字计算机的出现,FEA已经实现为FEA软件。基本上地,FEA软件提供有几何描述的模型和该模型内每个点处的相关材料属性。在这一模型内,分析中的系统的几何形状由被称为元素的各种大小的实体(solid)、壳(shell)和线束(beam)来代表。各元素的顶点称作为节点。模型包括有限数量的元素,所述元素分配有材料名称从而使各元素与材料属性相关联。该模型因此代表由分析中的物体连同其直接周边(immediate surrounding)所占据的物理空间。FEA软件接下来涉及一表格,在所述表格中将每个材料类型的属性(例如,应力-应变本构方程、杨格模量、泊松比、热导率)制成表。此外指定了该物体边界处的条件(即,负载、物理约束等)。照这样创建物体及其环境的模型。
FEA具有两个求解技术(solution technique):隐式(implicit)有限元分析(“隐式法”)和显式有限元分析(“显式法”)。两个方法均用于求解瞬态动力学方程且因此获得该方程的均衡解。各方法经由离散时间间隔或时间间隔Δt从时间(t)推进至时间((t+Δt)。上述方法有时被称为时间推进模拟,其包含许多连续时间步长或求解循环(solution cycle)。
本发明涉及显式法,只要时间步长非常小所述显式法则是稳定的—特别地,时间间隔必须小于弹性波从一元素的一侧传播到另一侧所花费的时间。在显式法中维持稳定解的最大时间步长称作为临界时间步长Δtcr。弹性波的速度是材料质量和结构的刚度的函数,有限元和元素大小或元素尺寸代表所述结构的刚度。对具有大致(substantially)相似材料的FEA模型而言,最小的元素通常控制临界时间步长。
FEA模型中,甚至一个大致较小的元素可以引起该FEA模型内大多数元素的临界时间步长不必要地变小。往往采用更精细的网格对结构的特定部分建模,以捕获更详细的结构响应。例如,与用于建模车身的壳元素形成对比的是,可能采用非常小的实体元素建模方向盘。两种大小间的比率可以是100或更大,从而使得采用100或小更多倍的时间步长分析整个结构(即轿车)。在另一情况下,当任何元素在时间推进模拟的中途变形为变得太小时,该问题也有可能发生。如果存在泡沫材料的话这会经常发生,这是由于在增大弹性波速度和降低最短尺寸时它们显著压缩。结果,对剩余模拟而言需要的时间步长将非常小。模拟不仅确实变得非常耗费时间,而且也变得不切实际。
针对这一问题存在许多现有技术的方法。然而没有一种方法提供令人满意的求解。它们的其中一种称作为子循环,其允许在有限元模型的不同部分使用不同时间步长。然而,子循环需要周期性地将各元素分类/再分类为不同二进制文件(bin),所述二进制文件随着不同时间步长动态变化。另一方法称作为质量缩放,其通过人工增加其质量密度而使一元素的临界时间步长增加。然而,将质量缩放人工应用于整个FEA模型可能会以不需要的方式改变结构的动态性能(即,人工地较高的质量)。
因此,需要使用显式FEA、采用子循环和质量缩放的组合技术来数值模拟产品的结构性能的方法及系统,以便高效且有效地执行模拟。
发明内容
本发明公开了使用显式FEA、采用子循环和质量缩放的组合技术数值模拟产品的结构性能的系统、方法和软件产品。根据本发明的示范性实施例,在其上安装有有限元分析应用模块的计算机系统中定义和接收代表一产品(例如,汽车)的有限元分析(FEA)模型。FEA模型包括多个节点和有限元(例如,壳元素、实体元素、泡沫元素(foam element)等)。另外,还指定产品的时间推进模拟的最小时间步长(Δtmin)(例如,通过用户(即,设计或帮助设计该产品的工程师、科学家))
基于其属性(即,刚度、质量和大小)计算FEA模型中每个有限元的临界时间步长。随后通过其邻近有限元或相邻有限元的较小的临界时间步长更新每个有限元的临界时间步长(Δtcr)。使用在以下“具体实施方式”分段中列示的Eq.4计算基于元素属性的临界时间步长。
基于各个所需的最小时间步长将FEA模型的有限元分类为多个元素组。建立多个元素组,其中第一元素组与作为所需的最小时间步长Δtmin相关联,而其他元素组与整数倍数的Δtmin相关联。特别地,公式(2n-1Δtmin,n=1,2,..,N)用于各个元素组,N是大于1的完整数。
当使用显式FEA的上述FEM模型用于时间推进模拟时,根据该组所需的最小时间步长处理或求解每个元素组的有限元。这意味着,仅在与上一元素组所需的时间步长(2N-1Δtmin)对应的主求解循环中处理上一元素组中的有限元,而在相应的子求解循环(即,Δtmin,2Δtmin,4Δtmin等)中处理其他元素。
为了确保在所需的最小时间步长(Δtmin)下处理第一元素组的FEA模型中的所有有限元,将质量缩放应用于第一元素组中过小或已变得过小以至于不能维持稳定解的那些有限元。
此外,可指定将一些有限元从周期性再分类中去除,例如代表汽车的金属部件的有限元(例如,车身)。为了确保可在初始分配的组中维持这些元素的稳定解,在时间推进模拟的过程中无论何时需要均可应用质量缩放。
在每个主求解循环结束时,基于它们各自的新的临界时间步长对有限元进行再分类,以确保在下一系列的子求解循环中获得稳定解。
一旦仔细阅读对本发明实施例的以下详细描述、连同加以考虑各附图,本发明的目标、特征和优点将变得显而易见。
附图说明
结合以下描述、所附权利要求和附图将更好地理解本发明的所有特征、方面和优点,所述附图如下:
图1是根据本发明一实施例的、阐述连同子循环和质量缩放的组合技术使用显式有限元分析(FEA)数值模拟产品的结构性能的示范性过程的流程图;
图2A-2C是依照本发明一个实施例的、图形化地示出用于建立FEA模型的有限元的各个临界时间步长的示范性程序的二维示意图;
图3A-3B是依照本发明一个实施例的、示出示范性元素组和各个所需的最小时间步长的示意图;
图4是依照本发明一个实施例的、示出可用在代表产品(例如轿车)的FEA模型中的两个示范性有限元的示意图;
图5是依照本发明一个实施例的、示出了质量缩放应用于的显式有限元分析中使用的示范性四边形壳元素以及相应质量矩阵和节点加速度矢量的示意图;
图6是根据本发明一个实施例的、示出了使用显式有限元分析的轿车碰撞的示范性时间推进模拟结果的示意图;以及
图7是示出了示范性计算机的突出组件的功能图,在所述示范性计算机中可实现本发明的一个实施例。
具体实施方式
图1是根据本发明一实施例的、阐述使用显式有限元分析(FEA)、采用子循环和质量缩放的组合技术数值模拟产品的结构性能的示范性过程100的流程图。优选地,连同其他图理解过程100,且过程100以软件方式实现。
通过接收有限元分析(FEA)模型开始过程100,所述有限元分析模型在步骤102代表一产品或结构(例如,汽车、飞机等)。FEA模型包括定义了产品的、具有相关联材料属性的多个节点和多个有限元。有限元的类型可包含但不限于实体、壳、线束、构架(truss)等(图4示出了示范性壳元素401和示范性实体元素402)。另外也指定了该产品的时间推进模拟的最小时间步长(Δtmin)(例如,通过用户(即设计或帮助设计该产品的工程师、科学家))。例如,假设最小时间步长确保时间推进模拟可在合理时间帧内因实际用途/实际关联而实施,负责汽车设计的工程师/科学家则处于一时间约束下(例如,通宵的计算机模拟)。另外,可指定将一个或多个可选集合的有限元从周期性再分类中去除,例如代表汽车的金属部件的有限元。为了从周期性再分类中去除这些有限元,在时间推进模拟的过程中无论何时有必要都能应用质量缩放。
接下来在步骤104,建立N个元素组。第一元素组与作为所需的最小时间步长Δtmin相关联,而其他元素组与整数倍数的Δtmin相关联。特别地,每个元素组依照以下公式对应于所需的最小时间步长:
Δti=2i-1Δtmin,i=1,2,..,N
换言之,第一组与(Δtmin)相关联,第二组与(2Δtmin)相关联,第三组与(4Δtmin)相关联,第N组与(2N-1Δtmin)相关联,如此等等。
在步骤106,基于有限元的各个临界时间步长将FEA模型的有限元分类为N个元素组。N是大于1的完整数或整数。每个有限元的临界时间步长大于或等于该元素组的所需的最小时间步长。
为了证明如何建立每个有限元的临界时间步长,图2B-2C示出了一示范性程序。FEA模型200包含各种大小的有限元。首先,使用Eq.4计算用于维持显式FEA中的稳定解的每个有限元的元素时间步长。图2A中示出的是,需要较小元素时间步长的较小有限元202由“1”表示,而较大有限元204由“2”表示。接着,通过分配与每个节点连接的所有元素最小的元素时间步长来确定每个节点的节点临界时间步长。图2B示出了FEA模型200的节点临界时间步长,在各个节点处为具有“1”212和“2”214的圆圈。注意用“B”210表示的有限元,其是在两个不同大小的边界附近的有限元(即,两个不同的临界时间步长)。“B”或边界有限元210包含具有不同的节点临界时间步长的节点。当应用子循环技术时,需要在较小时间步长下处理这些边界有限元210。图2C中示出的元素组“1”222和元素组“2”224是在FEA模型200上应用子循环技术的结果。
在这一实例中,元素组“1”222中的有限元分类为第一元素组,而组“2”224中的有限元分类为第二元素组(即N=2)。将使用Δtmin处理或求解第一元素组、使用(2Δtmin)处理或求解第二元素组。
在有限元已分类为N个元素组后,过程100检查第一元素组内的任何有限元是否需要质量缩放。将适当数量的质量缩放应用于第一元素组中视为过小以至于不能维持稳定解的那些有限元。例如,当元素的临界时间步长(Δtcr)小于所需的最小时间步长(Δtmin)时,可使用(Δtmin-Δtcr)的差和以下的Eq.6计算增大的适当数量的质量密度。注意的是,质量缩放仅应用于第一元素组。
接下来在步骤110,连同子循环和质量缩放的组合技术、使用显式FEA实施FEA模型所代表的产品的时间推进模拟。在每个主求解循环内,时间推进模拟包含多个主求解循环和多个子求解循环。每个主求解循环包含一个或多个子求解循环。主求解循环与(2N-1Δtmin)相关联,而子求解循环与较小倍数的Δtmin相关联。在这些求解循环中获得模拟产品响应。
在预设时间间隔结束时(例如,每个求解循环、一个主求解循环),过程100首先将质量缩放应用于指定从周期性再分类中去除的那些有限元,以确保可在初始分配的组内维持稳定解。接着在步骤112,过程100基于所有有限元各自的新临界时间步长对所有有限元进行再分类。
在第一元素组中,适当数量的质量缩放应用于已变得过小以至于不能维持稳定解的那些有限元。人工增大视为过小的任何元素的质量密度,以确保元素临界时间步长至少等于Δtmin。
最后,当在决定114处达到结束条件时则结束过程100,例如,模拟时间已经超出预设总模拟时间(例如,用户指定的总模拟时间)。否则,过程100移回步骤110,并重复直至决定114变为真(“是”)。
在将有限元分类为各个不同元素组的更通常的表述中,图3A示出了依照本发明一个实施例的、用于元素组管理的数据结构。有限元303a的分类导致了N个元素组301,其中N是大于1的完整数或整数。每个元素组301包含一列表的有限元303a。例如,第一组包含有限元310a-e,第二组包含元素320a-b,第三组包含元素330a-d,且第N组包含元素350a-m。各个元素组301需要的最小时间步长302是彼此的整数倍数。根据一个实施例,第一组需要最小时间步长Δt1,第二元素组需要最小时间步长2Δt1,第三元素组需要4Δt1,且第N个元素组需要2N-1Δt1。
当使用子循环和质量缩放的组合技术时,在时间推进模拟的过程中对有限元进行周期性再分类。在图3B所示的再分类之后,有限元303b可与不同元素组301相关联。
许多不同类型的有限元可用于FEA模型,例如,图4中示出的4-节点四边形壳元素401和8-节点实体元素402。
现在参考图5,其示出了质量缩放技术。在显式有限元分析中,系统(例如,结构)的运动方程列示如下:
a=M-1(fext-fint) (Eq.1)
其中a是代表节点加速度的未知加速度矢量,M是非对角项(off-diagonalterm)中具有零的集中(lumped)对角矩阵,fext和fint分别是外作用力和内节点力。力项fext-fint被称作为残余力,其包含但不限于作用力、接触力、内应力。由于M仅具有对角项,M-1的计算可以简单地按照除法、而不是代价高的矩阵求逆来执行。图5示出了一示范性的4-节点壳元素500的对角质量矩阵501和相应的加速度矢量520。在这一示范性质量矩阵510中,四个节点的每个具有由对角项mo代表的相同数量的质量(例如,集中节点质量501、502、503或504)。所有的非对角项506为零。加速度矢量520包含在每个节点处(即,a1521、a2522、a3523和a4524)的加速度。采用对角项乘以相应加速度来计算每个节点的惯性力,其如下所示:
其中N是元素中节点的数目。对于四边形壳元素而言,N等于4;对于六面体实体元素而言,N等于8。依照如下所示的Eq.2计算元素500的四个节点处的惯性力。
f1=moa1+(0)a2+(0)a3+(0)a4=moa1
f2=(0)a1+moa2+(0)a3+(0)a4=moa2
f3=(0)a1+(0)a2+moa3+(0)a4=moa3
f4=(0)a1+(0)a2+(0)a3+moa4=moa4 (Eq.3)
为了简化阐述和描述,图5中仅示出了每个节点处的一个元件(component)。在笛卡尔坐标系中,每个节点处有三个平移质量,因此元件的总数目为12(即,三乘四)。
采用时间积分的显式法求解Eq.1,逐一(at a time)采用一个时间步长执行所述显式法。只要时间步长Δt小于临界值Δtcr(即,临界时间步长),显式法则仅可产生稳定解;所述临界值由以下公式确定:
ωmax是系统的最大本征频率。可通过求解以下本征值问题获得本征频率:
det[K-ωi 2M]=0 (Eq.5)
其中K和M分别是系统的m×m刚度和质量矩阵,m是系统的自由度的数目。最大本征频率ωmax通常是一较大数(例如,大约为106Hz),因此临界时间步长Δtcr是非常小的数目(例如,大约为数微秒)。这一小时间步长限制了在模拟准静态过程(例如金属成形)中显式法的有用性。
增大临界时间步长的一种方式是降低最大本征频率ωmax,其可通过增大系统的质量得以实现。本发明允许用户指定所需的时间步长Δt。随后,确保显式法的稳定性所需的质量可如下所示来计算:
m=αk(Δt)2 (Eq.6)
其中m是所需的质量,k是刚度,α是取决于特定积分方案的常量。在有限元分析模型中,每个元素具有取决于几何形状和材料属性的刚度。因此,可使用Eq.6为每个求解循环的每个元素计算所需的质量。临界时间步长是质量的平方根的函数,因此增大临界时间步长需要更大的质量。下文中所需的质量也称作为附加质量,质量缩放指采用附加质量调节系统的真实质量的程序。
可采用许多技术应用质量缩放。其中一个最基本的方法是以试错法方式改变整个结构的材料质量密度。通常可由FEA的用户来执行上述程序、而不需要任何自动化计算,在数次材料质量密度调整后,可采用合适的材料质量密度获得FEA求解。
在应用模块(例如,FEA软件)中实现自动化过程。所需的时间步长通常由FEA的用户指定,应用模块包括某些逻辑或子程序,其将计算自动增大质量密度的因素,以便所需的时间步长可用于工程模拟以获得一求解方法。可在元素-对-元素(element-by-element)的基础上实施自动化程序。对FEA模型的每个元素而言,可基于元素的质量和几何形状计算临界时间步长。基于Eq.4,可通过降低最大本征频率增大临界时间步长Δtcr,这等同于增大了在其他方面未改变的元素的质量。
如果特定元素的临界时间步长小于所需的时间步长,可计算出附加质量并增加到特定元素,以确保该元素的临界时间步长变得至少与所需的时间步长一样大。当那些元素的临界时间步长大于或等于可由用户指定的所需的时间步长(例如,输入文件中的参数)时,一些元素可能不需要任何附加质量。
现在参考图6,其是示出了使用显式有限元分析的轿车碰撞的示范性时间推进模拟结果的示意图。换言之,数值模拟的图形显示由数值模拟的撞击事件的有限元分析造成。模拟中包括两个结构:汽车和障碍物。在有限元分析模型中可表示两个该结构。可采用超过1,000,000个有限元(例如,壳元素、实体元素和线束元素)、采用相当数量的节点对汽车建模。根据本发明的一个方面,采用子循环和质量缩放的组合技术增强显式有限元分析。
根据一个方面,本发明涉及能够实施本文所描述的功能性的一个或多个计算机系统。图7中示出了计算机系统700的实例。计算机系统700包括一个或多个处理器,例如处理器704。处理器704与计算机系统内通信总线702连接。依据这一示范性计算机系统描述各个软件实施例。在阅读这一描述之后,对本领域技术人员而言,如何使用其他计算机系统和/或计算机体系结构来实现本发明将变得显而易见。
计算机系统700还包括主存储器708、优选为随机存取存储器(randomaccess memory,RAM),其也可包括辅助存储器710。例如,辅助存储器710可包括一个或多个硬盘驱动器712和/或一个或多个卸除式存储驱动器714,其代表为软盘驱动器、磁带驱动器、光盘驱动器等。卸除式存储驱动器714以公知方式读取和/或写入卸除式存储单元718。卸除式存储单元718的代表为由卸除式存储驱动器714读取和写入的软盘、磁带、光盘等。应该理解的是,卸除式存储单元718包含具有存储其上的计算机软件和/或数据的计算机可用存储介质。
在替代性实施例中,辅助存储器710可包含其他相似装置,从而允许将计算机程序或其他指令载入计算机系统700中。例如,这种装置可能包含卸除式存储单元722和接口720。上述装置的实例可包含(例如在视频游戏设备中存在的)程序盒式存储器和盒式接口、卸除式存储芯片(例如可擦可编程只读存储器(Erasable Programmable Read-Only Memory,EPROM)、通用串行总线(Universal Serial Bus,USB)闪存、或PROM)和相关插口、和其他允许将软件和数据由卸除式存储单元722传递到计算机系统700的卸除式存储单元722和接口720。一般而言,由操作系统(operating system,OS)软件来控制和协调计算机系统700,所述操作系统软件执行诸如进程调度、内存管理、联网和I/O服务的任务。
还可能存在与总线702连接的通信接口724。通信接口724允许将软件和数据在计算机系统700和外部设备间进行传递。通信接口724的实例可包含调制解调器、网络接口(例如以太网卡)、通信端口、个人计算机存储卡国际协会(Personal Computer Memory Card International Association,PCMCIA)插槽和卡等。计算机700通过数据网络、基于一组特殊规则(即协议)与其他计算设备通信。其中一种常见协议是在因特网中普遍使用的TCP/IP(传输控制协议/互联网协议)。一般而言,通信接口724对于以下过程实施管理:将数据文件组装为在数据网络上传输的较小数据包、以及将接收的数据包重组为原始数据文件。另外,通信接口724处理每个数据包的地址码(address part),以便其到达正确目标或拦截目标是计算机700的数据包。在这一文件中,用语“计算机程序介质”和“计算机可用介质”通常用于指以下媒介:例如,卸除式存储驱动器714和/或安装在硬盘驱动器712中的硬盘。这些计算机程序产品是为计算机系统700提供软件的装置。本发明涉及这些计算机程序产品。
计算机系统700还可包含输入/输出(I/O)接口730,其为计算机系统700提供存取监控器、键盘、鼠标、打印机、扫描仪、绘图仪和类似物。
在主存储器708和/或辅助存储器710中,计算机程序(也称为计算机控制逻辑)存储为应用模块706。还可通过通信接口724接收计算机程序。当执行时,这些计算机程序使得计算机系统700能够执行本文所讨论的本发明的特征。尤其地,当执行时,计算机程序使得处理器704执行本发明的特征。因此,上述计算机程序代表计算机系统700的控制器。
在使用软件实现本发明的实施例中,软件可存储在计算机程序产品中,并使用卸除式存储驱动器714、硬驱动器712或通信接口724将其载入计算机系统700中。当处理器704执行应用模块706时,应用模块706使处理器704执行本文所描述的本发明的各功能。
主存储器708可载有一个或多个应用模块706,在有或没有通过I/O接口730的用户输入的情况下,一个或多个处理器704可执行一个或多个应用模块,从而实现希望实现的任务。在操作中,当至少一个处理器704执行其中一个应用模块706时,结果得到计算并保存在辅助存储器710中(即硬盘驱动器712)。以文本或图形表征的方式、通过I/O接口730向用户报导有限元分析的状态。
在一个实施例中,应用模块706用于促进代表产品的有限元分析模型的创建。应用模块706进一步允许连同显式有限元分析使用子循环和质量缩放的组合技术。在另一实施例中,应用模块706用于促进某些有限元的质量缩放,以维持稳定解。
尽管结合其特定实施例对本发明进行描述,但这些实施例仅仅是对本发明的阐释而非限制。本文向本领域技术人员暗示了对特别公开的示范性实施例的各个调整或变化。例如,虽然为如何实现子循环和质量缩放的组合技术显示和描述了四边形的壳有限元,但可替代使用其他类型的有限元,例如8-节点六面体元素、4-节点四面体元素、3-节点三角形元素等。进一步地,虽然具有两个元素组的简单FEA模型200已示出和描述了示出分类方案,但其他更复杂的FEA模型也可用于实现相同目的。总之,本发明的范围不应受限于本文所公开的特定示范性实施例,所有已提示过本领域技术人员的所有调整均应包含在这一申请的精神和追求的范围内,并且也包含在所附权利要求的范围内。
Claims (14)
1.一种使用显式有限元分析(FEA)、采用子循环和质量缩放的组合技术数值模拟产品的结构性能的方法,其特征在于,所述方法包括:
在计算机系统中接收FEA模型和最小时间步长(Δtmin),所述计算机系统具有安装其上的显式FEA应用模块,所述FEA模型代表由具有相关联材料属性的多个节点和多个有限元定义的产品,所述最小时间步长用于所述产品的时间推进模拟;
建立N个元素组,其中第一元素组与所述最小时间步长相关联,剩余的所述N个元素组与整数倍数的所述最小时间步长相关联,其中N是大于1的完整数;
基于所述有限元的各个临界时间步长将所述有限元分类为所述N个元素组;
将适当数量的质量缩放应用于在所述第一元素组中过小以至于不能维持稳定解的那些有限元;以及
通过使用FEA模型、采用在所述计算机系统中执行的显式FEA应用模块实施时间推进模拟,在许多主求解循环和每个主求解循环的一个或多个子求解循环中获得模拟结构性能;所述显式FEA应用模块基于所述有限元的各个新的元素临界时间步长、通过所述有限元的周期性再分类应用子循环和质量缩放的组合技术,以及所述显式FEA应用模块将适当数量的质量缩放应用于在所述第一元素组中已变得过小以至于不能维持稳定解的那些有限元。
2.如权利要求1所述的方法,其特征在于,所述方法还包括指定将所述有限元的一个或多个从所述周期性再分类中去除。
3.如权利要求1所述的方法,其特征在于,在所述每个主求解循环的所述周期性再分类之前,所述方法还包括无论何时有必要都将质量缩放应用于指定从所述周期性再分类中去除的所述有限元的一个或多个上。
4.如权利要求1所述的方法,其特征在于,基于公式Δti=2i-1Δtmin计算所述整数倍数的最小时间步长,其中i=1,2,...,N。
5.如权利要求1所述的方法,其特征在于,所述方法还包括基于所述每个有限元的属性计算每个有限元的元素时间步长;通过分配与所述每个节点连接的那些有限元的最小的元素时间步长建立每个节点的节点临界时间步长;以及将所述有限元的节点的最小的节点临界时间步长指定为所述每个有限元的临界时间步长。
6.如权利要求1所述的方法,其特征在于,所述组合技术还包括在那些有限元的所述每个主求解循环结束时使获得的模拟结构性能同步,所述那些有限元在所述子求解循环中需要求解。
7.如权利要求1所述的方法,其特征在于,所述方法还包括当达到预设结束条件时结束所述时间推进模拟。
8.如权利要求7所述的方法,其特征在于,所述预设结束条件包括用户指定的总模拟时间。
9.一种使用显式有限元分析(FEA)、采用子循环和质量缩放的组合技术数值模拟产品的结构性能的系统,其特征在于,所述系统包括:
主存储器,用于为显式有限元分析(FEA)应用模块存储计算机可读代码;
与所述主存储器连接的至少一个处理器,所述至少一个处理器执行所述主存储器中的计算机可读代码,以使得所述显式FEA应用模块通过以下方法执行各操作:
接收FEA模型和最小时间步长(Δtmin),所述FEA模型代表由具有相关联材料属性的多个节点和多个有限元定义的产品,所述最小时间步长用于所述产品的时间推进模拟;
建立N个元素组,其中第一元素组与所述最小时间步长相关联,剩余的所述N个元素组与整数倍数的所述最小时间步长相关联,其中N是大于1的完整数;
基于所述有限元的各个临界时间步长将所述有限元分类为所述N个元素组;
将适当数量的质量缩放应用于在所述第一元素组中过小以至于不能维持稳定解的那些有限元;以及
通过使用FEA模型、采用在所述系统中执行的显式FEA应用模块实施时间推进模拟,在许多主求解循环和每个主求解循环的一个或多个子求解循环中获得模拟结构性能;所述显式FEA应用模块基于所述有限元的各个新的元素临界时间步长、通过所述有限元的周期性再分类应用子循环和质量缩放的组合技术,以及所述显式FEA应用模块将适当数量的质量缩放应用于在所述第一元素组中已变得过小以至于不能维持稳定解的那些有限元。
10.如权利要求9所述的系统,其特征在于,还包括指定将所述有限元的一个或多个从所述周期性再分类中去除。
11.如权利要求9所述的系统,其特征在于,还包括在所述每个主求解循环的所述周期性再分类之前,无论何时有必要都将质量缩放应用于指定从所述周期性再分类中去除的所述有限元的一个或多个上。
12.如权利要求9所述的系统,其特征在于,基于公式Δti=2i-1Δtmin计算所述整数倍数的最小时间步长,其中i=1,2,...,N。
13.如权利要求9所述的系统,其特征在于,还包括基于所述每个有限元的属性计算每个有限元的元素时间步长;通过分配与所述每个节点连接的那些有限元的最小的元素时间步长建立每个节点的节点临界时间步长;以及将所述有限元的节点的最小的节点临界时间步长指定为所述每个有限元的临界时间步长。
14.如权利要求9所述的方法,其特征在于,所述组合技术还包括在那些有限元的所述每个主求解循环结束时使获得的模拟结构性能同步,所述那些有限元在所述子求解循环中需要求解。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/351,945 | 2012-01-17 | ||
US13/351,945 US8855977B2 (en) | 2012-01-17 | 2012-01-17 | Numerically simulating structural behaviors of a product using explicit finite element analysis with a combined technique of mass scaling and subcycling |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103207926A true CN103207926A (zh) | 2013-07-17 |
CN103207926B CN103207926B (zh) | 2017-03-01 |
Family
ID=48755147
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210422902.5A Active CN103207926B (zh) | 2012-01-17 | 2012-10-30 | 使用组合技术数值模拟产品的结构性能的方法和系统 |
Country Status (2)
Country | Link |
---|---|
US (1) | US8855977B2 (zh) |
CN (1) | CN103207926B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104636534A (zh) * | 2014-12-08 | 2015-05-20 | 浙江工业大学 | 一种提高管材弯曲成型过程仿真效率的方法及仿真方法 |
CN106611078A (zh) * | 2015-10-27 | 2017-05-03 | 利弗莫尔软件技术公司 | 产品的具有时间步长控制方案的高效显式有限元分析 |
CN110427668A (zh) * | 2019-07-19 | 2019-11-08 | 大连理工大学 | Cfrp三维细观切削仿真质量缩放系数的快速选取方法 |
CN110779795A (zh) * | 2019-11-04 | 2020-02-11 | 中国石油大学(华东) | 裂缝性储层地质力学建模网格单元大小确定方法 |
CN111199117A (zh) * | 2018-10-24 | 2020-05-26 | 通用汽车环球科技运作有限责任公司 | 用于加快显式有限元分析的自适应时间缩放 |
CN117709129A (zh) * | 2024-02-05 | 2024-03-15 | 国家超级计算天津中心 | 多尺度仿真方法、装置、设备及存储介质 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9098657B2 (en) * | 2013-02-05 | 2015-08-04 | Livermore Software Technology Corp. | Contact surface definition creation involving low order and quadratic finite elements in a numerical simulation of an impact event |
US10227145B2 (en) | 2015-02-27 | 2019-03-12 | Space Systems/Loral, Llc | Truss structure |
US10407189B1 (en) | 2016-10-27 | 2019-09-10 | Space Systems/Loral, Llc | Spacecraft exoskeleton truss structure |
US11286062B1 (en) | 2016-10-27 | 2022-03-29 | Space Systems/Loral, Llc | Spacecraft exoskeleton truss structure |
US10633123B2 (en) | 2017-04-05 | 2020-04-28 | Space Systems/Loral, Llc | Exoskeletal launch support structure |
EP3537316A1 (en) * | 2018-03-09 | 2019-09-11 | Ningbo Geely Automobile Research & Development Co., Ltd. | A method for optimizing the construction of a car body |
CN111079224B (zh) * | 2019-11-27 | 2023-06-13 | 广汽零部件有限公司 | 一种基于cae和显式动力学技术的汽车座椅刚度分析方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090228246A1 (en) * | 2008-03-05 | 2009-09-10 | Livermore Software Technology Corporation | Methods and systems of engineering analysis using a hybrid approach with FEM and adaptive SPH |
US20090299702A1 (en) * | 2008-05-27 | 2009-12-03 | Livermore Software Technology Corporation | Systems and Methods of Limiting Contact Penetration in Numerical Simulation of Non-linear Structure Response |
CN101685474A (zh) * | 2008-09-25 | 2010-03-31 | 利弗莫尔软件技术公司 | 有限元分析中初始化螺栓预紧力的改进方法和系统 |
CN101739490A (zh) * | 2008-11-26 | 2010-06-16 | 利弗莫尔软件技术公司 | 有限元分析中点焊焊点故障的确定方法和系统 |
CN101858713A (zh) * | 2009-04-01 | 2010-10-13 | 利弗莫尔软件技术公司 | 有限元分析中的爆炸模拟 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070244676A1 (en) * | 2006-03-03 | 2007-10-18 | Li Shang | Adaptive analysis methods |
JP5062249B2 (ja) * | 2007-03-07 | 2012-10-31 | 富士通株式会社 | 解析装置、解析方法及び解析プログラム |
DE102007012633A1 (de) * | 2007-03-16 | 2008-09-18 | Bayerische Motoren Werke Aktiengesellschaft | Automatisches Erzeugen einer Vernetzung eines Komponentenmodells |
US8401827B2 (en) * | 2008-04-14 | 2013-03-19 | Daa Draexlmaier Automotive Of America Llc | Processing device and method for structure data representing a physical structure |
US8374833B2 (en) * | 2009-09-29 | 2013-02-12 | Livermore Software Technology Corporation | Hybrid element enabling solid/SPH coupling effect |
US8494819B2 (en) * | 2010-08-25 | 2013-07-23 | Livermore Software Technology Corp. | Efficient data management for shell finite elements representing layered composite materials |
US20140049534A1 (en) * | 2012-08-14 | 2014-02-20 | Livermore Software Technology Corp | Efficient Method Of Rendering A Computerized Model To Be Displayed On A Computer Monitor |
-
2012
- 2012-01-17 US US13/351,945 patent/US8855977B2/en active Active
- 2012-10-30 CN CN201210422902.5A patent/CN103207926B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090228246A1 (en) * | 2008-03-05 | 2009-09-10 | Livermore Software Technology Corporation | Methods and systems of engineering analysis using a hybrid approach with FEM and adaptive SPH |
US20090299702A1 (en) * | 2008-05-27 | 2009-12-03 | Livermore Software Technology Corporation | Systems and Methods of Limiting Contact Penetration in Numerical Simulation of Non-linear Structure Response |
CN101685474A (zh) * | 2008-09-25 | 2010-03-31 | 利弗莫尔软件技术公司 | 有限元分析中初始化螺栓预紧力的改进方法和系统 |
CN101739490A (zh) * | 2008-11-26 | 2010-06-16 | 利弗莫尔软件技术公司 | 有限元分析中点焊焊点故障的确定方法和系统 |
CN101858713A (zh) * | 2009-04-01 | 2010-10-13 | 利弗莫尔软件技术公司 | 有限元分析中的爆炸模拟 |
Non-Patent Citations (1)
Title |
---|
马晓春等: "有限元数值模拟技术在汽车冲压件成形中的应用", 《浙江工业大学学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104636534A (zh) * | 2014-12-08 | 2015-05-20 | 浙江工业大学 | 一种提高管材弯曲成型过程仿真效率的方法及仿真方法 |
CN104636534B (zh) * | 2014-12-08 | 2017-12-15 | 浙江工业大学 | 一种提高管材弯曲成型过程仿真效率的方法及仿真方法 |
CN106611078A (zh) * | 2015-10-27 | 2017-05-03 | 利弗莫尔软件技术公司 | 产品的具有时间步长控制方案的高效显式有限元分析 |
CN106611078B (zh) * | 2015-10-27 | 2021-03-26 | 利弗莫尔软件技术公司 | 产品的具有时间步长控制方案的高效显式有限元分析 |
CN111199117A (zh) * | 2018-10-24 | 2020-05-26 | 通用汽车环球科技运作有限责任公司 | 用于加快显式有限元分析的自适应时间缩放 |
CN111199117B (zh) * | 2018-10-24 | 2023-09-05 | 通用汽车环球科技运作有限责任公司 | 用于加快显式有限元分析的自适应时间缩放 |
CN110427668A (zh) * | 2019-07-19 | 2019-11-08 | 大连理工大学 | Cfrp三维细观切削仿真质量缩放系数的快速选取方法 |
CN110427668B (zh) * | 2019-07-19 | 2021-01-19 | 大连理工大学 | Cfrp三维细观切削仿真质量缩放系数的快速选取方法 |
CN110779795A (zh) * | 2019-11-04 | 2020-02-11 | 中国石油大学(华东) | 裂缝性储层地质力学建模网格单元大小确定方法 |
CN110779795B (zh) * | 2019-11-04 | 2022-05-10 | 中国地质大学(武汉) | 裂缝性储层地质力学建模网格单元大小确定方法 |
CN117709129A (zh) * | 2024-02-05 | 2024-03-15 | 国家超级计算天津中心 | 多尺度仿真方法、装置、设备及存储介质 |
CN117709129B (zh) * | 2024-02-05 | 2024-05-14 | 国家超级计算天津中心 | 多尺度仿真方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
US20130185030A1 (en) | 2013-07-18 |
US8855977B2 (en) | 2014-10-07 |
CN103207926B (zh) | 2017-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103207926A (zh) | 使用显式有限元分析、采用子循环和质量缩放的组合技术数值模拟产品的结构性能 | |
US11520956B2 (en) | Systems and methods for automatically realizing models for co-simulation | |
Hartmann et al. | Model order reduction a key technology for digital twins | |
Jansson et al. | CUBE: A scalable framework for large-scale industrial simulations | |
Poole et al. | High-fidelity aerodynamic shape optimization using efficient orthogonal modal design variables with a constrained global optimizer | |
CN103207927A (zh) | 使用显式有限元分析、采用质量缩放增强型子循环技术数值模拟产品的结构性能 | |
Farhat et al. | Simulation of compressible viscous flows on a variety of MPPs: computational algorithms for unstructured dynamic meshes and performance results | |
JP2009289260A (ja) | 非線形の構造的応答の数値シミュレーションにおける接触貫通を制限するシステムおよび方法 | |
CN104573166A (zh) | 用于在固体力学中模拟物理行为的光滑粒子伽辽金公式 | |
Röck | Hardware in the loop simulation of production systems dynamics | |
CN101763448B (zh) | 在多处理器计算机系统中模拟撞击事件的改进方法 | |
EP2196926A2 (en) | Solid finite elements suitable for simulating large deformations and/or rotations of a structure | |
Papoutsis-Kiachagias et al. | Aerodynamic optimization of car shapes using the continuous adjoint method and an RBF morpher | |
EP2784704A1 (en) | Multi-component computational fluid dynamics simulations | |
EP2357578A1 (en) | Multiscale substructures in finite element analysis | |
CN106611078B (zh) | 产品的具有时间步长控制方案的高效显式有限元分析 | |
Noor | Computational structures technology: leap frogging into the twenty-first century | |
JP7042562B2 (ja) | 非圧縮過渡および定常状態ナビエ-ストークス方程式のための最適圧力射影法 | |
US7386428B1 (en) | Method and system for guided cable contact in finite element analysis | |
US8676551B1 (en) | Multi-solver simulation of dynamic systems in a modeling environment | |
US20150213165A1 (en) | Methods And Systems For Conducting Design Sensitivity Analysis | |
US20230034245A1 (en) | Lbm-based fluid analysis simulation device, method, and computer program | |
CN111199085A (zh) | 抽象表达模型及其状态的方法 | |
Pappalardo et al. | A General Method for Performing an Integrated CAD-MBD-FEM Analysis | |
Bungartz et al. | Fluid-acoustics interaction on massively parallel systems |
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 |