Bilag [1,7 MB] - Morten Christiansen

72 APPENDIKS G. PROGRAMLISTE
K_QuadIso_LinElas.m
Beregning af elementstivhedsmatrix for isoparametrisk 4-knuders element.
shape_3d_beam.m
Opstilling af formfunktioner for tredimensionelt bjælkeelement.
transformation.m
Indføring af matrix, der omregner et elements stivhedsmatrix i lokale koordinater
til globale. Filen anvendes ved tredimensionelle bjælker.
K_3d_beam.m
Beregning af elementstivhedsmatrix for tredimensionelt bjælkeelement.
Assemblering.m
Elementernes stivhedsmatricer samles til en global for det samlede legeme.
GaussRect.m
Bestemmelse af Gausspunkter til anvendelse ved numerisk integration ved
isoparametriske elementer.
solvequations.m
Løsning af ligningssystemet K a = f.
Ovenstående hovedprogram og filer er undervejs i projektet videreudviklet,
så udvalgte geometrier kan inddeles i et vilkårligt antal elementer. Udover
de nævnte filer er der lavet nye filer, hvor de filer, der bruges generelt, er
forklaret nedenfor, mens der for de øvrige henvises til næste afsnit.
K_QuadIso8_LinElas.m
Beregning af elementstivhedsmatrix for isoparametrisk 8-knuders element.
stressstrain.m
Bestemmelse af spændinger og tøjninger i elementerne.
G.2 Øvrige programmer
I forbindelse med projektet er der udover de førnævnte filer lavet en række
øvrige programmer og filer, der henvises til i rapporten.
1. FEM-bjælkemodeller
Mappe med FEM-filer anvendt til heksagonal og kvadratisk enhedscelle.
Udover de førnævnte filer til løsning af elementmetodeproblemer
indeholder denne mappe følgende: