D.1 Programs for first-order analysis

Program D.1 (spESTframeLaheWFI.m)   D.1 [*] - is used to compute the displacements and internal forces of a plane frame.
Called functions:
-       LaheFrameDFIm.m;
-       SisejoudPunktism.m.

Function D.1   (LaheFrameDFIm(baasi0,Ntoerkts,esQkoormus,esFjoud, sSolmF,tsolm,tSiire,krdn,selem)) D.2 [*], [*], [*], [*], [*], [*], [*] - is used to assemble and solve the boundary problem equations of a plane frame.
Called functions:
-       yzhqzm(baasi0,x,a,qx,qz,EA,EJ) D.3;
-       InsertBtoA(A,I,J,IM,JN,B,M,N) D.4;
-       spInsertBtoA(spA,IIv,IJv,spvF) D.5;
-       spSisestaArv(spA,iv,jv,sv) D.6;
-       SpTeisendusMaatriks2x2(NSARV,NEARV,VarrasN,krdn,selem) D.7;
-       SpTeisendusMaatriks(NSARV,NEARV,VarrasN,krdn,selem) D.8;
-       SpTeisendusUhikMaatriks2x2(VarrasN) D.9;
-       SpTeisendusUhikMaatriks0x1v(VarrasN) D.10;
-       SpTeisendusUhikMaatriks(VarrasN) D.11;
-       SpToeReaktsioonZvektor(NSARV,NEARV,VarrasN,krdn,selem) D.12;
-       SpToeReaktsioonXvektor(NSARV,NEARV,VarrasN,krdn,selem) D.13;
-       SpToeSiirdeFiVektor(VarrasN) D.14;
-       SpToeSiirdeUvektor(NSARV,NEARV,VarrasN,krdn,selem) D.15;
-       SpToeSiirdeWvektor(NSARV,NEARV,VarrasN,krdn,selem) D.16;
-       VardadSolmes(NSARV,NEARV,Solm,AB,ABB) D.17;
-       VardaPikkus(NSARV,NEARV,krdn,selem) D.18;
-       ylfhlin(baasi0,x,EA,GAr,EJ) D.19;
-       ysplfhlin(baasi0,x,EA,GAr,EJ)) D.20;
-       ysplvfmhvI(baasi0,x,l,EA,GAr,EJ) D.21;
-       yzfzv(baasi0,x,a,Fx,Fz,EA,EJ) D.22;
-       yzhqzm(baasi0,x,a,qx,qz,EA,EJ) D.23.

Function D.2   (SisejoudPunktism(VardaNr,X,AlgPar,lvarras,selem,
esQkoormus,esFjoud,suurused))
D.24 [*], [*], [*] - is used to compute the displace-ments and internal forces of the element 'VardaNr' at x = X.
Called functions:
-       ESTFrKrmus(baasi0,xx,Li,Fjoud,qkoormus,EA,EI) D.25.

Program D.2 (spESTframe93LaheWFI.m)   D.26 [*], [*] - is used to compute the displacements and internal forces of a plane frame.
Called functions:
-       LaheFrameDFIm.m;
-       SisejoudPunktism.m.

Program D.3 (spESTframe77LaheWFI.m)   D.27 [*] - is used to compute the displace-ments and internal forces of a plane frame.
Called functions:
-       LaheFrameDFIm.m;
-       SisejoudPunktism.m.

Program D.4 (spESTframe3hingeLaheWFI.m)   D.28 [*] - is used to compute the displacements and internal forces of a plane frame.
Called functions:
-       LaheFrameDFIm.m;
-       SisejoudPunktism.m.

Program D.5 (spESTframe3hingeLaheNQM.m)   D.29 [*] - is used to compute the internal forces of a plane frame.
Called functions:
-       LaheFrame3hingeNQM.m;
-       Sisejoud3LraamiPnktism.m.

Function D.3   (LaheFrame3hingeNQM(Ntoerkts,esQkoormus,esFjoud, sSolmF,tsolm,krdn,selem)) D.30 [*], [*] - is used to assemble and solve the boundary problem equations of a statically determinate plane frame.
Called functions:
-       ESTSKrmus(xx,Li,Fjoud,qkoormus) D.31;
-       InsertBtoA(A,I,J,IM,JN,B,M,N) D.32;
-       spInsertBtoA(spA,IIv,IJv,spvF) D.33;
-       spSisestaArv(spA,iv,jv,sv) D.34;
-       SpTeisendusMaatriks2x2D(NSARV,NEARV,VarrasN,krdn,selem) D.35;
-       SpTeisendusMaatriksD(NSARV,NEARV,VarrasN,krdn,selem) D.36;
-       SpTeisendusUhikMaatriks0x1v(VarrasN) D.37;
-       SpTeisendusUhikMaatriks2x2(VarrasN) D.38;
-       SpTeisendusUhikMaatriks(VarrasN) D.39;
-       VardadSolmesD(NSARV,NEARV,Solm,AB,ABB) D.40;
-       VardaPikkusD(NSARV,NEARV,krdn,selem) D.41;
-       ylSfhlin(x) D.42;
-       yspSlfhlin(x) D.43;
-       yspSlvfmhvI(x) D.44;
-       yzSfzv(x,a,Fx,Fz) D.45;
-       yzShqz(x,qx,qz) D.46.

Function D.4   (Sisejoud3LraamiPnktism(VardaNr,X,AlgPar,lvarras,
esFjoud,esQkoormus,suurused))
D.47 [*] - is used to compute the displacements and internal forces of the element 'VardaNr' at x = X.
Called function:
-       ESTSKrmus(xx,Li,Fjoud,qkoormus) D.48.

Program D.6 (spESTGerberBeamQM.m)   D.49 [*] - is used to compute the internal forces of a Gerber beam.
Called functions:
-       LaheGerberBeamQM.m;
-       SsjoudGrbrTalaPnktis.m.

Function D.5   (LaheGerberBeamQM(Ntoerkts,esQkoormus,esFjoud, sSolmF,tsolm,krdn,selem)) D.50 [*], [*] - is used to assemble and solve the boundary problem equations of a statically determinate beam.
Called functions:
-       VardaPikkusDT(NSARV,NEARV,krdn,selem) D.51;
-       yspSTlvfmhvI(x) D.52;
-       yspSTlfhlin(x) D.53;
-       ESTSTKrmus(xx,Li,Fjoud,qkoormus) D.54;
-       yzSTfzv(x,a,Fz) D.55;
-       yzSThqz(x,qz) D.56;
-       VardadSolmesDT(NSARV,NEARV,Solm,AB,ABB D.57;
-       SpTeisendusUhikMaatriks2x2(VarrasN) D.58;
-       SpTeisendusUhikMaatriks1x0(VarrasN) D.59;
-       ylSTfhlin(x) D.60;
-       SsjoudGrbrTalaPnktis(VardaNr,X,AlgPar,lvarras,esFjoud,esQkoormus; suurused) D.61;
-       spInsertBtoA(spA,IIv,IJv,spvF) D.62;
-       spSisestaArv(spA,iv,jv,sv) D.63;
-       InsertBtoA(A,I,J,IM,JN,B,M,N) D.64.

Function D.6   (SsjoudGrbrTalaPnktis(VardaNr,X,AlgPar,lvarras,
esFjoud,esQkoormus,suurused))
D.65 [*] - is used to compute the displacements and internal forces of the element 'VardaNr' at x = X.
Called functions:
-       ESTSTKrmus.m.

Program D.7 (spESTGerberBeamWFI.m)   D.66 [*] - is used to compute the displacements and internal forces of a Gerber beam.
Called functions:
-       LaheBeamDFI.m;
-       SisejoudTalaPunktis.m.

Program D.8 (spESTbeamLaheWFI.m)   D.67 [*] - is used to compute the displacements and internal forces of a beam.
Called functions:
-       LaheBeamDFI.m;
-       SisejoudTalaPunktis.m.

Function D.7   (LaheBeamDFI(baasi0,Ntoerkts,esQkoormus,esFjoud, sSolmF,tsolm,tSiire,krdn,selem)) D.68 [*] [*] - is used to assemble and solve the boundary problem equations of a beam.
Called functions:
-       VardaPikkusT(NSARV,NEARV,krdn,selem) D.69;
-       yspTlvfmhvI(baasi0,x,l,GAr,EJ) D.70;
-       yspTlfhlin(baasi0,x,GAr,EJ) D.71;
-       ESTtalaKrmus(baasi0,xx,Li,Fjoud,qkoormus,EI) D.72;
-       yzTfzv(baasi0,x,a,Fz,EJ) D.73;
-       yzThqz(baasi0,x,qz,EJ) D.74;
-       VardadSolmesT(NSARV,NEARV,Solm,AB,ABB) D.75;
-       SpTeisendusUhikMaatriks1x0(VarrasN) D.76;
-       SpTeisendusUhikMaatriks2x2(VarrasN) D.77;
-       SpToeSiirdeWvektorT(VarrasN) D.78;
-       SpToeSiirdeFiVektorT(VarrasN) D.79;
-       ylTfhlin(baasi0,x,GAr,EJ) D.80;
-       SisejoudTalaPunktis(VardaNr,X,AlgPar,lvarras,selem,esFjoud,esQkoormus; suurused) D.81;
-       spInsertBtoA(spA,IIv,IJv,spvF) D.82;
-       spSisestaArv(spA,iv,jv,sv) D.83;
-       InsertBtoA(A,I,J,IM,JN,B,M,N) D.84.

Function D.8   (SisejoudTalaPunktis(VardaNr,X,AlgPar,lvarras,selem,
esFjoud,esQkoormus,suurused))
D.85 [*] - is used to compute the displacements and internal forces of the element 'VardaNr' at x = X.
Called function:
-       ESTtalaKrmus.m.

Function D.9 (ESTtalaKrmus(baasi0,xx,Li,Fjoud,qkoormus,EI))   D.86 [*] - is used to compute the loading vector (q + F) for a continuous beam.
Called functions:
-       yzThqz(baasi0,x,qz,EJ) D.87;
-       yzTfzv(baasi0,x,a,Fz,EJ) D.88.

Program D.9 (spESTtrussLaheWFI.m)   D.89 [*] - is used to compute the displacements and internal forces of a plane truss.
Called function:
-       LaheTrussDFI.m.

Function D.10   (LaheTrussDFI(baasi0,Ntoerkts,sSolmF,tsolm,tSiire,
krdn,selem))
D.90 [*], [*] - is used to assemble and solve the boundary problem equations for a truss.
Called functions:
-       VardaPikkusTr(NSARV,NEARV,krdn,selem) D.91;
-       yspSRmhvI(baasi0,x,EA) D.92;
-       yspSRhlin(baasi0,x,EA) D.93;
-       VardadSolmesTr(NSARV,NEARV,Solm,AB,ABB) D.94;
-       SpTeisendusMaatriksTr2x2(NSARV,NEARV,VarrasN,krdn,selem) D.95;
-       SpTeisendusMaatriksTr2x1(NSARV,NEARV,VarrasN,krdn,selem) D.96;
-       SpTeisendusUhikMaatriks0x1v(VarrasN) D.97;
-       SpTeisendusUhikMaatriks2x2(VarrasN) D.98;
-       SpToeSiirdeUvektorTr(NSARV,NEARV,VarrasN,krdn,selem) D.99;
-       SpToeSiirdeWvektorTr(NSARV,NEARV,VarrasN,krdn,selem) D.100;
-       spInsertBtoAvect(spA,IM,JN,spB) D.101;
-       spInsertBtoA(spA,IIv,IJv,spvF) D.102;
-       spSisestaArv(spA,iv,jv,sv) D.103;
-       InsertBtoA(A,I,J,IM,JN,B,M,N) D.104.

Program D.10 (spESTtrussN2.m)   D.105 [*] - is used to compute the internal forces of a plane truss.
Called functions:
-       spSisestaArv(spA,iv,jv,sv) D.106;
-       spInsertBtoA(spA,IIv,IJv,spvF) D.107.

Function D.11 (spInsertBtoA(spA,IM,JN,spB))   D.108 [*], [*], [*], [*] - inserts a sparse matrix spB into a sparse matrix spA, starting at row index IM and column index JN. The overlapping elements of the matrices spA and spB are added together.

Function D.12 (InsertBtoA(A,I,J,IM,JN,B,M,N))   D.109 [*] - inserts matrix B (dimensions M, N) into matrix A (dimensions I, J), starting at row index IM and column index JN.

Function D.13 (spSisestaArv(spA,iv,jv,sv))   D.110 [*], [*] - inserts the number sv into sparse matrix spA, starting at row index iv and column index jv.

Function D.14   (SpToeSiirdeUvektor(NSARV,NEARV,VarrasS,krdn,
selem))
D.111 [*] - transforms a displacement of an element in local coordinates into a displacement Ux at the support in global coordinates.

Function D.15   (SpToeSiirdeWVektor(NSARV,NEARV,VarrasS,krdn,
selem))
D.112 [*] - transforms a displacement of an element in local coordinates into a displacement Wx at the support in global coordinates.

Function D.16 (SpToeSiirdeFiVektor(VarrasN))   D.113 [*] - transforms a displacement of an element in local coordinates into a displacement Fi at the support in global coordinates.

Function D.17 (ylfhlin(baasi0,x,EA,GAr,EJ)   D.114 [*] - is used to compute the transfer matrix for a frame (Sign Convention 2). Here,
EA - axial stiffness of the element,
GAr - shear stiffness of the element,
EJ - bending stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Function D.18 (yzhqz(baasi0,x,qx,qz,EA,EJ))   D.115 [*], [*], [*] - is used to compute the loading vector of a uniformly distributed load for a frame. Here,
qx - uniformly distributed load along the x-axis,
qz - uniformly distributed load along the z-axis,
EA - axial stiffness of the element,
EJ - bending stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Function D.19 (yzThqz(baasi0,x,qz,EJ)))   D.116 [*] - is used to compute the loading vector of a uniformly distributed load for a beam. Here,
qz - uniformly distributed load along the z-axis,
EJ - bending stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Function D.20 (yzShqz(x,qx,qz))   D.117 [*] - is usedis used to compute the loading vector of a uniformly distributed load for a three-hinged frame. Here,
qx - uniformly distributed load along the x-axis,
qz - uniformly distributed load along the z-axis.

Function D.21 (yzSThqz(x,qz))   D.118 [*] - is used to compute the loading vector of a uniformly distributed load for a Gerber beam. Here,
qz - uniformly distributed load along the z-axis.

Function D.22 (yzfzv(baasi0,x,a,Fx,Fz,EA,EJ))   D.119 [*], [*], [*] - is used to compute the loading vector of the point load for a frame. Here,
Fx - point load acting in the x direction,
Fz - point load acting in the z direction,
a - distance of a load point from a start point,
EA - axial stiffness of the element,
EJ - bending stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Function D.23 (yzmyv(baasi0,x,a,My,EJ))   D.120 [*] - is used to compute the loading vector of a moment load. Here,
My - load along the y-axis.

Function D.24 (yzTfzv(baasi0,x,a,Fz,EJ))   D.121 [*] - is used to compute the loading vector of a point load for a beam. Here,
Fz - point load acting in the z direction,
a - distance of a load point from a start point,
EJ - bending stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Function D.25 (yzSfzv(x,a,Fx,Fz))   D.122 [*] - is used to compute the loading vector of a point load for a three-hinged frame. Here,
Fx - point load acting in the x direction,
Fz - point load acting in the z direction,
a - distance of a load point from a start point.

Function D.26 (yzSTfzv(x,a,Fz))   D.123 [*] - is used to compute the loading vector of a point load for a Gerber beam. Here,
Fz - point load acting in the z direction,
a - distance of a load point from a start point.

Function D.27 (ysplfhlin(baasi0,x,EA,GAr,EJ))   D.124 [*] - is used to compute the transfer matrix for a frame (Sign Convention 2). Here,
EA - axial stiffness of the element,
GAr - shear stiffness of the element,
EJ - bending stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Function D.28 (yspTlfhlin(baasi0,x,GAr,EJ))   D.125 [*] - is used to compute the transfer matrix for a beam (Sign Convention 2). Here,
GAr - shear stiffness of the element,
EJ - bending stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Function D.29 (yspSlfhlin(x))   D.126 [*] - is used to compute the transfer matrix for a three-hinged frame (Sign Convention 2).

Function D.30 (yspSTlfhlin(x))   D.127 [*] - is used to compute the transfer matrix for a Gerber beam (Sign Convention 2).

Function D.31 (ysplvfmhvI(baasi0,x,l,EA,GAr,EJ))   D.128 [*] - is used to compute the sparse transfer matrix for a frame (Sign Convention 2). Here,
EA - axial stiffness of the element,
GAr - shear stiffness of the element,
EJ - bending stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Called functions:
-       ysplfhlin(baasi0,x,EA,GAr,EJ) D.129;
-       spInsertBtoA(spA,IM,JN,spB) D.130.

Function D.32 (yspTlvfmhvI(baasi0,x,l,GAr,EJ))   D.131 [*] - is used to compute the sparse transfer matrix for a beam (Sign Convention 2). Here,
GAr - shear stiffness of the element,
EJ - bending stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Called functions:
-       yspTlfhlin(baasi0,x,GAr,EJ) D.132;
-       spInsertBtoA(spA,IM,JN,spB) D.133.

Function D.33 (yspSlvfmhvI(x))   D.134 [*] - is used to compute the transfer matrix for a three-hinged frame (Sign Convention 2).
Called functions:
-       yspSlfhlin(x) D.135;
-       spInsertBtoA(spA,IM,JN,spB) D.136.

Function D.34 (yspSTlvfmhvI(x))   D.137 [*] - is used to compute the sparse transfer matrix for a Gerber beam (Sign Convention 2) .
Called functions:
-       yspSTlfhlin(x) D.138;
-       spInsertBtoA(spA,IM,JN,spB) D.139.

Function D.35 (yspSRmhvI(baasi0,Li,EA))   D.140 [*] - is used to compute the sparse transfer matrix for a truss.
Called functions:
-       yspSRhlin(baasi0,x,EA) D.141;
-       spInsertBtoA(spA,IM,JN,spB) D.142.

Function D.36 (yspSRhlin(baasi0,x,EA))   D.143 [*] - is used to compute the sparse transfer matrix for a truss. Here,
EA - axial stiffness of the element,
baasi0 - scaling multiplier for the displacements (i = EJ/l).

Function D.37 (ESTtalaKrmus(baasi0,xx,Li,Fjoud,qkoormus,EI))   D.144 [*] - is used to compute the loading vector (q + F) for a continuous beam.
Called functions:
-       yzThqz(baasi0,x,qz,EJ) D.145;
-       yzTfzv(baasi0,x,a,Fz,EJ) D.146.

Function D.38 (ESTSKrmus(xx,Li,Fjoud,qkoormus))   D.147 [*] - is used to compute the loading vector (q + F) for a three-hinged frame.
Called functions:
-       yzShqz(x,qx,qz) D.148;
-       yzSfzv(x,a,Fx,Fz) D.149.

Function D.39 (ESTSTKrmus(xx,Li,Fjoud,qkoormus))   D.150 [*] - is used to compute the loading vector (q + F) for a Gerber beam.
Called functions:
-       yzSThqz(x,qz) D.151;
-       yzSTfzv(x,a,Fz) D.152.

Program D.11 (spESTframe50WFI.m)   D.153 - is used to compute the displacements and internal forces of a plane frame.
Called function:
-       LaheFrameDFIm.m D.154.

Program D.12 (spESTframe3hinge1WFIm)   D.155 - is used to compute the displacements and internal forces of a plane frame.
Called function:
-       LaheFrameDFIm.m D.156.

Program D.13 (spESTframe3hinge1NQM.m)   D.157 - is used to compute the internal forces of a plane frame.
Called function:
-       LaheFrame3hingeNQM.m D.158.

andres
2014-09-09