fs.m

来自「基于力密度法的索杆张力结构找形程序」· M 代码 · 共 4 行

function y=fs(l,r,uz,s0,E,A,w,t)
c=abs(uz);
y=1/8*(-2*(w^2*r^2+4*c^2*t^2-4*c*w*r*t+4*l^2*t^2)^(1/2)*c*t+(w^2*r^2+4*c^2*t^2-4*c*w*r*t+4*l^2*t^2)^(1/2)*w*r+4*log(4*(-2*c*w*r*t+w^2*r^2+((w^2*r^2+4*c^2*t^2-4*c*w*r*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2))*l^2*t^2+2*(w^2*r^2+4*c^2*t^2+4*c*w*r*t+4*l^2*t^2)^(1/2)*c*t+(w^2*r^2+4*c^2*t^2+4*c*w*r*t+4*l^2*t^2)^(1/2)*w*r-4*log(4*(-2*c*w*r*t-w^2*r^2+((w^2*r^2+4*c^2*t^2+4*c*w*r*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2))*l^2*t^2)/t/r/w-s0-t/E/A*(l^2/r+c^2/r+1/12*w^2*r/t^2);