Euler-Bernoullibeamtheroy(欧拉伯努利梁理论与有限元计算).docVIP

Euler-Bernoulli beam 一、理论部分 Euler-Bernoulli beam 假设 由(1.1)式可得 虚位移原理 其中 令，为单元长度，则上式成为 单元节点位移取为 令 其中形函数 分别对式(1.7)、(1.8)求一阶和两阶导数得 由，可得 将(1.6)、(1.11)式代入(1.4)式，可得 刚度矩阵 二、算例 本文所举算例为悬臂梁端部受竖向载荷，挠度理论曲线为 具体数据如下： ，，，，， 划分10个单元，11个结点。 由上图可以看出，数值模拟结构跟理论解吻合非常好。 源代码(Matlab) % program_main initial_statistic; le=pretreatment(ND,NE,lamda,x,y); [Dis,K]=Euler_Bernoulli_beam(ND,NE,lamda,EA,EI,le,P,Node_c,Dis_c); post(Dis,x,y,ND,NE); % initial statistic; ND=11;NE=10; EA=2.1e8*ones(NE,1);EI=4.375e4*ones(NE,1); x=(0:0.1:1);y=zeros(11,1); lamda=[1,2;2,3;3,4;4,5;