弹性力学 上级报告供参习.docVIP

弹性力学上机报告 姓名： 董永飞 学号： 311107090113 专业： 工程力学 班级： 工程力学2011-1班 时间： 2013-2014年度第一学期 河南理工大学土木工程学院 实验一 实验名称：主应力的数值计算 实验目的： 1 进一步了解主应力、主应变、主方向。 2 掌握用MATLAB计算主应力、主应变、主方向的方法 基本理论： 主应力 在某点的应力状态中，如果在某个斜面上的切应力等于零，则该斜面上的正应力称为该点的一个主应力。 主应变 在某点的应变状态中，如果在某个斜面上的切应变等于零，则该斜面上的线应变称为该点的一个主应变。 主方向 对应于主应力或者主应变的斜面的法线方向称为应力主方向或者应变主方向。对于同一点，应力主方向与应变主方向是一致的。 二向状态下的主应力： 主方向： 三向状态下的主应力： 主应力满足方向： 式中： 其主方向需进一步讨论得到。 通过研究发现，如果把应力分量写成应力矩阵： 则主应力及其相对应的主方向就是应力矩阵的特征值及其对应的特征向量。 MATLAB中提供了计算矩阵特征值及特征向量的专门函数：eig()，该函数使用格式如下： [V,D] = EIG(X) 其中D和V分别为X矩阵的特征值及特征向量。学生需要再编程将特征向量换算成主方向的方向余弦。 力学模型： 程序命令： 运行结果： 结论与误差分析：以上计算结果告诉我们：该单元体