Mut FEMBase.jl:ään ja siitä sitten.
@jvaara:
struct Mamo
params
stress
other
end
function Mamo()
params = zeros(12)
stress = reshape(view(params, 1:9), 3, 3)
other = view(params, 10:12)
return Mamo(params, stress, other)
endTämä muuten itse asiassa toimii. Eli muistia käytetään vain tuo 12 alkiota, mutta sulla on stress 3x3 matriisi joka ottaa parametrivektorista ensimmäiset 9 lukua ja other ottaa kolme viimeistä vektoriin. Kaikista voi päivittää ja kaikkia voi lukea. Tuo params voi sitten mennä epälineaariseen iteraatioon parametrivektorina. Tässä näin esimerkkinä mitenkä saa koodia sievemmäksi.
tai toisinpäin
Todennäköisesti voi rakentaa jopa sellaisen viewin 6-pituinen vektori näkyy 3x3 symmetrisenä tensorina.
ajattelin refaktoroida sen
chaboche.jl:n kunhan tajuan että kuinka mun pitäs toimia ton gitin kans :)
vektorinotaatioon laskenta jne turhat tensorit pois kun on harjoteltu ja on testi etc
ihan mielenkiintosta nähdä kun benchmarkkaa että minkälaista nopeutusta tulee esim. siihen simulaattori testiin
niin sitten kun haluaisin alkaa täyttämään sinne tavaraa niin kun laitan
simulator.stresses[1][:] .= copy(material.stress) niin tulee ERROR: UndefRefError: access to undefined refence
set materiaalimalli pitänee määritellä problemissa
nyt on solverissa mutta sehän on täysin hassua.
vai pitäiskö määritellä malli elementissä. element.material_model = :IdealPlastic jne.
Tai pienin yksikkö missä sitä myös kutsutaan
@ahojukka5
@ovainola voit sitten käyttää sitä testeissä kunhan saan sen valmiiksi
JuliaFEM/Materials.jl#29 tuolla se on jo tulossa. On vähän kesken vielä. Jostakin syystä ne tallennetut jutut on dtime:n verran väärässä ajassa.
Muutoin tuo on jo valmis.
Saa korjata.
on aina pitäny opetella tollasia inkrementaalisia plastisuusmalleja mutta ei oo oikein löytyny hyvää lähdeteosta tai tarpeeksi aikaa
toisaalta kiinnostaisi myös ne variaatioepäyhtälöt plastisuuteen liittyen
@ovainola onko sulla tehtynä joku koodi jolla voi generoida sitä materiaalin kiderakennetta.
^ disclaimer: tuo oli status meidän projektista. Yliopiston projektiin voi olla jotain 2D koodia.
olen tehnyt siten, että on valmiiksi määritelty tietorakenne, jossa on "supermesh", joka on jonkinlainen yhdistetty verkko molempien rajapintojen pintaverkoista, ja jokaiselle supermeshin elementille kaks indeksiä, jotka kertoo sitä vastaavat pintaverkon elementit eri rajapinnoilla. tuollaisen tietorakenteen avulla assembly on tosi nopea ja kätevä tehdä, mutta tällaisen tietorakenteen tekeminen tyhjästä yleiselle geometrialle on varmaankin haastavaa enkä ole sitä miettinyt. myös suurten siirtymien tapauksessa tietorakennetta pitää vähän rukata jokaisella askeleella, joten voi olla että se ei ole fiksu tapa
prof. hannukainen on tehnyt matlab-koodin, joka generoi tällaisen supermesh-rakenteen kahdelle pintaverkolle ja olen käyttänyt sitä omissa testilaskuissa. jos muistan oikein niin siinäkin piti antaa jonkinlainen bijektiivinen kuvaus pinnalta toiselle arvauksena, joten sekään ei ole riittävän yleinen mielivaltaisiin suurten siirtymien laskuihin
https://github.com/JuliaInterop/CxxWrap.jl tässä pitää osata tehdä jaettu kirjasto oikein, hello world esimerkki tässä: https://github.com/JuliaInterop/CxxWrap.jl/blob/master/test/hello.jl
Ensin kaivetaan kirjastojen sijainti:
julia> using CxxWrap
julia> dirname(dirname(CxxWrap.libcxxwrap_julia))
"/home/jukka/.julia/packages/CxxWrap/sarOk/deps/usr"Sitten compilessa -fPIC päälle ja include-polkuun CxxWrap ja Julia:
jukka@jukka-XPS-13-9380:/tmp/cxxwrap$ g++ -c -I/home/jukka/.julia/packages/CxxWrap/sarOk/deps/usr/include -I/opt/julia/include/julia -fPIC hello.cppLinkkauksessa jotain tämän tapaista:
jukka@jukka-XPS-13-9380:/tmp/cxxwrap$ g++ hello.o -shared -o libhello.so -L/opt/julia/lib -Wl,-rpath,/home/jukka/.julia/packages/CxxWrap/sarOk/deps/usr/lib: /home/jukka/.julia/packages/CxxWrap/sarOk/deps/usr/lib/libcxxwrap_julia.so.0 -ljuliajukka@jukka-XPS-13-9380:/tmp/cxxwrap$ julia sayhello.jl
CppHello.greet() = "hello, world"