GS_program_FDM_(v1.1)_prazni_predlozak

2074 days ago by pgidak@grad.hr


Metoda gustoća sila

(Force Density Method)


[Ćelije u kojima je prvi redak %auto ne trebate izvoditi ( evaluate ); one se izvode „same” pri otvaranju radnoga lista.  Međutim, ako radni list ponovno pokrenete—ne izlazeći iz njega—s  Restart worksheet  (u izborniku  Action... ), i njih ćete morati izvesti.]


       

Proračun metodom gustoća sila:

ncoords = FDM (nodes, elems, supports, qs)

parametri:

nodes — lista koordinata čvorova; zadane moraju biti koordinate ležajnih čvorova

elems — lista štapova određenih indeksima početnoga i krajnjeg čvora

supports — lista indeksa ležajnih čvorova

qs — lista gustoća sila u štapovima

rezultat:

ncoords — matrica tipa $n_{\check{c}}\!\times 3$ koja sadrži koordinate čvorova upisane po redcima

Čvorovi:

       
       

Elementi:

       
       

Ležajevi:

       
       


 

Gustoće sila:

       
       
       
       
       
       
       
       
       
       
       
       

 

Zadavanje vrijednostî sila u kabelima:

       
       

 

Iteracijska primjena metode gustoća sila odabrani broj puta:

ncoords, forces, f_dens = multistepFDM (nodes, elems, supports, qs, fcs = [], lcs = [], l0cs = [], steps = 250)

parametri:

nodes — lista koordinata čvorova; zadane moraju biti koordinate ležajnih čvorova

elems — lista štapova određenih indeksima početnoga i krajnjeg čvora

supports — lista indeksa ležajnih čvorova

qs — lista gustoća sila u štapovima za prvi korak

fcs — lista zadanih vrijednosti sila u štapovima

komponente liste: parovi (idx, f); idx – indeks štapa, f – vrijednost sile

lcs — lista zadanih duljina štapova

komponente liste: parovi (idx, l); idx – indeks štapa, l – duljina

l0cs — lista zadanih nerastegnutih duljina štapova

komponente liste: parovi (idx, (l0, AE))idx – indeks štapa, l0 – duljina, AE – uzdužna krutost 

steps — broj koraka

rezultati:

ncoords — matrica tipa $n_{\check{c}}\!\times 3$ koja sadrži po redcima svrstane koordinate čvorova

forces — lista konačnih vrijednosti sila u štapovima

f_dens — lista gustoća sila u štapovima izračunanih u posljednjem koraku

       
       

 

• iteracijska primjena metode gustoća sila do postizanja zadane točnosti:

ncoords, forces = multistepFDM_wtol (nodes, elems, supports, qs, fcs = [], lcs = [], l0cs = [], tol_f = 1e-3, tol_l = 1e-3, abs_error = True, steps = 250)

parametri:

nodes — lista koordinata čvorova; zadane moraju biti koordinate ležajnih čvorova

elems — lista štapova određenih indeksima početnoga i krajnjeg čvora

supports — lista indeksa ležajnih čvorova

qs — lista gustoća sila u štapovima za prvi korak

fcs — lista zadanih vrijednosti sila u štapovima

komponente liste: parovi (idx, f); idx – indeks štapa, f – vrijednost sile

lcs — lista zadanih duljina štapova

komponente liste: parovi (idx, l); idx – indeks štapa, l – duljina

l0cs — lista zadanih nerastegnutih duljina štapova

komponente liste: parovi (idx, (l0, AE))idx – indeks štapa, l0 – duljina, AE – uzdužna krutost 

tol_f — točnost za vrijednosti sila

tol_l — točnost za duljine

abs_error — ako je True, izračunava se apsolutna, a ako je False, relativna pogreška

primjerice, za vrijednosti sila:  $e_{\mathrm{abs}} = \Big\| S_{\{i,j\}}^{(k)} - \bar{S}_{\{i,j\}} \Big\|_{\mathrm{max}}$  i  $e_{\mathrm{rel}} = \dfrac{\Big\|S_{\{i,j\}}^{(k)} - \bar{S}_{\{i,j\}} \Big\|_{\mathrm{max}}}{\bar{S}_{\{i,j\}} }$

steps — broj koraka

rezultati:

ncoords — matrica tipa $n_{\check{c}}\!\times 3$ koja sadrži po redcima svrstane koordinate čvorova

forces — lista konačnih vrijednosti sila u štapovima

       
       


Povećanje vrijednosti sile u jednom kabelu: