Scribd
Upload
41 views

Volume Optimization Problem: 1) - Muhammad Zargham Ali 2) - Hamad Khan Afridi

The document describes a Matlab code that uses symbolic optimization to find the optimal volumetric dimensions of a power converter design problem. The code defines objective and constraint functions in terms of the design variables of width (wc), width (ww), height (hw), and distance between coils (dc). It then solves the system of equations from setting the derivatives of the Lagrangian equal to zero to find the optimal symbolic solution, which is then evaluated numerically to find the optimal real positive values of the design variables and minimum volume.

Uploaded by

UmairRashid
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
41 views

Volume Optimization Problem: 1) - Muhammad Zargham Ali 2) - Hamad Khan Afridi

The document describes a Matlab code that uses symbolic optimization to find the optimal volumetric dimensions of a power converter design problem. The code defines objective and constraint functions in terms of the design variables of width (wc), width (ww), height (hw), and distance between coils (dc). It then solves the system of equations from setting the derivatives of the Lagrangian equal to zero to find the optimal symbolic solution, which is then evaluated numerically to find the optimal real positive values of the design variables and minimum volume.

Uploaded by

UmairRashid
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

MS Systems Engineering 2017-19

Volume Optimization
Problem
EE-529 Power Converter Design

1)_Muhammad Zargham Ali


2)_Hamad Khan Afridi
8-13-2018
Matlab Code:
clc;clear all;close all;
syms wc ww hw dc lmda;
VA_rating=1000; %in Volt Ampers
Jrms=3.2e6; %in Ampers per squared meter
Bmax=1.2; %in Tesla
Kcu=0.8; %Copper fill factor
f=50; %in Hertz
AcAw=sqrt(2)*VA_rating/(Kcu*pi*f*Bmax*Jrms)*1e8;%in cm4
lang=obj_func(wc,ww,hw,dc)-lmda*(const_func(wc,ww,hw,dc,AcAw));
eq_wc= diff(lang,wc)==0;
eq_ww= diff(lang,ww)==0;
eq_dc= diff(lang,dc)==0;
eq_hw= diff(lang,hw)==0;
eq_lm= diff(lang,lmda)==0;
s=solve(eq_wc,eq_ww,eq_dc,eq_hw,eq_lm);
sym_solution=[s.wc s.ww s.hw s.dc];
solution=find_real_positive(eval(sym_solution));%solution in cm
Wc_optimal=solution(1,1);
Ww_optimal=solution(1,2);
Hw_optimal=solution(1,3);
Dc_optimal=solution(1,4);
Volume_optimal=obj_func(Wc_optimal,Ww_optimal,Hw_optimal,Dc_optimal);

Matlab Functions Codes:


find_real_positive:
function out=find_real_positive(ip)
[rows,cols]=size(ip);
op=zeros(rows,cols);
ireal=1;
for i=1:rows
for j=1:cols
if(imag(ip(i,j))==0 && real(ip(i,j))>0)
op(i,j)=real(ip(i,j));
ireal=i;
end
end
out=op(ireal,:);
end
end

obj_func:
function vol=obj_func(wc,ww,hw,dc)
vol=(2.*wc.*ww.*hw+2.*dc.*ww.*hw+3.14157.*ww.*ww.*hw)+(2.*wc.*dc.*hw+wc.*dc
.*2.*(wc+ww));
end

const_func:
function rem=const_func(wc,ww,hw,dc,AcAw)
rem=AcAw-hw.*dc.*ww.*wc;
end
Output:
Wc_optimal = 2.9153 cm

Ww_optimal = 2.6060 cm

Hw_optimal = 5.8307 cm

Dc_optimal = 6.6161 cm

Volume_optimal = 851.9576 cm3

You might also like