wikipaom2015:lez20
Differenze
Queste sono le differenze tra la revisione selezionata e la versione attuale della pagina.
Entrambe le parti precedenti la revisioneRevisione precedenteProssima revisione | Revisione precedente | ||
wikipaom2015:lez20 [2015/05/19 20:29] – 208041 | wikipaom2015:lez20 [2015/09/30 11:22] (versione attuale) – 182648 | ||
---|---|---|---|
Linea 1: | Linea 1: | ||
+ | Nella presente lezione, si va ad analizzare un approccio semplificato per la soluzione di problemi assial simmetrici, quali forzamenti albero mozzo, tubi etc. con caricamenti uniformi anch' | ||
+ | Si consideri un elemento 4 nodi assialsimmetrico piano, quindi in condizione di tensione e deformazione piana, | ||
+ | |||
+ | |||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | Nel manipolatore Marc, il quadrilatero deve essere considerato accoppiato, all' | ||
+ | Fisicamente si andrà a disegnare l' | ||
+ | Gli spostamenti nodali sono 2, uniformandosi al sistema cilindrico standard, si avrà: | ||
+ | - u spostamento radiale | ||
+ | - w spostamento assiale | ||
+ | Spostamenti e forze applicate sui singoli nodi, sono da considerarsi distribuiti e applicati sulla circonferenza ottenuta dalla rivoluzione, | ||
+ | pari ad una forza per unità di lunghezza moltiplicata per la lunghezza della circonferenza 2πR, di conseguenza carichi uniformi sono rappresentabili con forze nodali linearmente crescenti in R (img. 1 e 2) | ||
+ | Valutando in via esemplificativa un tubo, sempre in ipotesi di caricamento assial simmetrico, si può dimostrare che ogni piano diametrale rappresenti un piano, di simmetria, risulta dunque lecito valutare porzioni sempre più piccole, ottenute mediante i suddetti piani, che possono degenerare fino all' | ||
+ | {{ : | ||
+ | | ||
+ | Il numero di gradi di libertà dell' | ||
+ | {{ : | ||
+ | |||
+ | La definizione di** __caricamento elastoplastico__** e l' | ||
+ | |||
+ | ======Modellazione elastica====== | ||
+ | |||
+ | |||
+ | Costruiamo la mesh di un tubo assialsimmetrico. | ||
+ | Creiamo la fettina da 1 mm e da 100 mm, i punti e nodi di riferimento. | ||
+ | Clicchiamo in sequenza i seguenti pulsanti: **mesh generation, pts, add,** scriviamo sulla riga di comando in basso le coordinate degli estremi del segmento associato all’__asse di simmetria__: | ||
+ | |||
+ | x,y,z | ||
+ | |||
+ | 0,0,0 | ||
+ | |||
+ | 1,0,0 | ||
+ | |||
+ | Ora mettiamo direttamente i nodi e non i punti: **Mesh generation, node, add,** poi inserisco da tastiera i valori che voglio: | ||
+ | |||
+ | x, y,z | ||
+ | |||
+ | 0,10,0 | ||
+ | |||
+ | 0,20,0 | ||
+ | |||
+ | 1,10,0 | ||
+ | |||
+ | 1,20,0 | ||
+ | |||
+ | n.b. non ha senso fisico caricare i nodi che giacciono sull' | ||
+ | |||
+ | Ora creiamo un cubo ad 1 mm di distanza, aggiungiamo un punto in 2,0,0 e proviamo a mettere un altro punto in 102,0,0 ma è troppo lontano, lo eliminiamo facendo undo, metto un punto in 42,0,0. | ||
+ | |||
+ | Ora aggiungiamo i nodi analoghi: **mesh generation nodes add**: | ||
+ | |||
+ | 2,10,0 | ||
+ | |||
+ | 2,20,0 | ||
+ | |||
+ | 42,10,0 | ||
+ | |||
+ | 42,20,0 | ||
+ | |||
+ | Stessa cosa con quota x =42. | ||
+ | |||
+ | Ora creiamo i due elementi quadrilateri, | ||
+ | |||
+ | Abbiamo il volume del materiale elastico. | ||
+ | |||
+ | {{: | ||
+ | |||
+ | Ora suddividiamo i nostri elementi: parto da quello più grande, **mesh generation subdivide**, | ||
+ | |||
+ | Definiti i paramentri della suddivisione li applichiamo all' | ||
+ | |||
+ | Con mesh del tutto identica abbiamo meshato il tubo più grande la fettina di tubo più piccola. Abbiamo 481 nodi. Suddividiamo ulteriormente? | ||
+ | |||
+ | Ora abbiamo la mesh ed impostiamo i **vincoli**: | ||
+ | |||
+ | A questo punto abbiamo definito completamente la geometria, andiamo | ||
+ | |||
+ | Procedimento: | ||
+ | |||
+ | {{: | ||
+ | |||
+ | Torniamo al menù principale e definiamo le **proprietà del materiale**. Material properties, creiamo una nuova scheda, materiale di tipo standard e lo chiamiamo alluminio (110 MPa), proprietà generali mettiamo la densità (2.7e-9 tonnellate al mm cubico), ma non verrà usata nei calcoli. Poi proprietà strutturali young moduls 70000, Poisson' | ||
+ | |||
+ | **Material properties, spunto material properties, new, name mettiamo alluminio, invio, entriamo in data categories, general, inseriamo il valore di mass density 2.7e-9, invio, ok, entriamo in structural, inseriamo il valore di Young' | ||
+ | |||
+ | Ora andiamo ad impostare le **boundary condition**: | ||
+ | |||
+ | |||
+ | **Boundary conditions, new, structural, edge load, poi entriamo in properties, spuntiamo pressure, mettiamo 1 MPa, invio, ok, selezioniamo con il mouse i nodi, end list.** | ||
+ | |||
+ | {{: | ||
+ | |||
+ | Gli diamo una pressione e lui si calcola i carichi nodali che sono funzione del raggio. | ||
+ | |||
+ | **Vincoli**: | ||
+ | |||
+ | Creiamo un vincolo new, strutturale, | ||
+ | **Boundary conditions, new, structural, fixed displacement, | ||
+ | |||
+ | {{: | ||
+ | |||
+ | Ora il problema è pronto per essere lanciato. | ||
+ | |||
+ | **Menu principale, jobs, new, structural, properties: initial load, pressione interna, simmetria e simmetriaeposizionamento attivati.** | ||
+ | |||
+ | JOBS.PROPERTIES.JOB RESULTS. e quindi selezioniamo stress, total strain, elastic strain, plastic strain (=0), equivalent von mises stress, ok, poi " | ||
+ | Lanciamo il calcolo elastico, | ||
+ | Tutti gli elementi della struttura sono di tipo 10. | ||
+ | |||
+ | {{: | ||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | Qual è la discretizzazione? | ||
+ | |||
+ | Ora facciamo **run, submit,** apriamo il file dei risultati. | ||
+ | |||
+ | {{: | ||
+ | |||
+ | Per andare ad analizzare i risultati: | ||
+ | < | ||
+ | POSTPROCESSING RESULTS -> OPEN POST FILE | ||
+ | scelgo come SCALAR PLOT :CONTOUR BANDS | ||
+ | </ | ||
+ | Vado quindi ad analizzare i diversi tipi di SCALAR: | ||
+ | Comp 33 of stress | ||
+ | |||
+ | **NOTA BENE**: | ||
+ | Comp 11 of stress: indica una tensione assiale i ìn direzione x. | ||
+ | Comp 22 of stress: indica una tensione radiale | ||
+ | Comp 33 of stress: indica una tensione circonferenziale | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Notiamo che, essendo stato applicato 1Mpa di pressione al bordo interno, si ottiene 1.749 MPa di tensione circonferenziale al bordo interno. | ||
+ | Per controllare se questa tensione che risulta sia corretta si apre un foglio Maxima e si confronta il valore ottenuto (1.749) con quello calcolato in Maxima (1.66666); quindi il valore ottenuto con Marc Mentat è leggermente eccessivo ma rientra nell’ordine di errore. | ||
+ | Si noti di come la singola fetta di tubo è rappresentativa dell’oggetto complesso e le curve di livello sono esattamente le stesse; questo è dovuto al fatto che questa fetta è libera di contrarsi assialmente perchè si è solo messo un vincolo di posizionamento ma tutto il resto dei due bordi sono liberi di muoversi; questa fetta è quindi in uno stato di tensione piana (la tensione in z dovrebbe essere nulla). | ||
+ | In questo caso nessuno vieta alla $\sigma_z$ di essere nulla , mentre nel modello della lezione precedente questa era un’ ipotesi di partenza e quindi non era verificabile, | ||
+ | |||
+ | **Comp 11 of stress** | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Per verificare se l’ ipotesi di stato di tensione piana sia applicabile al tubo vado a vedere la Comp 11 of stress che non è nulla in senso assoluto ma è comunque dell’ordine di $10^(-1)$ ovvero al di sotto di uno o due ordini di grandezza rispetto la tensione circonferenziale; | ||
+ | Quando non impongo delle ipotesi di partenza di stato di tensione piana (come la lezione precedente) ma la vado a verificare, non è verificata in maniera esatta e precisa ma possiede un ordine di grandezza al di sotto di 1 o 2 ordini di grandezza delle altre componenti tensionali. | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Si noti che c’è una piccola deviazione nell’intorno dell’area finale, che è dovuta da un errore numerico e che in teoria dovrebbe calare al calare della mesh. | ||
+ | |||
+ | |||
+ | **Comp 11 of Total Strain e Comp 11 of Plastic Strain** | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Si controlla ora se effettivamente se la struttura deforma uniformemente lungo l’asse ovvero se tutte le fette si compattano in maniera uniforme. | ||
+ | Ciò lo si va a vedere andando a guardare la Comp 11 of Total strain (la deformazione lungo x, in direzione assiale) e si nota che è compresa tra una forchetta che va da $3.4e^(-6)$ a $2.7e^(-6)$ mentre nel tubo finale è molto più informe. | ||
+ | Nella teoria della tensione piana questa quantità, ovvero la deformazione z, dovrebbe essere una quantità libera cioè se impongo una tensione in direzione assiale nulla lascio libera questa tensione. | ||
+ | In realtà per una particolare proprietà del tubo pur imponendo nulla la tensione in direzione fuori piano la deformazione rimane costante lungo la struttura. | ||
+ | Qui nel tubo appare costante con uno sporco che cala al calare delle dimensioni dell’elemento. | ||
+ | |||
+ | QUINDI quello che si sta vedendo è la componente x dello strain (allungamento lungo la direzione x) e si nota di come sia puramente dovuta alla quota elastica e non dalla elastica e plastica; infatti se si va a controllare la deformazione plastica (Comp 11 of Plastic strain) questa è omogeneamente nulla. | ||
+ | |||
+ | |||
+ | **Equivalent Von Mises strain** | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | L’ Equivalent Von Mises strain è omogenea e dà valori di 2.33 che indicano un piccolissimo snervamento. | ||
+ | Displacement x | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Notiamo come gli spostamenti in direzione x sono omogenei (con un piccolo sporco finale) lungo la sezione, ciò indica che tutti i punti di questa sezione si muovono in direzione assiale insieme. | ||
+ | |||
+ | __OSSERVAZIONI__ | ||
+ | |||
+ | La particolarità del tubo analizzato è che quest’ultimo è una struttura spessa che possiamo però modellare in tensione piana ; in questo caso non si è imposta una deformazione omogenea lungo la fettta ma, viceversa, una tensione assiale nulla ottenendo come risultato una deformazione omogenea in elasticità esatta. Nel metodo ad elementi finiti queste ipotesi valgono con i limiti della discretizzazione. Quindi il metodo di modellazione in tensione piana del tubo elastico è sostanzialmente valido e legato all’entità dell’errore. | ||
+ | |||
+ | |||
+ | ====Modellazione plastica=== | ||
+ | |||
+ | Proviamo a vedere cosa succede se usiamo una modellazione plastica: | ||
+ | Vado a modificare il modello quanto, per applicare una modellazione plastica, basta impostare un modello non lineare elastoplastico. | ||
+ | |||
+ | Per andare in stato plastico, vado a prendere dal libro di Progettazione assistita di organi meccanici (del Prof. Strozzi) ricavo pressione da applicare al bordo interno di tubo in funzione di quanto è la porzione di tubo che voglio plasticizzare. | ||
+ | Se applico 68.66Mpa come pressione interna al mio bordo non lineare, al bordo interno questo tubo verrà snervato fino a metà sezione ed avrà un raggio di transizione tra zona plastica e elastica dell’ordine di 15. | ||
+ | Nel tubo quindi ci sarà una zona plastica dal bordo interno fino ad un certo raggio di tranziozione e poi una zona elastica residua dal raggio di transizione al bordo esterno. | ||
+ | Questa formula calcola la pressione da un voluto raggio di transizione e non verrà esattamente uguale nel FEM un po per l’errore di discretizzazione ma anche perchè è vìcalcolato con tresca o $tao$ max mentre il marc usa la Von Mises. | ||
+ | Dobbiamo cambiare il nostro modello: | ||
+ | |||
+ | < | ||
+ | MAIN-> | ||
+ | e si cambiano le proprietà della p da 1 Mpa a 68.66 (dovrebbe dare snervamento fino a metà sezione). | ||
+ | </ | ||
+ | |||
+ | Se lancio un’analisi adesso : | ||
+ | < | ||
+ | JOB-> RUN-> | ||
+ | </ | ||
+ | La tensione secondo Von Mises risulta 153.3 Mpa nel bordo interno quindi ho una tensione secondo Von Mises che eccede il limite di elasticità. | ||
+ | |||
+ | Si può anche vedere quanto è l’**area al di sopra della plasticità** ,ovvero l’area che sfonda il muro dei 110, utilizzando una visualizzazione grafica: | ||
+ | |||
+ | OPEN POST FILE ed apriamo il file dei ridultati con la pressione maggiorata | ||
+ | metto la SCALAR: Equivalent Von Mises Stress | ||
+ | |||
+ | e guardiamo la mappa della Von Mises e vediamo che nel suo punto più elevato supera i 110MPa (che abbiamo supposto essere lo snervamento del materiale. | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | < | ||
+ | SCALARPOLOT-> | ||
+ | </ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Applicando una scala manuale che va da 0 a 110, possiamo isolare un fuori scala notando che tutti gli elementi che sono oltre i 110MPa in Equivalent Von Mises Stress sono evidenziati di grigio. | ||
+ | Nel nostro caso si nota che due file di elementi superano lo snervamento. | ||
+ | |||
+ | Per sapere qual’è il valore su cui è stato effettuato lo sweech elastico plastico da parte del solutore (ovvero sui punti di integrazione) vado nel menu: | ||
+ | < | ||
+ | EXTRAPOLATION-> | ||
+ | e disattivo il NODAL AVERAGING (la media nodale) cliccando su OFF | ||
+ | </ | ||
+ | |||
+ | Con il comando Translate prendo i quattro valori e li considero direttamente come valori nodali, quindi le mappe colorate son costituite coi valori dei punti di integrazione. | ||
+ | Togliendo la media nodale vedo esattamente i punti di integrazione dell’elemento. | ||
+ | Quindi se dovessi stimare un’area interessata dal fenomeno di plasticizzazione da questo calcolo lineare potrei dire che la plasticizzazione avviene solo sui due primi strati dell’ elemento in cui ho punti di integrazione che superano il valore dello snervamento, | ||
+ | |||
+ | Chiudo il file dei risultati ed imposto un vero e proprio **calcolo non lineare**. | ||
+ | Prima di tutto salvo con nome: | ||
+ | < | ||
+ | SAVE AS-> TUBO_ELASTICO | ||
+ | |||
+ | in seguito salvo nuovamente col nome: | ||
+ | |||
+ | SAVE AS-> TUBO_ELASTOPLASTICO | ||
+ | </ | ||
+ | |||
+ | ==TUBO ELASTOPLASTICO== | ||
+ | |||
+ | Attiviamo la plasticità della struttura così il sistema diventa elastoplastico e quindi disssipativo: | ||
+ | |||
+ | < | ||
+ | MAIN MENU-> | ||
+ | -> | ||
+ | </ | ||
+ | |||
+ | Così facendo applico uno snervamento di 110Mpa, mettendo solo uno snervamento ottengo solo un elastoplastico. | ||
+ | Se ci fosse incrudimento questo sarebbe isotropo e non cinematico (VEDI Appendice 1). | ||
+ | Adesso il modello è diventato non lineare, quindi il fenomeno è diventato non dissipativo e quindi lo stato della struttura è funzione della storia di carico. | ||
+ | |||
+ | Bisogna definire quindi una **storia di carico**: | ||
+ | < | ||
+ | MAIN MENU-> | ||
+ | </ | ||
+ | una table è in generale una funzione che ci permette di definire la storia del tempo; questa funzione è funzione del tempo e va a modulare il valore di 68.66 Mpa. | ||
+ | La funzione va tipicamente da 0 ad 1. | ||
+ | Per creare questa funzione modulante lo si fa andando dal menu TABLES creiamo una nuova tabella: | ||
+ | (si hanno fino a 4 variabili indipendenti) la funzione la chiamiamo “modula pressione”, | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Definiamo una storia di: carico scarico e ricarico ;ovvero definiamo una legge nella variazione del tempo della mia pressione inserendo 4 punti, in particolare quelli in cui la funzione modulante f(t): | ||
+ | -quando t = 0 la modulante vale 0; | ||
+ | -quando t=1 la modulante vale 1 (la modulante cresce fino a raggiungere un valore unitario al tempo 1); | ||
+ | -quando t=2 la modulante vale 0; | ||
+ | -quando t=12 la modulante vale 10; | ||
+ | Questa funzione si può esplicitare per “formula” o “data points”, noi scegliamo data points ed andiamo ad aggiungere | ||
+ | < | ||
+ | clicco DATA POINTS-> | ||
+ | (0,0) | ||
+ | (1,1) | ||
+ | (2,0) | ||
+ | (12,10) | ||
+ | </ | ||
+ | |||
+ | |||
+ | **NOTA BENE (1)**: per scalare il grafico e fargli contenere tutti i punti FILL | ||
+ | **NOTA BENE(2)**: l’ordine con cui si danno i punti non è importante. | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Dal grafico ottenuto notiamoche abbiamo fatto un caricamento, | ||
+ | |||
+ | Dopo aver fatto la tabella: | ||
+ | |||
+ | < | ||
+ | BOUNDARY CONDITION-> | ||
+ | </ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Per passare dalla visualizzazione della tabella a quella del modello : | ||
+ | |||
+ | < | ||
+ | SHORTCUTS-> | ||
+ | </ | ||
+ | E mi torna a mostrare il modello. | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Per adesso abbiamo definito la curva di carico, scarico ed esplosione ma non abbiamo definito il tempo della nostra analisi, in particolare le storie di carico si vanno a definire nel menu LOADCASES (storia di carico) e definiamo 3 pezzi di storia di carico: | ||
+ | -carico | ||
+ | -scarico | ||
+ | -continua a caricare | ||
+ | |||
+ | |||
+ | **Lodcase di carico** | ||
+ | < | ||
+ | LODCASES-> | ||
+ | </ | ||
+ | Le proprietà del pezzo di storia del carico dicono quali carichi e quali vincoli sono attivi in questo troncone di storia di carico, come si può vedere cliccando: | ||
+ | < | ||
+ | PROPERTIES-> | ||
+ | TOTAL LODCASE TIME->1 | ||
+ | </ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Il TOTAL LODCASE TIME indica quanto dura questo troncone di carico, nel nostro caso dura un secondo (46) | ||
+ | |||
+ | In seguito verranno analizzati i singoli comandi nella sezione STRUCTURAL STATIC LODCASE PROPERTIES e verranno fatte alcune piccole modifiche: | ||
+ | |||
+ | __SOLUTION CONTROL__ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Sono i parametri di Newton Rapson (NR), il NR è un procedimento iterativo in questo caso si vede nel numero massimo di iterazione che è 10 ed il minimo che è 0. | ||
+ | NOTA BENE con il termine RECICLE Marc intende “iterazioni di Newton Rapson” | ||
+ | |||
+ | -FULL NEWTON-RAPSON: | ||
+ | -MODIFIED NEWTON-RAPSON: | ||
+ | -N-R WITH STRAIN CORRECTION: non lo useremo | ||
+ | -NON POSITIVE DEFINITE: se attivato dichiariamo che la matrice è singolare, attiva anche una procedura automatica | ||
+ | -PROCEED WHEN NOT CONVERGED: Se dopo 10 iterazioni non è andato ad andare avanti continua fino a che non arriva a convergenza: | ||
+ | -ASSEMBLY EACH ITERATION: non lo useremo | ||
+ | |||
+ | < | ||
+ | Inseriamo: | ||
+ | MAX # RECYCLES -> 10 | ||
+ | MIN # RECYCLES -> 0 | ||
+ | </ | ||
+ | |||
+ | |||
+ | __CONVERGENCE TESTING:__ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Indica la convergenza di newton rapson quando i residui sono in norma più piccoli in una certa soglia: | ||
+ | -RELATIVE: | ||
+ | | ||
+ | | ||
+ | | ||
+ | -ABSOLUTE: | ||
+ | -RESIDUALS: il massimo valore assoluto in Newton che può avere residuo | ||
+ | |||
+ | < | ||
+ | RELATIVE-> | ||
+ | </ | ||
+ | |||
+ | __STEPPING PROCEDURE__ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | CONSTANT TIME STEP | ||
+ | |||
+ | -#STEPS | ||
+ | |||
+ | Quando si ha un sistema di carico non lineare e la storia di carco va percorsa a piccoli passettini | ||
+ | se metto 10 passettini vuol dire che la storia di carico viene divisa in 10 passi del 10%, se invece metto 40 vuol dire che applico un passo ogni 1/40 della storia di carico. | ||
+ | Ogni passo lancio un Newton Rapson ovvero da un istante i ad un istante i+1 comorta un’iterazione di newton rapson. | ||
+ | Questi istanti temporali sono chiamati incrementi, ogni passaggio da incremento i ad incremento i+1 comporta un certo numero di ricicli che sono le iterazioni del newton rapson. | ||
+ | Più sono gli istanti temporali e più i salti son piccoli e più la soluzione dell’incremento precedente è vicina alla soluzione dell’incremento successivo e quindi meno fatica fa Newton Rapson a convergere; viceversa se divido in 10000 incrementi, ho almeno 10000x2 ricicli di Newton Rapson (2 per ogni incremento) il che vuol dire che devo risolvere il sistema lineare associato a quella struttura 20000 volte; per questo noi usiamo pochi nodi, perchè ci danno un oggetto computazionalmente snello. | ||
+ | Bisogna fare attenzione perchè ogni riciclo di un calcolo non lineare costa come un calcolo lineare, quindi non devo mettere nè troppi nè pochi incrementi. | ||
+ | Nel problema dissipativo avere un passo fine sugli icrementi temporali vuol dire definire in modo preciso la storia di carico, essendo la soluzione in funzione della storia di carico, avere un’alta risoluzione nella risoluzione della storia di carico vuol dire avere una soluzione precisa. | ||
+ | |||
+ | **NOTA BENE** se si ha un dubbio si lancia a 40 steps e poi si lancia ad 80, se non cambia niente allora 40 va bene | ||
+ | |||
+ | Se Newton Rapson non converge dopo 10 tentativi riprovo Newton Rapson su un salto che è la metà (attraverso un cut back), l’intervallo di calcolo può essere 2^10 volte più piccolo dell’intervallo iniziale, se non converge esco con l’errore . | ||
+ | |||
+ | | ||
+ | |||
+ | -CONSTANT TIME STEP | ||
+ | Indica il tempo che si impone per percorrere un passettino | ||
+ | |||
+ | ARC LENGHT | ||
+ | Non ho step temporali imposti ma procedo per spostamenti costanti, ovvero ogni step voglio andare avanti per 1 mm e non mi interessa quanto ci vuole per fare 1mm, si usa per problemi istabili che non affrontiamo nel nostro corso. | ||
+ | |||
+ | < | ||
+ | Imponiamo: | ||
+ | # STEPS-> | ||
+ | </ | ||
+ | |||
+ | **LODCASE DI SCARICO E RICARICO FINO A ESPLOSIONE (PRESSURIZZAZIONE**) | ||
+ | |||
+ | < | ||
+ | copio due volte (ciccando copy) il lodcase di carico. | ||
+ | |||
+ | Il carico e scarico sono due lodcases gemelli, l’unica cosa che modifico è nel terzo impongo | ||
+ | # STEPS: 40 | ||
+ | </ | ||
+ | |||
+ | In seguito: | ||
+ | |||
+ | < | ||
+ | MAIN MENU-> JOBS-> | ||
+ | </ | ||
+ | D | ||
+ | isattivando la pressione interna sottolineo il fatto che all’istante 0 non applico carico e parto da scarico. | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | < | ||
+ | Clicco su AVAILABLE-> | ||
+ | </ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Selezionando la sequenza dei lodcases che voglio applicare nel mio job (carico, scarico | ||
+ | |||
+ | Infine clicco: | ||
+ | |||
+ | < | ||
+ | RUN->JOB | ||
+ | </ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Aspetto qualche secondo (Fino a che lo status non è “Complete”, | ||
+ | CURRENT INCREMENT: | ||
+ | |||
+ | ANALYSIS TIME:3.37 è il tempo su questa scala associate allo step 135, si può andare a vedere qual’e l’unità del tempo al quale si è impiantato il calcolo. | ||
+ | |||
+ | EXIT NUMBER: non è 3004 ma è 3009, la taglia del time step è diventata troppo piccola. | ||
+ | Tutti i risultati fino all’incremento 135 escluso sono stati scritti sul nostro t16. | ||
+ | |||
+ | Analizzando i vari risultati attraverso: | ||
+ | < | ||
+ | OPEN POST FILE | ||
+ | </ | ||
+ | |||
+ | Eequivalent Von Mises Strain | ||
+ | |||
+ | Per vedere l’area effettivamente snervata, lo faccio sogliando al di sotto dello snervamento, | ||
+ | < | ||
+ | PLOT-> | ||
+ | </ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | vediamo che 3 file di elementi sono effettivamente snervate (dal calcolo lineare ne avevamo 2 di file snervate) , la teoria mi dice che lo snervamento deve avvenire a metà del tubo. | ||
+ | |||
+ | **Displacement x** | ||
+ | |||
+ | < | ||
+ | DEFORMATION SCALING-> | ||
+ | SPOSTAMENTI ASSIALI (displacement x) | ||
+ | </ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Noto che in effetti l’espansione non è lineare, ho una parte plastica che tende a spanciare (incremento temporale di 76). | ||
+ | La fetta non mantiene espansioni costanti come il caso elastico, quindi nel caso elastoplastico non è più valida l’ipotesi che il tubo può essere modellato in tensione piana ma va modellato tenendo conto di una dilatazione differenziale in direzione assiale. | ||
+ | nell’asse y non è nulla ma ho uno spostamento radiale | ||
+ | |||
+ | **Component 11 of Stress(incremento 76)** | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | le componenti assiali di sforzo non sono nulle come nel tubo elastico. | ||
+ | |||
+ | **Von Mises** | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | all’incremento 134, il fronte plastico aumenta fino a che (in assenza di incrudimento) il materiale è interamente plasticizzato non ho più rigidezza residua. | ||
+ | |||
+ | **NOTA BENE** marc non lavora mai elastoplastico, | ||
+ | |||
+ | |||
+ | ====Appendice 1 ==== | ||
+ | |||
+ | **Caricamenti non assialsimmetrici** | ||
+ | |||
+ | Esistono infinite variazioni di questo elemento, per esempio si può assumere un caricamento assial simmetrico più torsione, a cui sono associato due gdl, sul piano, spostamento assiale e spostamento radiale, nonché una torsione fuori piano. | ||
+ | Nelle varianti possibili, si può deviare dall' | ||
+ | Si considerano 5 gdl per ogni nodo, di elemento piano. | ||
+ | - u spostamento radiale | ||
+ | - w spostamento assiale | ||
+ | inoltre | ||
+ | - u', w', spostamenti, | ||
+ | - v spostamento tangenziale valutato a 90° dal punto di valutazione dello spostamento radiale, necessario per definire spostamenti di tipo tagliante associati alla flessione. | ||
+ | Gli spostamenti di punti intermedi, sono scalati con legge cosinusoidale più costante, del tipo A+cosθ (che proiettato sul piano diviene una retta) sulla congiungente degli spostamenti assiali w, w', ciò consente di rappresentare i contributi flessionali, | ||
+ | |||
+ | |||
+ | {{: | ||
+ | |||
+ | Dunque elementi piani, in funzione della costruzione, | ||
+ | |||
+ | In questo elemento in particolare gli spostametni sono definiti da una quota costante più un contributo sinusoidale (nell' | ||
+ | Per esempio per una pressione P(X; | ||
+ | {{ : | ||
+ | Nel complesso, considerando un generico caricamento, | ||
+ | |||
+ | **Caricamento elastoplastico** | ||
+ | |||
+ | Lo stato tensionale di una struttura non risulta definito, esclusivamente dal carico istantaneo, ma sarà necessario conoscere la storia del caricamento, | ||
+ | |||
+ | |||
+ | **Incrudimento isotropo** | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Sebbene meno vicino al comportamento reale, questo modello è molto usato per la sua semplicità. | ||
+ | Nella fase di scarico si percorre il tratto elastico fino a un punto caratterizzato da un valore | ||
+ | simmetrico dello sforzo. Da lì in poi si segue la curva monotonica ribaltata. | ||
+ | |||
+ | Io parto da una funzione di scarico e comincio a caricare, parto dall' | ||
+ | |||
+ | **Incrudimento cinematico** | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | Sebbene meno vicino al comportamento reale, questo modello è molto usato per la sua semplicità. | ||
+ | Nella fase di scarico si percorre il tratto elastico fino a un punto caratterizzato da un valore | ||
+ | simmetrico dello sforzo. Da lì in poi si segue la curva monotonica ribaltata. | ||
+ | Tale modello riesce a descrivere l’effetto Bauschinger. | ||
+ | Avevo una forchetta dello stato tensionale da +Rs a – Rs e quindi ho un certo 2Rs del materiale elastico snervo in trazione poi incomincio a comprimere | ||
+ | |||
+ | **la teoria di Bauschinger** | ||
+ | |||
+ | Si definisce effetto Bauschinger (dal nome dell' | ||
+ | Il confine fra le due fasi è segnato da un valore di tensione di snervamento a trazione fy(t), ed uno a compressione fy( c ); supponendo di sottoporre un corpo ad uno stato di sforzo di trazione f'y(t) superiore a questo valore, dal momento in cui il materiale viene soggetto a scarico e successivamente ad una sollecitazione di compressione, | ||
+ | |||
+ | ====Esercizi svolti==== | ||
+ | |||
+ | {{: | ||
+ | |||
+ | {{: | ||
wikipaom2015/lez20.txt · Ultima modifica: 2015/09/30 11:22 da 182648