函数 [Material_State2,D_crco]=Fixed_Crack_Model(Material,Material_State,e) 输入： -------- 材料：包含材料属性 Material.E（弹性模量）的变量， Material.v（泊松比）、Material.f_t（抗拉强度）、Material.g_f（断裂能除以单元尺寸）和 Material.beta（剪切保持因子） Material_State：包含先前增量或迭代的材料状态变量的历史变量。 它包括 Material_State.s（应力向量）、Material_State.e（应变向量）和 Material.e_cr_l（裂纹应变向量）、Material.e_cr_n_max（最大达到正常裂纹应变）和 Material.theta（垂直于裂纹的角度） 这些历史变量的初始值应该为零，除了 theta 它是

解决了平面应力问题，即边缘受均匀张力的板。 使用等参 Q4 元素对板进行离散化。 预处理是使用标准 FEM 软件完成的。 二阶高斯积分用于获得刚度矩阵。 获得的结果与标准 FEM 软件进行了比较，两者的结果非常吻合。

二维平面应力问题中的弹性材料模型 [s，D] = Elastic_Model（材料，e） 输入： -------- 材料：包含材料属性 Material.E（弹性模量）和 Material.v（泊松比）的变量e: (3*1) 应变向量 输出： ----------- s: (3*1) 应力向量D: (3*3) 本构矩阵

函数 [Material_State2,D_crco]=Isotropic_Damage_Model(Material,Material_State,e) 输入： -------- 材料：包含材料属性 Material.E（弹性模量）、Material.v（泊松比）、Material.f_t（拉伸强度）和 Material.g_f（断裂能除以单元尺寸）的变量 Material_State：包含先前增量或迭代的材料状态变量的历史变量。 它包括 Material_State.s（应力向量）、Material_State.e（应变向量）、Material_State.e_eq（等效应变）、Material_State.k（硬化参数）和 Material_State.d（损坏） 这些历史变量的初始值应该为零 重要提示：对于这个模型，这个变量应该从材料的最后一个收敛状态（即前一个增量的结束而不

根据有限元中三角形常应变单元的理论知识，利用Matlab工具，编制了本求解器，用以求解平面应力问题，并将所求结果与Ansys软件所求结果及弹性理学中的理论值进行了比较，效果较好。