应用MATLAB进行有限元分析资料.docxVIP

 PAGE \* MERGEFORMAT 19  PAGE \* MERGEFORMAT 14 应用MATLAB进行有限元分析 姓名：苏红梅 学号班级：土建二班 电话目 录  TOC \o 1-2 \h \z \u  HYPERLINK \l _Toc294599959 一、MATLAB程序说明  PAGEREF _Toc294599959 \h 1  HYPERLINK \l _Toc294599960 （一）M文件源程序相关说明  PAGEREF _Toc294599960 \h 1  HYPERLINK \l _Toc294599961 （二）矩阵及输入量说明说明  PAGEREF _Toc294599961 \h 1  HYPERLINK \l _Toc294599962 二、习题计算  PAGEREF _Toc294599962 \h 1  HYPERLINK \l _Toc294599963 （一）单元及节点编号  PAGEREF _Toc294599963 \h 2  HYPERLINK \l _Toc294599964 （二）使用命令  PAGEREF _Toc294599964 \h 2  HYPERLINK \l _Toc294599965 （三）计算结果  PAGEREF _Toc294599965 \h 4  HYPERLINK \l _Toc294599966 三、源程序  PAGEREF _Toc294599966 \h 13 一、MATLAB程序说明 （一）M文件源程序相关说明 表1 源程序输入量输出量作用TDDAE、n、MTYD计算[T]、 [D]BSKINEk计算[B]、[S]、[k]AssembleelementK单元定位、计算[K]boundUV无化简[K]SLOVEPdelta使用LU分解求解[delta]STRESS无Sigma_e求解单元应力[Sigma_e]（二）矩阵及输入量说明说明 表2 名称含义名称含义AE材料类型m节点总数N单元总数K总体刚度矩阵MTY计算类型：1-平面应力； 2-平面应变UV控制节点约束矩阵：1-可动； 0-不动INE单元坐标矩阵P荷载矩阵T坐标转换矩阵delta节点位移矩阵D单元弹性矩阵Sigma_e单元应力矩阵K单元刚度矩阵二、习题计算 图1 原图 图2 计算图 E=1.0 t=1.0 μ=0.0 γ=0.0 G=1.0 (一)单元及节点编号 节点坐标 表3 节点号坐标节点号坐标1（0,2）10（3,2）2（0,1）11（3,1）3（0,0）12（3,0）4（1,2）13（4,2）5（1,1）14（4,1）6（1,0）15（4,0）7（2,2）16（5,2）8（2,1）17（5,1）9（2,0）18（5,0） 单元节点编号 表4 单元号（i,j,m）单元号（i,j,m）11,5,4118,9,1121,2,5129,12,1132,3,51310,14,1343,6,51410,11,1454,8,71511,12,1464,5,81612,15,1475,6,81713,17,1686,9,81813,14,1797,11,101914,15,17107,8,112015,18,17（二）使用命令 clear AE=[1.0 0.0 1.0 0.0 0.5 0.0 1.0 0.0;1.0 0.0 1.0 0.0 0.5 0.0 1.0 0.0;1.0 0.0 1.0 0.0 0.5 0.0 1.0 0.0;1.0 0.0 1.0 0.0 0.5 0.0 1.0 0.0;1.0 0.0 1.0 0.0 0.5 0.0 1.0 0.0;1.0 0.0 1.0 0.0 0.5 0.0 1.0 0.0;1.0 0.0 1.0 0.0 0.5 0.0 1.0 0.0;1.0 0.0 1.0 0.0 0.5 0.0