水平集拓扑优化经典程序.docVIP

function TOPLSM(DomainWidth, DomainHight, EleNumPerRow, EleNumPerCol, LV, LCur,FEAInterval, PlotInterval, TotalItNum) %=================================================================% % TOPLSM, a 199-line Matlab program, is developed and presented here for the % mean compliance optimization of structures in 2D, with the classical level set method. % % Developed by: Michael Yu WANG, Shikui CHEN and Qi XIA % First Version : July 10, 2004 % Second Version: September 27, 2004 % Last Modification:October 31, 2005, optimize the code. % % The code can be downloaded from the webpage: % .hk/~cmdl/download.htm % % Department of Automation and Computer-Aided Engineering, % The Chinese University of Hong Kong % Email: yuwang@.hk % %Main references: % (1.)M.Y. Wang, X. M. Wang, and D. M. Guo,A level set method for structural topology optimization, % Computer Methods in Applied Mechanics and Engineering, 192(1-2), 227-246, January 2003 % %(2.) M. Y. Wang and X. M. Wang, PDE-driven level sets, shape sensitivity, and curvature flow for structural topology optimization, % CMES: Computer Modeling in Engineering Sciences, 6(4), 373-395, October 2004. % %(3.) G. Allaire, F. Jouve, A.-M. Toader, Structural optimization using sensitivity analysis and a level-set method , % J. Comp. Phys. Vol 194/1, pp.363-393 ,2004. %Parameters: % DomainWidth : the width of the design domain; % DomainHight : the hight of the design domain; % EleNumPerRow : the number of finite elements in horizontal direction; % EleNumPerCol : the number of finite elements in vertical direction; % LV : Lagrange multiplier for volume constraint; % LCur : Lagrange multiplier for perimeter constraint whose shape sensitivity is curvature; % FEAInterval : parameters to specify the frequency of finite element % analysis; % PlotInterval : parameters to specify the frequency of plotting; % TotalItNum : total iteration number. %=================================================================% % Step 1: Data initializati