tcable.m

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

function t=Tcable(ux,uy,uz,s0,E,A,w)
%导入数据求出水平力。
%导入的数据有节点的坐标，原长，索面积

% syms x %ll rr cc EE AA S0 ww tt
l=sqrt(ux^2+uy^2);
r=sqrt(l^2+uz^2);
% siag1=subs(siag1)
% fs=simple(fs);

% fs=subs(FS,{ll,rr,cc,EE,AA,S0,ww,tt},{l,r,uz,E,A,s0,w,x});

t0=E*A*(r-s0)/s0;
if t0<0
    t0=5000;
end
x1=10;
x2=t0;
gent=0;
while abs(x2-x1)>100
    y1=fs(l,r,uz,s0,E,A,w,x1);
    y2=fs(l,r,uz,s0,E,A,w,x2);
    if y1>0&y2>0
        x1=x2;
        x2=2*x2;
    elseif y1>0&y2<0
        x3=(x1+x2)/2;
        y3=fs(l,r,uz,s0,E,A,w,x3);
        if y3>0
            x1=x3;
        elseif y3<0
            x2=x3;
        end
    elseif y1<0&y2<0
        x1=x2;
        x2=2*x2;
    elseif y1<0&y2>0
        x3=(x1+x2)/2;
        y3=fs(l,r,uz,s0,E,A,w,x3);
        if y3>0
            x2=x3;
        elseif y3<0
            x1=x3;
        end
    elseif y1==0
        t=x1;
        gent=1
    elseif y2==0
        t=x2;
        gent=1;
    end
end
if gent==0
    t=(x1+x2)/2;
end