好文档 - 专业文书写作范文服务资料分享网站

abaqus建模流程 - 学习笔记

天下 分享 时间: 加入收藏 我要投稿 点赞

行单元划分。当然,可以采用四面体单元(tetrahedral)利用自由网格技术对任何形状的模型区域进行单元划分。 ? ? ? ?

自由划分:可以应用到任意平面和曲面

结构划分:结构化的网格划分通常给出了对网格的最大的控制 扫掠划分:在源面中,相邻面之间的二面角不能和180°相差太远。

虚拟拓扑:在某些情况下,装配件的部件实例可能包含一些小的细节,比如表面和边。虚拟拓扑可以忽略这些不重要的细节。

如果部件实例中包含虚拟拓扑,那么它只能使用以下单元通过自由网格技术划分网格: ?自由网格:三角形和四面体单元、用波前法划分的四边形为四边形为主单元网格 ?扫略网格:六面体或者楔形单元

?映射网格:四边形, 三角形, 或者六面体单元 (2) 分解模型(partition)

可以对模型中的边(edge)、面(face)和体(cell)进行分解。用来将边、面、体分解成更小部分的点、边、面都成为模型中的特征体,这些特征体和其他特征体一样可以在特征体管理器中查看。(如:将一个体分解成两部分需要用一个面将体切割成两部分,这个面就成了模型中一个新的特征体。)。

有五种分割特征体可以将一个特征体分解:定义切割面(define cutting plane)、使用基准平面(use datum plane)、延伸平面(extended face)、挤压或旋转边(extrude/sweep edges)、N-sided patch。

一次分解操作仅仅只是将被分解的对象分解成两部分,并不能改变被分解对象所在特征体(部件实例)的整体性。不能对刚体进行单元划分。 (3)设置网格种子 选择Seed part或SeedEdge ? ?

设置全局种子 设置边上的种子 边上的种子无约束圆圈 边上的种子受部分约束三角形 边上的种子受完全约束方形

对于可以用扫掠方法进行网格划分的区域,边的网格种子从选定的边到匹配边自动传播,由分区创建的新边自动继承总体网格种子。 (4)单元划分控制

在 Assign mesh control中指定单元划分方法(结构划分法(structured)、自由划分法(free)、扫掠法(sweep)等等)。操作的对象是被分解后的边、面、和体,可以对同一实例(装配模型的特征体)分解后产生的不同边、面、体分别采用不同的单元划分方法,指定不同的单元类型。

在Assign element type中指定单元类型(六面体单元、四面体单元等等),选择单元库(standard、explicit)、确定线性单元(linear)或者二次单元(quadratic)、确定这两种单元的特性:杂交元(hybrid formulation)、缩减积分(reduced integration)、非协调单元模式(incompatible modes)。 1)

单元形状 二维问题

a) Quad:网格中完全使用四边形单元。

b) Quad-dominated:网格中主要使用四边形单元,但在过渡区域允许出现三角形单元,更容易实现从粗网

格向细网格过渡。

c) Tri: 网格中完全使用三角形单元。 三维问题

a) Hex:网格中完全使用六面体单元。

b) Hex-dominated:网格中主要使用六面体单元,但在过渡区域允许出现锲形(三棱柱)单元。

c) Tet:网格中完全使用四面体单元。 d) Wedge:网格中完全使用楔形单元。

Quad单元(二维)和Hex单元(三维)可以用较小的计算代价得到较高的精度。自由网格采用Tri单元(二维)和Tet单元(三维),一般应选择带内部节点的二次单元来保证精度。结构化网格和扫掠网格一般采用quad单元(二维单元)和hex单元(三维),分析精度较高,因此在划分单元时选择后两种。 2)

三维实体单元类型 ? ? ?

线性(linear)单元又称一阶单元,仅在单元的角点处布置节点,在各方向都采用线性插值; 二次(quadratic)单元又称二阶单元,在每条边上有中间节点,采用二次插值;

修正的(modified)二次单元只有Tri或Tet单元才有这种类型,即在每条边上有中间节点,并采用修正的二次插值

b) 连续体单元

ABAQUS/Standard 的连续体单元库包括二维和三维的线性单元和二次单元,分别可以采用完全积分或减缩积分,另外还有修正的二次Tri 和Tet单元,以及非协调模式单元和杂交单元。

ABAQUS/Explicit 的连续体单元库包括二维和三维的线性减缩积分单元,以及修正的二次Tri和Tet单元。ABAQUS/Explicit 中没有二次完全积分的连续体单元。 c) 线性完全积分( linear full-integration ) 单元

保持默认的Linear参数,取消对Reduced integration(减缩积分)的选择,例如CPS4单元(4节点四边形双线形平面应力完全积分单元)和C3D8单元(8节点六面体线性完全积分单元)。所谓“完全积分”是指当单元具有规则形状时,所用的离斯积分点的数目足以对单元刚度矩阵中的多项式进行精确积分。承受弯曲荷载时,线性完全积分单元会出现剪切自锁(shear locking)问题,造成单元过于刚硬,即使划分很细的网格,计算精度仍然很差。

d) 二次完全积分(quadratic-full-integration ) 单元

在Element Type对话框选择 Quadratic 参数,取消对Reduced integration的选择,例如CPS8 (8节点四边形二次平面应力完全积分单元)和C3D20 (20节点六面体二次完全积分单元)。

二次完全积分单元的优点如下。

? 对应力的计算结果很精确,适于模拟应力集中问题; ? —般情况下没有剪切自锁。

但使用这种单元时要注意以下问题 ? 不能用于接触分析;

? 对于弹塑性分析,如果材料是不可压缩性的(例如金属材料),则容易产生体积自锁 (volumetric locking); ? 当单元发生扭曲或弯曲应力有梯度时,有可能出现某种程度的自锁。 e) 线性减縮积分(linear reduced-integration) 单元

对于Quad单元和Hex单元ABAQUS/CAE 默认的单元类型是线性减缩积分单元,例如CPS4R (4节点四边形双线形平面应力减缩积分单元)和C3D8R(8节点六面体线性减缩积分单元)。减缩积分单元比通常的完全积分单元在每个方向少用一个积分点。线性减缩积分单元只在单元的中心有一个积分点,由于存在所谓“沙漏”(hourglass) 数值问题而过于柔软。ABAQUS在线性减缩积分单元中引人了“沙漏刚度”以限制沙漏模式的扩展。模型中的单元越多,这种刚度对沙漏模式的限制越有效。可以选择沙漏控制参数为Enhanced , RelaxStiffness ,Stiffness ,Viscous或combined。采用线性减缩积分单元模拟承受弯曲荷载的结构时,沿厚度方向上至少应划分四个单元。

线性缩减积分单元有以下优点: ? 对位移的求解结果较精确。

? 网格存在扭曲变形时(例如 Quad单元的角度远远大于或小于 90°),分析精度不会受到大的影响。 ? 在弯曲荷载下不容易发生剪切自锁。 a) 节点数目和插值阶数

其缺点如下:

? 需要划分较细的网格来克服沙漏问题。

? 如果希望以应力集中部位的节点应力作为分析指标,则尽量不要使用线性减缩积分单元,而应使用二次单元,因为线性缩减积分单元只在单元的中心有一个积分点,相当于常应力单元,它在积分点上的应力结果是相对精确的,而经过外插值和平均后得到的节点应力则不精确。如果在应力集中部位进行了网格细化,使用二次减缩积分单元与二次完全积分单元得到的应力结果相差不大,而二次减缩积分单元的计算时间相对较短。

Note:在查看模型的应力结果时有两种选择:

? 查看节点上的应力:这是通常用的方法,其优点是简单方便,但事实上,后处理中得到的节点应力是对单元积分点上的应力进行外插值和平均后得到的并不精确;

? 查看单元积分点上的应力:这是ABAQUS所推荐的方法。线性缩减积分单元只有一个积分点,可以很方便地查看积分点上的分析结果,但其他类型的单元有多个积分点,就需要详细了解节点的编号顺序,并根据模型的实际情况来决定查看哪个积分点,这一过程很烦琐。另外要注意,单元积分点上的应力值往往不是应力集中区城的最大应力。

(用户可以自己在上述两种方法中做出选择,需要注意的是,如果希望查看节点上的应力,就尽量不要使用线性缩减积分单元;如果使用了线性缩减积分单元,就应该查看单元积分点上的分析结果,并且要在应力变化剧烈的部位划分足够细的网格。如果外推应力值与积分点应力值差别很大,说明单元间应力变化剧烈,单元网格过于粗糙,计算的应力不够精确。这种外推应力误差会因为单元网格细化而减小,但总是存在。所以在使用单元变量的节点值时要谨慎。默认的应力不变量的计算方法是computer scalars before averaging(先不变量再平均),得到的节点应力偏大,作为工程分析的结果会更安全;也可以在result options中设置为computer scalars after averaging(先平均再不变量),得到的节点应力偏小。事实上只有单元积分点上的应力结果是相对精确的,一般在所关心的高应力部位细化网格,使用默认设置即可。更广义地说,有限元作为一种数值计算方法,其本身就是不精确的。 f) 二次缩减积分(quadratic-reduced-integration)单元

对于Quad单元或Hex单元,可以在 Element Type 对话框中将单元类型设置为二次减缩积分单元,如 CPS8R (8节点四边形二次平面应力减缩积分单元)和C3D20R(20节点六面体二次缩减积分中元),这种单元不但保待了前面介绍的线性缩减积分单元的优点,而且还具有以下特性: ? 即使不划分很细的网格也不会出现严重的沙漏问题 ? 即使在复杂应力状态下,对自锁问题也不敏感 但使用这种单元时需要注意以下问题: ? 不能在接触分析中使用 ? 不适于大应变问题

? 存在与线性减缩积分单元相类似的问题,由于积分点少,得到的节点应力的精度往往低于二次完全积分单元。

g) 非协调模式( incompatible modes) 单元

对于Quad单元或Hex单元,可以在Element Type 对话框中将单元类型设为非协调模式单元。例如CPS4(4节点四边形双线形平曲应力非协调单元)和C3D8I(8节点六面体线性非协调模式单元)。仅在 ABAQUS/Standard中有非协调模式单元。其目的是克服在线性完全积分单元中的剪切自锁问题。

非协调摸式单元的优点如下:

? 克服了剪切自锁问题,在单元扭曲比较小的情况下,得到的位移和应力结果很精确;

? 在弯曲问题中,在厚度方向上只需很少的单元,就可以得到与二次单元相当的结果,而计算成本却明显降低;

? 使用了增强变形梯度的非协调模式,单元交界处不会重叠和开洞,因此很容易扩展到非线性、有限应变的位移。

但使用这种单元时需注意,如果所关心部位的单元扭曲比较大,尤其是出现交错扭曲时,分析精度会降低。

h) Tri单元和Tet单元

使用Tri单元和Tet单元时应注意以下问题:

? 线性Tri单元和Tet单元的精度很差,所以不要在模型中所关心的部位及其附近区域使用。

? 二次 Tri 单元和 Trt 单元精度较高,而且能模拟任意的几何形状,但计算代价比Quad单元或Hex单元大,因此如果模型中能够使用 Quad单元或Hex单元,尽量不要使用Tri 单元或Trt单元。 ? 二次Trt单元(C3D10) 适于ABAQUS/standard中的小位移无接触问题;修正的二次Tet单元(C3D10M) 适于ABAQUS/Explicit,以及 ABAQUS/standard中的大变形和接触问题。 ? 使用自由网格不易通过布置种子来控制实体内部的单元大小。 i) 杂交(hybrid) 单元

在ABAQUS/standard中,每一种实体单元(包括所有缩减积分和非协调模式单元)都有其相应的杂交单元,用于不可压縮材料(泊松比为0.5) 或近似不可压缩材料(泊松比大于0.475),橡胶就是一种典型的不可压缩材料。除了平面应力问题之外,不能用普通单元来模拟不可压缩材料的响应,因为此时单元中的压应力是不确定的。杂交单元在它的名字中含有字母H。ABAQUS/Explicit 中没有杂交单元。 j) 混合使用不同类型的单元

当三维实体几何形状较复杂时,无法在整个实体上使用结构化网格或扫掠网格划分技术,得到Hex单元网格,这时一种常用的做法是对于实体不重要的部分使用自由网格划分技术,生成Tet单元网格,而对于所关心的部分采用结构化网格或扫掠网格生成 Hex 单元网格,在生成这样的网格时,ABAQUS会给出提示将生成非协调的网格,在不同单元类型的交界处将自动创建绑定(tie)约束。

需要注意的是,在不同单元类型网格的交界处.,即使单元角部节点是重合的,仍然有可能出现不连续的应力场,而且在交界处的应力可能大幅度地增大。如果在实体中混合使用线性和二次单元,也会出现类似的问题。因此在混合使用不同类型的单元时,应确保其交界处远离所关心的区域,并仔细检查分析结果 是否正确。

对于无法完全采用Hex单元网格的实体,还可以使用以下方法:

? 对整个实体划分Tet 单元网格,使用二次单元C3D10或修正的二次单元 C3D10M,同样可以达到所需的精度,只是计算时间较长。

? 改变实体中不重要部位的几何形状,然后对整个实体采用Hex单元网格。 k)

选择三维实体单元类型的基本原则

计算代价,提高计算精度。当几何形状复杂时,也可以在不重要的区域使用少量楔形(Wedge) 单元。 ? 如果使用了自由网格划分技术,Tet单元的类型应选择二次单元。在ABAQUS/Explicit中应选择修正的 Tet单元C3D10M,在,ABAQUS/Standard中可以选择 C3D10,但如果有大的塑性变形,或模型中存在接触,而且使用的是默认的“硬”接触关系,则也应选择修正的Tet单元 C3D10M。

? ABAQUS的所有单元均可用于动态分析,选取单元的一般原则与静力分析相同。但在使用 ABAQUS/Explicit 模拟冲击或爆炸荷载时,应选用线性单元,因为它们具有集中质量公式,模拟应力波的效果优于二次单元所采用的一致质量公式。

如果使用的求解器是 ABAQUS/Standard, 在选择单元类型时还应注意以下方面:

? 对于应力集中问题,尽量不要使用线性减缩积分单元,可使用二次单元来提高精度。如果在应力集中部位进行了网格细化,使用二次减缩积分单元与二次完全积分单元得到的应力结果相差不大,而二次减缩积分单元的计算时间相对较短。

? 对于弹塑性分析,如果材料是不可压缩性的(例如金属材料),则不能使用二次完全积分单元(C3D20),否则会出现体积自锁问题,也不要使用二次Tri 单元或Tet单元。推荐使用的是修正的二次Tri单元或Tet单元(C3D10M)、非协调单元(C3D8I),以及线性减缩积分单元(C3D8R)。如果使用二次减缩积分? 对于三维区域,尽可能采用结构化网格划分技术或扫掠网格划分技术,从而得到Hex 单元网格,减小

单元,当应变超过20%-40%时要划分足够密的网格。

? 塑性材料和接触面上都不能用 C3D20R 和 C3D20 单元,这可能是你收敛问题的主要原因。如果需要得到应力,可以使用 C3D8I (在所关心的部位要让单元角度尽量近 90 度),如果只关心应变和位移,可以使用 C3D8R, 几何形状复杂时,可以使用 C3D10M。

? 如果模型中存在接触或大的扭曲变形,则应使用线性Quad或Hex单元,以及修正的二次Tri单元或Tet 单元,而不能使用其他的二次单元。

? 对于以弯曲为主的问题,如果能够保证在所关心部位的单元扭曲较小,使用非协调单元(例如C3D8I单元)可以得到非常精确的结果。

? 除了平面应力问题之外,如果材料是完全不可压缩的(例如橡胶材料),则应使用杂交单元;在某些情况下,对于近似不可压缩材料也应使用杂交单元。 3)

选择壳单元的类型

如果一个薄壁构件的厚度远小于其典型整体结构尺寸(一般为小于 1/10),并且可以忽略厚度方向的应力,就可以用壳单元来模拟此结构。壳体问题可以分为两类:薄壳问题(忽略横向剪切变形)和厚壳问题(考虑横向剪切变形)。对于单一各向同性材料,一般当厚度和跨度的比值小于1/15 时,可以认为是薄壳;大于 1/15 时,则可以认为是厚壳。对于复合材料,这个比值需要更小一些。 ABAQUS 的壳单元可以有多种分类方法,按照薄壳和厚壳可划分为: ? 通用目的(general-purpose) 壳单元:此类单元对薄壳和厚壳问题均有效。

? 特殊用途(special-purpose) 壳单元:包括纯薄壳 (thin-only) 单元和纯厚壳 ( thick-only) 单元。 根据单元的定义方式,还可以将 ABAQUS 壳单元划分为:

? 常规(conventional) 壳单元:通过定义单元的平面尺寸、表面法向和初始曲率来对参考面进行离散,只能在截面厲性中定义壳的厚度,而不能通过节点来定义壳的厚度。

? 连续体(continuum) 壳单元:类似于三维实体单元,对整个三维结构进行离散。 选择壳单元的类型时可以遵循以下原则:

? 对于薄壳问题,常规壳单元的性能优于连续体壳单元;而对于接触问题,连续体壳单元的计算结果更加精确,因为它能在双面接触中考虑厚度的变化。

? 如果需要考虑薄膜模式或弯曲模式的沙漏问题,或模型中有面内弯曲,在 ABAQUS/Standard 中使用S4单元(4节点四边形有限薄膜应变线性完全积分壳单元)可以获得很高的精度。 ? S4R单元(4节点四边形有限薄膜应变线性减缩积分壳单元)性能稳定,适用范围很广。

? S3/S3R 单元(3节点三角形有限薄膜应变线性壳单元)可以作为通用壳单元使用。由于单元中的常应变近似,需要划分较细的网格来模拟弯曲变形或高应变梯度。

? 对于复合材料,为模拟剪切变形的影响,应使用适于厚壳的单元(例如 S4、S4R、S3、S3R、S8R),并要注意检查截面是否保持平面。

? 四边形或三角形的二次売单元对剪切自锁或薄膜自锁都不敏感,适用于一般的小应变薄壳。

? 在接触模拟中,如果必须使用二次单元,不要选择 STRI65单元(三角形二次壳单元),而应使用 S9R5 单元(9节点四边形壳单元)。

? 如果模型规模很大且只表现几何线性,使用 S4R5单元(线性薄壳单元)比通用壳单元更节约计算成本。 ? 在ABAQUS/Explicit中,如果包含任意大转动和小薄膜应变,应选用小薄膜应变单元。 4)

选择梁单元的类型

如果一个构件横截面的尺寸远小于其轴向尺度(一般的判据为小于1/10),并且沿长度方向的应力是最重要的因素,就可以用梁单元来模拟此结构。ABAQUS中的所有梁单元都是梁柱类单元,即可以产生轴向变形、弯曲变形和扭转变形。Timoshenko 梁单元还考虑了横向剪切变形的影响。B21 和 B31 单元(线性梁单元)以及B22和 B32单元(二次梁单元)是考虑剪切变形的Timoshenko梁单元,它们既适用于模拟剪切变形起重要作用的深梁,又适用于模拟剪切变形不太重要的细长梁。这些单元的横截面特性与厚壳单元的横截面特性相同。

abaqus建模流程 - 学习笔记

行单元划分。当然,可以采用四面体单元(tetrahedral)利用自由网格技术对任何形状的模型区域进行单元划分。????自由划分:可以应用到任意平面和曲面结构划分:结构化的网格划分通常给出了对网格的最大的控制扫掠划分:在源面中,相邻面之间的二面角不能和180°相差太远。虚拟拓扑:在某些情况下,装配件的部件实例可能包含一些小的细节,比如
推荐度:
点击下载文档文档为doc格式
7dbtz4plgr10e619mkhh
领取福利

微信扫码领取福利

微信扫码分享