有限元法15年作业综述.docx

有限元法作业 班级： 学号： 姓名： 作业提交至：femhomework@126.com, 标题：班级+学号+姓名+有限元作业（不含+） 1、已知弹性模量E=29.5x104N/mm2，A=100mm2，求结点位移、单元应力及支座反力。 解： 1）计算程序 1.主程序 prep; kw=fun_kw(XY,ELB,EA); P=fun_load(XY,NPQ); [kw,P]=fun_restrain1(kw,P,SU); U=fun_mguass(kw,P); stress=fun_interf(XY,ELB,U,EA); posp; 2.前处理程序 clear;clc; XY=[1,0,0 2,400e-3,0 3,400e-3,300e-3 4,0,300e-3]; EA=[295e9,1e-4]; ELB=[1,1,2,1 2,2,3,1 3,1,3,1 4,3,4,1]; NPQ=[3,20000 6,-25000]; SU=[1,0 2,0 4,0 7,0 8,0]; 3.后处理程序 [N,m]=size(XY); U=reshape(U,2,N); U=U; U1=1:N; U=[U1,U] [MB,m]=size(ELB); stress1=1:MB; stress=[stre