有限元MATLAB_第1页
有限元MATLAB_第2页
有限元MATLAB_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流有限元MATLAB.精品文档.MATLAB报告Matlab程序求解简要过程如下：（1） 求取单元节点位移提取矩阵T单元节点位移提取矩阵T本质上是置换矩阵群中的一个，结果可将任意杂乱的节点顺序置换成统一的顺序。另一方面其作用是对单元刚度矩阵进行“升维操作”，将单元刚度矩阵统筹到整体刚度矩阵上来，便于对总体节点位移矩阵和支座反力进行求取。本程序分析过程中对单元1的节点提取是按顺序编号1-2-3，对单元2的节点提取是按顺序编号2-3-4。单元1的节点位移提取矩阵如下：单元2的节点位移提取矩阵如下：（2） 求取单元几何矩阵B单元1的节点按编号顺序1-

2、2-3分别进行对几何函数矩阵或算子矩阵的bi逆时针操作，对ci顺时针操作；单元2的节点按编号顺序2-3-4分别进行对几何函数矩阵的bi顺时针操作，对ci逆时针操作.在MATLAB程序中通过mod()取模函数来达到对节点的顺时针或逆时针循环操作。单元1的几何矩阵如下：单元2的几何矩阵如下：（3） 求取应力矩阵S单元应力矩阵满足S=D*B,其中D为弹性矩阵，B为单元几何矩阵各单元的弹性矩阵如下：单元1的应力矩阵如下：单元2的应力矩阵如下：（4） 求取单元刚度矩阵K单元刚度矩阵K满足公式K=B*D*B*t*A,其中t为平面板的厚度，A为单元面积，且单元刚度矩阵为对称矩阵。单元1的刚度矩阵如下：单元2

3、的刚度矩阵如下：（5） 求取总体刚度矩阵sumKK由上述步骤求得的单元刚度矩阵K利用单元虚功原理和刚度方程可导出K*=f，其中为单元节点位移列阵，f为单元等效节点载荷列阵，为了能将各个单元刚度方程统一到一个整体，便需要步骤（1）的单元节点提取矩阵对单元刚度方程进行变换，将两个变换结果联立便得到总体刚度方程，其中也可得到总体刚度矩阵sumKK，且总体刚度矩阵可由sumKK= T*K*T求得。总体刚度矩阵如下：（6） 求取总体节点位移矩阵和支座反力利用上述步骤提到的总体刚度方程sumKK*delta=F，其中delta为总体节点位移矩阵，F为总体等效节点载荷列阵。利用总体边界条件可将总体刚度方程分成两部分，一部分先求取总体节点位移矩阵delta，最后代入总体刚度方程求取总体等效节点载荷列阵F，由此便得到支座反力。总体节点位移矩阵如下：支座反力如：（7） 求取单元应力各单元应力满足=S

人人文库> 全部分类> 教育资料 > 备课教案

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论