材料力学之应力分析算法：塑性变形分析：塑性变形分析的实验验证方法.docx

import matplotlib.pyplot as plt

strain = [0.00, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.10, 0.11, 0.12, 0.13, 0.14, 0.15]

stress = [0.00, 10.00, 20.00, 30.00, 40.00, 50.00, 60.00, 70.00, 80.00, 90.00, 100.00, 105.00, 110.0

0, 115.00, 120.00, 125.00]

plt.figure(figsize=(10, 6))

import numpy as np

from scipy.optimize import fsolve

E = 200e9 #

s += p * np.eye(3)

return s

def update_plastic_strain(s, p, ep):

f = np.sqrt(3 / 2 * np.dot(s - p * np.eye(3), s - p * np.eye(3))) - sigma_y - H * ep

return f

def strain_increment(s, p, ep, d):

ep_new = fsolve(update_plastic_strain, ep, args=(s, p))

d_ep = ep_new - ep