《有限元法》实验报告
专业班级
姓名
力学(实验)
田诗豪
1601
学号
10
提交日期
实验编号实验一实验二实验三总分
得分
实 验 一(30 分)
一、实验内容
编写一个计算平面
3 结点三角形单元的应变矩阵、应力矩阵和单元刚度矩阵的
MATLAB
函数文件
[B3,S3,K3]
= ele_mat_tri3(xy3,mat)
,其中:输入变量
xy3
为结点坐标数组,
mat
为材料参数矩阵;输出变量B3 为应变矩阵, S3 为应力矩阵,K3 为单元刚度矩阵。 (要求给
出 3 个不同算例进行验证,并绘制出单元形状和结点号)
二、程序代码
通用函数
function
[B3,S3,K3] = ele_mat_tri3(xy3,mat)
%生成平面3 结点三角形单元的应变矩阵、应力矩阵和单元刚度矩阵的功能函数
%*********变量说明 ****************
%xy3------------------结点坐标数组
%mat------------------材料参数矩阵(弹性模量,泊松比,壁厚)
%B3-------------------应变矩阵
%S3-------------------应力矩阵
%K3-------------------单元刚度矩阵
%********************