Ciao Guglielmo, innanzitutto ti ringrazio per l'interessamento.
I calcoli di cui sopra la prima volta li ho incontrati nello studio della
scienza delle costruzioni, ma si ritrovano anche in corsi successivi più settoriali inerenti l'
ingegneria strutturale e ovviamente poi anche nel mondo lavorativo. Il problema di fondo è che, a livello lavorativo, nessuno calcola direttamente tali quantità a meno di casi particolarmente cruciali, bensì si trovano alla vecchia maniera consultando dei prontuari o alla nuova maniera affidandosi ad uno dei millemila software tutto fare che li calcola tramite
Finite Element Method.
D'altro canto, a me è sempre piaciuto automatizzare i calcoli affidandomi alle mie limitate conoscenze di matematica, quindi mi sono creato un semplice foglio Excel in cui in input fornisco le coordinate dei vertici di una
regione piana il cui bordo è costituito da segmenti rettilinei (perlomeno questo è il mio target prefissato).
Il caso più semplice che permette un minimo di discussione è quello di una
regione rettangolare, diciamo di vertici \((0,0)\), \((b,0)\), \((b,h)\), \((0,h)\) (il sistema di riferimento iniziale è del tutto arbitrario, non influisce molto). Quindi, il codice di calcolo che sto scrivendo triangolarizza la regione di piano diciamo in un primo triangolo di vertici \((0,0)\), \((b,0)\), \((0,h)\) e in un secondo triangolo di vertici \((b,0)\), \((b,h)\), \((0,h)\). E fin qui tutto regolare.
A questo punto, noti i vertici di tali triangoli, per ciascuno il codice è in grado di calcolare in modo automatico area, momenti del primo ordine, momenti del secondo ordine e così via, tutte quantità che preventivamente ho determinato simbolicamente su carta e poi trasferito in VBA di Excel. Sommando tra loro tali quantità si ottengono quelle rispettive della regione rettangolare in input, che è quanto desiderato.
I grattacapi nascono quando passo al calcolo dell'integrale \(I\) sopra esposto. Se da un lato riesco a calcolare analiticamente tale integrale considerando la regione piana nel suo complesso (a titolo d'esempio ho mostrato i passaggi per il caso di regione circolare e rettangolare), dall'altro lato sono in difficoltà nel capire come poter giungere a tali risultati considerando i triangoli di cui è composta una regione del tipo qui considerato.
EDIT: non riuscendo a trovare una formulazione in forma chiusa, ho risolto tramite
integrazione numerica. Sostanzialmente ho spazzolato la
regione piana da \(y=y_{\min}\) ad \(y=y_{\max}\) e per ogni corda ne ho calcolato la lunghezza grazie alla conoscenza dei vertici dei triangoli che compongono la regione. Ciò fatto, ho calcolato il primo integrale tramite la
formula di Simpson composita e successivamente ho calcolato il secondo integrale, quello che effettivamente mi interessa, applicando nuovamente la formula di Simpson. Buona notte, ciao!