Page 4 - Demo
P. 4
Planar Truss Example for Comrel Add-on RCP Consult, 2023-2026 Page 4//* Global stiffness matrix double K[ndof][ndof]; for(int i=0; i=0; i--) { ac[i]=fc[i]; for(int j=ndof-1; j>i; j--) ac[i]=ac[i]-K[i][j]*ac[j]; ac[i]=ac[i]/K[i][i]; }//* Axial loads in bars double NN[nfe]; for(int e=0; e
