programme pour calculer cette matrice

[waloudu93]

voila j'ai trouvé sa sur un forum c'est un programme en pascal pour calculer la matrice de rigidité ( peut on faire la même chose sur ti89)


uses crt,graph;
var i2,j2,npro,nfn,nelt,nnt,nnel,ndim,ndl,nnb,nnapp,i,j,pilote,mode,x1,y1,x2,y2,x1f,y1f,x2f,y2f:integer;
pp,y11f,x11f,i3,k,nt,n,p,q,n1,n2,nfen:integer;
tga,qt,ql,a,b,c,pivot,cosa,sina,lx,ly,l,x21,x11,y11,x22,y22,y21,xmax,ymax:real;
alp,ro,x1c,x2c,y1c,y2c,x11c,y11c:real;
rox,roy:integer;
case1:boolean;
f1,f2,f3,f4:text;
fd,fr,fdon,fint,nature:string;
type co=array[1..6]of integer;
di=array[1..3] of real;
rest=array[1..6]of integer;
rig=array[1..6]of real;
gr=array[1..50]of real;
var corg:array[1..50]of di;
qa,qqa:array[1..6] of real;
fen:array[1..50]of rig;
f:array[1..50]of real;
ff:array[1..50]of real;
fff:array[1..50]of real;
r:array[1..50]of real;
u:array[1..50] of real;
uu:array[1..50] of real;
connec:array[1..50] of co;
res:array[1..50]of rest;
fn:array[1..50]of di;
prop:array[1..50] of di;
tr:array[1..6] of rig;
rigid:array[1..6] of rig;
grigid:array[1..6] of rig;
w:array[1..6]of integer;
grig:array[1..50]of gr;
procedure don;
begin
{Lecture des donnes generales}
readln(f1,nnt,nelt,nnel,ndim,nnapp,nfn,nfen,npro);
{Lecture des coordonnees des noeuds}
for i:=1 to nnt do
begin
for j:=1 to ndim do
begin
read(f1,corg[i,j]);
end;
end;
{Lecture de la connectivit‚}
for i:=1 to nelt do
begin
for j:=1 to nnel +1 do
begin
read(f1,connec[i,j]);
end;
end;
{Lecture des conditions d'appuis}
ndl:=3;
for i:=1 to nnapp do
begin
read(f1,nnb);
for j:=1 to ndl do
begin
read(f1,res[nnb,j]);
end;
for q:=ndl downto 1 do
begin
res[nnb,q+1]:=res[nnb,q];
end;
res[nnb,1]:=nnb;
end;
for i:=1 to nnt do
begin
res[i,1]:=i;
end;

{Lecture des charges nodales}
for i:=1 to nfn do
begin
readln(f1,fn[i,1],fn[i,2],fn[i,3]);
n1:=trunc(fn[i,1]);
n2:=trunc(fn[i,2]);
f[3*(n1-1)+n2]:= fn[i,3];
end;
{Lecture des charges entre noeuds}
for i:=1 to nfen do
begin
for j:= 1 to 5 do
begin
read(f1,fen[i,j] );
end;
end;
{Lecture des proprietes}
for i:=1 to npro do
begin
readln(f1,prop[i,1],prop[i,2],prop[i,3]);
prop[i,1]:=prop[i,1]*prop[i,2];
prop[i,2]:=prop[i,1]*sqr(prop[i,2])/12;
end;
end;
procedure impdon;
begin
writeln(f3,'Nombre de noeuds = ',nnt:2);
writeln(f3,'Nombre d''elements = ',nelt:2);
writeln(f3,'Nombre de noeuds par element = ',nnel:2);
writeln(f3,'Dimension du probleme = ',ndim:2);
Writeln(f3);
writeln(f3);
writeln(f3,'Coordonnees des noeuds');
writeln(f3,'**********************');
case ndim of
3:writeln(f3,'Noeud x y z');
2:writeln(f3,'Noeud x y' );
end;
for i:=1 to nnt do
begin
write(f3,' ',i:2);
for j:=1 to ndim do
begin
write(f3,' ',corg[i,j]:12:3);
end;
writeln(f3);
end;
writeln(f3);
writeln(f3);
writeln(f3);
writeln(f3,'Connectivite des elements');
writeln(f3,'*************************');
writeln(f3,'Element noeud 1 noeud 2 Prop');
for i:=1 to nelt do
begin
write(f3,' ',i:2);
for j:=1 to nnel+1 do
begin
write(f3,' ',connec[i,j]:10);
end;
writeln(f3) ;
end;
writeln(f3);
writeln(f3);
writeln(f3,'degres de libertes bloquees');
writeln(f3,'***************************');
writeln(f3,'noeud ux uy rz');
for i:= 1 to nnt do
begin
for j:=1 to ndl+1 do
begin
write ( f3,' ',res[i,j]:5);
end;
writeln(f3);
end;
writeln(f3);
writeln(f3);
writeln(f3,'charges nodales');
writeln(f3,'***************');
writeln(f3,'noeud dir valeur');
for i:= 1 to nfn do
begin
writeln ( f3,trunc(fn[i,1]):2,' ',fn[i,2]:0:0,' ',fn[i,3]:9:3);
end;
writeln(f3,'charges entre noeuds');
writeln(f3,'********************');
writeln(f3,'noeud valeur nature');
for i:= 1 to nfen do
begin
if fen[i,4]=2 then
nature:='CONCENTREE'
else
nature:='REPARTIE';
writeln ( f3,trunc(fen[i,1]):2,' ',fen[i,3]:9:3,' ', nature);
end;
writeln(f3);
writeln(f3);
writeln(f3,'proprietes des elements');
writeln(f3,'***********************');
writeln(f3,'elem Sec iner E');
for i:= 1 to nelt do
begin
write( f3,i:2,' ',prop[connec[i,nnel+1],1]:5:4,' ',prop[connec[i,nnel+1],2]:6:4,' ');
writeln(f3,prop[connec[i,nnel+1],3]:9:4);
end;
end;
begin
clrscr;
writeln;
writeln;
writeln;
writeln('=====================================================================');
writeln('= *************************************************************** =');
writeln('= =');
writeln('= CALCUL DES STRUCTURES PAR LA METHODE DES ELEMENTS FINIS =');
writeln('= =');
writeln('= *************************************************************** =');
writeln('= =');
writeln('= =');
writeln('= Programme de Mr: =');
writeln('= ********************************* =');
writeln('= * SIMO2309 * =');
writeln('= * ETUDIENT * =');
writeln('= * DEPARTEMENT DE GENIE CIVIL * =');
writeln('= * UNIVERSITE MOHAMED BOUDIAF * =');
writeln('= ********************************* =');
writeln('= =');
writeln('=====================================================================');
writeln;
writeln;
writeln;
writeln(' APPUIYEZ SUR POUR CONTINUER');
readln;
clrscr;
WRITELN;
writeln;
writeln('DONNEZ LE NOM DU FICHIER DE DONNEES EN INDIQUANT SON CHEMIN D''ACCES');
readln(fd);
fdon:=fd+'.don';
assign(f1,fd);
assign(f3,fdon);
reset(f1);
don;
close(f1);
nt:=3*nnt;
rewrite(f3);
impdon;
close(f3);
writeln('Consultez les fichiers:');
writeln(' ',fdon,' pour les donn‚es');
writeln(' ');
writeln(' ' );
writeln(' Appyuez pour continuer');
readln;
end.

[Loulou 54]

Poua ya du boulot ! Je ne connais pas le pascal, perso... Je vais regarder "l'exemple très détaillé". Il est en train de charger depuis tt à l'heure... [Oui je n'ai pas l'adsl ]

= Bon, je ne pense pas avoir ni le tps ni le courage de m'y coller , je risquerai de me tromper, ne sachant pas quels doivent être les résultats et tout...
Demande ptèt aux gens de ta classe, (il y en a bien un qui se débrouille en TI-basic) comme ça il y'aura un intérêt commun à la réalisation de ce programme.
Parce-que faire un prog sans savoir exactement à quoi il sert et à quels paramêtres faire attention... c'est pas évident.

(Si tu parviens à en avoir un, mets-le sur TI-Bank ! Je suis en S-SI, ça pourrait m'être utile plus tard... Et à d'autres ! )

[waloudu93]

mercii

[Loulou 54]

Ah tu as déjà rep !^^ De rien ben regarde mon message au dessus que j'ai édité...

[Bisam]

En fait, la fonction "read" en Pascal sert à récupérer des données dans un fichier (ou éventuellement par le clavier ou n'importe quel autre entrée).
Ici, le programme récupère des montagnes de données dans un fichier f1 puis les redispose dans un fichier f3 sous une forme plus claire (probablement une matrice mais aussi plein d'autres choses.

Ce n'est pas du tout ce que tu veux pour 2 raisons :
- Tu n'as pas envie d'écrire tout ce qu'il faut mettre dans la matrice, tu voudrais que la calculatrice le fasse toute seule.
- Il y a des tas de choses qui ne correspondent pas à ce que tu nous as demandé.

Pour conclure, en tant que prof d'algorithmie (à très faible niveau, je l'admets), je me permets de te signaler une remarque que je fais souvent à mes élèves :

"Il est inutile de chercher pendant des heures, voire des jours, à faire faire à un ordinateur (ou une calculatrice) un calcul simple qui vous prend moins de 5 minutes... à moins que vous ne le fassiez TRES souvent".

[waloudu93]

En fait, la fonction "read" en Pascal sert à récupérer des données dans un fichier (ou éventuellement par le clavier ou n'importe quel autre entrée).
Ici, le programme récupère des montagnes de données dans un fichier f1 puis les redispose dans un fichier f3 sous une forme plus claire (probablement une matrice mais aussi plein d'autres choses.

Ce n'est pas du tout ce que tu veux pour 2 raisons :
- Tu n'as pas envie d'écrire tout ce qu'il faut mettre dans la matrice, tu voudrais que la calculatrice le fasse toute seule.
- Il y a des tas de choses qui ne correspondent pas à ce que tu nous as demandé.

Pour conclure, en tant que prof d'algorithmie (à très faible niveau, je l'admets), je me permets de te signaler une remarque que je fais souvent à mes élèves :

"Il est inutile de chercher pendant des heures, voire des jours, à faire faire à un ordinateur (ou une calculatrice) un calcul simple qui vous prend moins de 5 minutes... à moins que vous ne le fassiez TRES souvent".

Mais le problème c'est que le calcul demande ou moins 10 minutes et en plus la durée de nos examen sont de 1h15 (et il y a 4 exo avec plein de truck intégral triple, équation diff, calcul matricielle ect.) (je suis en 4eme année d'école d'ingénieurs à L'ENPC)
moi je maitrise bien la matrice mais je chercher à être plus rapide dans les calcul j'ai juste à comprendre le portique et mettre les bonnes valeurs dans la matrice.

On n'oublie alors le programme que j'ai trouver et on part à zéro, qu'est ce que vous en pensez??

[Bisam]

A vrai dire, je veux bien t'aider (je le fais même) mais il va falloir que tu fasses un cahier des charges plus précis.

1) Quelles sont les données ?
2) Que veux-tu rentrer dans la calculatrice ?
3) Que doit renvoyer la calculatrice ?
4) Sous quelle forme ?

J'ai l'impression que pour remplir ta matrice il y a tellement de contraintes que donner ces contraintes à la calculatrice sera au moins aussi long que d'écrire toi-même la matrice.
Je ne comprends pas (encore) pourquoi il faut 10 minutes pour trouver cette matrice à partir du moment où toutes les données (contraintes) sont connues.

[waloudu93]

bonjour

et bien je pense que je vais laisser tomber ce programme car j'ai oublier de précisé qu'il me le faut pour le lundi qui vient
donc je pense que je vais me concentrer de le refaire plusieurs fois afin que je devienne rapide

merci pour ton aide bissam

[Bisam]

je pense que je vais me concentrer de le refaire plusieurs fois afin que je devienne rapide

Voilà une décision qui me semble très sage !
Tu iras loin.

merci pour ton aide bissam

Mais, je t'en prie, le forum est là pour ça.