【布局优化】基于改进粒子群算法实现充电桩选址优化问题附matlab代码

1 简介

当前世界环境污染和能源危机问题凸显,电动汽车以零排放和低耗能的优势得到各国的大力关注和支持。以电动汽车为代表新能源汽车产业,成为国家七大战略性新兴产业之一。电动汽车具有良好的发展前景,市场规模也在逐步扩大。但是作为电动汽车的配套设施中的充电设施网络的不完善在很大程度上阻碍了电动汽车的使用和普及。充电设施对于电动汽车作用相当于加油站对于燃油车的作用,充电设施的选址和设施内充电桩的数量直接影响到了电动汽车的发展,也间接影响了充电设施的服务质量和运营效益。因此,本文在一定预算费用限制下,需要对电动汽车充电设施的选址位置进行确定的同时,还要确定每个选址点要建设多少个快速充电桩。

2 部分代码

clear; clc; close all;

%Input
%--------------------------------------------------------------------------
S=24; %Span
N=12; %Number of Bays (even)
H1=1; %Height at edge
H2=3.7; %Height at midspan
DL=-4000; %Dead Load on Top Chord
LL=-4000; %Live Load on Top Chord

%Nodes Coordinates and Degrees of Freedoms
%------------------------------------------
for i=1:1:N+1
Node{i}.x=S*(i-1)/N; Node{i}.y=0;
Node{i}.DOFx=2*i-1; Node{i}.DOFy=2*i;
end

for i=N+2:1:1.5*N+2
Node{i}.x=S*(i-N-2)/N; Node{i}.y=H1+(H2-H1)*(i-N-2)/(N/2);
Node{i}.DOFx=2*i-1; Node{i}.DOFy=2*i;
end

for i=1.5*N+3:1:2*N+2
Node{i}.x=S*(i-N-2)/N; Node{i}.y=H1+(H2-H1)*(2*N+2-i)/(N/2);
Node{i}.DOFx=2*i-1; Node{i}.DOFy=2*i;
end

%Elements Connectivity, Section, and Modulus of Elasticity, and Yield Strength
%--------------------------------------------------------------------------
for i=1:1:N %Bottom Chord
Element{i}.Nodes=[i i+1];
Element{i}.A=1472e-6;
Element{i}.rmin=0.0237;
Element{i}.E=200e9;
Element{i}.Fy=250e6;
end

for i=1:1:N %Top Chord
Element{N+i}.Nodes=[N+1+i N+2+i];
Element{N+i}.A=1472e-6;
Element{N+i}.rmin=0.0237;
Element{N+i}.E=200e9;
Element{N+i}.Fy=250e6;
end

for i=1:1:N/2 %Diagonal Members
Element{2*N+i}.Nodes=[i+1 N+1+i];
Element{2*N+i}.A=1472e-6;
Element{2*N+i}.rmin=0.0237;
Element{2*N+i}.E=200e9;
Element{2*N+i}.Fy=250e6;
end

for i=N/2+1:1:N %Diagonal Members
Element{2*N+i}.Nodes=[i N+2+i];
Element{2*N+i}.A=1472e-6;
Element{2*N+i}.rmin=0.0237;
Element{2*N+i}.E=200e9;
Element{2*N+i}.Fy=250e6;
end

for i=1:1:N+1 %Vertical Members
Element{3*N+i}.Nodes=[i N+1+i];
Element{3*N+i}.A=1472e-6;
Element{3*N+i}.rmin=0.0237;
Element{3*N+i}.E=200e9;
Element{3*N+i}.Fy=250e6;
end

%Support, Support Displacement, and Applied Force
%-------------------------------------------------
Support=[1 2 2*N+2];
U_s=[0 0 0]';
Force=zeros(4*N+4,1);
Force(2*N+4)=(DL+LL)/2;
Force(4*N+4)=(DL+LL)/2;
Force(2*N+6:2:4*N+2)=DL+LL;

%Processing
%--------------------------------------------------------------------------
[Node,Element,F_s]=Analysis(Node,Element,Support,U_s,Force);

%Check AISC code
%--------------------------------------------------------------------------
Element=Check_Code(Element);

%Output
%--------------------------------------------------------------------------
Mag_Factor=20;
Plot_Results(Node,Element,Mag_Factor,0);
Plot_Results(Node,Element,Mag_Factor,1); title('Strain');
Plot_Results(Node,Element,Mag_Factor,2); title('Stress');
Plot_Results(Node,Element,Mag_Factor,3); title('Axial Force');
Plot_Results(Node,Element,Mag_Factor,4); title('Design Ratio');

3 仿真结果

4 参考文献

[1]齐琳, 姚俭, 王心月. 基于改进粒子群算法的电动汽车充电站布局优化[J]. 公路交通科技, 2017, 34(6):8.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。