matlab有限单元法计算桁架算例代码
matlab有限单元法计算桁架算例代码
E=1.0e7;
A=1.5;
node_number=6;
element_number=8;
nc=[0,0;0,40;40,0;40,40;80,0;80,40];
en=[1,3;1,4;2,4;3,4;3,5;4,5;4,6;5,6];
ed(1:node_number,1:2)=1;
constraint=[1,1;1,2;2,1;2,2];
for loopi=1:length(constraint);
ed(constraint(loopi,1),constraint(loopi,2))=0;
end
dof=0;
for loopi=1:node_number
for loopj=1:2
dof=dof+1;
ed (loopi,loopj)=dof
end
end
ek=EA[1,0,-1,0;0,0,0,0;-1,0,1,0;0,0,0,0];
k(1:dof,1:dof)=0;
thea(1:8)=0;
el(1:8)=0;
e2s(1:4)=0;
for loopi=1:element_number
for zi=1:2
e2s((zi-1)*2+1)=ed(en(loopi,zi),1);
e2s((zi-1)*2+2)=ed(en(loopi,zi),2);
end
el(loopi)=sqrt((nc(en(loopi,1),1)-nc(en(loopi,2),1))2+(nc(en(loopi,1),2)-nc(en(loopi,2),2))2);
theta(loopi)=asin((nc(en(loopi,1),2)-nc(en(loopi,2),2))/el(loopi)) ;
lmd=[cos(theta(loopi)) sin(theta(loopi));-sin(theta(loopi)) cos(theta(loopi))];
t=[lmd zeros(2);zeros(2) lmd] ;
dk=t’ekt/el(loopi);
for jx=1:4
for jy=1:4
if (e2s(jx)*e2s(jy)~=0);
k(e2s(jx),e2s(jy))=k(e2s(jx),e2s(jy))+dk(jx,jy);
end
end
end
end
num=[5 6 7 8 9 10 11 12]
kk=k(num,num)
F=[0 1;-2000 2;0 3;0 4;2000 5;0 6;4000 7;6000 8];
F1=zeros(8,1);
for i=1:8
F1(i)=F(i,1)
end
v=inv(kk)F1;
num1=[1 2 3 4];
kk1=k(num1,1:12);
u1=zeros(12,1);
for i=5:12
u1(i)=v(i-4)
end
F2=kk1u1
strain=zeros(8,1)
for loopi=1:element_number
uc1=zeros(2,1)
uc2=zeros(2,1)
el(loopi)=sqrt((nc(en(loopi,1),1)-nc(en(loopi,2),1))2+(nc(en(loopi,1),2)-nc(en(loopi,2),2))2);
theta(loopi)=asin((nc(en(loopi,1),2)-nc(en(loopi,2),2))/el(loopi)) ;
lmd1=[cos(theta(loopi)),sin(theta(loopi))];
uc1(1)= u1(2en(loopi,1)-1)
uc1(2)= u1(2 en(loopi,1))
uc2(1)= u1(2en(loopi,2)-1)
uc2(2)= u1(2 en(loopi,2))
ub=lmd1uc1;
ue=lmd1uc2;
strain (loopi)=(ub-ue)/el(loopi)
end
stress=E*strain
matlab有限单元法计算桁架算例代码相关推荐
- 用Matlab求解一维非稳态周期性导热问题(有限单元法+隐式离散+高斯赛德尔迭代法)
本次求解不一定对,请先看最后说明 一.问题描述与分析 本次问题条件如下: 计算模拟如下一维常物性无内热源非稳态导热的温度场,以及内外壁面的热流密度,并进行温度场和热流的特点分析,相关参数如下. 室内温 ...
- 二维有限体积 matlab,二维有限体积法计算热传导及源码.pdf
二维有限体积法计算热传导及源码 //#include "stdafx.h" #include #include #include #include #include using n ...
- 二阶偏微分方程组 龙格库塔法_有限单元法(Finite Element Method)实现声波方程模拟(Part 2)...
2.1 前言 承接上一篇文章,前面我们已经介绍了一维声波方程有限元求解: 蓝不是蓝:有限单元法(Finite Element Method)实现声波方程模拟(Part 1)zhuanlan.zhih ...
- 分布式电源接入对配电网的影响matlab程序(IEEE9节点系统算例)
分布式电源接入对配电网的影响matlab程序(IEEE9节点系统算例) 摘 要:分布式电源的接入使得配电系统从放射状无源网络变为分布有中小型电源的有源网络.带来了使单向流动的电流方向具有了不确定性等等 ...
- 岩土工程渗流问题之有限单元法:理论、模块化编程实现、开源程序手把手实操技术
有限单元法在岩土工程问题中应用非常广泛,很多商业软件如Plaxis/Abaqus/Comsol等都采用有限单元解法.尽管各类商业软件使用方便,但其使用对用户来说往往是一个"黑箱子" ...
- Plaxis Python 命令流自动化处理、岩土工程渗流问题之有限单元法
目录 岩土工程渗流问题之有限单元法:理论.模块化编程实现.开源程序手把手实操应用 基于python命令流及代码的Plaxis自动化建模与典型案例实践应用 岩土工程渗流问题之有限单元法:理论.模块化编程 ...
- matlab求刚度,matlab直接刚度法计算结构频率
matlab直接刚度法计算结构频率 syms E I K r l x cja1=-E*I*K*r^3*(cos(x)*sinh(x)+sin(x)*cosh(x))jc1=-E*I*K*r^2*sin ...
- matlab用jc法计算可靠度,用MATLAB实现JC法计算结构可靠度程序
软件设计开发本栏目责任编辑 :谢媛媛 ComputerKnowledgeand Technology 电脑知识与技术第 5 卷第 29 期 (2009 年 10 月) 用 MATLAB 实现 JC 法 ...
- matlab多自由度弹簧响应,用matlab编程实现法计算多自由度体系的动力响应..doc
用matlab编程实现法计算多自由度体系的动力响应. 用matlab编程实现法计算多自由度体系的动力响应 姓名:王 辉 学号 班级:硕1205班 专业:桥梁与隧道工程 2012年11月27日 用mat ...
- matlab用jc法计算可靠度,基于MATLAB的截尾分布下JC法计算可靠度
摘 要: 在水工结构可靠度分析中,随机变量的分布形式常因几何尺寸.物理环境等条件限制,传统JC法已经不适用,因此需要对部分变量进行截尾分布处理.在此借助MATLAB丰富的函数资源,编制出截尾分布处理后 ...
最新文章
- 十分钟上手2022latex安装与入门
- Contest Leaderboard(分组、聚合、连接)
- “中文四六级”考试来了!《国际中文教育中文水平等级标准》发布
- 后端Coder如何做好代码设计?
- VRRP协议简介与实例配置
- 利用 Angular Directive 和 @HostBinding 实现输入文本框随着键盘输入自动变色效果
- ASP.NET Core 6.0对热重载的支持
- CString常用方法
- 下次激活策略10_(好文重读)低风险之可转债4大卖出策略
- POJ 2289 Jamie's Contact Groups 【二分】+【多重匹配】(模板题)
- 微信小程序转发,发送朋友圈(uni-app)
- oracle startup open ora 03113,oracle宕机,startup报错ora03113
- 风口来了第二期——电子科学与技术专业现状和前景介绍分享
- Linux下套接字详解(八)----select模式下服务器(非阻塞,单进程+多进程+多线程)
- Goby - 新一代安全工具
- corrupted double-linked list: 0x086a6c50(内存问题)
- 神经架构搜索(NAS)2020最新综述:挑战与解决方案
- 帆软设计器连接远程服务器,帆软设计器远程连接数据库问题
- 招聘牛人的五步骤蔡文贤_招聘技术人才的5个步骤
- 域名Goldrush第6部分-追溯到未来?