cdm.ing.unimore.it

Questa pagina è in sola lettura. Puoi visualizzare il sorgente, ma non puoi modificarlo. Contatta l'amministratore se pensi che ci sia un errore.
====== Risposta dinamica delle strutture ======
riprendiamo da {{:wikitelaio2015:forcella_su_profilo_v5.mfd|modello lezione scorsa}}.

Il modello era sollecitato da forzante di entità $$98.70 \left(\frac{f}{25\mathrm{Hz}}\right)^2$$ sull'intera struttura, ridotto a $$49.35 \left(\frac{f}{25\mathrm{Hz}}\right)^2$$ sulla porzione modellata per simmetria.

Utilizzo smorzamento in forma di Rayleigh, da cui 
$$
\zeta_i = \frac{1}{2}\left(\frac{\alpha}{\omega_i} + \beta \omega_i\right)
$$
Suppongo $\zeta=0.01$ per ogni modo proprio, e $\alpha=0$, da cui $$\beta=\frac{2\zeta}{\omega}=\frac{\zeta}{\pi f}$$

...

Controllo risposta in risonanza:
$$
\left|\xi_i\right| = \frac{\left|q_i\right|}{\left(2 \pi f_i\right)^2}\frac{1}{2\zeta}
$$
{{:wikitelaio2015:controllo_risonanza.ods|}}

====== modellazione pannello honeycomb ======
{{:wikitelaio2015:honeycomb_sandwich_design_technology.pdf|}}

====== instabilita ======
{{:wikitelaio2015:piramide_buckling_base_v2012.mfd|}}

====== Costruzione del grafico “History Plot” ======
Si tratta di un grafico che mostra l'evoluzione del caso nel tempo. Voglio plottare lo spostamento in y del punto centrale del corpo rigido (simulacro di compressore) al variare della frequenza dell'eccitante.

Set location > Scelgo il punto centrale del corpo rigido
Increment range > 0:1	0:61		1 (Questo comando memorizza la “storia” del punto selezionato dall'istante 1 al 61 con passo 1.

Add curves > Single location > Variabile x: Frequency; Variabile y: Displacement y

In caso ci sia una risonanza nel range di lavoro, uso uno smorzatore (dal menu “Links”, oppure uso lo smorzamento strutturale:
Material properties > Structural > Damping
   * Stiffness matrix multiplier=0,01
   * Creo la tabella in cui inserisco la funzione che lega lo smorzamento alla frequenza: $ \beta=\frac{\zeta}{\pi f} $ 
new table > v1:frequency; formula:1/(v1*pi)

Loadcase > Restringo l'intervallo del loadcase a 21,39 – 21,69 Hz
Job > Properties > Attivo il “Complex Damping”

Analisi dei risultati: Ricordare di controllare sia la parte reale sia la parte immaginaria degli spostamenti (fasi 0° e 270°)

{{:wikitelaio2015:lezione_39.odt|{{:wikitelaio2015:lezione_39.odt|}}}}
====== Instabilità nei fenomeni non lineari ======

Considero un fenomeno di tipo softening; la curva carico–spostamento cala la sua pendenza gradualmente, fino diventare orizzontale. Da lì in poi, si hanno infinite condizioni di equilibrio; in formule: se det(K)=0 → K*Δδ=0 , e ciò comporta la possibilità di avere variazione di spostamento non nulla a fronte di variazione di forza nulla.
Dati i punti del grafico di coordinate (δ1,F1) e (δ2,F2), applicando il metodo di estrapolazione lineare posso trovare un altro punto approssimato della curva F-δ: Fλ = F1+λ*(F2-F1)
Faccio l'ipotesi che la matrice di rigidezza vari linearmente con  λ (formula analoga) 
Mi chiedo ora quando det(K λ)=0. Ciò equivale a det(K1+λ*(K2-K1)
Ottengo un problema della stessa forma di quello della ricerca dei modi propri (F1+λ*(F2-F1))*Δδ=0. Di default, Marc-Mentat usa come punto 1 la condizione di nessun caricamento.

===== Piramide buckling =====
Carico 1 =0, Carico 2 = 1000N
Job1: Attivo i carichi (lineari)
Loadcase: Buckle (Lascio il default)

Analisi dei risultati: 
   * Step 0: Pura compressione/Trazione
   * Step 1: Buckling per un carico uguale a 8.46 volte il carico di prova (fattore mostrato in alto nella schermata)
N.B.: Jobs > Options > Analysis options: Attivare “Large rotations”

La struttura geometricamente perfetta ha una discontinuità di comportamento a cavallo dell'instabilità: fino all'attimo precedente, presenta una pura compressione, poi improvvisamente collassa a flessione per instabilità. Per cogliere numericamente questo comportamento, bisogna mitigare la discontinuità, affinché il metodo di Newton-Rhapson converga a una soluzione. Introduciamo allora una piccola perturbazione geometrica, spostando un nodo di uno spigolo della piramide di 0,5 mm.
N.B.: Per giungere più velocemente ad un risultato in questo tipo di studi, procedere in controllo di spostamento piuttosto che di carico.