Differenze

Queste sono le differenze tra la revisione selezionata e la versione attuale della pagina.

--- wikitelaio2015:lez39 [2015/06/18 07:18] – [Costruzione del grafico “History Plot”] 168316
+++ wikitelaio2015:lez39 [2015/07/02 08:44] (versione attuale) – [Instabilità nei fenomeni non lineari] 168316
@@ Linea 1: / Linea 1: @@
+====== 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.