۴-۷ پیشنهادات
در ادامه پژوهش حاضر موضوعات زیر پیشنهاد می شود.
بررسی تأثیر تعداد المانهای مرز خارجی برروی عملکرد الگوریتم یکی از مسائلی میباشد که می تواند از جمله عواملی باشد که برروی عملکرد الگوریتم تأثیر گذار باشد. از آنجا که اطلاعات روی مرز خارجی جزء اطلاعات معلوم مسئله به شمار میروند، ممکن است با افزایش تعداد آنها عملکرد الگوریتم بهتر شود.
بررسی کلی روشهای بهینه سازی و استفاده از آنها و مقایسه نتایج به دست آمده از آنها با نتایج حاصل از این پژوهش به منظور پیدا کردن بهترین روش بهینه سازی برای حل مسائل معکوس.
( اینجا فقط تکه ای از متن فایل پایان نامه درج شده است. برای خرید متن کامل پایان نامه با فرمت ورد می توانید به سایت feko.ir مراجعه نمایید و کلمه کلیدی مورد نظرتان را جستجو نمایید. )
مطالعه برای پیدا کردن جنبه های کاربردی الگوریتم پیشنهاد شده تا حالت عملیاتیتری پیدا کند.
از آنجا که حل معکوس مسئله دو حفرهای در طی تحقیقات به عمل آمده بر روی کارهای گذشته به ندرت انجام گرفته است، تحقیق بر روی شرایط مرزی دیگر و قیود اعمالی در روی مرزهای جسم و همچنین بسط آنها به مسائل بیشتر از دو حفره می تواند مفید باشد.
بررسی عملکرد الگوریتم با اضافه کردن تعداد حفرهها و همچنین به کار بردن المانهای مرتبه بالاتر مثل المانهای مرتبه سه، چرا که این المانها میتوانند مرز را با تعداد کمتری گره شبیه سازی کنند.
پیوست
clear
clc
%% COORDINATES OF BOUNDARIES
xb1= %coordinate of nodes of cavity number one in x direction
yb1=%coordinate of nodes of cavity number one in y direction
x=coordinates of cavities
xb12=[x(1)-x(3) x(1)-x(4)*cos(pi/6) x(1)-x(5)*cos(pi/3) x(1) x(1)+x(7)*cos(pi/3) x(1)+x(8)*cos(pi/6) x(1)+x(9) x(1)+x(10)*cos(pi/6) x(1)+x(11)*cos(pi/3) x(1) x(1)-x(13)*cos(pi/3) x(1)-x(14)*cos(pi/6) x(1)-x(3)];
yb12=[x(2) x(2)+x(4)*sin(pi/6) x(2)+x(5)*sin(pi/3) x(2)+x(6) x(2)+x(7)*sin(pi/3) x(2)+x(8)*sin(pi/6) x(2) x(2)-x(10)*sin(pi/6) x(2)-x(11)*sin(pi/3) x(2)-x(12) x(2)-x(13)*sin(pi/3) x(2)-x(14)*sin(pi/6) x(2)];
xb13=[x(15)-x(17) x(15)-x(18)*cos(pi/6) x(15)-x(19)*cos(pi/3) x(15) x(15)+x(21)*cos(pi/3) x(15)+x(22)*cos(pi/6) x(15)+x(23) x(15)+x(24)*cos(pi/6) x(15)+x(25)*cos(pi/3) x(15) x(15)-x(27)*cos(pi/3) x(15)-x(28)*cos(pi/6) x(15)-x(17)];
yb13=[x(16) x(16)+x(18)*sin(pi/6) x(16)+x(19)*sin(pi/3) x(16)+x(20) x(16)+x(21)*sin(pi/3) x(16)+x(22)*sin(pi/6) x(16) x(16)-x(24)*sin(pi/6) x(16)-x(25)*sin(pi/3) x(16)-x(26) x(16)-x(27)*sin(pi/3) x(16)-x(28)*sin(pi/6) x(16)];
hold on
plot(xb1,yb1);
plot(xb12,yb12,’k’);
plot(xb13,yb13);
scatter(xb1,yb1,’DisplayName’,’yb1 vs xb1′,’XDataSource’,’xb1′,’YDataSource’,’yb1′);figure(gcf)
scatter(xb12,yb12,’DisplayName’,’yb12 vs xb12′,’XDataSource’,’xb12′,’YDataSource’,’yb12′);figure(gcf)
scatter(xb13,yb13,’DisplayName’,’yb13 vs xb13′,’XDataSource’,’xb13′,’YDataSource’,’yb13′);figure(gcf)
hold off
%% Coordinate of nodes
xc=[-2 -1 0 1 2 2 2 2 2 1 0 -1 -2 -2 -2 -2 xb12(1) xb12(2) xb12(3) xb12(4) xb12(5) xb12(6) xb12(7) xb12(8) xb12(9) xb12(10) xb12(11) xb12(12) xb13(1) xb13(2) xb13(3) xb13(4) xb13(5) xb13(6) xb13(7) xb13(8) xb13(9) xb13(10) xb13(11) xb13(12)];
yc=[-2 -2 -2 -2 -2 -1 0 1 2 2 2 2 2 1 0 -1 yb12(1) yb12(2) yb12(3) yb12(4) yb12(5) yb12(6) yb12(7) yb12(8) yb12(9) yb12(10) yb12(11) yb12(12) yb13(1) yb13(2) yb13(3) yb13(4) yb13(5) yb13(6) yb13(7) yb13(8) yb13(9) yb13(10) yb13(11) yb13(12)];
%% Mechanical properties on damaines
Gs1=79300;
nu1=0.292;
%% Number of nodes and elements
TNE1=8;
TNE2=6;
TNE3=6;
TN1=16;
TN2=12;
TN3=12;
TN=TN1+TN2+TN3;
TNE=TNE1+TNE2+TNE3;
%% Relation between nodes and Boundary number 1
H1=zeros(TN*2,4*TNE1);
G1=zeros(TN*2,6*TNE1);
w=[0.066671344308688 0.149451349150581 0.219086362515982 0.269266719309996 0.295524224714753 0.295524224714753 0.269266719309996 0.219086362515982 0.149451349150581 0.066671344308688];
xi=[-0.973906528517172 -0.865063366688985 -0.679409568299024 -0.433395394129247 -0.148874338981631 0.148874338981631 0.433395394129247 0.679409568299024 0.865063366688985 0.973906528517172];
for i = 1 : TN
for j = 1 :TNE1
A=xb1(2*j+1)-2*xb1(2*j)+xb1(2*j-1);
B=0.5*(xb1(2*j+1)-xb1(2*j-1));
C=yb1(2*j+1)-2*yb1(2*j)+yb1(2*j-1);
D=0.5*(yb1(2*j+1)-yb1(2*j-1));
DE=4*pi*(1-nu1);
for l = 1 : 10
F1=0.5*xi(l)*(xi(l)-1);
F2=1-xi(l)^2;
