提取单元刚度矩阵_第1页
提取单元刚度矩阵_第2页
提取单元刚度矩阵_第3页
全文预览已结束

下载本文档

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

文档简介

本文格式为Word版，下载可任意编辑——提取单元刚度矩阵单元刚度矩阵的提取

刚度矩阵在有限元求解过程中扮演者十分重要的角色，以最小位能原理求解过程为例最终越是转换为含有结构刚度矩阵的能量泛函的取值问题。有限元过程中涉及到三类刚度：单元刚度矩阵，组合结构刚度矩阵和最终求解刚度矩阵。

其中单元刚度矩阵：仅与单元的自身自由度有关，同一编号的单元矩阵的维数是固定。组合结构刚度：矩阵根据求解的初始变量个数决定刚度矩阵的维数，属于单元组装后的初始刚度，维数和整个单元初始变量个数相等。最终求解刚度矩阵：代入边界条件简化后的刚度。以《FiniteElementAnalysis-TheoryandApplicationWithANSYS》中的梁单元例子为例，解释刚度提取过程：

此模型的单元刚度矩阵：（学则beam3梁单元后，该单元包含两节点，每个节点具有三个自由度，因此对应单元刚度矩阵为6*6的方阵）

组合结构刚度矩阵：（该结构含有三个节点，每个节点具有三个原始自由度，因此组合结构刚度矩阵具有9*9阶的形式）

最终求解刚度矩阵：（由于边界条件的存在，该结构中，1,3点的自由度不存在，求解参数中有六个参数已知，因此对最终求解刚度矩阵为三阶方阵）

通过最终的刚度矩阵组成的方程，求解出2节点的位移解，再以这些原始解得出应力，应变，支反力的其他的解。ansys实现过程：

提取思路如下：通过/debug提取单元刚度矩阵，通过filname.full文件提取后两者的矩阵ansys实现过程如下：finish/clear/filname,k,1/prep7N,1N,2,120

N,3,120,-108et,1,beam3mp,ex,1,3.0e7mp,prxy,1,0.3R,1,7.65,204,10E,1,2E,2,3

/debug,-1,,,1,,,,,

/OUTPUT,k,MATRix,finish/soluallsel,allsolve

!生产单个刚度矩阵,查看k.MATRix文件看出单元1,2各次刚度矩阵如下

/AUX2

FILE,'k','full','

HBMAT,'k','txt','',ASCII,STIFF,NOFINISH

!读取k.txt文件，得出初始结构刚度刚度矩阵如下：是个9*9阶矩阵，含有非零参数23个，（以harwell-boeing文件格式生成的，只要该格式图区方式详见站内）结果如下：

!为参与边界条件的刚度矩阵finish/solu

SFBEAM,1,1,PRES,800/12

nsel,s,,,1,3,2/GO

D,all,,,,,,ALLallsel,all

!/debug,-3,,1,,,,,,!/OUTPUT,k4,MATRix,solve/AUX2

FILE,'k','full',''

HBMAT,'k','txt','',A

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

温馨提示

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

最新文档

评论

0/150

提交评论