clearvars close all InputFile='input_File_1'; %% Pre-Processing %DrawStructure? iDrawStructure=1; %Read Input File iReadFromExcel=0; switch iReadFromExcel case 0 %Directly From Matlab eval(InputFile); case 1 [anElement,Nodes,aStructure]=readFromExcel(InputFile); end %Draw Structure v=DrawStructure(iDrawStructure,0, anElement, Nodes); %% Processing %Indexing, Structural Matrix Allocations Process_Input; % Evaluate Elemental properties and finally the Structural Stiffness Matrix for i1=1:Nelements %% Evaluate Element Length %anElement{i1}.Length=GetLength(anElement{i1}.iEnd,Nodes); %% Evaluate Transformation Matrix %anElement{i1}.Lambda=Evaluate_Transformation_Matrix(anElement{i1}.Length,... % anElement{i1}.iEnd,Nodes); %% Evaluate Element Local Stiffness Matrix %anElement{i1}.k=Evaluate_Local_Stiffness(anElement{i1}); %% Evaluate Element Global Stiffness Matrix %anElement{i1}.kg=Evaluate_Global_Stiffness(anElement{i1}); %% Add Contribution of the Element Global Stiffness to the Structural %aStructure.Kg=aStructure.Kg+Evaluate_Elemental_Contribution(anElement{i1},Ndofs); end %Keep for output %% Rotate Nodes if required %% Boundary Conditions %Fetch rearrangement matrix %aStructure.V=Get_V_Mat(aStructure.Bdofs,Ndofs); %% External Loading %aStructure.Pext=Apply_External_Loading(Nodes,Ndofs); %% Analysis %%END OF PROCESSING %%BEGIN POST-PROCESSING %% Post-Processing %Get Element forces %Get Nodal Displacement %Draw Deformed Structure