六面体单元matlab后处理,《有限元基础教程》_【MATLAB算例】基于节点六面体单元的空间块体分析(HexahedralDNode).doc...
《有限元基础教程》_【MATLAB算例】基于节点六面体单元的空间块体分析(HexahedralDNode)
【MATLAB算例】4.8.2(1) 基于8节点面体单元。基于MATLAB平台,用一个空间8节点六面体单元计算各个节点位移、支座反力以及单元的应力。
(a) 问题描述 (b) 有限元分析模型
图4-23 右端部受集中力作用的空间块体
解答:对该问题进行有限元分析的过程如下。
(1)结构的离散化与编号
将结构离散为一个8节点六面体单元,节点编号如图4-23(b)所示,节点的几何坐标见表4-13。
表4-13 节点的坐标
节点节点坐标/mxyz10.20020.20.80300.80400050.200.660.20.80.6700.80.68000.6节点位移列阵
(4-194)
总的节点载荷列阵
(4-195)
其中,节点外载;支反力为,,,,,,,,,,,;其余节点载荷分量为零。
(2)计算单元的刚度矩阵(以国际标准单位)
首先在MATLAB环境下,输入弹性模量E和泊松比NU,然后针对题中单元节点坐标,调用函数Hexahedral3D8Node_Stiffness,就可以得到单元的刚度矩阵k1 (24×24)。
>>E=1.0e10;
>>NU=0.25;
>>lx=0.2;
>>ly=0.8;
>>lz=0.6;
>>k1=Hexahedral3D8Node_Stiffness(E,NU,lx,0,0,lx,ly,0,0,ly,0,0,0,0,lx,0,lz,lx,ly,lz,0,ly,lz,0,0,lz);
(3) 建立整体刚度方程
由于该结构共有8个节点,则总共的自由度数为24,因此,结构总的刚度矩阵为K K(24×24),先对KK清零,然后调用函数Hexahedral3D8Node_Assembly进行刚度矩阵的组装。由于本题中只用了一个单元,因此总体刚度矩阵KK与单元刚度k1相同,此处不再列出,调用函数的过程如下:
>>KK=zeros(24,24);
>>KK=Hexahedral3D8Node_Assembly(KK,k1,1,2,3,4,5,6,7,8);
(4) 边界条件的处理及刚度方程求解
由图4-23(b)可以看出,节点1,4,5和8的3个方向的位移将为零,即。因此,将针对节点2,3,6和7的位移进行求解,这4个节点的位移将对应KK矩阵中的4~9行、4~9列,4~9行、16~21列,16~21行、4~9列,以及16~21行、16~21列,则需从KK(24×24)中提取相应行和列的数据,置给k,然后生成对应的载荷列阵p,再采用高斯消去法进行求解,即MATLAB中的反斜线符号“\”求解。
>> k=[KK(4:9,4:9),KK(4:9,16:21);KK(16:21,4:9),KK(16:21,16:21)];
>> p=[0;0;0;0;0;0;0;0;-1e5;0;0;-1e5];
>> u=k\p
u = 1.0e-003 *
0.0223 -0.2769 -0.6728 -0.0223 -0.2769 -0.6728 [将列排成行排量 [] ]
-0.0129 0.3108 -0.7774 0.0129 0.3108 -0.7774 [将列排成行排量 [] ]
由此可以看出,所求得的结果为(单位为m):
(5)支反力的计算
由方程可知,在得到整个结构的节点位移后,由原整体刚度方程就可以计算出对应的支反力。先将上面得到的位移结果与位移边界条件的节点位移进行组合(注意位置关系),可以得到整体的位移列阵U(24×1),再代回原整体刚度方程,计算出所有的节点力P(24×1),按式的对应关系就可以找到对应的支反力。
>> U=[0;0;0;u(1:6);0;0;0;0;0;0;u(7:12);0;0;0];
>> P=KK*U
P = 1.0e+005 *
-0.2509 1.3333 0.6938 -0.0000
六面体单元matlab后处理,《有限元基础教程》_【MATLAB算例】基于节点六面体单元的空间块体分析(HexahedralDNode).doc...相关推荐
- matlab patch 六面体,《有限元基础教程》_【MATLAB算例】4.8.2(1) 基于8节点六面体单元的空间块体分析(Hexahedral3D8Node)...
[MATLAB 算例]4.8.2(1) 基于8节点六面体单元的空间块体分析(Hexahedral3D8Node) 如图4-23所示的一个空间块体,在右端部受两个集中力F 作用,其中的参数为: 1051 ...
- householder变换qr分解matlab_【基础教程】Matlab实现傅里叶变换
傅立叶变换 傅立叶变换是一种常见的分析方法,傅立叶变换将满足一定条件的函数表示为一些函数的加权和(或者积分).可以分为四个类别: 1. 非周期连续性信号 对应于傅里叶变换,频域连续非周期 2. 周期性 ...
- 黑马程序员全套Java教程_Java基础教程_异常(含扩展)(二十三)
黑马程序员全套Java教程_Java基础教程_异常(含扩展)(二十三) 1.1 异常概述与异常体系结构 1.2 JVM遇到异常时的默认处理方案 1.3 异常处理 1.4 异常处理之try--catch ...
- 【Axure基础教程】第19章 树节点
Axure10基础教程系列,适合小白打基础,基础扎实的同学可以在公众号查看其它案例教程. 本章就主要讲解树节点展开折叠的交互 [Axure基础教程]第19章 树节点
- 非线性有限元:基本理论与算法及基于Python、Fortran程序实现与案例分析实践技术
有限单元法在岩土工程问题中应用非常广泛,很多商业软件如Plaxis/Abaqus/Comsol等都采用有限单元解法.尽管各类商业软件使用方便,但其使用对用户来说往往是一个"黑箱子" ...
- 非线性有限元:基本理论与算法及基于Python、Fortran程序实现与案例分析
非线性有限元:基本理论与算法及基于Python.Fortran程序实现与案例分析 (qq.com) 有限单元法在岩土工程问题中应用非常广泛,很多商业软件如Plaxis/Abaqus/Comsol等都采 ...
- 岩土工程--非线性有限元:基本理论与算法及基于Python、Fortran程序实现与案例分析
非线性有限元:基本理论与算法及基于Python.Fortran程序实现与案例分析实践技术 有限单元法在岩土工程问题中应用非常广泛,很多商业软件如Plaxis/Abaqus/Comsol等都采用有限单元 ...
- matlab 简单算例,(简单算例)基于Matlab的电力系统潮流编程计算.pdf
(简单算例)基于Matlab的电力系统潮流编程计算 基于Matlab的电力系统潮流编程计算 口黄扬威吴喜春郭志峰张斯翔 (三峡大学电气与新能源学院湖北·宜昌443002) 摘要:通过介绍电力系统的实际 ...
- matlab中右三角形方向,《有限元基础教程》_【MATLAB算例】4.7.1(2) 基于3节点三角形单元的矩形薄板分析(Triangle2D3Node)...
[MATLAB 算例]4.7.1(2) 基于3节点三角形单元的矩形薄板分析(T riangle2D3Node) 如图4-20所示为一矩形薄平板,在右端部受集中力100 000F N =作用,材料常数为 ...
- 五杆桁架matlab有限元分析,《有限元基础教程》_【MATLAB算例】3.2.5(2)__四杆桁架结构的有限元分析(Bar2D2Node)...
[MATLAB算例]3.2.5(2)四杆桁架结构的有限元分析(Bar2D2Node) 如图3-8所示的结构,各个杆的弹性模量和横截面积都为42 E N mm =?, 29.510/ 2 =.试基于MA ...
最新文章
- android入门学习一 基本概念
- What is ORM ?
- 揭开并发包底层AQS的神秘面纱
- LocalDB 静默安装
- 聊一聊单机、集中式、分布式和云原生存储
- C# 站点IP访问频率限制 针对单个站点
- php 字符串的比较大小,php怎么比较两个字符串的大小
- android开发卡死代码,Android Studio编译卡死(示例代码)
- CentOS7 修改静态IP
- 华为专利全球第一:哪里跌倒,哪里爬起!
- android错误详细教程四
- 洛谷 P2383 狗哥玩木棒
- mybatis 多租户saas_SaaS 微服务脚手架
- 自定义RPC通讯框架,实现dubbo远程RPC服务治理功能
- 论文笔记:微表情识别综述1
- 高等数学(第七版)同济大学 习题2-3 个人解答
- java servlet面试题_JavaWeb之JSP和servlet面试题附答案
- 中控考勤仪IFace302多线程操作时无法订阅事件
- 第三章 基本粒子
- Android手机截图怎么做,手机截屏怎么弄,教您手机截图方法