ABAQUS二次开发-带孔薄板圆角处网格研究

带孔薄板在单向拉力作用下，圆孔上的圆弧各处材料都处于单向受力状态。为方便研究，取四分之一圆弧，其应力分布的理论公式为 σ=σ(1-2cos2θ)，应力分布如下图

 

我们先用ABAQUS建立初始模型，获得rpy程序；

将后缀rpy改成后缀为py，就可以用python进行编辑了

接着删去了不必要的session操作语句；

然后对程序关键部分进行修改；

运行不同网格下的工况，生成曲线。

当使用一阶完全积分时，分别绘制了圆弧边缘种子数目为4,6,8,16划分网格后的模拟情况。

 

    结果分析：可以看到当使用一阶完全积分的方法和圆弧边缘种子数目为16的时候，拟合效果最贴近理论曲线。 

程序如下：

# -*- coding: mbcs -*-
from math import *
from abaqus import *
from abaqusConstants import *
from caeModules import *
###################parameters#################################
L = 300.0
W = 300.0
R = 3.0
Press = 10.0
NumList = [2, 3, 4, 8]
MkList = [HOLLOW_CIRCLE, HOLLOW_SQUARE, HOLLOW_DIAMOND, HOLLOW_TRI,
    CROSS, XMARKER, POINT]
curveList = []
Mdb()
vp = session.viewports['Viewport: 1']
for (i, NR) in enumerate(NumList):
    dR = R
    NumR = NR
    Size = 4.0
    NumdR = int(dR/(pi/4.0*R/NumR))
    ModelName = 'HolePlate'+str(i)
    INPName = 'HolePlate'+str(i)
    PNGName = 'CPS4_'+str(NumR)+'Element'
    dataName = 'Num='+str(NumR*2)
    #constructe the sketch
    m = mdb.Model(name=ModelName)
    s = m.ConstrainedSketch(name='plate', sheetSize=200.0)
    g = s.geometry
    g1 = s.Line(point1=(0.0, 0.0), point2=(0.0, L/2.0))
    g2 = s.Line(point1