函数 [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（损坏） 这些历史变量的初始值应该为零 重要提示：对于这个模型，这个变量应该从材料的最后一个收敛状态（即前一个增量的结束而不

用谱随机有限元方法分析了弹性模量不确定并承受确定性分布载荷的二维平面应力实体。 此参考示例在 Ghanem 和 Spanos 的第 (5.3) 节：“随机有限元：谱方法”中进行了描述。 主函数“main_ssfem_plate.m”包含此问题的逐步解决方案。 为了更好地理解代码，编写了一些对方程的引用和有用的注释。 程序估计 Karhunen-Loève 展开式（变量 'KLterms'）和多项式混沌的给定阶数（变量 'p_order'）、均值和标准差位移响应场中的数项； 还有在给定的实体自由度下位移的相应 PDF 和 CDF。 文件夹“result_MCS”包含一个 .txt 文件，其中包含用于比较目的的蒙特卡罗模拟方法获得的解决方案。 该方法的完整描述在： http://www.bdigital.unal.edu.co/51037/ 值得一提的是，可以在FERUM软件http://

软件还不是很完善，希望对大家的学习有所帮助

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

有限元平面应力问题（应力集中），平面单元，采用c++编程求解。

包含一个桁架机构算例和平面应力算例，均较为基础，代码后大多有注释，可供初学者参考