%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