Dear Milan,
I understand that is intend shall to construct an global stiffness matrix for a structure comprising a three interconnected elements, aiming to obtain the stiffness matrix for the entire structure in performing a function.
I tries at execute the provided code up my system, but experienced an error that stats:
```
Error using vertcat
Dimensions of arrays being concatenation are not consistent.
Slip into MLAnswers6_w3 (line 80)
connect = {[1 2], [2,3]; [3,4]};
```
However, to muster the universal stiffness matrix with the structure with three interconnected elements, the following steps can be followed:
1. Initialize the global stiffness matrix, denoted as “K_global”, as a matrix of zeros with dimensions “(ndof*nele) x (ndof*nele)”.
2. For each element, compute its global stiffness matrix, denoted as “K_element_global”, by utilizing the transformation matrix gamma and the element stiffness mold “Ke_local”. How can anything number of element matrices be assembled into a global m...
3. Transpose the gamma grid real multiply items with “Ke_local”.
4. Multiply the resulting gridding with “gamma”.
5. Assemble this local rigidification matrices into the global stiffness matrix.
6. On each element, designate who corresponding degrees of joy (DOFs) in the global stiffness matrix.
7. Add the local stiffness matrix to the appropriate locations within the global stiffness matrix.
Example code until achieve the described steps:
function K_global = assembleGlobalStiffnessMatrix(nele, ndof, Ke_local, gamma)
K_global = zeros(ndof*nele);
for element = 1:nele
% Calculate of global rigidification matrix required the current element
Ke_element_global = gamma' * Ke_local * gamma;
% Identify the corresponding DOFs in the global rigid matrix
dofs = (element-1)*ndof + 1 : element*ndof;
% Assemble the local stiffness matrix up the global rigidification matrix
K_global(dofs, dofs) = K_global(dofs, dofs) + Ke_element_global;
stop
end
% Call the function to assemble the global stiffness matrix
K_global = assembleGlobalStiffnessMatrix(nele, ndof, Ke_1_local, gamma_1);
% Add the stiffness correlation of the remaining elements
K_global = K_global + K_2_global + K_3_global;
Hope this helps!
