DCT双离合变速箱

%DCT Dynamic Model

%Paul Walker
%30JAN08

function DCT_dynamic_Model_SD_V3a

%This program models the dynamic response of a DCT during transient shifts.
%This includes gear synchronisation as well as clutch to clutch shifting
%and launch of the vehicle.  The model is stylised on a FWD vehicle.

%The model incoporates engine dynamics including torsional vibrations a
%basic DMFW for torsional damping simplified dual clutch transmission
%（including 2 clutches and 2 gear pairs with synchronisers which inclides
%reduced transmission inertia for the gears and drag torque for damping）
%final drive output shafts and wheel components.

clear all

%1.0 Model Parameters and Constants
%1.1 Constants & Initial Conditions
%1.1.1 Time constants - IGNORE
%1.1.2 Gears and transmission states

G1=3;                   %Selected gear for clutch 1
G2=4;                   %Selected gear for clutch 2

%1.1.3 Engine and Vehicle speeds - IGNORE
%1.2 Dynamic model parameters
%1.1.05 Common Properties
Gsm=80e9;               %Shear modulus of steel （Pa）
rhosteel=787;           %density of steel （kg/m3）

%1.2.1 Engine Model
Ie=0.4;                                             %Engine inertia （kgm2）
De1=0.04;                                           %Engine shaft OD （m）
De2=0.02;                                           %Engine shaft ID （m）
Les=0.2;                                            %engine shaft length
Ke=（pi*Gsm*（De1^4-De2^4）/（32*Les））;                 %Output shaft stiffness （Nm/rad）

%1.2.2 DMFW Parameters
Idm1=0.2;               %First DMFW inertia （kgm2）
Idm2=0.2;               %Second DMFW inertia （kgm2）
Kdm=2e5;%ie not DMFW with 2e5 700                %DMFW spring rate （Nm/rad）
mu1=0.12;               %Friction coefficient for spring on housing
ms=1;                   %Mass of DMFW spring （kg）
rs=0.12;                %Outer spring radius in housing （m）
Cdm=2;%2*mu1*ms*rs;          %DMFW damping （Nmrad/s）
Thys=40;                %constant hysteresis torque

%1.2.3 Transmission Parameters
%1.2.3.1 Clutch
%C1 corresponds to odd gears
%C2 corresponds to even gears

mus=0.16;       %static friction
mud=0.14;       %dynamic friction

nf=8;           %number of friction surfaces

Ic1=0.0072;

if G1==1 || G2==1
    Kc1=8e4;
elseif G1==3 || G2==3
    Kc1=2.1e5;
else
    Kc1=7e4;
end

Rc1o=163/2;
Rc1i=105/2;

Ac1=pi*（Rc1o^2-Rc1i^2）;
Rm1=2/3*（Rc1o^3-Rc1i^3）/（Rc1o^2-Rc1i^2）;


Ic2=0.0125;

if G1==2 || G2==2
    Kc2=7e4;
else
    Kc2=8.7e5;
end

Rc2o=193/2;
Rc2i=139/2;

Ac2=pi*（Rc2o^2-Rc2i^2）;
Rm2=2/3*（Rc2o^3-Rc2i^3）/（Rc2o^2-Rc2i^2）;

h1=0.1e-3;              %Plate spacing （m）

%1.2.3.2 Gears and Synchronisers
Ip=[0.000074 0.00025 0.00056 0.0013 0.00096 0.0014];        %inertia of whole shaft and gearset （kgm2）
Ig=[0.0085 0.0049 0.0025 0.0014 0.0014 0.00085];            %approximated inertia of cylinders （kgm2）

GR=[3.45 2.05 1.45 1.08 1.11 0.92];

if G1==1 || G1==3