gamar_fct.m

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

function [gamarx,gamary,gamarz]=gamar_fct(R,L,Ux,Uz,Uy,Tforce,E,Area,weight,s0,cab_num,Etype)
%求各个的aifa的系数
for i=1:1:cab_num
    if Etype(i)==0
        aifax(i,1)=Bar_aifa(Ux(i,i),R(i,i),E(i),Area(i));
        aifay(i,1)=Bar_aifa(Uy(i,i),R(i,i),E(i),Area(i));
        aifaz(i,1)=Bar_aifa(Uz(i,i),R(i,i),E(i),Area(i));
    elseif Etype(i)==1
        aifax(i,1)=aifax_fct(R(i,i),Ux(i,i),Tforce(i,i),E(i),Area(i),L(i,i),Uz(i,i),weight(i),s0(i));
        aifay(i,1)=aifax_fct(R(i,i),Uy(i,i),Tforce(i,i),E(i),Area(i),L(i,i),Uz(i,i),weight(i),s0(i));
        aifaz(i,1)=aifaz_fct(R(i,i),Uz(i,i),Tforce(i,i),E(i),Area(i),L(i,i),weight(i),s0(i));
    end
end
[aifax,aifay,aifaz];
%求bia的系数
for i=1:1:cab_num
    if Etype(i)==0
        bita(i)=-s0(i);
    elseif Etype(i)==1
        bita(i)=bita_fct(R(i,i),E(i),Area(i),L(i,i),Tforce(i,i),Uz(i,i),s0(i),weight(i));
    end
end
% bita
%将aifa系数和bia系数组合
for i=1:1:cab_num
    pantpanux(i,i)=-aifax(i)/bita(i);
    pantpanuy(i,i)=-aifay(i)/bita(i);
    pantpanuz(i,i)=-aifaz(i)/bita(i);
end
gamarx=inv(R)*Tforce-Ux^2*inv(R^3)*Tforce+Ux*inv(R)*pantpanux;
gamary=inv(R)*Tforce-Uy^2*inv(R^3)*Tforce+Uy*inv(R)*pantpanuy;
gamarz=inv(R)*Tforce-Uz^2*inv(R^3)*Tforce+Uz*inv(R)*pantpanuz;