Metoda Sił
w ujęciu macierzowym
Belka statycznie
niewyznaczalna
Sprawdzenie
Wykonanie w programie Matlab
Zapisanie skryptu
Proponuję, aby skrypt sprawdzający został zapisany do
osobnego pliku, którego nazwa powstanie poprzez dodanie
przyrostka  spr do dotychczas używanej.
Przykładowo:
Nazwa pliku ze skryptem rozwiązującym belkę 
agat.m
Nazwa pliku ze skryptem sprawdzającym rozwiązanie 
agat_spr.m
Na początku skrypt sprawdzający będzie zatem
identyczny ze skryptem rozwiązującym.
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Nowy układ podstawowy
W celu kontroli poprawności uzyskanego rozwiązania
obliczamy przemieszczenie odpowiadające zwolnionemu
więzowi w dodatkowym, nowym UP i sprawdzamy czy to
przemieszczenie jest równe ZERO.
Nowy UP
X3 X4
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Dodatkowy wykres jednostkowy
Sporządzamy nowy wykres momentów dla jednostkowego
obciążenia odpowiadającego niewiadomej X3.
2.0
1 4.0
7.0
M3
10.0
Rzędne tego wykresu wprowadzamy do skryptu:
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Wyniki
Wyniki, które widać obok zostały
uzyskane jako rozwiązanie dla
przykładowej belki.
Teraz należy wprowadzić te rzędne do
skryptu sprawdzającego.
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Wyniki - c.d.
Wyniki, które są wyświetlane jako efekt
działania skryptu rozwiązującego belkę
można szybko przenieść do skryptu
sprawdzającego.
Zaznaczamy w tym celu myszką wektor
wynikowy tak, jak to widać obok
i wywołujemy w dowolny sposób polecenie
Kopiuj.
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Wyniki - c.d.
W edytorze, w którym pracujemy nad
skryptem sprawdzającym wpisujemy
nazwę zmiennej, znak równości i nawias
otwierający a następnie wydajemy
polecenie Wklej.
Pozostaje dopisanie średników po każdej
z wklejonych liczb oraz dodanie nawiasu
zamykającego.
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Edycja skryptu
Zmiany, które należy również wprowadzić w skrypcie to:
Usunięcie niewykorzystywanych danych czyli macierzy
M1, M2 oraz Mp.
Zmiana definicji macierzy momentów jednostkowych
M  teraz posługujemy się tylko jednym wykresem
jednostkowym, innym niż dotychczas:
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Edycja skryptu  c.d.
Zamiast fragmentu skryptu, w którym obliczane były
macierz podatności ", wektor wyrazów wolnych "p,
rozwiązywane było równanie MS oraz obliczane były
rzędne wykresu momentów zginających wpisujemy co
następuje:
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Uruchomienie skryptu
Po wpisaniu nazwy pliku, do którego został zapisany
skrypt sprawdzający otrzymamy na ekranie wynik. Jest to
obliczone przemieszczenie odpowiadające zwolnionemu
więzowi w dodatkowym, nowym UP.
Przemieszczenie to powinno być równe ZERO, gdyż
punkt, w którym zostało ono obliczone jest miejscem
występowania w rozważanej belce więzów czyli podpory.
ZERO
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Sprawozdanie
Sprawozdanie z tej części projektu powinno zawierać:
1) schemat belki
2) rysunek nowego układu podstawowego
3) dodatkowego wykres jednostkowy
4) obliczone przemieszczenie
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003
Politechnika Rzeszowska, Katedra Mechaniki Konstrukcji, 10'2003