Bar2D2Node曾攀有限元基础教程（3.2.5）

 bar2d2node

主命令文件

clc
%本文件提供所属调用函数中所需的参数，函数返回结果可作为后续调用函数的输入参数
E=2.95e11;
A=0.0001;
x1=0;
y1=0;
x2=0.4;
y2=0;
x3=0.4;
y3=0.3;
x4=0;
y4=0.3;
alpha1=0;
alpha2=90;
alpha3=atan(0.75)*180/pi;
Bar2D2Node=Bar2D2NodeFunc;%自建函数的集合
k1=Bar2D2Node.Stiffness(E,A,x1,y1,x2,y2,alpha1)
k2=Bar2D2Node.Stiffness(E,A,x2,y2,x3,y3,alpha2)
k3=Bar2D2Node.Stiffness(E,A,x1,y1,x3,y3,alpha3)
k4=Bar2D2Node.Stiffness(E,A,x4,y4,x3,y3,alpha1)
KK=zeros(8);%（8*8）初始KK==0提供位置
KK=Bar2D2Node.Assembly(KK,k1,1,2); %杆1
KK=Bar2D2Node.Assembly(KK,k2,2,3);%加杆2
KK=Bar2D2Node.Assembly(KK,k3,1,3);%加杆3
KK=Bar2D2Node.Assembly(KK,k4,4,3) %加杆4

k=KK([3 5 6],[3 5 6]) %加入边界条件采用消去法化简
p=[20000;0;-25000];
u=k\p  %高斯消去法（左除）

q=[0 0 u(1) 0 u(2) u(3) 0 0]';
P=KK*q

u1=[q(1);q(2);q(3);q(4)]
stress1=Bar2D2Node.Stress(E,x1,y1,x2,y2,alpha1,u1)
u2=[q(3);q(4);q(5);q(6)]
stress2=Bar2D2Node.Stress(E,x2,y2,x3,y3,alpha2,u2)
u3=[q(1);q(2);q(5);q(6)]
stress3=Bar2D2Node.Stress(E,x1,y1,x3,y3,alpha3,u3)
u4=[q(7);q(8);q(5);q(6)]
stress4=Bar2D2Node.Stress(E,x4,y4,x3,y3,alpha1,u4)

bar2d2nodefun

调用函数

%%% 4个函数放在一个文件中（总体函数）
function Bar2D2Node=Bar2D2NodeFunc
Bar2D2Node.Stiffness=@Bar2D2Node_Stiffness;
Bar2D2Node.Assembly=@Bar2D2Node_Assembly;
Bar2D2Node.Stress=@Bar2D2Node_Stress;
Bar2D2Node.Forces=@Bar2D2Node_Forces;
end
%%%%% Bar2D2Node %%begin%%%%
function k=Bar2D2Node_Stiffness(E,A,x1,y1