aifax_fct.m

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

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