bita_fct.m

基于力密度法的索杆张力结构找形程序

function bita=bita_fct(r,E,A,l,t,uz,s0,w)
c=abs(uz);
bita=1/8*(-1/(w^2*r^2+4*c^2*t^2-4*c*w*r*t+4*l^2*t^2)^(1/2)*c*t*(8*c^2*t-4*c*w*r+8*l^2*t)-2*(w^2*r^2+4*c^2*t^2-4*c*w*r*t+4*l^2*t^2)^(1/2)*c+1/2/(w^2*r^2+4*c^2*t^2-4*c*w*r*t+4*l^2*t^2)^(1/2)*w*r*(8*c^2*t-4*c*w*r+8*l^2*t)+(4*(-2*c*w*r+1/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*((8*c^2*t-4*c*w*r+8*l^2*t)/l^2/t^2-2*(w^2*r^2+4*c^2*t^2-4*c*w*r*t+4*l^2*t^2)/l^2/t^3)-((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/t*w^2*r^2+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)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2)-8*(-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^3/(w^2*r^2/t^2/l^4)^(1/2)+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^7/t^5/(w^2*r^2/t^2/l^4)^(3/2)*w^2*r^2)/(-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^5*t^4*(w^2*r^2/t^2/l^4)^(1/2)+8*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+1/(w^2*r^2+4*c^2*t^2+4*c*w*r*t+4*l^2*t^2)^(1/2)*c*t*(8*c^2*t+4*c*w*r+8*l^2*t)+2*(w^2*r^2+4*c^2*t^2+4*c*w*r*t+4*l^2*t^2)^(1/2)*c+1/2/(w^2*r^2+4*c^2*t^2+4*c*w*r*t+4*l^2*t^2)^(1/2)*w*r*(8*c^2*t+4*c*w*r+8*l^2*t)-(4*(-2*c*w*r+1/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*((8*c^2*t+4*c*w*r+8*l^2*t)/l^2/t^2-2*(w^2*r^2+4*c^2*t^2+4*c*w*r*t+4*l^2*t^2)/l^2/t^3)-((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/t*w^2*r^2+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)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2)-8*(-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^3/(w^2*r^2/t^2/l^4)^(1/2)+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^7/t^5/(w^2*r^2/t^2/l^4)^(3/2)*w^2*r^2)/(-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^5*t^4*(w^2*r^2/t^2/l^4)^(1/2)-8*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)/t/r/w-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^2/r/w-1/E/A*(l^2/r+c^2/r+1/12*w^2*r/t^2)+1/6/t^2/E/A*w^2*r;