====== Lezioni Master Aprile 2018 ====== mentat2013.1 -ogl -glflush ===== Telaietto simil-FSAE ===== {{ :wikitelaio2017:telaio_monocoque_2017_v008.mud | modello di partenza}} {{ :wikialtro2018:telaio_monocoque_2017_v010.mud | modello a tarda lezione}} {{ :wikifemfuchde2018:telaio_monocoque_2018_v001.mud |}} {{ :wikifemfuchde2018:telaio_monocoque_2018_solomasse.mfd | inertial elements}} {{ :wikialtro2018:telaio_monocoque_2018_v002.mfd |}} {{ :wikialtro2018:telaio_monocoque_2018_v003.mfd |}} :!: {{ :wikialtro2018:telaio_monocoque_2018_v004.mfd |}} Discussione: * elementi shell (elem. 75), beam(elem. 14), truss(elem. 9); * vincoli cinematici interni RBE2, RBE3, inserts; Loadcases: * Prova statica di rigidezza torsionale; * verifica preliminare modellazione cinematismo ({{:wikitelaio2016:telaio_monocoque_comeimportato_tutorial_optistruct.mud|versione originale tutorial Optistruct}} m( ); * discussione ruolo barra antirollio; * Bump su anteriore dx (inertia relief); * Risposta dinamica: analisi modale; * Risposta dinamica a sollecitazioni armoniche: forza verticale ant. dx, campionamento risposta su range 1-100 Hz in 99*4+1=397 steps (passo 1/4 Hz); [[https://www.youtube.com/watch?v=-LFLV47VAbI|back view]] [[https://www.youtube.com/watch?v=RihcJR0zvfM|side view]] effetto vincolamento improprio su comportamento dinamico di struttura. :!: {{:wikitelaio2016:inerzia_parallelepipedo_equivalente_motore.ods|}} :!: {{:wikitelaio2016:inerzia_parallelepipedo_equivalente_gruppi_mozzo.ods|}} :!: [[https://it.wikipedia.org/wiki/Inertanza]] :!: {{:wikitelaio2016:telaio_monocque_solomasse_ok.mfd|masse concentrate da inserire}} Foglio di calcolo di confronto tra i risultati in termini di inertanza di una risposta dinamica in bassa frequenza e la simulazione statica in //inertia relief//. :!: {{:wikitelaio2016:telaio_monocoque_inertanze_centroruota_antsx_z.ods|}} {{ :wikitelaio2017:plascore_crushlite.pdf |datasheet}} honeycomb per assorbimento impatti (idea [[https://www.google.com/patents/US3130819|non nuovissima]]) {{ :wikialtro2018:animacomp.proc | procedura per animare oscillazioni complesse}} ---- ===== Integrazione: risposta smorzata ===== Per inserire uno smorzamento strutturale (es. 1% del critico) occorre seguire i seguenti step: * entro nel menu menu ''MAIN -> MATERIAL PROPERTIES -> MATERIAL PROPERTIES''; * definisco preventivamente una table da menu ''TABLES'', ''NEW -> 1 INDIPENDENT VARIABLE'' * definisco ''NAME'' come ''modula_stiffmatmult'' * setto il ''TYPE'' della //Indipendent variable v1// a ''frequency'' * definisco la //table// per ''FORMULA'' e batto ''1/pi/v1'', ossia definisco una $g(f)=\frac{1}{\pi f}$ * torno nel menu menu ''MAIN -> MATERIAL PROPERTIES -> MATERIAL PROPERTIES'' con ''RETURN''; * seleziono qui il materiale "titanio", quindi entro nel menu ''STRUCTURAL -> DAMPING'' e attivo ''DAMPING''; * lascio a valore nullo lo ''MASS MATRIX MULTIPLIER'' * definisco uno ''STIFFNESS MATRIX MULTIPLIER'' pari allo smorzamento frazione del critico desiderato, nel caso specifico ''0.01'', * lo modulo per una //TABLE// cliccando sul menu ''TABLE'' che affianca //stiffness matrix multiplier// * scelgo la table ''modula_stiffmatmult'' appena definita, quindi do ''OK'' e ancora ''OK'' per tornare al menu material properties * in questo modo ho definito lo smorzamento fraz. del critico in funzione dei coefficienti $\alpha$ e $\beta$ del Rayleigh damping, supponendo nullo $\alpha$ e quindi il contributo della matrice massa alla matrice smorzamento. In pratica ho $\zeta = \frac{1}{2}(\frac{\alpha}{2 \pi f}+2 \pi f \beta)$ con $\alpha=0$ e $\beta= 0.01 \cdot g(f)=\frac{0.01}{\pi f}$, da cui $\zeta=0.01$ come desiderato. * passo quindi al menu ''MAIN -> JOBS'' e creo una copia del job di risposta in frequenza non smorzata //risposta// selezionandolo e utilizzando il comando ''COPY'' in alto a sx * rinomino in //rispostasmorzata// il nuovo job * entro nel menu ''PROPERTIES'', da cui seleziono ''ANALYSIS OPTIONS''; qui attivo ''COMPLEX DAMPING'' dal blocco //dynamic harmonic//, quindi esco con ''OK'' * vado quindi in ''JOB RESULTS'' e disattivo ''Stress'' e ''Equivalent von Mises stress'' * inserisco al loro posto da //AVAILABLE ELEMENT SCALARS// * '' Equivalent Real Harmonic Stress '', layers ''MAX & MIN'' * '' Equivalent Imag Harmonic Stress '', layers ''MAX & MIN'' * le //REAL HARMONIC// e //IMAG HARMONIC// equivalenti delle componenti di sollecitazione dell'elemento trave come da paragrafo successivo, , layers ''DEFAULT'', oltre al comune ''Beam Orientatio Vector'' * inserisco da //AVAILABLE ELEMENT TENSORS// * '' Real Harmonic Stress '', layers ''ALL'' * '' Imag Harmonic Stress '', layers ''ALL'' * Si procede quindi a lanciare il calcolo come di consueto da ''RUN -> SUBMIT'' e aprendo il file dei risultati con ''OPEN POST FILE (RESULTS MENU)'' * La deformata appare ora visualizzabile //con fase// entro il ciclo di oscillazione (vedere menu ''DEFORMED SHAPE SETTINGS''); nel caso senza smorzamento la fase poteva essere solo 0° o 180°, casistica rappresentabile mediante una variazione di segno delle componenti di spostamento o di tensione monitorate. Ricordo che la componente reale ha fase 0° (modulata in $\cos(\omega t)$) mentre la componente immaginaria ha fase 270° (modulata in $-\sin(\omega t)$). Ricordo inoltre che alla risonanza si amplifica fortemente la **componente immaginaria** della risposta, mentre si annulla quella reale (la risposta è infatti sfasata di ~90° rispetto all'eccitante). * Raccogliamo a titolo di esempio lo spostamento in direzione $z$ del nodo al centro impronta a terra della ruota (o equivalentemente al centro ruota): * dal menu POSTPROCESSING ''RESULTS'', con file dei risultati t16 aperto, procedere entro il menu ''HISTORY PLOT'' * definire il punti di campionamento con ''SET LOCATIONS'', fornendo quindi al prompt ''146'' [invio], seguito da un ''END LIST'' * definire il range di sottoincrementi di campionamento da ''INC RANGE'', fornendo quindi al prompt ''0:1'' [invio], ''0:397'' [invio], ''1'' [invio], come sottoincrementi di inizio, fine e passo di campionamento. * procedere a compilare diagrammi sulla base dei dati appena campionati dal menu ''ADD CURVES'', quindi ''ALL LOCATIONS'' (ho selezionato un solo punto di campionamento); richiedo la compilazione di un grafico che abbia come asse delle ascisse la global variabl ''Frequency'', e come asse delle ordinate ''Displacement Z Magnitude''. Con ''FIT'' adatto le scale del grafico alla curva. * Con ''RETURN'' torno al menu HISTORY PLOT, ove posso ridurre la frequenza delle etichette indicanti l'incremento portando ''SHOW IDS'' da '1' a '10'; inserendo valore '0' ometto la visualizzazione delle etichette. * I picchi di risposta alle risonanze appaiono ora finiti (non erano limitati in assenza di smorzamento), ed è ora visibile che alcuni modi propri, pur non essendo strettamente ortogonali all'eccitante, risultano scarsamente accoppiati e facilmente contenuti da un ridotto smorzamento strutturale. ==== Riferimenti per valutazione damping strutturale ==== {{ :wikipaom2017:structural_damping_values_jdstevenson.pdf |}} {{ :wikipaom2017:damping_cross-reference_and_material_properties.pdf |}} {{ :wikipaom2017:f_orban_damping_of_materials_and_members_in_structures.pdf |}} {{ :wikipaom2017:tom_irvine_damping_in_bolted_and_welded_joints.pdf |}} {{ :wikipaom2017:estratto_afwal-tr-84-3089v2_dampingvalues.pdf |estratto}} vol. 2, sezione 8 di Soovere, J., and M. L. Drake. Aerospace Structures Technology Damping Design Guide.LOCKHEED-CALIFORNIA CO BURBANK, 1985. ===== materiale didattico spiccio ===== {{ :dispensa_progettazione_assistita.pdf | Appunti per il corso di Progettazione Assistita}}; contiene note su //servo-links//, //RBE2//, //RBE3//, algoritmo di Newton-Raphson {{ ::dispensa_matrice_massa_bozzabozzabozza.pdf | Dispensa matrice massa}} {{ :wikitelaio2016:mia_dispensa_dinamica_updated.pdf | Dispensa dinamica updated}} {{ :wikifemfuchde2018:dispensa_2018_v2018-06-07.pdf | WIP dispensa FFCD2018}} ===== RBE3 (vs. RBE2) ===== [[https://cdm.ing.unimo.it/files/progettazione_assistita/corso_2012_2013/esempi_vari/dimostratore_RBE2_vs_RBE3.mud | Dimostratore RBE2 vs RBE3]] immagine cinematica RBE3. {{:wikipaom2015:cinematica_rbe3.png?300|}} {{:wikipaom2015:cinematica_rbe3_v1.pdf|sorgente ipe}} ===== Singolarità stato tensionale - cuscinetti ===== {{:wikipaom2016:sinclair_i.pdf|}} {{:wikipaom2016:sinclair_ii.pdf|}} {{:wikipaom2015:singolarita_cuscinetti_tensioni.gif?400|}} {{:wikipaom2015:singolarita_cuscinetti_cedimenti.gif?400|}} {{:wikipaom2016:singolarita_cuscinetti_dubbia_modellazione.mud| modello marc/mentat}} ===== Grattacapi ===== {{ :wikipaom2017:flangia_con_braccio_di_piatto_ovvero_trova_le_differenze.mfd |}} ===== Instabilità flesso-torsionale trave I ===== {{ :wikipaom2018:es_instabilita_flextors_raffinato_iniziolez.mud | solo mesh}} * acciaio S355JR, Rs=355 MPa; * campata: 1120 mm; * altezza profilo: 120 mm (da piano medio a piano medio ali); * larghezza profilo: 40 mm; * spessore anima: 2 mm; * spessore ali: 4 mm; * spessore fazzoletti di rinforzo ai supporti: 4mm; * carico: 50kN; :!: {{ :wikialtro2018:es_instabilita_flextors_raffinato_master.mud | mezzo svolto}} {{ :wikipaom2018:es_instabilita_flextors_raffinato.mud | svolto}} perturbazione: spostamento * x -> x * y -> y +x*z/40/1200*(1.0) * z -> z agli elementi dell'ala superiore (ala superiore lievemente torta, slope max 1.0/40) {{ :wikipaom2018:es_instabilita_flextors_raffinato_perturbato.mud | svolto, forma perturbata}} ===== Altro materiale ===== [[https://box.ing.unimo.it/owncloud/index.php/s/So3nMdJOCaQvdhC|Corso strutture per Maserati - febbraio 2017]]