simulink_sfcl-master

sfcl superconducting fault current limiter

/**
 * SFCL resistance characteristic
 *
 * Copyright （c） 2011 Steven Blair
 *
 * Permission is hereby granted free of charge to any person obtaining a copy
 * of this software and associated documentation files （the “Software“） to deal
 * in the Software without restriction including without limitation the rights
 * to use copy modify merge publish distribute sublicense and/or sell
 * copies of the Software and to permit persons to whom the Software is
 * furnished to do so subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included in
 * all copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED “AS IS“ WITHOUT WARRANTY OF ANY KIND EXPRESS OR
 * IMPLIED INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM DAMAGES OR OTHER
 * LIABILITY WHETHER IN AN ACTION OF CONTRACT TORT OR OTHERWISE ARISING FROM
 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 * THE SOFTWARE.
 */

function [Rsfcl E J Tsc state] = fcn（current Ta Ts）

% superconductor dimensions and thermal properties
sc_length = 50.0;
sc_diameter = 0.004;                            % superconductor diameter （m）
sc_area = pi * （sc_diameter / 2）^2;             % superconductor cross-sectional area （m^2）
Cv = 1.0e6;
k = 1.5e3;
Csc = sc_length * sc_area * Cv;
Rsc = 1 / （k * sc_length * sc_diameter * pi）;   % includes the total surface area （minus the “ends“）
%nominal_resistance = rho_Tc * sc_length / sc_area;  % SFCL resistance at temperature=Tc

% superconductor characteristics
Jc_77K=1.5e7;               % critical current density at 77K （A/m^2）
E0=0.1;                     % E-field for transition from superconducting state to flux-flow state （V/m）
Ec=1e-6 * 100;              % definition of E-field required for Jc （V/m）
alpha_77K=6.0;              % exponent value during superconducting state （dimensionless）
beta=3;                     % exponent value during flux-flow state （dimensionless）
Tc=95;                      % critical temperature （K）
rho_Tc=1.0e-6;              % superconductor resistivity at Tc （ohm-m^2）

% inputs from previous iterations
persistent T;               % the instantaneous temperature of the superconductor
persistent prevE;           % store the E-field value from the previous iteration
persistent RsfclNew;        % history of Rsfcl values for averaging
persistent Qremoved;
persistent Qsc;
persistent time;            % elapsed simulation time

if isempty（T）
    T = Ta;                 % Ta is the ambient temperature of the coolant; it is assumed that this does not change
end
if isempty（prevE）
    prevE = 0.0;
end
if isempty（time）
    time = 0.0;
end
if isempty（RsfclNew）
    RsfclNew = [1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6];
end
if isempty（Qremoved）
    Qremoved = 0.0;
end
if isempt

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-09-23 12:50  simulink_sfcl-master\
     文件         285  2011-09-23 12:50  simulink_sfcl-master\README
     文件       57008  2011-09-23 12:50  simulink_sfcl-master\SFCL.mdl
     文件        5578  2011-09-23 12:50  simulink_sfcl-master\sfcl_resistance_characteristic.m