PDE-Toolbox: Time dependent applied Temperature

12 ビュー (過去 30 日間)
Andreas Kohl
Andreas Kohl 2019 年 8 月 5 日
回答済み: Kavya Vuriti 2019 年 8 月 8 日
Hello,
for the simulation of a brake disk, I want to apply the temperature for 7.2 seconds. Then there should only be convection, emissivity and Thermal Conductivity. How do I tell the toolbox to stop heating the disk after 7.2 seconds?
We've tried to put dT into a for-loop but the toolbox just takes the last known value.
Below there is a extract from our current program
%Variables
mm=200; %mass of the motorbike in kg
v0ms=250000/9; %velocity at the start in mm/s
c=460e6; %specific heat capacity
mb=0.43; %mass of the brake disk
thermalmodel = createpde('thermal','transient');
importGeometry(thermalmodel,'Model.stl');
generateMesh(thermalmodel);
thermalIC(thermalmodel,20);
%Equations
Ekin=0.5*mm*(v0ms^2); %kinetic energy of the bike
dT=Ekin/(c*mb); %Temperature change
tlist=0:1:20;
thermalBC(thermalmodel,'Face',95,'Temperature',dT); %Temperature, that should be applied for 7.2 seconds
%thermal losses
thermalmodel.StefanBoltzmannConstant=5.670373E-8;
thermalBC(thermalmodel,'Face',1:94,'ConvectionCoefficient',0.05,'AmbientTemperature',20);
thermalBC(thermalmodel,'Face',1:94,'Emissivity',0.4,'AmbientTemperature',20);%Konvektionskoeffizient
thermalBC(thermalmodel,'Face',96:100,'ConvectionCoefficient',0.05,'AmbientTemperature',20);
thermalBC(thermalmodel,'Face',96:100,'Emissivity',0.4,'AmbientTemperature',20);%Konvektionskoeffizient
thermalProperties(thermalmodel,'ThermalConductivity',lam,'MassDensity',rho,'SpecificHeat',c);
R = solve(thermalmodel,tlist);

回答 (1 件)

Kavya Vuriti
Kavya Vuriti 2019 年 8 月 8 日
The temperature boundary condition can be specified as a function handle to specify temperature that depends on time and it can be passed to the function thermalBC.
dT= @(location,state) {code to set dT according to time}
thermalBC(thermalmodel,'Face',95,'Temperature',dT)
Where state is structure array and state.time field contains time at evaluation points.
For more information, refer the following link:

製品


リリース

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by