Strumenti Utente

Strumenti Sito


wikipaom2018:lez_2018-05-25

Differenze

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

Link a questa pagina di confronto

Entrambe le parti precedenti la revisione Revisione precedente
Prossima revisione
Revisione precedente
wikipaom2018:lez_2018-05-25 [2018/05/28 11:53]
244399
wikipaom2018:lez_2018-05-25 [2018/05/28 17:33] (versione attuale)
244399
Linea 1: Linea 1:
 +A cura di Giancarlo Degli Esposti
 +
 +======ESERCITAZIONE MARC.MENTAT: TUBO ELASTOPLASTICO======
 +
 +===== Introduzione =====
 +
 +Il problema di contatto oggetto della presente esercitazione è costituito da un tubo elastoplastico caricato in mezzeria di cui è stata modellata un quarto della struttura complessiva (oggetto dell'esercitazione precedente) sfruttando le proprietà di simmetria del sistema. Al seguente link: [[https://cdm.ing.unimo.it/files/index.php?dir=progettazione_assistita_2018/Tubo_elastoplastico]] è possibile scaricare il file //sigma_epsilon_AISI_4130.txt// che sarà necessario nell'impostazione del caso. 
 +
 +{{ :wikipaom2018:sistema_paom_25-5-18_1_.jpg?400 |}}
 +
 +===== Impostazione delle condizioni di contatto =====
 +
 +Si determinano i corpi che saranno a contatto: 
 +
 +  * il tubo impostato come //MESHED (DEFORMABLE)//
 +
 +  * i due cilindri impostati come //GEOMETRIC//
 +
 +Per il tubo si procede nel seguente modo:
 +
 +dal ''MAIN'' menu: ''CONTACT -> CONTACT BODIES -> NEW -> MESHED (DEFORMABLE) -> NAME: tubo -> ELEMENTS: ADD -> ALL: EXIST.''
 +
 +Per il primo cilindro:
 +
 +''CONTACT -> CONTACT BODIES -> NEW -> GEOMETRIC -> NAME: cilindro_mobile -> 3-D: SURFACES: ADD'' e si seleziona la superficie esterna del cilindro mobile che carica la struttura al centro; cliccare, infine, su ''END LIST (#)''
 +
 +
 +Per il secondo cilindro:
 +
 +
 +''CONTACT -> CONTACT BODIES -> NEW -> GEOMETRIC -> NAME: cilindro_fisso -> 3-D: SURFACES: ADD'' e si seleziona la superficie esterna dell'altro cilindro (quello all'estremità libera del tubo che simula l'appoggio); cliccare, infine, su ''END LIST (#)''
 +
 +Per verificare che quanto fatto sia corretto si possono evidenziare i corpi a contatto premendo su ''ID CONTACT'' dal menu ''CONTACT BODIES''
 +
 +
 +Per visualizzare meglio il sistema è consigliabile visualizzare le superfici esterne come solide: 
 +
 +
 +''PLOT -> SURFACES: SOLID -> REGEN''
 +
 +{{ :wikipaom2018:contact_bodies.png?600 |SURFACES: SOLID}}
 +
 +Tramite il comando ''CONTACT BODIES -> ID BACKFACES'' è possibile vedere che le faccia interna della superficie dei due cilindri è interpretata dal codice come //outside// e quella esterna come //inside//. Per ribaltare le facce cliccare dal menu ''CONTACT BODIES'' su ''FLIP SURFACES -> ALL: VISIB.''
 +
 +{{ :wikipaom2018:contact_id_backfaces.png?600 |FLIP SURFACES}}
 +
 +Ora è necessario definire nel codice quale dei cilindri è quello fisso e quale quello mobile, cioè quello che simulerà il caricamento del tubo: si torna a visualizzare il primo corpo di contatto (cilindro_mobile) cliccando su ''PREV'' e si seguono i passaggi descritti:
 +
 +
 +''PROPERTIES -> BODY CONTROL:'' al posto di '' VELOCITY'' selezionare ''POSITION'' e cliccare su ''PARAMETRES -> POSITION (CENTER OF ROTATION): y:'' e digitare  ''-15''
 +
 +{{ :wikipaom2018:screen3.png?600 |BODY CONTROL: POSITION}}
 +{{ :wikipaom2018:screen4.png?600 |POSITION (CENTER OF ROTATION): y: -15}}
 +
 +A questo punto è necessario definire come il cilindro mobile si avvicina al tubo poiché, se non specificassimo nessuna legge di moto, il contatto sarebbe istantaneo e produrrebbe errori nei risultati. Si impone, quindi, una legge di moto, mediante una tabella, secondo la quale lo spostamento del cilindro superiore avvenga in modo graduale così come il contatto.
 +
 +Dal menu ''CONTACT BODIES'' (deve essere visualizzato ''cilindro_mobile'') cliccare su 
 +
 +''TABLES -> NEW -> 1 INDEPENDENT VARIABLE -> TYPE -> time -> ADD'' digitare quindi le coordinate dell'origine del diagramma che rappresenta la legge di moto lineare: ''0 (invio) 0 (invio)'' e successivamente quelle del punto finale del diagramma: ''1 (invio) 1 (invio)''
 +
 +{{ :wikipaom2018:screen5.png?600 |1 INDEPENDENT VARIABLE -> TYPE -> TIME}}
 +
 +Si carica, infine, la tabella appena creata nelle proprietà del cilindro mobile:
 +
 +Dal menu ''CONTACT BODIES'' (deve essere visualizzato ''cilindro_mobile'') cliccare su 
 +
 +''PROPERTIES -> PARAMETRES -> POSITION (CENTER OF ROTATION): y: TABLE ->'' selezionare la tabella creata in precedenza e premere ''OK'' due volte.
 +
 +Si impostano ora le //CONTACT INTERACTIONS// tra i due cilindri e il tubo:
 +
 +Dal menu ''MAIN -> CONTACT'' cliccare su: ''CONTACT INTERACTIONS -> NEW -> MESHED (DEFORMABLE) vs. GEOMETRIC -> NAME: tubo_vs_fisso -> PROPERTIES -> CONTACT TOLERANCE -> VALUE 0.01 -> FRICTION -> FRICTION COEFFICIENT: 0.2''
 +
 +{{ :wikipaom2018:screen6.png?600 |CONTACT TOLERANCE -> VALUE 0.01}}
 +{{ :wikipaom2018:screen6.1.png?600 |FRICTION -> FRICTION COEFFICIENT: 0.2}}
 +
 +Copiare la stessa CONTACT INTERACTION premendo su ''COPY'' e rinominarla in //tubo_vs_mobile//
 +
 +Successivamente dal menu ''CONTACT'' cliccare su: ''CONTACT TABLES -> NEW -> PROPERTIES'' e in corrispondenza della riga 1 (tubo) cliccare sul pulsante ''2'' cliccare, quindi, su ''CONTACT INTERACTION -> ACTIVE -> tubo_vs_mobile -> OK''
 +
 +{{ :wikipaom2018:screen7.png?600 |in corrispondenza della riga 1 (tubo) cliccare sul pulsante ''2''}}
 +{{ :wikipaom2018:screen8.png?600 |CONTACT INTERACTION -> ACTIVE -> tubo_vs_mobile -> OK}}
 +
 +In corrispondenza della riga 1 (tubo) cliccare sul pulsante ''3'' cliccare, quindi, su ''CONTACT INTERACTION -> ACTIVE -> tubo_vs_fisso -> OK''
 +
 +
 +==== Diagrammi $\sigma-\epsilon$ ====
 +
 +
 +Si rende ora necessario impostare nel codice il comportamento elastoplastico del materiale.
 +I codici FEM non utilizzano la curva σ-ε ingegneristica perché non tiene conto della strizione del materiale all'aumentare della tensione; vengono invece implementate le curve cosiddette //true// di cui si riporta un esempio:
 +
 +{{:wikipaom2018:sigma_epsilon_1_.jpg?600 |Diagramma ingegneristico vs. diagramma vero}}
 +
 +
 +Il diagramma vero, tiene conto della strizione del materiale, cioè della riduzione dell'area della sezione su cui agisce la forza che genera la $\sigma = F/A$.
 +La deformazione elastica $\epsilon_e$ è sempre calcolabile come:
 +
 +$$\epsilon_e=\frac{\sigma}{E}$$
 +
 +mentre, la deformazione plastica è ricavabile dalla precedente:
 +
 +$$\epsilon_p=\epsilon_{tot}-\epsilon_e=\frac{\sigma_{e}}{E}$$
 +
 +Il diagramma che usano i software FEM per modellare il comportamento plastico è, dunque, il seguente:
 +
 +{{:wikipaom2018:sigma_epsilon_2_1_.jpg?600 |Diagramma usato dal codice}}
 +
 +Tale diagramma viene caricato nel codice tramite una tabella, motivo per cui viene fornito il file di testo di cui sopra.
 +
 +Per caricare la tabella si impostano le proprietà del materiale:
 +
 +Dal ''MAIN'' menu, si entra in ''MATERIAL PROPERTIES -> MATERIAL PROPERTIES -> TABLES -> READ -> RAW'' si seleziona il file //sigma_epsilon_AISI_4130.txt// ''TYPE -> eq_plastic_strain''
 +
 +{{:wikipaom2018:screen8.5.png?600 |TABLES -> READ -> RAW}}
 +
 +Dal menu ''MATERIAL PROPERTIES'' (deve essere visualizzato ''acciaio'') si clicca su: ''STRUCTURAL -> PLASTICITY -> PLASTICITY -> YIELD STRESS: 1 -> TABLE -> sigma_epsilon_plast''
 +
 +Dal grafico si può notare che la tensione parte da 320 MPa, cioè dallo snervamento. Per uno stato di tensione triassiale è necessario impostare uno //Yield criterion//:
 +
 +{{:wikipaom2018:screen9.png?600 |la tensione parte da 320 MPa}}
 +
 +Dal menu delle proprietà plastiche: ''YIELD CRITERION -> VON MISES'' e come criterio di incrudimento si seleziona: ''ISOTROPIC''.
 +
 +Dal menu ''MATERIAL PROPERTIES'' selezionare tutti gli elementi esistenti come elementi a cui applicare le proprietà appena definite: ''ELEMENTS -> ADD -> ALL: EXIST.''
 +
 +==== Incrudimento ====
 +
 +
 +Si ricorda che l'incrudimento è una trasformazione (indurimento) che il materiale subisce non appena entra in campo plastico.
 +In assenza di evidente snervamento si adotta come carico di snervamento convenzionale $R_{p0.2}$, cioè quel carico che induce una deformazione percentuale dello 0.2%.
 +
 +{{ :wikipaom2018:rp02.jpg?400 |deformazione percentuale dello 0.2%}}
 +
 +Esistono due tipi principali di incrudimento: isotropico e cinematico.
 +
 +  * Incrudimento isotropico: una volta plasticizzato il materiale fino al punto 1 e dopo che viene scaricato, questo presenterà una deformazione permanente corrispondente al punto 2. Se caricato nuovamente, il materiale seguirà la deformazione elastica dal punto 2 al punto 1 e, aumentando ulteriormente il carico, proseguirà la deformazione plastica. Se il materiale, una volta raggiunto il punto 2, viene, invece, caricato in direzione opposta (compressione) il materiale plasticizza quando il tratto elastico in compressione è uguale a quello in trazione.
 +
 +  * Incrudimento cinematico: il materiale, quando viene caricato in compressione, snerva molto prima del punto 1', segue una curva di plasticizzazione speculare a quella che si avrebbe in trazione. Se viene interrotta la deformazione prima della rottura e il provino viene caricato in direzione opposta, il diagramma presenterà un andamento di tipo isteretico come in figura.
 +
 +{{ :wikipaom2018:sigma_epsilon_isotropico_vs_cinematico.jpg?600 |DIAGRAMMA INCRUDIMENTO ISOTROPICO VS. CINEMATICO}}
 +
 +Gestione da parte del codice dell'incrudimento isotropico e cinematico: 
 +
 +  * Incrudimento isotropico: fino a quando la σ sta nell'intervallo $(-R_{p0.2}, R_{p0.2})$, il materiale non snerva. Nel momento in cui il la tensione dovesse superare il carico di snervamento questo viene "agganciato" e aumentato in trazione e, specularmente, anche in compressione.
 +
 +  * Incrudimento cinematico: quando viene superato il carico di snervamento, il valore di $R_{p0.2}$ viene aggiornato in aumento come nel caso dell'incrudimento isotropico, ma il carico di snervamento in compressione viene abbassato in modo da mantenere costante l'ampiezza dell'intervallo $(-R_{p0.2}, R_{p0.2})$
 +
 +==== Ellissoide di Von Mises e diagramma di Mohr ====
 +
 +{{:wikipaom2018:ellissoidi_isotropico_vs_cinematico.png?400 |Ellissoide di Von Mises}}
 +
 +Nei codici numerici, viene implementato l'ellissoide di Von Mises poiché il bordo dell'ellisse è continuo e derivabile in ogni punto al contrario del diagramma di Mohr.
 +In caso di incrudimento isotropico, l'ellissoide si modifica in espansione, ma non in posizione; in caso di incrudimento cinematico l'ellissoide mantiene le stesse dimensioni ma trasla nel piano $\sigma_1 - \sigma_2$.
 +
 +===== Caricamento struttura e calcolo =====
 +
 +Si procede con la definizione delle condizioni di carico e di contatto ai fini del calcolo:
 +Dal ''MAIN'' menu nella sezione ''ANALYSIS'', cliccare su ''LOADCASES -> NEW -> STATIC -> PROPERTIES -> CONTACT -> CONTACT TABLE -> ctable1 -> OK''
 +
 +{{ :wikipaom2018:screen10.png?600 |CONTACT -> CONTACT TABLE -> ctable1 -> OK}}
 +
 +Dal sotto menu ''PROPERTIES'' del loadcase cliccare su ''CONVERGENCE TESTING -> DISPLACEMENT -> RELATIVE'' e successivamente su ''DISPLACEMENT'' impostare la ''RELATIVE DISPLACEMENT TOLERANCE'' a 0.01 e premere ''OK''.
 +
 +{{ :wikipaom2018:screen11.png?600 |CONVERGENCE TESTING -> DISPLACEMENT -> RELATIVE}}
 +
 +Dal sotto menu ''PROPERTIES'' del loadcase cliccare su ''SOLUTION CONTROL -> MAX # RECYCLES: 25 -> OK'' tornare nel menu ''PROPERTIES'' e impostare ''STEPS'' a 50.
 +
 +{{ :wikipaom2018:screen12.png?600 |SOLUTION CONTROL}}
 +
 +{{ :wikipaom2018:screen12.1.png?600 |MAX # RECYCLES: 25 -> OK}}
 +
 +{{ :wikipaom2018:screen12.2.png?600 |STEPS -> 50}}
 +
 +Si inseriscono ora i dati per il //JOB//:
 +
 +Dal ''MAIN'' menu: ''JOBS -> NEW -> STRUCTURAL -> PROPERTIES -> CONTACT CONTROL -> INITIAL CONTACT -> CONTACT TABLE -> ctable1 -> OK -> OK''
 +
 +{{ :wikipaom2018:screen12.51.png?600 |CONTACT CONTROL -> INITIAL CONTACT -> CONTACT TABLE -> ctable1 -> OK -> OK}}
 +
 +Dal menu ''STRUCTURAL JOB PROPERTIES'' cliccare su ''JOB RESULTS'' e selezionare le quantità da calcolare:
 +
 +Dall'elenco //AVAILABLE ELEMENT TENSORS//: ''Stress, Total Strain, Elastic Strain, Plastic Strain''
 +
 +Dall'elenco //AVAILABLE ELEMENT SCALARS//: ''Equivalent Von Mises Stress''
 +
 +Dall'elenco //AVAILABLE NODAL QUANTITIES//: ''Contact Normal Stress, Contact Normal Force, Contact Friction Stress, Contact Friction Force, Contact Status, Displacement, External Force, Reaction Force''
 +
 +{{ :wikipaom2018:screen13.png?600 |JOB RESULTS}}
 +
 +Dal menu ''STRUCTURAL JOB PROPERTIES'' nei //LOADCASES//, alla voce //AVAILABLE//, selezionare il loadcase creato in precedenza.
 +
 +Infine viene impostato l'attrito come Coulombiano: 
 +
 +Dal menu ''STRUCTURAL JOB PROPERTIES'' cliccare su ''CONTACT CONTROL -> FRICTION: TYPE: COULOMB BILINEAR (DISPLACEMENT) -> OK''
 +
 +{{ :wikipaom2018:screen12.5.png?600 |CONTACT CONTROL -> FRICTION: TYPE: COULOMB BILINEAR (DISPLACEMENT) -> OK}}
 +
 +infine lanciare il calcolo tramite: ''JOBS -> RUN -> SUBMIT''