20 lines
878 B
Matlab
20 lines
878 B
Matlab
clc,clear,close all
|
||
tic
|
||
d1=30;d2=2*d1;d3=3*d1;d4=4*d1;d5=5*d1;d6=6*d1;d7=7*d1;d8=8*d1;
|
||
h1=80;h2=2*h1;h3=3*h1;h4=4*h1;
|
||
x=[d4 d3 d5 d2 d4 d6 d1 d3 d5 d7 d4 0 d8]; % 节点x轴方向坐标
|
||
y=[h4 h3 h3 h2 h2 h2 h1 h1 h1 h1 0 0 0]; % 节点y轴方向坐标
|
||
A=3;E=2.1E005; % 定义截面面积和弹性模量
|
||
%单元信息：编号，节点1编号，节点2编号，截面面积，弹性模量
|
||
ele=[1 1 2 A E;2 1 3 A E;3 2 3 A E;4 2 4 A E;5 2 5 A E;6 3 5 A E;7 3 6 A E;...
|
||
8 4 5 A E;9 5 6 A E;10 4 7 A E;11 4 8 A E;12 5 8 A E;13 5 9 A E;14 6 9 A E;...
|
||
15 6 10 A E;16 7 8 A E;17 8 9 A E;18 9 10 A E;19 7 12 A E;20 8 12 A E;...
|
||
21 7 11 A E;22 8 11 A E;23 9 11 A E;24 10 11 A E;25 9 13 A E;26 10 13 A E];
|
||
%载荷信息：节点编号，x方向力，y方向力，z方向力
|
||
Load=[1 200 100]
|
||
%约束：节点编号，x方向约束，y方向约束，z方向约束
|
||
Constr=[11 0 0;12 0 0;13 0 0];
|
||
[U, Strain, Stress, AxialForce]=Truss2DFEA(x,y,ele,Load,Constr,100 )
|
||
toc
|
||
|