python - Abaqus/Python 修复了播种约束

Question

我正在使用 Abaqus/Python 进行有限元分析。我遇到了以下问题:我有一个 3D 域，它由包含边和圆段的 2D 域组成，并且在三维中拉伸(stretch)。现在我想对该域进行网格划分，约束条件是只能使用沿厚度方向的一个元素。以下是 Python 代码中拉伸(stretch)四分之一圆的最小示例:

from abaqus import *
from abaqusConstants import *
import __main__

# Sketch quarter circle
s = mdb.models['Model-1'].ConstrainedSketch(name='__profile__', 
    sheetSize=50.0)
g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraints
s.setPrimaryObject(option=STANDALONE)

s.ArcByCenterEnds(center=(0.0, 0.0), point1=(0.0, 5.0), point2=(5.0, 0.0),
    direction=CLOCKWISE)
s.Line(point1=(0.0, 5.0), point2=(0.0, 0.0))
s.Line(point1=(0.0, 0.0), point2=(5.0, 0.0))


# Create part from sketch and extrusion in the depth
p = mdb.models['Model-1'].Part(name='Part-1', dimensionality=THREE_D,
    type=DEFORMABLE_BODY)
p.BaseSolidExtrude(sketch=s, depth=0.1)
s.unsetPrimaryObject()
session.viewports['Viewport: 1'].setValues(displayedObject=p)
del mdb.models['Model-1'].sketches['__profile__']

# Mesh
# Global seeds
p.seedPart(size=0.35, deviationFactor=0.1, minSizeFactor=0.1)
# 'fixed' seeding constraint: only one element in thickness direction
e = p.edges
pickedEdges = e.findAt(((5.0, 0.0, 0.025), ), ((0.0, 5.0, 0.025), ), 
    ((0.0, 0.0, 0.025), ))
p.seedEdgeByNumber(edges=pickedEdges, number=1, constraint=FIXED)

# Generate mesh
p.generateMesh()

在最后一行中，当我想要生成网格时，无法对域进行网格划分，因为无法遵守“固定”播种约束。如果我使用

p.generateMesh(seedConstraintOverride=ON)

为了生成网格，可以对域进行网格划分，但只能在厚度方向上使用两个单元。有没有一种方法可以在厚度方向上仅使用一个单元对域进行网格划分？如果不是，为什么？

Answer 1

有几种方法可以做到这一点。一种选择是简单地将网格技术更改为扫描。例如，假设您的零件由单个几何单元组成(如示例代码中所示)，您可以使用以下内容:

part_cells = p.cells()
p.setMeshControls(regions=(part_cells[0],), technique=SWEEP)
p.generateMesh()

您必须将 region 关键字参数指定为 Region 或 Set 对象的序列。请注意，如果扫描方向的长度小于全局种子大小，则无需指定任何边缘种子约束，并且只会生成一个元素。