2021年ABAQUS空气弹簧(基于表面的流体模型)

1、*欧阳光明*创编2021.03.07*欧阳光明*创编2021.03.0711.5基于表面的流体模型欧阳光明(2021.03.07)“基于流体腔表面：概述，“第11.5.1“液腔的定义，“第11.5.2“流体交换的定义，“部分11.5.3“打气筒的定义，”部分11.5.4surface-based流体腔：概述概述基于表面的充液腔是由：用标准有限元法对充液结构进行建模；利用表面定义来提供充液结构的变形和流体在结构腔边界上施加的压力之间的耦合；定义流体行为；使用流体交换定义来模拟流体在腔和环境之间或多个腔之间的流动；使用充气器定义将气体混合物注入流体腔以模拟汽车安全气囊的膨胀。基于表面的流体腔能力可

2、以用来模拟液体或充气结构。它取代了基于元件的静液腔能力的功能，不需要用户定义流体或流体链接元素。介绍在某些应用中，可能需要预测充液或A的机械响应。充气结构。例如压力容器，液压或气动驱动机构，汽车安全气囊。解决此类应用程序的一个主要困难是结构变形与含流体对结构的压力。图1说明了一个简单的-1151对系统外部进行流体填充结构的例子*欧阳光明*创编2021.03.07*欧阳光明*创编2021.03.07荷载.结构的响应不仅取决于外部载荷，还取决于压力。由流体施加的，而流体又受结构变形的影响。基于表面流体腔能力提供了分析腔体可以使用的情况所需的耦合。假设完全由具有相同性质和状态的流体填充。应用与意义空

3、腔内的空间变化不能用这个特性建模。例如，考虑流体-结构相互作用和耦合EulerianLagrangian能力的应用涉及晃动和波传播通过流体（见“欧拉分析，“第14.1.1条；“流体结构互动”中的“协同仿真：概述，“第17.1.1；和“流体结构的协同仿真和共轭传热，”部分17.3.2）。离散液腔流体腔的边界由一个基于单元的表面定义，表面法线指向腔内。底层元素可以是标准的固体或结构元素，也可以是表面元素。表面元素可以用来模拟结构中的空洞或填充刚性区域。在刚性或其他承载元素不存在（见“表面元素，”部分32.7.1）。护理当使用表面元素，使节点完全被唯一的表面包围时，必须使用该元素。元素有适当的边界条

4、件。considertheexamplepresentedinfigurell.5.1-1。solidelementsaredefinedonthetopandsideof图2表明，腔-11.5.1。在底部刚性边界上定义了一个曲面单元。没有标准元素存在的空洞。节点位于对称轴的交点处。腔的下刚性边界必须被约束在r和z方向是因为它仅与表面元素相连接。定义空腔的表面是基于底层固体和表面元素。在Abaqus中，可以将额外的用户定义卷添加到实际或几何体中。空腔体积。如果腔的边界不是由一个基于单元的表面定义的，流体腔被假定有一个固定的体积等于附加的体积。定义空腔参考节点的位置一个被称为空腔参考节点的单一节

5、点与流体腔相关联。这一腔的参考nodehasasingledegreeoffreedomrepresentingthepressureinsidethefluidcavity。thecavityreference节点也用于计算腔体体积。如果空腔不受对称平面的约束，则定义空腔的表面必须完全。封闭空腔以确保其容积的正确计算。在这种情况下，空腔的位置参考节点是任意的，不必在空腔内。如果由于对称性，只有空腔边界的一部分用标准元素建模，腔参考节点必须位于对称平面或轴线（图1151-2）。如果多个对称平面存在，腔参考节点必须位于对称的交点上。飞机（图11.5.1-3）。对于轴对称分析，空腔参考节点必须位于

6、对称轴。这些要求是流体腔未完全封闭的结果。定义空腔的表面。有限元计算采用体积元法对基于表面的空腔进行有限元计算。在“流体静压计算描述，“ABAQUS理论指导第3.8.1。体积腔的元素是由Abaqus利用表面小面几何和你定义的空洞参考节点。在Abaqus/标准中，表面刻面用以下元素类型：fax2和F2D2（这是线性，两级，轴对称和平面元素，分别）和f3d3和f3d4（这是线性的，分别为3-4节点三维元素，）。Abaqus的二阶刻面进一步细分为多个线性面或元素.液腔的行为流体填充腔中的流体的行为可以基于液压或气动模型。水力模型可以模拟几乎不可压缩的流体行为。Abaqus/标准中的不可压缩行为。通过

7、定义一个体积来引入压缩性。模量。气动模型是以理想气体为基础的。气体可以由多个物种来定义。在ABAQUS/显式中，你可以指定气体的温度，或者根据绝热行为假设。具有绝热温度更新的多物种理想气体是汽车安全气囊的合适模型。模腔流入或流出的模型ABAQUS中有许多方法来模拟流体进出腔。流量可以是specifiedasaprescribedmassorvolumefluxhistoryorcanmodelphysicalmechanismsduetoapressuredifferentialsuchasventingthroughanexhaustorificeorleakagethroughaporou

8、sfabric。fluidexchange定义用于此目的，可以模拟流体腔与其周围环境之间的流动。两液腔之间（见“流体交换的定义，“部分11.5.3，详情）。此外,ABAQUS/Explicit有能力模型充气机用于汽车安全气囊的展开。可以直接指定充气器上的条件，也可以使用罐试验数据（见“充气机定义”）。第11.5.4，详情）。模拟多室许多充液系统，如气囊，有多个腔室，流体在腔室之间通过孔或织物泄漏流动。在其他情况下，将一个单独的物理室分割成多个腔室，并用假想的壁来模拟整个物理室的压力梯度是有利的。一些虚拟泄漏机制通过室间墙可以被定义以获得合理的行为。当模拟气囊的复杂展开时，这可能是一种有用的建模

9、技术。为了模拟多个腔室，为每个腔室定义一个流体腔，并将流体腔与适当的流体交换定义联系起来。平均性能的多室模型可以输出如果要求（见“液腔的定义，“第11.5.2，详情）。在动态过程中定义流体惯性流体腔内的流体惯性或空穴之间流体交换的惯性不会自动考虑在内。为了增加惯性效应，在腔体边界上使用质量单元。你应该确保总的附加质量相当于腔中流体的质量，质量元素的分布是结构承受荷载的分布流体质量的合理表示。只有流体惯性的整体效应可以被建模；腔中的均匀压力假设使得不可能模拟任何压力梯度驱动的流体运动。因此，该方法假定激励的时间尺度与流体的典型响应时间相比非常长。涉及空腔边界的接触模拟如果从空腔中取出大量流体或围

10、绕空腔的材料是非常灵活的，空腔可以部分塌陷，并且空腔壁的部分可以彼此接触。通过使用Abaqus现有的标准技术来模拟接触，可以有效地处理空腔壁的自接触和与周围结构的接触。ABAQUS/Explicit也能说明堵塞流出由于接触表面的空腔（见“占堵塞由于接触界面”中的“体液交换的定义，“部分11.5.3）。在某些应用程序中，在求解过程中会遇到负本征值。这些负的本征值并不一定意味着超过了分叉或屈曲载荷。如果预测响应似乎是合理的，这些消息可以忽略。详细描述如何负特征值可以开发静压流体元问题的解决在“流体静压的计算方法中，“ABAQUS理论指导第3.8.1。程序基于表面的流体腔的能力可以用在除了耦合孔隙流

11、体扩散/应力分析程序以外的任何程序（见“耦合孔隙流体扩散和应力分析，“部分6.8.1）。初始条件初始流体压力和温度可指定（见“初始条件在Abaqus/标准和ABAQUS/Explicit”部分34.2.1）o对于理想气体，初始压力表示高于或超过环境压力的压力表。初始温度应在所使用的温标中给出。绝对零度的温度范围，对于理想气体另行规定（见“液腔的定义，“第11.5.2）o如果膜元件作为流体腔的基本元素，参考网格（IMM）也可以指定（见“初始条件在Abaqus/标准和ABAQUS/Explicit”部分34.2.1）o边界条件空腔参考节点的自由度（自由度8）是问题的主要变量。因此，可以通过定义边界

12、条件的规定（见“边界条件在Abaqus/标准和ABAQUS/Explicit”部分34.3.1），其方式类似于结构节点位移可以规定。处方的压力腔参考节点相当于施加均匀压力腔边界采用分布式负载定义（见“分布载荷，”部分34.4.3）o如果压力是用边界条件规定的，则流体体积自动调节以填充空腔（即，流体被假定为进入并离开所需的腔以维持规定的压力）。这种行为对于在引入流体作用之前变形的空腔是有用的。在随后的步骤中您可以删除在自由的压力程度的边界条件（参见拆卸边界条件”在“边界条件在Abaqus/CFD,”部分34.3.2），因此“封”与当前的流体体积的腔。荷载分布式压力和身体的力量，以及集中的节点力，

13、可应用于充满液体的结构，如在集中荷载描述”部分34.4.2，和“分布载荷，”部分34.4.3o预定义场预定义的温度场和用户自定义场变量可以为充满液体的结构和封闭的流体的定义，如“预定义场描述”部分34.6.1o温度流体温度可在所有的指定参考节点腔预定义的字段（见“预定义温度”在“预定义字段，”部分34.6.1），除非一个绝热过程是指定或耦合温度位移程序是用。应用和初始温度之间的任何差异都会引起气动流体的热膨胀，如果给定热膨胀系数，则会对液压流体产生热膨胀。一个特定的温度场也可以影响温度相关的材料性质，如果有的话，对于流体填充的结构和封闭流体。场变量用户定义的字段变量的值可以在所有腔参考节点指定

14、（见“预定义场变量”中的“预先定义的领域，”部分34.6.1）。这些值将影响封闭流体的场变量相关的材料性质。输出腔内流体状态可利用节点输出变量度和cvol历史输出，代表测量流体的压力和容积，分别。在稳态的动态程序的大小和流体压力的相位角可以得到节点变差。ABAQUS/Explicit还提供了把腔温度、腔的表面积，和流体质量（节点输出变量ctemp,csarea,和cmass，分别）。输出变量ctemp仅当理想气体模型是在绝热条件下使用。如果节点设置为输出的要求是包含一个以上的液腔,历史的平均流体压力、总体积、流体平均温度,和所有的外腔的表面区域,和这些空腔的总质量也提出了利用节点输出变量apc

15、av,tcvol,actemp，tcsarea,和tcmass，分别。在Abaqus/明确,当模型包括流体交换的定义,利用节点输出变量CMFL和cmflt获得的总流量和累计流量的历史输出一腔创面肉芽和ceflt获得的总热流量和累计流量历史热能输出腔。如果为一个空腔定义了不止一种流体交换，那么质量或热能流量和累积的质量或热能流的时间历程就会从每个流体交换腔也将被输出。如果流体腔是由一个混合理想气体模型，分子的质量分数各流体的物种内流体腔时程可以利用节点输出变量CMF获得。如果使用充气器，利用节点输出变量minfl，minflt，和tinfl获得质量流量的历史，积累了大量的流量，和每一个打气筒充气

16、温度定义（见“ABAQUS/Explicit输出变量标识符，“422节）。输入文件模板流体静力学分析：*HEADINGFLUIDCAVITY,NAME=cavity_name,BEHAVIOR=behavior_name,REFNODE=cavity_reference_node,SURFACE=surface_nameFLUIDBEHAVIOR,NAME=behavior_nameFLUIDDENSITYDatalinetodefinedensityFLUIDBULKMODULUSDatalinetodefinebulkmodulusFLUIDEXPANSIONDatalinetodefin

17、ethermalexpansion*FLUIDEXCHANGE,NAME=exchange_name,PROPERTY=exchange_property_namecavity_reference_node*FLUIDEXCHANGEPROPERTY,NAME=exchange_property_name,TYPE=MASSFLUXDatalinetodefinemassflowrateperunitarea*INITIALCONDITIONS,TYPE=TEMPERATUREDatalinetodefineinitialtemperature*INITIALCONDITIONS,TYPE=F

18、LUIDPRESSUREDatalinetodefineinitialpressure*STEP*TEMPERATUREDatalinetodefinetemperature*FLUIDEXCHANGEACTIVATIONexchange_name*ENDSTEP*标题液腔，名称=cavity_name，行为=behavior_name,参考节点=surface_namecavity_reference_node，表面*流体的行为，behavior_name名称=*流体密度定义密度的数据线流体体积模量定义体积模量的数据线*流体膨胀定义热膨胀的数据线*体液交换，名字=exchange_name，

19、产权=exchange_property_namecavity_reference_node*流体交换性能、名称=exchange_property_name，类型=质量通量确定单位面积质量流量的数据线*初始条件，类型=温度确定初始温度的数据线*初始条件，类型=流体压力确定初始压力的数据线*步*温度定义温度的数据线*流体交换活化exchange_name*结束步骤用理想气体混合的气囊分析：*HEADING*FLUIDCAVITY,NAME=chamber_1,MIXTURE=MOLARFRACTION,ADIABATIC,REFNODE=chamber_1_reference_node,SUR

20、FACE=surface_name_1blanklineOxygen,0.2Nitrogen,0.75Carbon_dioxide,0.05*FLUIDCAVITY,NAME=chamber_2,BEHAVIOR=Air,ADIABATIC,REFNODE=chamber_2_reference_node,SURFACE=surface_name_2blankline*FLUIDBEHAVIOR,NAME=Air*CAPACITY,TYPE=POLYNOMIALDatalinetodefineheatcapacitycoefficient*MOLECULARWEIGHTDatalinetode

21、finemolecularweight*FLUIDBEHAVIOR,NAME=Oxygen*CAPACITY,TYPE=POLYNOMIALDatalinetodefineheatcapacitycoefficient*MOLECULARWEIGHTDatalinetodefinemolecularweight*FLUIDBEHAVIOR,NAME=Nitrogen*CAPACITY,TYPE=POLYNOMIALDatalinetodefineheatcapacitycoefficient*MOLECULARWEIGHTDatalinetodefinemolecularweight*FLUI

22、DBEHAVIOR,NAME=Carbon_dioxide*CAPACITY,TYPE=POLYNOMIALDatalinetodefineheatcapacitycoefficient*MOLECULARWEIGHTDatalinetodefinemolecularweight*FLUIDINFLATOR,NAME=inflator,PROPERTY=inflator_propertychamber_1_reference_node*FLUIDINFLATORPROPERTY,NAME=inflator_property,TYPE=TEMPERATUREANDMASSDatalinestod

23、efinemassflowrateandgastemperature*FLUIDINFLATORMIXTURE,TYPE=MOLARFRACTION,NUMBERSPECIES=2Carbon_dioxide,NitrogenTabletodefinemolecularmassfraction*FLUIDEXCHANGE,NAME=exhaust,PROPERTY=exhaust_behaviorchamber_1_reference_node*FLUIDEXCHANGEPROPERTY,NAME=exhaust_behavior,TYPE=ORIFICEDatalinetospecifyor

24、ificebehavior*FLUIDEXCHANGE,NAME=leakage_1,PROPERTY=fabric_behaviorchamber_1_reference_node*FLUIDEXCHANGE,NAME=leakage_2,PROPERTY=fabric_behaviorchamber_2_reference_node*FLUIDEXCHANGEPROPERTY,NAME=fabric_behavior,TYPE=FABRICLEAKAGEDatalinetospecifyfabricleakagebehavior*FLUIDEXCHANGE,NAME=chamber_wal

25、l,PROPERTY=wall_behavior,EFFECTIVEAREA=chamber_1_reference_node,chamber_2_reference_node*FLUIDEXCHANGEPROPERTY,NAME=wall_behavior,TYPE=ORIFICEDatalinetospecifyorificebehavior*AMPLITUDE,NAME=amplitude_nameDatalinetodefineamplitudevariations*PHYSICALCONSTANTS,UNIVERSALGASCONSTANT=*INITIALCONDITIONS,TY

26、PE=FLUIDPRESSUREDatalinetodefineinitialpressure*INITIALCONDITIONS,TYPE=TEMPERATUREDatalinetodefineinitialtemperature*STEP*TIME*FLUIDEXCHANGEACTIVATIONexhaust,leakage_1,leakage_2,chamber_wall*FLUIDINFLATORACTIVATION,INFLATIONAMPLITUDE=amplitude_nameinflator*ENDSTEP*标题*液腔，名称=chamber，混合=摩尔分数，绝热，参考节点=su

27、rface_name_lchamber_l_reference_node，表面空白行氧气，0.2氮，0.75carbon_dioxide，0.05*液腔，名称=chamber_2，行为=空气绝热，参考节点=surface_name_2chamber_2_reference_node，表面空白行*流体行为，名称=空气*容量，类型=多项式定义热容量系数的数据线*分子量确定分子量的数据线*流体行为，名称=氧气*容量，类型=多项式定义热容量系数的数据线*分子量确定分子量的数据线*流体行为，名称=氮气*容量，类型=多项式定义热容量系数的数据线*分子量确定分子量的数据线*流体的行为，carbon_diox

28、ide名称=*容量，类型=多项式定义热容量系数的数据线*分子量确定分子量的数据线*液充气，充气inflator_property属性名称=,=chamber_1_reference_node*流体的充气性能，名称=inflator_property，类型=温度和质量确定质量流量和气体温度的数据线*流体充气混合物，类型=摩尔分数，种数=2carbon_dioxide,氮分子质量分数表*体液交换，名称=排气，产权=exhaust_behaviorchamber_1_reference_node*流体交换性能、名称=exhaust_behavior,类型=指定管口行为的数据线*体液交换，名字=lea

29、kage，产权=fabric_behaviorchamber_1_reference_node*体液交换，名字=leakage_2，产权=fabric_behaviorchamber_2_reference_node*流体交换性能、名称=fabric_behavior，类型=织物泄漏指定织物泄漏行为的数据线*体液交换，名字=chamber_wall，产权=wall_behavior，有效面积=chamber_1_reference_node，chamber_2_reference_node*流体交换性能、名称=wall_behavior，类型=指定管行为的数据线*振幅，名称=amplitude

30、_name定义振幅变化的数据线*物理常数，通用气体常数=*初始条件，类型=流体压力确定初始压力的数据线*初始条件，类型=温度确定初始温度的数据线*步*流体交换活化排气，leakage，leakage_2,chamber_wall*流体充激活，通胀幅度=amplitude_name打气筒*结束步骤11.5.2液腔的定义产品：ABAQUS/标准ABAQUS/显式ABAQUS/CAE参考资料“基于流体腔表面：概述，“第11.5.1“流体交换的定义，“部分11.5.3*能力*流体行为流体体积模量*流体腔流体密度分子量“定义流体腔相互作用，对ABAQUS软件的用户指南部分15.13.11，本指南中的HT

31、ML版本“定义流体腔相互作用特性，对ABAQUS软件的用户指南部分15.14.4，本指南中的HTML版本概述一种基于表面的流体腔：可以用来模拟充液或充气的结构；与一个称为空腔参考节点的节点相关联；通过指定一个完全包围空腔的表面来定义；仅适用于在任何时间点内流体在特定空腔内的压力和温度均匀的情况；可以用理想气体混合物在绝热条件下的假设来模拟安全气囊；有一个名称，可以用来标识与该腔相关联的历史输出。1定义流体腔必须将名称与每个流体腔相关联。使用输入文件：*FLUIDCAVITY,NAME=name使用ABAQUS/CAE：相互作用模块：CreateInteraction:Fluidcavity,N

32、ame:name指定空腔参考节点每个流体腔都必须有一个相关的腔参考节点。与流体腔名称一起，参考节点用于识别流体腔。此外，还可以通过流体交换和充气器定义来参考。参考节点不应连接到模型中的任何元素。InputFileUsage:FLUIDCAVITY,REFNODE=nAbaqus/CAEUsage:Interactionmodule:CreateInteraction:Fluidcavity:selectthefluidcavityreferencenode1.2指定流体腔的边界流体腔必须完全封闭的有限元建模除非对称面（见“基于表面的流体腔：概述，“第11.5.1）。表面元素可用于非结构的空腔表

33、面的部分。空腔的边界是用一个基于单元的表面来描述的，该表面覆盖包围空腔的元素，表面法线指向内部。默认情况下，如果表面的底层元素没有一致的法线，就会发出错误消息。或者，也可以跳过对表面法线的一致性检查。InputFileUsage:使用下面的选项定义具有一致正常检查的表面：FLUIDCAVITY,SURFACE=surface_name,CHECKNORMALS=YES使用下面的选项定义表面没有一致的法线检查:*FLUIDCAVITY,SURFACE=surface_name,CHECKNORMALS=NOAbaqus/CAEUsage:Interactionmodule:CreateInter

34、action:Fluidcavity:选择流体腔边界面；切换或关闭检查表面法线。1.3在流体腔中指定附加容积可以为Abaqus/显式流体腔指定附加容积。当空腔的边界由指定的表面定义时，附加体积将被加到实际的体积中。如果没有指定形成流体腔边界的表面，则流体腔被假定为具有与所添加体积相等的固定体积。InputFileUsage:*FLUIDCAVITY,ADDEDVOLUME=rAbaqus/CAEUsage:在AbaqusCAE中不支持附加体积规范.1.4指定最小体积当流体腔的体积非常小时，显式动态过程中的瞬变会使体积变为零甚至负，使有效空腔刚度值趋于无穷大。为了避免数值问题，您可以在Abaqu

35、s中指定流体的最小体积。如果腔的体积（等于实际体积加上体积）低于最小值，则用最小值来计算流体压力。您可以直接指定最小体积，也可以指定流体腔的初始体积。如果使用后一种方法，流体腔的初始体积为负值，则最小体积等于t=l=零。输入文件用法：使用下面的选项直接指定最小体积：*FLUIDCAVITY,MINIMUMVOLUME=minimumvolume使用以下选项指定最小音量等于初始体积：*FLUIDCAVITY,MINIMUMVOLUME=INITIALVOLUMEAbaqus/CAEUsage:最小体积的规范不支持。2流体腔行为的定义流体腔的行为支配腔压力、体积和温度之间的关系。Abaqus/标准

36、中的流体腔只能包含一种流体。在Abaqus/Explicit，一个空洞可以包含一种液体或理想气体的混合物。均质流体的流体行为为了定义由单流体构成的流体腔行为，指定一种流体行为来定义流体性质。必须将流体行为与名称关联起来。然后，这个名称可以用来将某种行为与流体腔定义联系起来。输入文件的使用：使用以下选项：*FLUIDCAVITY,NAME=fluid_cavity_name,BEHAVIOR=behavior_name*FLUIDBEHAVIOR,NAME=behavior_nameAbaqus/CAEUsage:Interactionmodule:CreateInteractionProperty:Fluidcavity,Name:behavior_nameAbaqus/Explicit理想气体混合的流体行为在Abaqus/Explicit中，你可以定义一个由多种气体组成的流体腔行为。为了定义由多种气体组成的流体腔行为，你必须指定多种流体行为来定义流体性质。指定流体行为的名称和定义混合物的初始质量或摩尔分数，以便将某一组行为与流体腔定义