/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/ /* [ Created with wxMaxima version 0.8.7 ] */ /* [wxMaxima: input start ] */ kill(all); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] pulisco memoria [wxMaxima: comment end ] */ /* [wxMaxima: comment start ] caratteristiche di sollecitazione [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Mf1 : C - P/2 * r * sin(phi) $ /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ N1 : P/2 * sin(phi) $ /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ T1 : P/2 * cos(phi) $ /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ Mf2 : C - P/2 * r $ /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ T2 : 0 $ /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ N2 : P/2 $ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] energia interna al tratto 1 [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ U1 : integrate ( ( Mf1^2/2/E/A/delta/r + N1^2/2/E/A + xi * T1^2/2/G/A + N1*Mf1/ E/A/r ) *r, phi, 0, %pi/2 ); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] energia interna al tratto 2 [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ U2 : integrate( Mf2^2/2/E/J + N2^2/2/E/A + xi* T2^2/2/G/A , x, 0, l ); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] energia potenziale del quarto di maglia [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ U: U1 + U2$ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] scrivo equazione di compatibilità col vincolo bipendolo al punto di applicazione di C uso Castigliano [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ diff ( U , C ) = 0; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ linsolve ( % , C ) , globalsolve=true; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] controllando, trovo che i contributi a U di P introdotti mediante la coppia iperstatica C non vengono esplicitati forzo la sostituzione di C con la sua espressione in funzione di P per avere espliciti tutti i contributi energetici di P [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ U : ev(U); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] calcolo con castigliano lo spostamento del punto di applicazione della forza P riconducendomi a meta' struttura superiore (due quarti di maglia) per avere un'applicazione di Castigliano formalmente "standard" [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Delta : diff ( 2 * U , P ); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] definisco allungamento di riferimento sulla base di un tondino pari materiale e pari sezione [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Delta_rif : P*(l + r - d/2)/E/A; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ ratio : Delta / Delta_rif; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] suppongo materiale omogeneo e isotropo, nu = 3/10 [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ ip_materiale : [ G = E / 2 / (1+3/10) ]; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] suppongo sezione circolare [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ ip_sezione : [ A = %pi*d^2/4 , J = %pi * d^4/64, rn = (re - ri)^2/8/((ri+re)/2-sqrt(ri*re)), delta = r - rn, xi = 111/100 ]; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ ip_sezione[3] : ev(ip_sezione[3], ri = r-d/2, re=r+d/2); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ ip_sezione[4] : ev( ip_sezione[4] , ip_sezione[3] ); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] forma finale delle ipotesi sezione [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ ip_sezione ; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ ev ( ratio , ip_materiale,ip_sezione); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] ipotizzo specifico dimensionamento [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ DIN764_w6 : [ r=7, d=6, l=6.5 ,E=210000]; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ ev ( ratio , ip_materiale, ip_sezione, DIN764_w6, infeval, numer ); /* [wxMaxima: input end ] */ /* Maxima can't load/batch files which end with a comment! */ "Created with wxMaxima"$