嘿,亲!知识可是无价之宝呢,但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下,绝对物超所值哦!如有下载和支付问题,请联系我们QQ(微信同号):813200300
本次赞助数额为: 4 元微信扫码支付:4 元
请留下您的邮箱,我们将在2小时内将文件发到您的邮箱
平板结构中有一个洞,已知条件为:E=70Mpa,泊松比u=0.25,板材厚度t=0.02,P=20KN,试进行该问题的有限元分析
clc;
clear;
E=70E6;%弹性模量
pu=0.25;%泊松比
t=0.02;%厚度
ly=0.9;%y方向长度
lx=0.9;%x方向长度
ls=30;%x y方向划分单元的数目
o=ly/ls;%单元的长度
%划分网格数为800个单元
jds=880;%结点总数
zyd=880*2;%总自由度数
%应力
%%单元刚度矩阵的计算及组装
K=zeros(zyd,zyd);
for m=1:(11-1)
for i=1:(31-1)
a=[o*(i-1),o*(m-1)];
b=[o*(i),o*(m-1)];
c=[o*(i),o*(m)];
d=[o*(i-1),o*(m)];
B=[a,b,c,d];
k=SBGJ(E,pu,t,B(1),B(2),B(3),B(4),B(5),B(6),B(7),B(8),1);
C=[i 31*(m-1),i 1 31*(m-1),31 i 1 31*(m-1),31 i 31*(m-1)];
K=QJGDJZ(K,k,C(1),C(2),C(3),C(4));
end
end
for i0=(1:11-1)
a=[o*(i0-1),0.3];
b=[o*(i0),0.3];
c=[o*(i0),0.3 o];
d=[o*(i0-1),0.3 o];
B=[a,b,c,d];
k=SBGJ(E,pu,t,B(1),B(2),B(3),B(4),B(5),B(6),B(7),B(8),1);
C=[310 i0,310 i0 1,310 31 i0 1,310 31 i0];
K=QJGDJZ(K,k,C(1),C(2),C(3),C(4));
end
【文件目录】
信号处理作业
├── five.m
├── four.m
├── one.m
├── seven.m
├── six.m
├── three.m
├── two.m
├── 信号处理作业1.docx
├── 信号处理作业2.docx
├── 信号处理作业3.docx
├── 信号处理作业4.docx
├── 信号处理作业5.docx
├── 信号处理作业6.docx
├── 信号处理作业7.docx
├── 信号处理作业8.docx
└── 研数字信号处理考查题.doc
0 directories, 16 files