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);
%-------------------------------------------------------------------------%