多约束拓扑优化MATLAB,最大固有频率为目标，体积比为约束拓扑优化matlab程序

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

while (l2-l1 > 1e-4)

lmid = 0.5*(l2+l1);

x(:,:,m+1) = max(0.001,max(x(:,:,m)-move,min(1.,min(x(:,:,m)+move,x(:,:,m).*((-dc./lmid).^0.5)))));%%%%sqrt=0.5 is damping coefficient

if sum(sum(x(:,:,m+1))) - volfrac*nelx*nely > 0

l1 = lmid;

else

l2 = lmid;

end

end

%%%%%%%%%% MESH-INDEPENDENCY FILTER %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [dcn]=check(nelx,nely,rmin,x,dc,m)

dcn=zeros(nely,nelx);

for i = 1:nelx

for j = 1:nely

sum=0.0;

for k = max(i-floor(rmin),1):min(i+floor(rmin),nelx)

for l = max(j-floor(rmin),1):min(j+floor(rmin),nely)

fac = rmin-sqrt((i-k)^2+(j-l)^2);           %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%rimn is filter radius

sum = sum+max(0,fac);

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！