带孔平板的应力分布及应力集中系数的计算.docx

带孔平板的应力分布及应力集中系数的计算.docx

《带孔平板的应力分布及应力集中系数的计算.docx》由会员分享，可在线阅读，更多相关《带孔平板的应力分布及应力集中系数的计算.docx（18页珍藏版）》请在冰豆网上搜索。

带孔平板的应力分布及应力集中系数的计算

带孔平板的应力分布及应力集中系数的计算

、问题重述

计算带孔平板的应力分布及应力集中系数。

、模型的建立与计算

在ANSY防建立模型，材料的设置属性如下

分析类型为结构（structural），材料为线弹性（LinearElastic），各向同

性（Isotropic，。

弹性模量、泊松比的设定均按照题目要求设定，以N、cm为标

准单位，实常数设置中设板厚为1

采用solid4node42板单元，ElementBehavior设置为Planestrsw/thk。

/\PLANE42eementtypecations

OptionsforPLANE42,ElementTypeRer

ElementcaardsystemdefinedKI

EsrtradisplacementshapesK2

HelpI

Elementbehavior

Ejctrastressoutput

Extrasurfaceoutput

OK

K3

K5

K6

Cancel

建立模型时先建立完整模型，分别用单元尺度为5cm左右的粗网格和单元尺

度为2cm左右的细网格计算

然后取四分之一模型计算比较精度，为了使粗细网格单元数与完整模型接近,

四分之一模型分别用单元尺度为2.5cm左右的粗网格和单元尺度为1cm左右的细网格计算

（1）完整模型的计算

①粗网格

单元网格的划分及约束荷载的施加如图（单元尺度为5cm）

AIM

h«R30201414:

19:

23

Stress

约束施加时在模型左侧边界所有节点上只施加x方向的约束，即令UX=0，在

左下角节点上施加x、y两个方向的约束，即U=0、UY=0o荷载施加在右侧边界

上，大小为100

对模型进行分析求解得到:

节点应力云图（最大值222.112）

单元应力云图（最大值256.408）

可看出在孔周围有应力集中现象，其余地方应力分布较为均匀，孔上部出现最大应力。

②细网格

单元网格的划分及约束荷载的施加如图（单元尺度为2cm

AN

MK302014

twrrg1

约束及荷载的施加方法如前，对模型进行分析求解得到:

str^s仏tHiLution.azialy止

节点应力云图（最大值272.484）

E_ii-iL.rrj（.u7rict:

SI3B-1

TZMEXL

直；刘J-

SI-K-

右加站.盟9

3?

.me

Stl吕3吕diwtrlLutiH」创-jalywis

（2）取1/4模型的计算

①粗网格

单元网格的划分及约束荷载的施加如图（单元尺度为2.5cm）

ELEMlcrS

AN

2.-.uir；4

15:

04:

49PLOTQ1

约束施加时在模型左侧边界所有节点上只施加x方向的约束，即UX=0,在下

侧边界所有节点上只施加y方向的约束，即LY=0o荷载施加在右侧边界上，大小

为100o

对模型进行分析求解得到:

NCCfiL90UUTICH

SIEE^l

5UE=1

TM>1

5ECV腔刪

EK<.024958

3bNYL254swt=£51,333

瞬30201415:

20147PLOTNO*1

7亠心64止列"'W111

4'523157.9651°-'-（7722^^1251+^

str湖ganddl^trlLutiorjanalysis

节点应力云图（最大值251.333）

EMHEWT90UJTICN

STTEP-1

SCE=1

TM=1

:

却-NOA^

MX=.n74?

，5F

拠=55.147

139.0^1

jS.147a?

.□&9

^1J18113.06

Less如distiibatioaanalysis

165.

190.

②细网格

AN

bfiR30201415:

21:

02ttnrbio,1

单元网格的划分及约束荷载的施加如图（单元尺度为1cm）

Sti^ss如dist£Lb.ition.

AN

MK为201415:

34;ID

PLJIg

约束及荷载的施加方法如前，对模型进行分析求解得到:

NODAL沁LUT丄工

S7EP-1

SUB=1

TIME-15ECV（AVG）CMK=.025023

廉刑-^S.327

SMK=29O^7a

57.455115.711

Stress知di5txijxiti

NPlR铀201415:

36:

58ZTJ3TNO.1

173.966232,222右乩4%

EJlIE'TSCiETrCJ]

SUE=1

TIME^lS37JUHK

S>H日⑼

1OI兀

=+025Q23-26,512-.97.137

■・■

■■

如如阪时的遡ne.72倔小讹•遡

二g疝

PLOTNO.

stressauddistriLuti^aiialysi^

AN

MHR30201415捋了；361

节点应力云图（最大值290.478）

（3）计算结果比较

下面按照弹性力学理论求解带孔平板的应力集中系数。

（参考曾攀.有限元

分析及应用.清华大学出版社）

设在无限大薄板中有一个半径为R0的圆孔，该无限大薄板在x方向受有

匚xx的均匀荷载如图

对于无限大板宽的孔边应力集中问题，基于以上平面极坐标下的三大类基本

方程，可以得到以下弹性状态下的解析解:

最大的环向应力为

口（八90,「=Ro）=3q°

从以上推导可知，对此类带孔平板应力集中系数的弹性力学精确解为3。

计算应力集中系数，对结果进行汇总，与弹性力学精确解进行比较如下：

（应力集中系数采用节点应力云图数据/100计算）

项目

粗网格

细网格

弹性力学

精确解

误差（%

粗网格

细网格

完整模型

2.221

2.725

3.000

26.0

9.2

四分之一模型

2.513

2.905

3.000

16.2

3.2

结论：

此种以单元大小为基准划分网格的方式使得计算结果偏小。

模型计算时，单元划分的越精细，结果越精确。

计算报告中由于完整模型划分的单元数与四分之一模型划分的单元数

接近，使得四分之一模型的计算精度较高，但也可看出利用对称性取完

整模型的四分之一计算也能得到较好的结果。