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.
====== Marc/Mentat - Lezione 01/04/2015 ======



In questa lezione vediamo la distorsione di un elemento isoparametrico e verifichiamo se il software riesce a risolvere il problema. (controllo che lo Jacobiano converta bene l'elemento finito)

Teoria della trave su Marc Mentat; modello della trave 2D.

Realizzeremo in Marc un elemento quadrangolare e lo vincoliamo isostaticamente con una cerniera e un carrello.

Applicheremo un carico di 1N verticale sul nodo (1;1;0).

===== Case Study 1: elemento distorto =====

==== Creazione Mesh ====

//Inseriamo le coordinate dei nodi dell'elemento finito
//

''Mesh generation -> Nodes -> add : (0;0;0), (0,3;0,6;0), (1;1;0), (0;1;0)''

//Controllo che su "Element class" sia attivo "Quad(4)"
//

''Mesh generation -> elements -> add'' : seleziono i nodi in senso antiorario per avere circuitazione positiva.

{{:wikitelaio2015:schermata_2015-04-17_alle_15.23.51.png|}}

//Definiamo la natura dello stato tensionale (tensione piana)
//

==== Proprietà geometriche e del materiale ====


''Geometric properties -> new -> structural -> planar -> plane stress''
 
(Diamo il nome "tensione_piana")

{{:wikitelaio2015:screenshot-2.png|}}

//Definiamo il materiale
//

''Material properties  -> material properties  -> new -> standard -> structural -> (E = 70000 MPa  Poisson ratio = 0.33 )''

(Diamo il nome "Alu")

{{:wikitelaio2015:screenshot-3.png|}}

N.B. ulteriori informazioni sui materiali sul sito http://www.matweb.com/index.aspx 

==== Boundary conditions ====


**Vincolamento e caricamento dell'elemento finito.
**

Cerniera

''Boundary conditions -> new -> structural -> fixed displacement -> properties -> (displacementX=0 displacementY=0) -> nodes -> add: (0;0;0)''

(Diamo il nome "cerniera") 

{{:wikitelaio2015:screenshot-4.png|}}

Carrello

''Boundary conditions -> new -> structural -> fixed displacement -> properties -> (displacementX=0) -> nodes -> add: (0;1;0)''

(Diamo il nome "carrelloX") 


Forza

''Boundary conditions -> new -> structural -> point load -> properties -> (forceY=-1N) -> nodes -> add: (1;1;0)''

(Diamo il nome "forza") 

{{:wikitelaio2015:screenshot-7.png|}}

**Controllo della circuitazione dell'elemento
**

''Mesh generation -> check -> upside-down'' (risulta nessun elemento invertito)

**Controllo della distorsione dell'elemento (verifica della Jacobiano)
**

''Mesh generation -> check -> distorted'' (risulta un elemento distorto)

{{:wikitelaio2015:screenshot-9.png|}}

**Controllo "aspect ratio"**: quanto l'elemento è aghiforme.

Aspect ratio: rapporto perimetro-area.

''Mesh generation -> check -> aspect ratio'' (risulta zero elementi aghiformi) (threshold=0.5)

{{:wikitelaio2015:screenshot-10.png|}}


==== Job ====


//Lanciamo la simulazione
//

''Jobs -> new -> structural -> properties (facciamo attenzione che in "analysis dimension" sia selezionato "plane stress") ->  job result -> plane stress (seleziono: "stress", "elastic strain", "equivalent Von Mises stress")''

{{:wikitelaio2015:screenshot-11.png|}}

In "output file" seleziono "element matrices" (chiedo al solutore di inserire nel file di output le matrici di rigidezza degli elementi, così possiamo controllarle)

{{:wikitelaio2015:screenshot-12.png|}}

''Run -> submit -> monitor''

In "exit number" invece di 3004 (calcolo corretto), ottengo 13 (non riesce a trasformare l'elemento in isoparametrico).

{{:wikitelaio2015:screenshot-14.png|}}

''Plot -> solid element -> regen''

{{:wikitelaio2015:screenshot-15.png|}}

Fatto ciò possiamo tramite il manipolatore Maxima ricavare lo Jacobiano dell'elemento finito analizzato.


===== Case Study 2: elemento trave =====



//Analizziamo ora l'elemento trave e confrontiamo la soluzione ottenuta tramite il software con la soluzione analitica (teorema di Castigliano).
//

**Teorema di Castigliano**: In un sistema elastico, la derivata parziale dell'energia totale elastica rispetto ad una forza (o coppia) definisce lo spostamento (o rotazione) del punto di applicazione nella direzione di applicazione della forza (o coppia).

{{:wikitelaio2015:fullsizerender.jpg|}}


Materiale: alluminio

Carico: P=5N

Lunghezza: L=400mm

Sezione rettangolare b * l = 30mm * 6mm

{{:wikitelaio2015:schermata_2015-04-17_alle_16.56.12.png|}}

Verifichiamo il risultato in Mentat

==== Creazione Mesh ====


//Realizziamo l'elemento trave
//

Mesh generation -> nodes -> add: (400;0;0),(0;0;0).

Controllo che su "Element class" sia attivo "Line(2)"

''Mesh generation -> elements -> add'' : seleziono i nodi

''Subdivide -> divisions (100,1,1) bias factor (1,1,1)''

{{:wikitelaio2015:screenshot-19.png|}}

''Mesh generation -> sweep -> sweep all'' risulta una cancellazione di due nodi agli estremi.

Con gli elementi shell dobbiamo fare attenzione alle proprietà geometriche.

In questo caso non è soltanto planar ma va selezionato "straight beam" per rappresentare la teoria della trave dritta (nel software l'elemento è "element5").

Ogni nodo della trave ha 3 DOF: traslazioni in x e y, rotazione in z.

Teoricamente potrei dare un incastro. (nel Tria(3) non potevo)

Dimensioni della matrice di rigidezza 3*n X 3*n (con n numero di nodi = 101, 3 rappresenta il numero di DOF).

{{:wikitelaio2015:screenshot-20.png|}}

{{:wikitelaio2015:screenshot-26.png|}}

==== Proprietà geometriche e del materiale ====


Definiamo le dimensioni della sezione.

Il programma richiede altezza e area. 

''Geometrical properties -> new -> structural -> planar -> straight beam -> properties'' (inserisco altezza e area)

{{:wikitelaio2015:screenshot-21.png|}}

Materiale: alluminio (alu)

''Material properties  -> material properties  -> new -> standard -> structural -> (E = 70000 MPa  Poisson ratio = 0.33 )''

==== Boundary conditions ====


**Vincolamento e caricamento della trave
**

Incastro

''Boundary conditions -> new -> structural -> fixed displacement -> properties -> (displacementX=0 displacementY=0 displacementZ=0) -> nodes -> add: (0;0;0)''

N.B. Per bloccare la rotazione in Z seleziono "displacementZ" al posto di "rotationZ" in quanto avrei una deformata sbagliata a causa di un errore interno al software.

Forza

''Boundary conditions -> new -> structural -> point load -> properties -> (forceY=-5N) -> nodes -> add: (400;0;0)''

(Diamo il nome "forzaY") 

{{:wikitelaio2015:screenshot-25.png|}}


==== Job ====


Lancio la simulazione

''Jobs -> new -> structural -> properties (facciamo attenzione che in "analysis dimension" sia selezionato "2D") ->  job result -> plane stress (seleziono: "stress", "elastic strain", "equivalent Von Mises stress", "Beam axial force", "Beam bending xy", "Beam shear xy", "Beam torsional", "Beam bimoment")''

{{:wikitelaio2015:screenshot-29.png|}}

Accanto a "Stress" c'è selezionato "default"; seleziono "out & mid" in quanto indago lo stato tensione agli estremi e sull'asse neutro.

Lo applico solo allo stato tensionale essendo indifferente per gli altri parametri.

{{:wikitelaio2015:screenshot-30.png|}}

''Renumber-all'': rende bandata la matrice di rigidezza.

''Run -> submit -> monitor -> open post file -> more -> maximum principal value stress -> on'' 

{{:wikitelaio2015:screenshot-33.png|}}

{{:wikitelaio2015:screenshot-32.png|}}

I valori numerici calcolati dal software (anche la deformata della trave all'estremità libera) corrispondono ai valori calcolati con la teoria di Castigliano con le formule sopra riportate.