% Multiphysics fem=multiphysics(fem); % Extend mesh fem.xmesh=meshextend(fem); % Solve problem fem.sol=femnlin(fem, ... 'solcomp',{'T'}, ... 'outcomp',{'T'}); % Save current fem structure for restart purposes fem0=fem; % Integrate I1=postint(fem,'T', ... 'dl',[1]); % Integrate I2=postint(fem,'1', ... 'dl',[1]); kN(i)=I1/I2; disp(I1/I2); %end comsol end; tN=kN; %end function b=checkCoords(x,r,bitsNV,varLim,varN,popN,bestCh) for i=1:popN uklad=1; for j=1:2:varN p1r=(x(i,j)/2^bitsNV(j))*(varLim(j,2)-varLim(j,1))+varLim(j,1); p1f=(x(i,j+1)/2^bitsNV(j+1))*(varLim(j+1,2)-varLim(j+1,1))+varLim(j+1,1); for k=(j+2):2:varN p2r=(x(i,k)/2^bitsNV(k))*(varLim(k,2)-varLim(k,1))+varLim(k,1); p2f=(x(i,k+1)/2^bitsNV(k+1))*(varLim(k+1,2)-varLim(k+1,1))+varLim(k+1,1); d=sqrt(p1r^2+p2r^2-2*p1r*p2r*cos(p2f-p1f)); %disp(d); if(d
Załącznik 2 – kod programu - rozwiązanie zagadnienia Motza metodą rozwiązań podstawowych Program napisano w języku Matlab/Comsol Script function v=motz2(pkt,NZ,saveF) %function v=motz2() S=0.2; M=12; %saveF=1; %NZ=12; %[XZ,YZ]=zrodla(NZ,S); XZ=pkt(1,:); YZ=pkt(2,:); N0=6*M-3-NZ; N0=3; % kollokacja z rozwiązaniami podstawowymi %NZ=24; %podzielne przez 6; [XAB,YAB,XCB,YCB,XDC,YDC,XED,YED,XEO,YEO,XOA,YOA,E,F,P,O,H,G,T,Z]=wkolMFS(M); [A,B]=muklMFS(M,N0,NZ,XZ,YZ,XAB,YAB,XCB,YCB,XDC,YDC,XED,YED,XEO,YEO,XOA,YOA,E,F,P,O,H,G, T,Z); X=linsolve(A,B'); bd1=0; bd2=0; bd3=0; bd4=0; bd5=0; bd6=0; xb1=1; xb2=0; xb3=1; xb4=0; yb1=0; yb2=0; num=201; for i=1:num wsp=(i-1)/num; [dfu,xdf,ydf]=rozprzy(1,wsp,XZ,YZ,X,NZ,N0); pd1=abs(dfu-1); if(bd1
- Page 1 and 2:
Politechnika Poznańska Wydział Bu
- Page 3 and 4:
5 Optymalizacja położenia źróde
- Page 5 and 6:
Streszczenie Praca poświęcona jes
- Page 7 and 8:
zegowych równań całkowych, co zm
- Page 9 and 10:
wyliczane są współczynniki wagow
- Page 11 and 12:
W pracy analizowane są zagadnienia
- Page 13 and 14:
włókien, oraz przypadek odwrotny.
- Page 15 and 16:
2 Optymalizacja 2.1 Wprowadzenie Dl
- Page 17 and 18:
W powyższym wzorze f oznacza maksy
- Page 19 and 20:
z punktów w przestrzeni poszukiwa
- Page 21 and 22:
osobniki co oznacza, że w wyniku r
- Page 23 and 24:
manipulowanie wartościami poszczeg
- Page 25 and 26:
twórcę programu znającego charak
- Page 27 and 28:
metody jest liczność populacji ty
- Page 29 and 30:
3 Metoda rozwiązań podstawowych 3
- Page 31 and 32:
a ij = Bφ i x j , j = 1, … , N,
- Page 33 and 34:
Tabela 3.1 Przykładowe zbiory funk
- Page 35 and 36:
najnowszych opublikowanych badań z
- Page 37 and 38:
do obliczeń współrzędnych punkt
- Page 39 and 40:
gdzie L i B oznaczają operator lin
- Page 41 and 42:
4.3 Rozwiązanie dwuwymiarowego ust
- Page 43 and 44:
Uwzględniając warunki brzegowe ot
- Page 45 and 46:
5 Optymalizacja położenia źróde
- Page 47 and 48:
Na brzegu Γ zadana jest stała tem
- Page 49 and 50:
Ponieważ istnieje tak duża liczba
- Page 51 and 52:
5.5 Wyniki numeryczne Do obliczeń
- Page 53 and 54:
Jak wspomniano wcześniej, w algory
- Page 55 and 56:
W drugim cyklu obliczeń populacja
- Page 57 and 58:
Rysunek 5.12 Konfiguracja 4 element
- Page 59 and 60:
Zestawienie uzyskanych wyników dla
- Page 61 and 62:
Rysunek 5.20 Rozwiązanie optymalne
- Page 63 and 64: W związku z tym, że położenie p
- Page 65 and 66: gdzie ∂ ∂n oznacza pochodną w
- Page 67 and 68: Tabela 6.1 Zestawienie wyników obl
- Page 69 and 70: Wykresy ilustrujące rozwiązanie i
- Page 71 and 72: Występująca na brzegu osobliwoś
- Page 73 and 74: odsunięta jest o 2 od brzegu obsza
- Page 75 and 76: 7 Wyznaczanie efektywnego współcz
- Page 77 and 78: 7.3 Równanie rządzące i warunki
- Page 79 and 80: T i = c ij ln r j 2 j dla i = M, F
- Page 81 and 82: 7.7 Wyniki numeryczne 7.7.1 Komórk
- Page 83 and 84: Rysunek 7.5 Układ włókien dla os
- Page 85 and 86: włókien. Parametry materiału prz
- Page 87 and 88: Rysunek 7.8 Układ włókien dla os
- Page 89 and 90: przedstawiono w tabeli 7.18. Wykorz
- Page 91 and 92: Rysunek 7.11 Układ włókien dla o
- Page 93 and 94: 7.7.5 Komórka zawierająca 3 włó
- Page 95 and 96: Rysunek 7.14 Układ włókien dla o
- Page 97 and 98: Tabela 7.32 Parametry metody rozwi
- Page 99 and 100: 8 Wnioski i podsumowanie Niniejsza
- Page 101 and 102: Literatura [Ake2000] Akella M. R.,
- Page 103 and 104: [Fai1998] Fairweather G., Karageorg
- Page 105 and 106: wykorzystaniem metody kollokacji br
- Page 107 and 108: [Pal2008] Paluch B., Grediac M., Fa
- Page 109 and 110: Załączniki W załącznikach zawie
- Page 111 and 112: wM(nw1,bw)=wM(nw1,bw)-2^(k); wM(nw2
- Page 113: end Temp=I1/I2; newline=NaN; plik='
- Page 117 and 118: dfu=dfu+X(i)*fln(xw,yw,XZ(j),YZ(j))
- Page 119 and 120: Załącznik 3 - kod programu - wyzn
- Page 121 and 122: AP( PK_z+i, 4*PZ_z+KW_cnt*PZ_ww+ j
- Page 123 and 124: colormap(jet(128)) % surf(iX,iY,iU)