【正文】 參考文獻(xiàn) 1 唐文艷 .顧元憲 .趙國忠復(fù) 合材料層合板鋪層順序優(yōu)化遺傳算法 [期刊論文 ]大連理工大 學(xué)學(xué)報(bào) 20xx,44(2) 2 周兵 .楊靜 .李偉平主副簧鋼板彈簧離散梁動(dòng)力學(xué)參數(shù)反求 [期刊論文 ]計(jì)算機(jī)仿 20xx(1) 3 高占遠(yuǎn) .鹿曉陽 .陳世英 .賈虎基于單純形模擬退火算法的網(wǎng)架結(jié)構(gòu)優(yōu)化設(shè)計(jì) [期刊論文 ] 山東建筑大學(xué)學(xué)報(bào) 20xx(1) 4 。 SA算法進(jìn)行結(jié)構(gòu)優(yōu)化，結(jié)構(gòu)重分析的次數(shù)太多，這是阻礙其在大型結(jié)構(gòu)優(yōu)化中應(yīng)用的主要問題。 3，本文主要是 利用屈曲、震動(dòng) 間題來說明 SA算法的特性。 2， SA算法的應(yīng)用非常廣泛 .但是它的理論還不夠完善，這也阻礙了 它的發(fā)展。對(duì)于這些參數(shù)的選擇還需要進(jìn)一步的研究。 % Equality constraints geq=[]。\n%39。E:\MAT_ANSYS\39。s function % z=(1u(1))^2+(u(2)u(1)^2)^2+(1u(3))^2+(1u(4))^2 % All constraints function [g,geq]=constraints(x) % Inequality constraints g=[]。新設(shè)計(jì)變量 39。)。,39。 %z=fscanf(39。)。,39。 end % Objective functions function z=fobj(u) % Welded Beam Design Optimization % K. Ragsdell and D. Phillips, Optimal design of a class of welded % strucures using geometric programming, % J. Eng. Ind., 98 (3):10211025, (1976). % Best solution found in literature % [, , , ) with the objective % by Cagnina et al., Solving engineering optimization problems with the % simple constrained particle swarm optimizer, Informatica, 32 (20xx)319326 z=textread(39。 end % Test if equalities hold function H=geteqH(g) if g==0, H=0。 end % Test if inequalities hold function H=getH(g) if g=0, H=0。 % Z=sum(lameq*geq(k1).^2.*geteqH(geq(k1)))。 % Inequality constraints for k=1:length(g), Z=Z+ lam*g(k)^2*getH(g(k))。 lameq=10^15。 function Z=getnonlinear(u) Z=0。 %調(diào)用 ANSYS程序 !C:\Program Files\Ansys Inc\v100\ANSYS\bin\intel\ b p ane3fl i E:\MAT_ANSYS\ o E:\MAT_ANSYS\ % Objective z=fobj(u)。,u)。 fprintf(fid,39。w39。E:\MAT_ANSYS\39。 else ns=ns。 ns_tmp(J)=Ub(J)。 ns_tmp(I)=Lb(I)。 %% function ns=bounds(ns,Lb,Ub) if length(Lb)0, % Apply the lower bound ns_tmp=ns。 end s=bounds(s,Lb,Ub)。 init_flag==1, s=Lb+(UbLb).*rand(size(u0))。 end % Update the estimated optimal solution fmin=E_old。 accept=accept+1。 best = ns。amp。 j = 0。 E_old = E_new。 % Decide to accept the new solution DeltaE=E_newE_old。 end totaleval=totaleval+1。 else init_flag=0。 hold on end % Function evaluations at new locations if totalevalinitial_search, init_flag=1。kx39。)。 ylabel(39。迭代次數(shù) 39。 diedai=diedai+1。The best found so far =39。 % Cooling according to a cooling schedule T = cooling(alpha,T)。 % Check if max numbers of run/accept are met if (i = max_run) || (accept = max_accept) % reset the counters i = 1。 % initially guessed values diedai=1。 E_new=E_old。 E_init = Fun(u0)。 totaleval = 0。 j = 0。 % Boltzmann constant Enorm=1e5。 % Maximum number of accept initial_search=20xx。 % Maximum number of rejections max_run=500。 % Finial stopping temperature F_min = 1e+100。 % Dimension of the problem % Initializing parameters and settings T_init = 。)。 if length(Lb) ~=length(Ub), disp(39。 Ub=[90 90]。)。 end % Display usage disp(39。 利用 ANSYS中 shell99建立模型得到的 ANSYS分析文件如下： finish /clear *dim,aa,1,2 *vread,aa(1,1),e:\mat_ansys\sa_shuju,dat,JIK,2,1 () x1=aa(1,1) x2=aa(1,2) t0= !單層板厚度 a0=55e2 b0=a0 geshu=10 !長(zhǎng)度為 b的邊網(wǎng)格個(gè)數(shù) ! /PREP7 ET,1,SHELL99 !ET,1,SHELL91 ! KEYOPT,1,8,1 R,1,4,1 rmore rmore,1,x1,t0,1,x2,t0, ! mp,ex,1, ! mp,EY,1, MP,EZ,1, MP,PRXY,1, MP,PRYZ,1, !有待確認(rèn) MP,PRXZ,1, MP,GXY,1, MP,GYZ,1, !有待確認(rèn) MP,GXZ,1, mp,dens,1, !密度， 有待確認(rèn) ! RECTNG,a0/2,a0/2,b0/2,b0/2, TYPE,1 MAT,1 REAL,1 ESYS,0 !ESIZE,0,10 lsel,s,loc,x,a0/2 lesize,all,geshu !劃分 9等份 allsel lsel,s,loc,y,b0/2 lesize,all,a0/b0*geshu !劃分 12等份 allsel mshape,0,2D mshkey,1 ! amesh,all finish ! /SOL antype,modal !模態(tài)分析 modopt,lanb,1 mxpand,1 !施加約束 nsel,s,loc,x,a0/2 d,all,uz allsel nsel,s,loc,y,b0/2 d,all,uz allsel nsel,s,loc,x,a0/2 d,all,uz allsel nsel,s,loc,y,b0/2 d,all,uz allsel nsel,s,loc,x,0 d,all,ux allsel nsel,s,loc,y,0 d,all,uy allsel solve finish ! /post1 set,first *get,plv1,mode,1,freq *dim,jieguo,1,1 jieguo(1,1)=plv1 *mwrite,jieguo,e:\mat_ansys\ans_meici,dat,JIK,1,1 () MATLAB調(diào)用 ANSYS文件程序 程序如下： % % % Matlab Programs included the Appendix B in the book: % % XinShe Yang, Engineering Optimization: An Introduction % % with Metaheuristic Applications % % Published by John Wiley amp。對(duì)于以上模型，將每層又分成纖維鋪設(shè)角與 x軸成 θ 角與 θ 角的兩層進(jìn)行計(jì)算。材料常數(shù) : 彈性模 E1= ， E2=. 剪切模量G,2=, 泊松比 v12= 。 利用 ANSYS中 shell99建立模型得到的 ANSYS分析文件如下： finish /clear *dim,aa,1,2 *vread,aa(1,1),e:\mat_ansys\sa_shuju,dat,JIK,2,1 () x1=aa(1,1) x2=aa(1,2) t0= !單層板厚度 a0=55e2 b0=a0 geshu=10 !長(zhǎng)度為 b的邊網(wǎng)格個(gè)數(shù) ! /PREP7 ET,1,SHELL99 !ET,1,SHELL91 ! KEYOPT,1,8,1 R,1,4,1 rmore rmore,1,x1,t0,1,x2,t0, ! mp,ex,1, ! mp,EY,1, MP,EZ,1, MP,PRXY,1, MP,PRYZ,1, !有待確認(rèn) MP,PRXZ,1, MP,GXY,1, MP,GYZ,1, !有待確認(rèn) MP,GXZ,1, ! RECTNG,a0/2,a0/2,b0/2,b0/2, TYPE,1 MAT,1 REAL,1 ESYS,0 !ESIZE,0,10 lsel,s,loc,x,a0/2 lesize,a