Strutture reticolari
File MATLAB
-
Funzione per il calcolo della matrice di rigidezza:
kbar.m
. -
Funzione ausiliaria per il calcolo di una base delle rotazioni infinitesime in $ \mathbb{R}^D $:
rots.m
. -
Script per la soluzione guidata di un problema in $ \mathbb{R}^2 $:
barND.m
. Questo script ha validità generale e può essere facilmente modificato per risolvere problemi in $ \mathbb{R}^3 $. -
Esempio di una esecuzione dello script
barND.m
.
Questo primo insieme di file MATLAB è stato scritto senza curarsi dell’ottimizzazione del codice, ma cercando di tradurre in maniera più semplice possibile le formulazioni matematiche di base. Molte operazioni potrebbero essere riscritte in una forma più concisa ed efficiente.
Il codice ha validità generale in $ \mathbb{R}^D $ con $ D \geq 1 $. Questo è possibile perché la formulazione utilizzata (basata semplicemente sul calcolo dell’energia elastica e del potenziale dei carichi) è indipendente dalle dimensioni dello spazio in cui viene definita la struttura reticolare. Ovviamente soluzioni con $ D \geq 4 $, pur essendo formalmente corrette, non hanno significato fisico.