www.gusucode.com > 基于多边形网格的流动问题拓扑优化设计元代码 > MatIntFnc.m
%------------------------------ PolyTop ----------------------------------% % Ref: A Pereira, C Talischi, GH Paulino, IFM Menezes, MS Carvalho % % "Fluid Flow Topology Optimization in PolyTop: Stability and % % Computational Implementation", Struct Multidisc Optim, % % DOI 10.1007/s00158-014-1182-z % %-------------------------------------------------------------------------% function [E,dEdy,V,dVdy] = MatIntFnc(y,param) mu0 = param(1); q = param(2); epsilon = 4e-5; E = (mu0/epsilon)*q.*(1-y)./(y+q); dEdy = -(mu0/epsilon)*(1+q)*q./((y+q).^2); V = y; dVdy = ones(size(y,1),1); %-------------------------------------------------------------------------%