====== PROBLEMA LASTRA FORATA ====== Importiamo la geometria di una lastra da portare a trazione con tensione nominale 1MPa. Raggio di curvatura del foro 1 mm, larghezza del quarto di piastra 8 mm. Nella geometria importata è presente una suddivisione non legata alla geometria del pezzo, si tratta di un elemento AUX per costruire la mesh. Ci sono due modi per fare una mesh: si fa una mesh molto fine, l’operatore impiega meno tempo per impostare il problema, ma il calcolatore lo risolve più lentamente; oppure l’operatore impiega tempo per impostare il problema più precisamente e il calcolatore svolge il calcolo più velocemente. Utilizzeremo taglie non uniforme in base alle zone di interesse della geometra: per la zona nei pressi del foro useremo taglie differenziate, via via più fini man mano che ci avvicina al foro; mentre per le zone remote una taglia uniforme. Le curve della geometria importata non sono circonferenze esatte ma poligonali, quindi vanno rese più curvilinee. PLOT> CURVES-SETTINGS> HIGH>REGEN DIRECTION>REGEN LABELS>REGEN abbiamo selezionato High per definire il livello di resa della curv elevato, direction affinchè sia resa visibile la direzione delle curve, labels per assegnare un nome numerico progressivo ad ogni curva. {{:wikipaom2017:cattura.png?600|}} Dato che al punto non vengono associati gdl mentre ai nodi si, devo creare nodi ed elementi sulla geometria MESH GENERATION> AUTOMESH> CURVE DIVISION> APPLY CURVE DIVISION> seleziono la curva 14 END LIST# in curve division inserisco il numero di divisioni, in questo caso 6; {{:wikipaom2017:cattura3.png?600|}} Dopo aver creato i nodi creiamo gli elementi tipo linea sulle curve. MESH GENERATION> AUTOMESH> CURVE MESHING> MESH CURVES> seleziono la curva 14 END LIST# devo espandere gli elementi creati sulla curva 14 fino alla curva 21, nel menu expande, che trasforma elementi linea in elementi quadrilateri, ed elementi quadrilateri in elementi solidi (esaedri). è una espansione in scalatura, quindi devo definire parametri di espansione. {{:wikipaom2017:cattura4.png?600|}} devo espandere gli elementi per scalatura attorno al centro del foro, inserisco il numero di divisioni all'interno del foglio di calcolo dato al lezione, che ci restituisce i parametri n numero dei passi e lambda fattore di scala di ogni passo. ciò è necessario per far sembrare i quadrilateri più simili possibile a dei quadati ideali. MESH GENERATION> EXPANDE> REPETITION=6 CENTROIDE=0 0 0 SCALE FACTOR=X Y Z MODE-REMOVE ELEMENTS> seleziono gli elementi END LIST# con repetition decido il numero dei passi, centroide 0 0 0 perchè faccio partire la misura dell'origine degli assi, cioè il centro del foro scale factor-> X=Y=a^(1/b) Z=1 remove fa in modo che gli elementi generatori verranno rimossi successivamente. {{:wikipaom2017:cattura5.png?400|}} {{:wikipaom2017:cattura7.png?400|}} Devo asegnare delle proprietà ai lati, i lati dei quadrilateri esterni giacciono sulla curva 21 e 19. faccio ciò con il comando attach che associa i lati alle curve (con detach elimino le associazioni) MESH GENERATION> ATTACH> EDGES->CURVES seleziono la cuva 21 e i lati da associarle END LIST# {{:wikipaom2017:cattura8.png?600|}} Itero con le curve 14 19 28. associare curve ai lati sarà utile per associare le boundary condition e affinare la mesh. riprendo il foglio di calcolo facendo attenzione alla taglia dell'elemento a raggio esterno; la formula usata è specificata nello screen del foglio di calcolo che segue. {{:wikipaom2017:cattura9.png?600|}} risulta taglia 1.047, caratteristica della curva 34 e 30 nell'intorno dei punti di intersezione delle suddette curve con l'elemento AUX. MESH GENERATION> AUTOMESH> CURVE DIVISIONS> TYPE-VARIABLE L1->L2 inserisco i valori L1 e L2 APPLY CURVE DIVISIONS END LIST# {{:wikipaom2017:cattura10.png?600|}} la taglia L1=1.047*1.414=1.48 e L2=1.047 è caratteristica delle curve 30 e 34 la taglia L1=L2=1.48 è caratteristica delle curve 11 e 7; applico quest'ultima taglia a tutte le altre curve. andaimo a suddividere la curva 28 MESH GENERATION> AUTOMESH> CURVE DIVISIONS> TYPE-UNIFORM #DIVISIONS=6 APPLY CURVE DIVISION seleziono la curva 28 END LIST# creiamo gli elementi sull'area delimitata dalle curve 30 28 34 7 14 MESH GENERATION> AUTOMESH> 2-D PLANAR MESHING> QUAD MESH! seleziono il perimetro (curve 30 28 34 7 14) END LIST# {{:wikipaom2017:cattura13.png?600|}} Se il numero di suddivisioni è dispari ènormale che nella mesh compaiano elementi triangolari. itero il procedimento per gli altri perimetri chiusi. {{:wikipaom2017:cattura14.png?600|}} i nodi a cavallo tra due domini sono duplicati.devono collassare in un unico nodo se no quando faccio la deformata ottengo risultati incoerenti. MESH GENERATION> SWEEP> TOLERANCE=0.01 ALL devo ricordare di rimettere la tolleranza di default a 0.0001 GEOMETRIC PROPERTIES> NEW> STRUCTURAL> PLANAR> PLANE STRESS> name: pstress PROPERTIES> THICKNESS=1 OK ELEMENTS-ADD ALL:EXIST. con thickness associo uno spessore unitario per modellare una lastra unitaria {{:wikipaom2017:cattura16.2.png?600|}} MATERIAL PROPERTIES> MATERIAL PROPETRIES> NEW> STANDARD nome: alluminio STRUCTURAL> YOUNG'S MODULUS=70000 POISSON'S RATIO=0.3 ELEMENTS-ADD ALL:EXIST. Nel menu genera c'è la possibilità di indicare la densità, 2.9e-09. nel nostro caso non ci interessa metterla dato che non agiscono forze gravitazionali, inerziali o centrifughe; quindi lascio 0 come di default. {{:wikipaom2017:cattura18.png?600|}} inseriamo le condizioni al contorno, di tre tipi: -simmetria sul piano yz, normale a x, -simmetria sul piano xz, normale a y, -trazione nominale di 1 MPa nel lato sotto, ottenuta applicando una tensione remota di 0.875 di MPa BOUNDARY CONDITION> NEW> STRUCTURAL> FIXED DISPLACEMENT name:symm_nx_pyz PROPERTIES> DISPLACEMENT X OK CURVES-ADD seleziono le curve 2 3 4 30 21 END LIST# {{:wikipaom2017:cattura19.png?400|}} {{:wikipaom2017:cattura20.png?400|}} BOUNDARY CONDITION> NEW> STRUCTURAL> FIXED DISPLACEMENT name:symm_ny_pzx PROPERTIES> DISPLACEMENT Y OK CURVES-ADD seleziono le curve 19 34 END LIST# {{:wikipaom2017:cattura21.png?400|}} {{:wikipaom2017:cattura22.png?400|}} BOUNDARY CONDITION> NEW> STRUCTURAL> EDGE LOAD name:trazione_sigman_1MPa PROPERTIES> PRESSURE= -7/8 OK CURVES-ADD seleziono la curva 1 END LIST# {{:wikipaom2017:cattura23.png?400|}} {{:wikipaom2017:cattura24.png?400|}} lanciamo il calcolo JOBS> NEW> STRUCTURAL PROPERTIES> JOB RESULTS> AVAIABLE ELEMENT TENSORS-STRESS AVAIABLE ELEMENT SCALARS-Equivalent Von Mises Stress OK OK ELEMENT TYPE> ANALYSIS DIMENTIONS> PLANAR SOLID> PLANE STRESS FULL INTEGRATION : 3 OK ALL EXIST RUN> SUBMIT(1) OPEN POST FILE (RESULTS MENU)> SETTING> DEFORMATION SCALING:10000 SHOW FACTOR RETURN attenzione: ricordare di chiudere il file dei risultati con CLOSE {{:wikipaom2017:cattura25.png?400|}} {{:wikipaom2017:cattura26.png?400|}} {{:wikipaom2017:cattura27.png?400|}} {{:wikipaom2017:cattura28.png?400|}} {{:wikipaom2017:cattura29.png?400|}} {{:wikipaom2017:cattura30.png?400|}} negli ultimi screen si può notare che la deformata si presenta più lunga di quella originale e che il foro ovalizza. ora posso raffinare modello e mesh, lo faccio duplicando il modello MESH GENERATION> DUPLICATE> TRANSLATIONS 10 0 0 REPETITIONS 2 DUPLICATE ALL EXIST. {{:wikipaom2017:cattura31.png?600|}} con la duplicazione non ho duplicato solo le geometrie, ma anche le boundary conditions. ho quindi tre lastre: la prima che lascio come riferimento: bilineare a tratti; la seconda verrà raffinata per suddivisione degli elementi; la terza verrà raffinata per grado polinomiale dell'elemento, quindi passo ad una biquadratica a tratti; avrò gli stessi domini, ma auento i gradi del polinomio sul dominio stesso. MESH GENERATION> SUBDIVIDE> DIVISION 2 2 1 ELEMENTS seleziono tutto il secondo elemento END LIST# {{:wikipaom2017:cattura32.png?600|}} ora ci occupiamo del terzo elemento MESH GENERATION> CHANGE CLASS> TO QUADRILATIC ELEMENT seleziono il terzo elemento END LIST# {{:wikipaom2017:cattura33.png?600|}} MESH GENERATION> SWEEP> ALL SAVE lo sweep ha collassato i nodi vicini: in totale nel nostro caso 832 {{:wikipaom2017:cattura34.png?600|}} JOBS> RUN> SUBMIT(1) OPEN POST FILE (RESULTS MENU) {{:wikipaom2017:cattura35.png?400|}} {{:wikipaom2017:cattura36.png?400|}} nell'ultima immagine si può notare l'equivalent von mises stress, però non sappiamo a quale figura si riferiscono i valori della legenda; quindi procederemo per identificare i valori delle tre figure per valutare l'errore relativo. sulla seconda figura, nel post processing SELECT> ELEMENTS> seleziono l'elemento di interesse edalcuni intorno END LIST# MAKE VISIBLE la scala della Von Mises è riferita adesso a quei 4 elementi {{:wikipaom2017:cattura38.png?600|}} rifaccio la stessa cosa per la terza figura {{:wikipaom2017:cattura39.png?600|}} riporto i dati in un foglio di calcolo e riporto l'errore relativo secondo la formula [(BASE)-(href)]/(href) risulta un errore del 5%