材料力学之应力分析算法：应力强度因子（SIF）计算：SIF计算软件工具介绍与实践.docx

from fenics import *

import numpy as np

length = 1.0

def boundary(x, on_boundary):

return on_boundary

E = 1e5 #

nu = 0.3 #

mu = E / (2 * (1 + nu)) #

V = VectorFunctionSpace(mesh, 'Lagrange', 1)

bc = DirichletBC(V, Constant((0, 0)), boundary)

u = Function(V)

solve(a == L, u, bc)

path = np.linspace(crack_position - crack_length / 2, crack_position + crack_length / 2, 100)

SIF = 0

for i in range(len(path) - 1):

x1, x2 = path[i], path[i + 1]

SIF += sqrt(2 * pi * mu) * (u(x1) - u(x2)) / (x2 - x1)

sigma = 100 #

a = 0.01 #

K_I = sigma * math.sqrt(math.pi * a)