首页 / 专利库 / 物理 / 应力 / 平面应力 / 平面双材料圆环界面应力奇异性特征值的分析方法

平面双材料圆环界面应奇异性特征值的分析方法

阅读:506发布:2020-05-13

专利汇可以提供平面双材料圆环界面应奇异性特征值的分析方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种平面双材料圆环界面应 力 奇异性特征值的分析方法,其特征是首先判定近似奇异点,利用 选定 的可行方向和计算间隔选取测算点,记录测算点的正 应力 分量及各测算点到近似奇异点的距离,获得测算点各正应力分量和距离的双对数分布图,根据双对数分布图判断应力奇异性特征值为单解或为双解,在应力奇异性特征值为双解时继续选取测算点的切应力分量,并按相应的方法分析获得应力奇异性特征值。本发明方法能够方便地计算得到双材料圆环界面的应力奇异性特征值。,下面是平面双材料圆环界面应奇异性特征值的分析方法专利的具体信息内容。

1.平面双材料圆环界面应奇异性特征值的分析方法,其特征按如下步骤进行:
步骤1、在平面双材料圆环界面的应力分布图中取应力集中部分的中心为近似奇异点A;自近似奇异点A起向某一选定方向取一段线段AB作为计算间隔,在所述线段AB中平均选取Z个测算点,设定θ为所述选定方向的极,设定k为1到Z之间的整数、k=1,
2,3……Z,分别记录各测算点到所述近似奇异点A的距离rk以及各测算点的正应力分量σθk,对于所述极角θ以及线段AB的设定,要求所述正应力分量σθk和距离rk满足线性关系,令所述选定方向为可行方向;
步骤2、利用式(1)获得正应力分量σθk和距离rk的双对数分布图P1:
logσθk=C0-(1-λ)logrk (1)
式(1)中:C0为常数,λ为平面双材料圆环界面的应力奇异性特征值;
若所述双对数分布图P1呈直线性,则是只有一个λ值,即为具有单一应力奇异性,求取所述双对数分布图P1的斜率即为特征值λ;反之则是λ为双解,即为具有二重应力奇异性,对于λ为双解继续如下步骤3;
步骤3、对于λ为双解:在线段AB中选取所述Z个测算点的切应力分量τrθk,并要求所述切应力分量τrθk与距离rk同样具有线性关系,否则返回步骤1,选定不同的可行方向和计算间隔重新操作,直至切应力分量τrθk与距离rk具有线性关系,则有:
式(2)中:B1和B2为常数,λ1和λ2为平面双材料圆环界面的两个应力奇异性特征值,系数A1和A2的值由步骤4确定;
步骤4、在Z个测算点中选取距离rk成等比数列的三个点,分别为M点、N点和O点,令M点、N点和O点与近似奇异点A的距离分别为r0,ρr0,ρ2r0,其中ρ为比例系数;令M点的正应力分量和切应力分量分别为σθM和τrθM;N点的正应力分量和切应力分量分别为σθN和τrθN;O点的正应力分量和切应力分量分别为σθO和τrθO,则将式(2)转换为式(3):
式(3)中:
1 2 2
当b=4ac时,A1和A2为重根;当b >4ac时,A1和A2为两个实根;当b <4ac时,A1和A2为一对共轭复根;
至少取10组不同的r0与ρ值,分别计算各组A1和A2,要求求得的各组A1值之间相差不超过0.05,并且各组A2值之间亦相差不超过0.05,否则返回步骤1,选定不同的可行方向和计算间隔重新操作,直至各组A1值之间相差不超过0.05,并且各组A2值之间亦相差不超过0.05;取各组A1的平均数为终值 取各组A2的平均数为终值
步骤5、所述终值 和 具有如下三种不同的形式
形式一: 和 为重根,表明近似奇异点A仅具有单一应力奇异性,即为步骤2中具有单一应力奇异性的情形;
形式二: 和 为两个实根,则有式(4):
利用式(4)获得测算点各组合应力分量 n=1,2和距离rk的双对数分
布图P2;若是所述双对数分布图P2中两个图像均呈直线性,求取所述双对数分布图P2中两个图像的斜率即分别为实数特征值λ1和λ2,同时表明该圆环界面具有二重实应力奇异性;
形式三: 和 为一对共轭复根,特征值λ1和λ2也为一对共轭复数;令:
λ1=λR+λMi,λ2=λR-λMi;其中AR为 和 的
实部,AM为 和 的虚部;λR为λ1和λ2的实部,λM为λ1和λ2的虚部;i为虚数单位;则有:
式 (5) 中 :
e为自然常数,D1和D2为常
数;
利用式(5)分别获得测算点各组合应力分量σ0和距离rk以及ψ0和距离rk的双对数分布图P3;若是所述双对数分布图P3中两个图像均呈直线性,求解所述双对数分布图P3中两个图像的斜率分别为λR和λM,最终得到两个共轭复数特征值λ1和λ2,同时表明该圆环界面具有二重振荡应力奇异性。

说明书全文

平面双材料圆环界面应奇异性特征值的分析方法

技术领域

[0001] 本发明涉及应力奇异性分析方法,更具体地说是一种平面双材料圆环界面应力奇异性特征值的分析方法。

背景技术

[0002] 双材料圆环结构作为一种新型复合材料结构,已被广泛地应用于海底石油输送管道等领域。由于两侧材料性能差异性,存在应力奇异场,承受载荷时容易在界面附近形成严重的应力不均,是引起界面裂纹及其脱粘的主要原因。
[0003] 很多学者基于二维弹性理论,针对平面双材料正交楔形体点附近应力场的奇异性问题,通过Airy函数得出双材料界面端具有应力奇异性。借助界面端奇异性Dundurs参数,采用Mellin变换可以推导任意接合角组合条件下界面端附近应力奇异性特征方程,确定双材料的材料特性、界面端的几何形状和界面端应力场强度三者之间的相互关系。用界面端特征值λ表征界面端应力场强度,λ最多有6个解,1-λ为应力奇异性次数。应用Goursat复变应力函数,可以推导具有任意几何结合形状的界面端的应力和位移场解析解。已有研究对开式界面应力奇异性次数进行了研究,但对于圆环状的闭式界面并未给出相关的结论。

发明内容

[0004] 本发明是为避免上述现有技术所存在的不足之处,提供一种平面双材料圆环界面应力奇异性特征值的分析方法,以期能够方便地计算得到双材料圆环界面的应力奇异性特征值,为评价平面双材料圆环状闭式界面结构的力学性能提供依据。
[0005] 本发明为解决技术问题采用如下技术方案:
[0006] 本发明平面双材料圆环界面应力奇异性特征值的分析方法的特点是按如下步骤进行:
[0007] 步骤1、在平面双材料圆环界面的应力分布图中取应力集中部分的中心为近似奇异点A;自近似奇异点A起向某一选定方向取一段线段AB作为计算间隔,在所述线段AB中平均选取Z个测算点,设定θ为所述选定方向的极角,设定k为1到Z之间的整数、k=1,2,3......Z,分别记录各测算点到所述近似奇异点A的距离rk以及各测算点的正应力分量σθk,对于所述极角θ以及线段AB的设定,要求所述正应力分量σθk和距离rk满足线性关系,令所述选定方向为可行方向;
[0008] 步骤2、利用式(1)获得正应力分量σθk和距离rk的双对数分布图P1:
[0009] logσθk=C0-(1-λ)logrk (1)式(1)中:C0为常数,λ为平面双材料圆环界面的应力奇异性特征值;
[0010] 若所述双对数分布图P1呈直线性,则是只有一个λ值,即为具有单一应力奇异性,求取所述双对数分布图P1的斜率即为特征值λ;反之则是λ为双解,即为具有二重应力奇异性,对于λ为双解继续如下步骤3;
[0011] 步骤3、对于λ为双解:在线段AB中选取所述Z个测算点的切应力分量τrθk,并要求所述切应力分量τrθk与距离rk同样具有线性关系,否则返回步骤1,选定不同的可行方向和计算间隔重新操作,直至切应力分量τrθk与距离rk具有线性关系,则有:
[0012]
[0013] 式(2)中:B1和B2为常数,λ1和λ2为平面双材料圆环界面的两个应力奇异性特征值,系数A1和A2的值由步骤4确定;
[0014] 步骤4、在Z个测算点中选取距离rk成等比数列的三个点,分别为M点、N点和O2
点,令M点、N点和O点与近似奇异点A的距离分别为r0,ρr0,ρr0,其中ρ为比例系数;
令M点的正应力分量和切应力分量分别为σθM和τrθM;N点的正应力分量和切应力分量分别为σθN和τrθN;O点的正应力分量和切应力分量分别为σθO和τrθO,则将式(2)转换为式(3):
[0015]
[0016] 式 (3) 中: b=σθMτrθO+σθOτrθM-2σθNτrθN,2 2
当b=4ac时,A1和A2为重根;当b >4ac时,A1和A2为两个实根;当
2
b<4ac时,A1和A2为一对共轭复根;
[0017] 至少取10组不同的r0与ρ值,分别计算各组A1和A2,要求求得的各组A1值之间相差不超过0.05,并且各组A2值之间亦相差不超过0.05,否则返回步骤1,选定不同的可行方向和计算间隔重新操作,直至各组A1值之间相差不超过0.05,并且各组A2值之间亦相差不超过0.05;取各组A1的平均数为终值 取各组A2的平均数为终值
[0018] 步骤5、所述终值 和 具有如下三种不同的形式
[0019] 形式一: 和 为重根,表明近似奇异点A仅具有单一应力奇异性,即为步骤2中具有单一应力奇异性的情形;
[0020] 形式二: 和 为两个实根,则有式(4):
[0021]
[0022] 利用式(4)获得测算点各组合应力分量 n=1,2和距离rk的双对数分布图P2;若是所述双对数分布图P2中两个图像均呈直线性,求取所述双对数分布图P2中两个图像的斜率即分别为实数特征值λ1和λ2,同时表明该圆环界面具有二重实应力奇异性;
[0023] 形式三: 和 为一对共轭复根,特征值λ1和λ2也为一对共轭复数;令:
[0024] λ1=λR+λMi,λ2=λR-λMi;其中AR为 和的实部,AM为 和 的虚部;λR为λ1和λ2的实部,λM为λ1和λ2的虚部;i为虚数单位;则有:
[0025]
[0026] 式(5)中: e为自然常数,D1和D2为常数;
[0027] 利用式(5)分别获得测算点各组合应力分量σ0和距离rk以及ψ0和距离rk的双对数分布图P3;若是所述双对数分布图P3中两个图像均呈直线性,求解所述双对数分布图P3中两个图像的斜率分别为λR和λM,最终得到两个共轭复数特征值λ1和λ2,同时表明该圆环界面具有二重振荡应力奇异性。
[0028] 与已有技术相比,本发明有益效果体现在:
[0029] 1、本发明能够方便地计算得到双材料圆环界面的应力奇异性特征值,为评价平面双材料圆环状闭式界面结构的力学性能提供了依据,也可以为双材料圆环结构材料组合的选择和设计提供依据。
[0030] 2、本发明方法中采用数值分析方法,引入了测算点,其过程简单、易于操作;采用近似奇异点代替奇异点,显著简化了计算过程。
[0031] 3、本发明方法对于可行方向以及对应的计算间隔的选取进行了多重限定,提高精确性。
[0032] 4、本发明方法中选取测算点中多组成等比数列的三个点分别进行A1值和A2值的计算,并检验各组A1值和A2值之间的差值,通过求取平均数的方法得到终值,避免了偶然性,进一步提高了精确性。
[0033] 5、本发明方法中采用图像法,利用双对数分布图判断并计算应力奇异性特征值,直观可视,进一步简化了计算过程。附图说明
[0034] 图1为本发明方法中涉及的平面双材料圆环结构的应力分布图;
[0035] 图2为本发明方法中平面双材料圆环结构的应力集中图;
[0036] 图3为本发明方法中单一应力分量的双对数分布图P1;
[0037] 图4为本发明方法中组合应力分量的双对数分布图P2;
[0038] 表1为本发明方法中等比数列点组及其相对应的A1值和A2值。

具体实施方式

[0039] 本实施例中平面双材料圆环界面应力奇异性特征值的分析方法是按如下过程进行:
[0040] 步骤1:在平面双材料圆环界面的应力分布图中取应力集中部分的中心为近似奇异点A,圆环界面的应力奇异点很难精确地确定,可以确定的是应力奇异点一定在应力集中部分中心的附近,因此本实施例中用近似奇异点代替奇异点,以此显著简化计算过程,并且对计算结果不产生影响;自近似奇异点A起向某一选定方向取一段线段AB作为计算间隔,在线段AB中平均选取Z个测算点,设定θ为选定方向的极角,设定k为1到Z之间的整数、k=1,2,3......Z,分别记录各测算点到近似奇异点A的距离rk以及各测算点的正应力分量σθk,对于极角θ以及线段AB的设定,要求正应力分量σθk和距离rk满足线性关系,令选定方向为可行方向;若是正应力分量σθk和距离rk不满足线性关系,则表明取定的可行方向和计算间隔远离应力奇异性所影响的区域,此时选定的测算点不能用来计算应力奇异性特征值,故必须重新选定可行方向和计算间隔;测算点、各测算点到近似奇异点A的距离rk以及各测算点的正应力分量σθk的选取可以借助于通用的有限元软件
[0041] 步骤2:利用式(1)获得正应力分量σθk和距离rk的双对数分布图P1:
[0042] logσθk=C0-(1-λ)logrk (1)
[0043] 式(1)中:C0为常数,λ为平面双材料圆环界面的应力奇异性特征值;
[0044] 经大量计算发现,平面圆环结构普遍只具有单解或双解;若是双对数分布图P1呈直线性,则是只有一个λ值,即为具有单一应力奇异性,求取双对数分布图P1的斜率即为特征值λ;反之则是λ为双解,即为具有二重应力奇异性;对于λ为双解继续如下步骤3;
[0045] 步骤3:对于λ为双解:在线段AB中选取Z个测算点的切应力分量τrθk,并要求切应力分量τrθk与距离rk同样具有线性关系,否则返回步骤1,选定不同的可行方向和计算间隔重新操作,直至切应力分量τrθk与距离rk具有线性关系,这一方式是为了保证可行方向和计算间隔在应力奇异性所影响的区域内,从而进一步保证测算点的有效性,则有:
[0046]
[0047] 式(2)中:B1和B2为常数,λ1和λ2为平面双材料圆环界面的两个应力奇异性特征值,系数A1和A2的值由步骤4确定;
[0048] 步骤4、在Z个测算点中选取距离rk成等比数列的三个点,分别为M点、N点和O2
点,令M点、N点和O点与近似奇异点A的距离分别为r0,ρr0,ρr0,其中ρ为比例系数;
令M点的正应力分量和切应力分量分别为σθM和τrθM;N点的正应力分量和切应力分量分别为σθN和τrθN;O点的正应力分量和切应力分量分别为σθO和τrθO,则将式(2)转换为式(3):
[0049]
[0050] 式 (3) 中: b=σθMτrθO+σθOτrθM-2σθNτrθN,
[0051] 当b2=4ac时,A1和A2为重根;当b2>4ac时,A1和A2为两个实根;当b2<4ac时,A1和A2为一对共轭复根;
[0052] 至少取10组不同的r0与ρ值,分别计算各组A1和A2,要求求得的各组A1值之间相差不超过0.05,并且各组A2值之间亦相差不超过0.05,否则返回步骤1,选定不同的可行方向和计算间隔重新操作,直至各组A1值之间相差不超过0.05,并且各组A2值之间亦相差不超过0.05;取多组不同的r0与ρ值,并要求各组A1和A2值的相关性,这样既可以避免计算中的偶然性,又可以减小通用有限元软件计算模拟的误差;一旦各组A1和A2值的相关性不满足要求,则同样表明取定的可行方向和计算间隔远离应力奇异性所影响的区域,此时选定的测算点不能用来计算应力奇异性特征值,故必须返回步骤1重新操作;取各组A1的平均数为终值 取各组A2的平均数为终值
[0053] 步骤5、终值 和 具有如下三种不同的形式
[0054] 形式一: 和 为重根,表明近似奇异点A仅具有单一应力奇异性,即为步骤2中具有单一应力奇异性的情形;
[0055] 形式二: 和 为两个实根,则有式(4):
[0056]
[0057] 利用式(4)获得测算点各组合应力分量 n=1,2和距离rk的双对数分布图P2;若是双对数分布图P2中两个图像均呈直线性,求取双对数分布图P2中两个图像的斜率即分别为实数特征值λ1和λ2,同时表明该圆环界面具有二重实应力奇异性;
[0058] 形式三: 和 为一对共轭复根,特征值λ1和λ2也为一对共轭复数;令:
[0059] λ1=λR+λMi,λ2=λR-λMi;其中AR为 和的实部,AM为 和 的虚部;λR为λ1和λ2的实部,λM为λ1和λ2的虚部;i为虚数单位;则有:
[0060]
[0061] 式(5)中: e为自然常数,D1和D2为常数;
[0062] 利用式(5)分别获得测算点各组合应力分量σ0和距离rk以及ψ0和距离rk的双对数分布图P3;若是双对数分布图P3中两个图像均呈直线性,求解双对数分布图P3中两个图像的斜率分别为λR和λM,最终得到两个共轭复数特征值λ1和λ2,同时表明该圆环界面具有二重振荡应力奇异性。
[0063] 参见图1,本实施例中在内圆环1和外圆环2之间形成有平面双材料圆环界面3,其中内圆环1为6061-T651合金、外圆环2为AZ91D镁合金,材料参数分别为:镁合金剪切模量G1=16.67GPa,泊松比ν1=0.35;铝合金剪切模量G2=26.69GPa,泊松比ν2=0.3。如图1所示,圆环界面的应力奇异点具有对称性,故只计算圆环界面上的一个奇异点,即图1所示的第二象限中的奇异点,图2为其局部放大图,具体实施步骤如下:
[0064] 第1步、在平面双材料圆环界面的应力分布图中取应力集中部分的中心为近似奇异点A;自近似奇异点A起向某一选定方向取一段线段AB作为计算间隔,线段AB长度取定为0.76mm;在线段AB中平均选取50个测算点,即Z取定为50,选定方向的极角θ取定为0.75π,设定k为1到50之间的整数、k=1,2,3......50,可行方向选定为由近似奇异点指向圆心;分别记录各测算点到近似奇异点A的距离rk以及各测算点的正应力分量σθk,经检验近似奇异点A的距离rk以及各测算点的正应力分量σθk满足线性关系;
[0065] 第2步、利用式(1)获得正应力分量σθk和距离rk的双对数分布图P1;如图3所示,双对数分布图P1不呈直线性,则表明该圆环界面具有二重应力奇异性,即λ有双解;
[0066] 第3步、在线段AB中选取50个测算点的切应力分量τrθk,经检验所述切应力分量τrθk与距离rk同样具有线性关系;
[0067] 第4步、在50个测算点中选取距离rk成等比数列的三个点,分别为M点、N点和O2
点,令M点、N点和O点与近似奇异点A的距离分别为r0,ρr0,ρr0,其中ρ为比例系数;
令M点的正应力分量和切应力分量分别为σθM和τrθM;N点的正应力分量和切应力分量分别为σθN和τrθN;O点的正应力分量和切应力分量分别为σθO和τrθO,本实施例选取其中11组成等比数列的测算点,并将各组M点、N点和O点的应力分量带入式(3),分别计算各组A1和A2值,具体r0和ρ的取值及相对应的计算所得的A1和A2值如表1所示。经检验求得各组A1值之间相差不超过0.05,并且各组A2值之间亦相差不超过0.05,取各组A1的平均数为终值 取各组A2的平均数为终值 分别为
[0068] 第5步、由于求得的 和 为两个实根,将其带入式(4),并利用式(4)获得测算点各组合应力分量 n=1,2和距离rk的双对数分布图P2;如图4所示,双对数分布图P2呈直线性,求取双对数分布图P2中两个图像的斜率,分别得到实数特征值λ1=0.9332、λ2=0.9455;同时表明本例中该圆环界面具有二重实应力奇异性。
[0069] 表1
[0070]
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈