/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/ /* [ Created with wxMaxima version 14.12.1 ] */ /* [wxMaxima: input start ] */ kill(all); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ FA:0; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Mf1 : C - P/2*r*sin(theta)+ FA*r*(1-cos(theta)); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ Mf2 : C -P/2*r+ FA*(r+x); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] ricaviamo energia potenziale elastica della struttura ricavo energia potenziale elastica sul tratto 1, arco di circonferenza [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ U1 : integrate( Mf1^2/2/E/J *r , theta , 0 , %pi/2 ) ; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] ricavo energia potenziale elastica nel tratto rettilineo [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ U2 : integrate( Mf2^2/2/E/J , x , 0 , l) ; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] ricavo energia potenziale elastica totale [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ U : U1+U2; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] utilizzo il teorema di Castigliano e ricavo la rotazione al punto A [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ rotA : diff ( U , C , 1 ) ; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] derivata parziale d( sin(x*y)*exp(y) ) -------------------- dx dy^2 diff( sin(x*y)*exp(y) , x , 1 , y , 2 ); [wxMaxima: comment end ] */ /* [wxMaxima: comment start ] definisco equazione di compatibilità con il vincolo doppiopendolo originario in A [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ eqcomp : rotA = 0; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] risolvo l'equazione lineare nell'incognita iperstatica C [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ linsolve ( eqcomp , C ) , globalsolve=true; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ C; /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] voglio trovare l'estensione dell'intera maglia sotto carico P ricavo l'energia potenziale elastica dell'intera maglia come 4 volte quella del quarto [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Utot : ev( 4*U , infeval ); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] applico Castigliano per ricavare l'allungamento della maglia di catena sotto carico P [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ delta : diff (Utot, P); /* [wxMaxima: input end ] */ /* Maxima can't load/batch files which end with a comment! */ "Created with wxMaxima"$