有限元上机实验报告.docxVIP

有限元法基础及应用上机报告 南京理工大学 2015年12月上机实验一1 实验题目设计一个采用减缩积分线性四边形等参元的有限元模型，通过数值试验来研究网格密度、位移约束条件与总刚度矩阵奇异性、沙漏扩展、求解精度的关系，并验证采用减缩积分时保证总刚度矩阵非奇异的必要条件。总结出你的研究结论，撰写实验报告。2 实验目的通过实验来研究减缩积分方案中网格密度和位移约束条件对总体刚度矩阵奇异性和求解精度的影响，以此加深对有限元减缩积分的理解，和对减缩积分中保证总体刚度矩阵非奇异性的认识。3建模概述先保持位移约束条件不变，研究网格密度对总体刚度矩阵奇异性和求解精度的影响，并验证采用减缩积分时保证总刚度矩阵非奇异的必要条件。如下图1所示，建立一个简支和链杆的约束条件，然后不断增加网格密度，通过ABAQUS来计算位移和应力的变化规律。 简支（两个约束） 链杆（一个约束）P 积分点（3个独立关系式） 节点（两个自由度）4 计算结果分析讨论与结论1）1*1单元四边形减缩积分实验 载荷 布种/单元 P应力云图2）2*1单元四边形减缩积分实验 载荷 单元应力云图3）4*4单元四边形减缩积分实验 载荷 布种 单元应力云图结果分析项目自由度数比较独立关系式结果1)4*2-3=51*3=3奇异2)6*2-3=92*3=6奇异3)25*2-3=4716*3=48非奇异5 实验体会与小结 单元刚度矩阵的特征：对称性奇异性主元恒正平面图形相似、弹性矩阵D、厚度t相同的单元，相同的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两列，其位置与结点位置对应。整体刚度矩阵的特征：对称性奇异性主元恒正稀疏性非零元素呈带状分布。[K]的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足力和力矩的平衡。为消除[K]的奇异性，需要引入边界条件，至少需给出能限制刚体位移的约束条件。对于一个给定形式的单元，如果采用精确积分，则插值函数中所有项次在|J|=常数的条件下能被精确积分，并能保证刚度矩阵的非奇异性。如果采用减缩积分，因为插值函数中只有完全多项式的项次能被精确积分，因此需要进行刚度矩阵非奇异必要条件的检查。上机实验二1 实验题目图示一个简支梁平面应力模型。梁截面为矩形，高度h=160mm，长度L=1000mm，厚度t=10mm。上边承受均布压力q =1N/mm2，材料E=206GPa，μ=0.29。X方向正应力弹性力学理论解为：分别应用3节点三角形单元、4节点线性等参元（完全积分、减缩积分、非协调模式）、8节点二次等参元完全积分进行下列数值实验：1）用较粗单元网格求解梁中部应力分量σx的最大值和上下边法向应力分量，并对各单元计算精度进行比较分析；2）对粗网格下梁中部最大位移进行对比和分析；3）通过网格加密对比试验3节点三角形单元和8节点二次等参元的收敛速度。总结出研究结论，撰写实验报告。2 实验目的（1）通过ABAQUS软件，用有限元法分析整个梁上的σ_x的分布规律，讨论σ_y的有限元解与材料力学解的区别；（2）用有限元法求梁顶边和底边中点正应力σx的最大值；（3）比较3节点三角形单元、4节点线性等参元（完全积分、减缩积分、非协调模式）、8节点二次等参元完全积分求得的参数与理论解的接近程度；（4）逐步加密单元网格，把有限元法求得的σx值与理论值进行对比，考察解的收敛性；（5）针对以上力学模型，对比分析3节点三角形平面单元和8节点四边形平面单元的求解精度和收敛性；（6）绘制σx的误差——计算次数曲线，并进行分析说明。3 建模概述（1）启用ABAQUS/CAE程序；（2）创建部件（Module：Part），选择2D Planar，Approximate size=2000，绘制长度为1000，高度为200的矩形；（3）创建材料和截面属性（Module：Property），弹性模量为E=210000 MPa，泊松比为0.3；（4）定义装配件（Module：Assembly），选择Independent；（5）设置分析步（Module：Step）；（6）定义边界条件和载荷（Module：Load），将左右两边界在中点分成两段，在梁的上侧面施加1Mpa的载荷，选取左边界中点设置U1=0和U2=0，选取右边界中点设置U2=0；（7）划分网格（Module：Mesh），全局尺寸设为50，采用3节点线性平面应力三角形单元CPS3划分网格；（8）提交分析作业（Module：Job）；（9）后处理（Module：Visualization）： ①显示的应力云图； ②查询底边中点的最大值，与理论最大应力值比较考察其精度； ③绘制底边上各点的应力的曲线；（10）细化网格验证收敛性，重设单元尺寸为20和10，给出和的应力云图，考察收敛于理论解的程度；（11）高阶单元分析，将单元类