XFEM 模拟裂纹扩展INP文件中裂纹定义

wrj001

各位虫子:XFEM，模拟裂纹扩展INP ，注意黑体部分最后两个数字的含义[照PDF文件做的例子生成的INP]

*Node #模型节点

……

*Element, type=C3D8R,elset=Elem #模型单元

……

*Solid Section, elset=Elem, material=Aluminium #单元赋材料

*Nset, nset=_PickedSet10, internal,instance=Solid-1, generate

22,  484,   22

*Elset, elset=_PickedSet10, internal,instance=Solid-1, generate

21,  441,   21

*Nset, nset=_PickedSet11, internal,instance=Solid-1, generate

  1,  463,   22

*Elset, elset=_PickedSet11, internal,instance=Solid-1, generate

  1,  421,   21

*Elset, elset=__PickedSurf8_S6, internal,instance=Solid-1, generate

   1,  9241,    21

*Surface, type=ELEMENT, name=_PickedSurf8,internal

__PickedSurf8_S6, S6

*Elset, elset=__PickedSurf9_S4, internal,instance=Solid-1, generate

  21,  9261,    21

*Surface, type=ELEMENT, name=_PickedSurf9,internal

__PickedSurf9_S4, S4

*Enrichment, name=EdgeCrack,type=PROPAGATION CRACK, elset=_PickedSet6 #定义Enrichment区域

*Material, name=Aluminium #定义材料

*Damage Initiation, criterion=MAXPS

500.,

*Damage Evolution, type=DISPLACEMENT

1.,

*Elastic

70000., 0.33

*Initial Conditions, type=ENRICHMENT #含裂纹单元定义

6626,1,EdgeCrack, -0.0952381, -0.0476191

6626,2,EdgeCrack, 0.0952381, -0.0476191

……

*Step, name=Loading #定义载荷步，载荷

*Static

1., 1., 1e-05, 1.

*Dsload #载荷

_PickedSurf9, P, -1.

_PickedSurf8, P, -1.

*Enrichment Activation, name=EdgeCrack,activate=ON

*Restart, write, frequency=0

*Enrichment Activation, nam

*Output, field #输出要求

*Node Output

CF, PHILSM, RF, U

*Element Output, directions=YES

LE, PE, PEEQ, PEMAG, S

*Contact Output

CDISP, CSTRESS

*Output, history, variable=PRESELECT

*EndStep

2022-4-17 附加说明：*Initial Conditions, type=ENRICHMENT
6626, 1,EDGECrack, -0.0952381, -0.0476191
数据行：
1 元素编号或元素集标签；
2 构成元素连接性的节点的相对位置；
3 在*ENRICHMENT 上指定的扩展功能的名称；
4 第一个有符号距离函数的值。
5 第二个有符号距离函数的值。如果只需要第一个有符号距离函数，请将此项留空。
根据需要重复此数据行，以在各种元素或元素集中定义初始有符号距离函数。必须在元素内的所有节点上定义有符号距离函数。
每个节点通常需要两个有符号距离函数来描述裂纹几何体中的裂纹位置，包括裂纹尖端的位置。第一个有符号距离函数描述裂纹表面，而第二个用于构造正交表面，以便两个表面的交点定义裂纹前缘。第一个有符号距离函数仅指定给与裂纹相交的单元节点，而第二个有符号距离函数仅指定给包含裂纹尖端的单元节点。不需要裂纹的显式表示，因为裂纹完全由节点数据描述。

wrj001

分析一下，
-0.0952381,这个数字是4/21/2， 模型尺寸为4，单元有二十一层，裂纹位于第十一层中间，两侧最近单元节点的坐标大致是这个数；
-0.0476191 是上个数的一半。为距离裂尖的距离。
[2022-4-17重新编辑]