巴西圆盘中心裂纹摩擦接触分析.docxVIP

巴西圆盘中心裂纹摩擦接触分析 随着材料和结构材料的应用，材料、结构和结构之间的关系已成为研究的重要内容。胶结是材料和结构连接中重要的一种连接方式,与其他的一般连接形式相比,胶结具有不需要连接孔和紧固件,载荷分布区域大、结构重量轻等优点。然而胶结带来界面两侧材料性质失配(如热胀系数、弹性模量等) 含界面裂纹的复合结构承载时,结构往往沿界面发生复合模态破坏,而界面断裂韧性G 实际上,当双材料Brazilian盘所加荷载的荷载角增加或裂纹长度较长时,裂尖附近的裂纹面可能产生相互嵌入的不合理现象。为此,本文提出了一种新的考虑摩擦的接触算法——逐点Lagrange乘子法来研究Brazilian盘中心裂缝问题。为避免对总刚度阵求逆,将Lagrange乘子逐点对转到局部坐标系下,采用Gauss-Seidel迭代法求解局部坐标系下的法向和切向乘子,注意在求解的过程中加上对切向乘子约束修正。之后代换到整体坐标系下采用拟高斯迭代法求解整体坐标系下的乘子,将乘子全部求解更新后代入方程总刚度阵中迭代求解位移U,如此反复迭代修正乘子uf06c和位移U直至达到收敛条件。由于采用逐点对求解,计算机存储量很小,对接触面上点的数目要求不高,因此可以采用精细网格进行计算,提高了接触应力的计算精度。基于本文的算法开发了相应的程序,对Brazilian盘中心裂纹和双材料Brazilian盘中心界面裂纹进行了计算,通过位移线性外插法得到了裂尖应力强度因子,并与已有参考文献进行了比较。对光滑接触时考虑荷载角和裂缝长度对位移,接触区和应力强度因子的影响,以及考虑摩擦时摩擦系数对位移和应力强度因子的影响。 1 局部法向乘子法 考虑有摩擦的接触问题,引入Lagrange乘子uf06c,根据库仑摩擦模型,接触摩擦问题等价于如下约束优化问题: 其中:[A]为总刚度阵;{λ}表示接触力;{U}是位移场;[B]表示位移与接触力的相关矩阵,对于每个接触点对;{g}为点对中的初始间隙;{F}为荷载;表示局部坐标系下的乘子,[T]为整体坐标与局部坐标之间的转换矩阵;表示局部法向乘子;表示局部切向乘子;为摩擦系数。uf06d 式(1)等价于: 将式(2)中第一式写成分量形式: 式(3)写成矩阵形式: 对二维情形,其中i表示第i个点对; 于是结合点对约束条件,类似摩擦接触的做法,由下式逐点求解Lagrange乘子 待式(5)求解完成后,再通过 2 算法的迭代算法 具体逐点Lagrange乘子算法如下: 步I提取[A 步II迭代步赋初值P=0,接触力赋初值 步III迭代步赋值P=P+1,逐点对循环 用Gauss-Seidel迭代法对式(6)中法向接触力和切向接触力进行交替迭代,在求解过程中采用式(7)对乘子进行修正: 并用类似Uzawa的算法格式 步IV将所有的乘子全部更新后,通过 步V如此反复,直到 3 局部充放电坐标系x轴 在图1所示坐标系下,由两种各向同性材料制成的界面裂纹的双材料弹性体,其裂尖附近的应力场为 其中:r和θ是极坐标;局部笛卡尔坐标系x轴是沿双材料界面的; 式中: 界面裂纹前沿(θ=0)的应力为: 式(10)中引入固定特征长度 于是式(10)变为: 沿着裂缝面(θ=π)的位移场表示为: 其中, 类似地,通过解式(12),单独的应力强度系数在裂尖能够由裂尖位移场确定 当裂尖的应力和位移场通过有限元计算得到,应力强度因子可通过外插法得到。 4 计算值的例和讨论 4.1 集中力p的荷载角 对Brazilian盘(如图2),取圆盘的半径R=30mm,半裂缝长度a=15mm。材料的弹性模量E=2.9GPa,n=0.29。集中力P在荷载角θ=5°、10°、15°、20°、25°、35°、45°、55°、65°、75°、85°下进行加载。有限元计算模型见图3。对θ=5°、10°、15°、20°时,沿整个裂缝面是张开的,因此K 通过计算得到的裂纹面法向间隙Δu 对θ=35 图6给出了对摩擦系数f=0.0,0.5,1.0时, 4.2 bragilin盘中心界面裂纹面法向间隙的应力强度因子 双材料Brazilian盘(如图7),几何尺寸同均质材料,圆盘的半径R(28)30mm,取半裂缝长度与半径比0.1≤a/R≤0.8。材料1玻璃弹性模量E 裂纹面法向间隙Δu 对裂缝长度比a/R=0.5在不同荷载角θ下的相对滑移Δu 当荷载角θ=15°时,不同裂缝长度比下的法向间隙 A端无量纲应力强度因子 双材料Brazilian盘中心界面裂纹法向间隙和切向相对位移与摩擦系数的关系如图16、图17。在θ=30°,裂缝长度比a/R=0.5下,随摩擦系数增大,法向间隙减小,接触面增大,切向相对位移减小。无量纲应力强度因子随摩擦系数变化见图18所示。可知在f=1时相对于 5 结果验证和应用 (1)采用了逐点Lagr