弹性力学上报告.docVIP

弹性力学上机报告 姓名： 董永飞 学号： 311107090113 专业： 工程力学 班级： 工程力学2011-1班 时间： 2013-2014年度第一学期 河南理工大学土木工程学院 实验一 实验名称：主应力的数值计算 实验目的： 1 进一步了解主应力、主应变、主方向。 2 掌握用MATLAB计算主应力、主应变、主方向的方法 基本理论： 主应力 在某点的应力状态中，如果在某个斜面上的切应力等于零，则该斜面上的正应力称为该点的一个主应力。 主应变 在某点的应变状态中，如果在某个斜面上的切应变等于零，则该斜面上的线应变称为该点的一个主应变。 主方向 对应于主应力或者主应变的斜面的法线方向称为应力主方向或者应变主方向。对于同一点，应力主方向与应变主方向是一致的。 二向状态下的主应力： 主方向： 三向状态下的主应力： 主应力满足方向： 式中： 其主方向需进一步讨论得到。 通过研究发现，如果把应力分量写成应力矩阵： 则主应力及其相对应的主方向就是应力矩阵的特征值及其对应的特征向量。 MATLAB中提供了计算矩阵特征值及特征向量的专门函数：eig()，该函数使用格式如下： [V,D] = EIG(X) 其中D和V分别为X矩阵的特征值及特征向量。学生需要再编程将特征向量换算成主方向的方向余弦。 力学模型： 程序命令： 运行结果： 结论与误差分析：以上计算结果告诉我们：该单元体的两个主应力从大到小为： σ1=170MPa，σ2=70Mpa 对应的主方向分别为：-71.5651度， 18.4349度 若按方向角描述：两个主方向的方向角分别是： α1= 108.4349 度， β1=18.4349度 α2=161.5651度， β2=108.4349度； 程序命令： 运行结果： 结论与误差分析：以上计算结果告诉我们：该单元体的三个主应力从大到小为： σ1=106.7216MPa，σ2=43.7072Mpa，σ3=-0.4288MPa 三个主应力对应的三个主方向的方向角分别是： α1= 36.6001度， β1=65.2562度，γ1=64.8741度； α2=100.1432度， β2=120.9251度，γ2=32.9054度； α3=124.7240度，β3=41.5137度，γ3=70.1920度。 实验二 实验名称： 弹性力学公式的辅助推算 实验目的： 1 初步了解MATLAB强大的符号推演功能。 2 用MATLAB的符号计算功能，导出弹性力学中的方程。 例1：导出用应变表示应力的物理方程： 程序命令： 运行结果： 例2极坐标解答中，若设定应力函数为，试用MATLAB导出各应力分量。 已知应力分量表示式：， 相容方程：。 程序命令： 运行结果： 实验三 实验名称：应力变分法在平面问题中应用 实验目的： 1 掌握平面问题的应力变分方法的原理及计算过程。 2 用MATLAB实现平面问题的应力变分方法的计算过程，并讨论形函数的收敛性。 基本理论：平面问题的应变余能可以表示为： 引入应力函数后： 应力变分方程是： 由于面力是给定的，所以上式右端等于零。即。 如果应力函数的构成是通过形函数的系数构成，则余能的变分方程就写成如下的偏微分形式： 力学模型： 例题（P284，11-2）：正方形薄板，边长为，如图所示在左右两边受有按抛物线分布的拉力，即。试用应力变分法按如下的应力函数求解： 。 程序命令： syms a x y FI sx sy sxy q a1 a2 a3 a4 FI=q*y^4/12/a^2+q*a^2*(1-x^2/a^2)^2*(1-y^2/a^2)^2*(a1); s1=diff(FI,2,y); s2=diff(FI,2,x); s3=diff(diff(FI,x),y); s4=s1^2+s2^2+2*s3^2; VC=int(s4,x,-a,a) ; VC=int(VC,y,-a,a); s1=diff(VC,a1) ; A1=solve(s1,a1); sx=diff(FI,2,y) ; sx=subs(sx,a1,A1) ; sx=subs(sx,x,0); sx=vpa(sx,2); sx1=simple(sx) %以下计算3阶 clear