utmosttorsion.m

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

% 截面塑性极限扭矩的计算

Z=[Cf -Cf]';

% 输入最大剪应力屈服极限
ty=input('Please input the maximum shearing strength(MPa): ');
Ty=ty*[(T*I1)' (T*I1)']';

% 截面回路数
Nc=Ne-Nv+1;
lb=zeros(Nc,1);

% 用线性规划的方法优化求出回路剪流值
Fpc=linprog(-Ac',Z,Ty,[],[],lb);

% 求出每一边的剪流值
Fpe=Cf'*Fpc;

% 求沿薄壁单元分布的扭矩值
Mps=((T*I1).*(T*I1)-(Fpe/ty).*(Fpe/ty)).*(L*I1)*ty/4;

% 求整个截面塑性极限扭矩值（包括塑性环流的扭矩值和沿薄壁单元分布的扭矩值）
Mp=2*Ac'*Fpc+sum(Mps);