bendanalysis.m

来自「薄壁结构有限元计算」· M 代码 · 共 21 行

% 闭口截面受弯所产生的剪流分析及弯曲中心判定

% 输入产生弯矩的外力荷载和截面的惯性矩
Ny=input('Please input the y-direction load: ');
Ix=input('Please input the inertia constant to x-axis: ');

% 计算s/t的沿程积分值
tj=sum((L*I1)./(T*I1));

% 输入用于计算各回路剪流值的方程系数矩阵和值列
Af=input('Please input the matrix to calculate shear current: ');
bf=input('Please input the array to calculate shear current: ');

% 计算各回路剪流值
f0=Af\(-Ny*bf/Ix);

% 输入剪流沿程对某点取矩的积分值
fh=input('Please input the integral value of shear current: ');

% 计算弯曲中心的位置
xa=-(fh+2*Ac'*f0)/Ny;