1988
年第
1
j
到
北方交通犬学学报
JOURNAL
OF
NORTHERN
JIAOTONG
UNIVERSITY
No.l
1988
求解应力强度因子的一种新方法
吕晓寅
孙燕君
(工程力学研究所〉
摘要
本文介绍了求解混合型应为强度因子的一种新方法,即求解非线性方程组
的随机方法一一随机投点浊。这种方法是在一个选定的区域内进行投点,每投
→个点,计算一次目标函数的值,将所有的函数值进行比较,选出其中的最小
值,再以这个函数最小值所对应的根均中心,选定→个区域,在此区域内投
点,重复以上过程,这样偌环住复、直到求出合适的值均止。这种方法形象直
观,并且证明,它和其它好的算法一样,能得到满意的结果。
关键词
z
应为强度因子,随机投点。
一、引
古口
要得到应力强度因予的数值解,有很多方法,如
z
理论估计法、
Smith
外插浩、二
分法、一般Ir
win
岳、
Sanford-Da
lly
的最小二乘法、
Th
eocaris-Gdoutos
法等,其
中理论估计法需要乘→个放大系数1.
1
,这个放大系数可以从
Smith
、
Peters
和.A
ndo
nian(l)
的结果中推出,但这种方法需要对薄板上的穿透裂纹进行预测,所以得到的应
力强度因子的值被认为是近似的。
Olaosebikan
(2)
对
Smith
的外插法进行了改进。
二分法需要获得大量的数据,而且仅在
8=θm
时的数据是可靠的,因此要借助
lrwin
的
条件。一般
lrwin
法有坚实的数学基础,因此是比较可靠的。
Sanford-D
抖
ly
方法需要
比较小心地选取初始佳,否则用这个方法进行计算,问题会变的发散,即使收敛,也可
能收敛到→个不可靠的值。
1981
年
A.
S.
Kobayashi
等人
(3)
提出了用超定值最小二乘
法由动态等色线求应力强度因子的方法,这种方法利用以常速扩展裂纹的近场应力状态
导出动态应力强度因子
Kr
、
Krr
和远场应力分量
σ
ox
与动态等色线之间的关系,这种关
系与超定值最小二乘法一起,构成了从运动裂纹周围记录的动态光弹条纹图中求
Kr
、
Krr
和
σ
ox
的数据处理方法。
Freund
等人
(4)
【
5)
通过在不同载荷下由动态裂纹扩展的
本文
1987
年
8
月
22
日收到
-
29
~