二维/三维零厚度单元批量插入_Abaqus-仿真秀干货文章

二维/三维零厚度单元批量插入

零

概述

零厚度单元室abaqus中的一种特殊单元，该单元在初始状态下为零厚度，但是同一个空间位置却有两个节点编号，当单元参与计算以后，对应的节点发生位移，零厚度改变为有厚度。

这种单元常常用于模拟界面问题，在规则的界面中，可以手动创建零厚度单元，在不规则界面，如混凝土内部的骨料-砂浆界面，需要插入数万计的零厚度单元，单靠手动创建是不可能的，这时需要编写程序创建插入零厚度单元。

不同的网格需要插入不同的零厚度单元，下面是一些示例。

壹

三角形网格中插入COH2D3单元

贰

四边形网格中插入COH2D4单元

叁

三角形与四边形混合网格中插入COH2D4单元

肆

四面体网格中插入COH3D6单元

伍

楔形体网格中插入COH3D6和COH3D8单元

陆

六面体网格中插入COH3D8

柒

四面体与楔形体混合网格中插入COH3D6和COH3D8单元

捌

六面体与楔形体混合网格中插入COH3D6和COH3D8单元

玖

六面体、楔形体与四面体混合网格中插入COH3D6和COH3D8单元

来源：有限元先生

abaqus dload子程序实现复杂移动荷载

零概述 abaqus给用于提供了dload子程序用以满足复杂的荷载需求，帖子采用该子程序模拟了四轮小车在平板上移动的运动过程。壹几何模型信息模型为四轮车，每个轮子与地面的接触区域为0.1x0.3（不考虑实际情况），车辆在平板居中放置，从右边界启动，向左移动采用C3D8单元划分网格网格横断面车轮经过的区域网格设置加密，车辆中间区域与外部区域网格尺寸设置渐变。其中，车轮经过区域的网格贰程序 ! 模拟四个车轮经过桥面，动荷载 SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT, 1 COORDS,JLTYP,SNAME)C INCLUDE 'ABA_PARAM.INC'C DIMENSION TIME(2), COORDS (3) CHARACTER*80 SNAME! FIXED INCRANTMENT TIME IS 0.01S DOUBLE PRECISION INC_TIME DOUBLE PRECISION F,DISP INC_TIM=0.01! 10s，增量步长为0.01s，总增量步数为1000 ! 荷载表达式：F=sin(2*3.14/10)! 桥梁尺寸：桥长40，桥宽6 X=COORDS(1) Y=COORDS(2) Z=COORDS(3) DISP=TIME(2)*(40/10) IF((Z>=2.2.AND.Z<=2.5).OR.(Z>=3.5.AND.Z<=3.8))THEN IF((X>DISP-0.1).AND.X<(DISP))THEN F=1 ELSEIF((X>DISP-0.1-2).AND.(X<DISP-2))THEN F=1 ELSE F=0 ENDIF ELSE F=0 ENDIF RETURN END 来源：有限元先生