matlab力学大作业答案,计算力学大作业报告.docx

计算力学大作业报告

计算力学基础大作业报告航21班 苏浩 尹思凡一、求解流程说明本次作业的求解过程主要分为三个步骤：1. 前处理网格划分——采用ABAQUS软件进行画网格，并输出前处理文件。形成前处理文件——利用matlab编程，形成inmesh文件。2. 有限元计算采用上机课上提供的FEATP程序计算，程序略有修改，使其能够计算单元和节点数目更多的算例。本次大作业所有的图中str1~3,分别代表x(r)、y(theta)方向的正应力，剪力。3. 后处理采用Tecplot软件画图，得到位移场和应力场的云图。二、所编程序及功能简介程序名称题目功能abs2pt.m1、2、3、4(两种情况)、7将ABAQUS的前处理文件进行处理，得到p、t矩阵semibw.m1、2、3、4(两种情况)、7优化总刚矩阵并计算半带宽(已修改为适合3、4、6、8节点的情况)viewmsh.m1、2、3、4(两种情况)、7显示网格(主要目的是检查边界条件是否施加正确，调bug用)readdata.m1、2、3、4直接输出featp程序可以读取的节点信息和单元信息，可以直接复制到inmesh中。readdata_7.m7直接输出featp程序可以读取的节点信息和单元信息，可以直接复制到inmesh中。bianjie_n.m(n为对应题号)1、2、3、4(两种情况4_1平面应变和4_2轴对称)、7(每道题的程序都不同)直接输出featp程序可以读取的位移边界条件条件和力边界条件，可以直接复制到inmesh中。xy2rtheta_n.m (n为对应题号)3、4(平面应变)将笛卡尔坐标系中求得的位移和应力转换到极坐标系中，可以用tecplot画图work8.m8读取abaqus的前处理文档，并计算出半带宽，直接输出featp程序可以读取的节点信息、单元信息、位移边界条件条件和力边界条件，可以直接复制到inmesh中。work5.m5得到5题域内温度场,可以直接在tecplot中画图output_stress.m1、2、4、7将不同情况算得的应力进行比较，以确定单元收敛性。get_road3.m3沿不同路径对应力进行比较，观察应力沿截面的分布。get_road4.m4沿不同路径对应力进行比较，观察应力沿截面的分布。get_road4_2.m4对第四题轴对称问题取得路径上的lilunjie4.m4得到第四题的理论解二、具体算例1.悬臂梁问题1.1 问题描述图1所示悬臂梁，一端固支，承受集中力 P=1000N，梁的长度为 100mm，高度为10mm，厚度为1mm，材料弹性模量为2×1011Pa，泊松比为0.32，假定在变形过程中界面始终保持完整，利用有限单元法计算位移场分布，并利用求得的位移场计算结点应力，讨论与理论解的误差。分别使用不同类型不同阶次单元进行计算，比较结果精确度，单元收敛性等。1.2 解析计算采用材料力学方法进行计算，建立主轴坐标系，根据材料力学的推导，有如下公式： (1.1)其中是z轴正方向的弯矩，是方向为y轴正方向的剪力，是梁对z方向的惯性矩，是横截面对z轴的静面矩，t为厚度。由题目中的数据，可以得到： (1.2)将(1.2)中的公式带入(1.1)中，可以得到: (1.3)材料力学中的挠度公式和转角公式： (1.4)而x方向的位移： (1.5)将(1.2)中的计算结果和材料参数带入，并对式(1.4)进行积分，得到： (1.6)位移边界条件： 得到位移场： (1.7)1.3 问题讨论(1)位移边界条件：固定位移：左边界——中点 u=0 v=0 其余点 u=0(2)载荷条件：右上角点1.4 用不同类型的单元求解位移场和应力场1.4.1 三角形三节点单元(1)网格划分采用三种不同疏密度的网格计算：(2)计算结果这里仅选取网格密度最大的第三种情况的结果位移场应力场 1.4.2 四边形四节点单元(1)网格划分采用四种不同疏密度的网格计算：(2)计算结果位移场应力场1.4.3 四边形八节点单元(1)网格划分采用三种疏密度不同的网格计算(2)计算结果位移场应力场1.5 单元精度与单元收敛性1)考虑平面应力问题，理论上v的挠度位移： (1.8)代入数据：P=1000N，L=100mm E=2×105MPa G=E/(2(1+))=75757.6MPa得到理论挠度：以(100,5)处的节点位移为例网格密度5*52.5*2.51*10.5*0.5三节点三角形10.966916.611819.496119.9965四节点四边形17.946819.562320.071520.1480八节点四边形20．165820．171120.17332)再对应力场进行检验以(50,10)处的正应力值作为收敛对象，理论5*52.5*2.51*10.5*0.5三节点三角形0.109718E+040.