?? oc.m
字號:
% optimality criteria update
function [xnew]=oc(nelx,nely,x,volfrac,dc)
l1=0;l2=100000;move=0.2;
while (l2-l1>1e-4)
lmid=0.5*(l2+l1);
xnew =max(0.001,max(x-move,min(1.,min(x+move,x.*sqrt(-dc./lmid)))));
if sum(sum(xnew))-volfrac*nelx*nely>0;
l1=lmid;
else
l2=lmid;
end
end
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -