Parametric Sweep for a Simscape Thermal Model

Set Up Variables

Set the required variables for the example.

Model = 'sltestProjectorFanSpeedExample';
Harness = 'FanSpeedTestHarness';
TestSuite = 'sltestProjectorFanSpeedTestSuite.mldatx';
open_system(Model);

Test Plan and System Requirements

This test demonstrates sweeping through several fan speeds to determine the optimal value. In short, the optimal fan speed results in the fastest response without damaging the system. In detail, the optimal fan speed:

The document sltestProjectorFanSpeedExampleRequirements.slreqx captures these detailed requirements and the test procedure.

Test-specific model items reside in the test harness, keeping main models free of unnecessary blocks, suitable for code generation, and suitable for integration with other models.

Open the Test File

Open the Test Manager to view the test suite controlling the parameter sweep.

open(TestSuite)

Description of the Test

The test investigates the transient and steady-state thermal characteristics of the system. The test sequence initializes the system to ambient temperature, then powers the projector lamp. When the system reaches a steady-state condition, the lamp switches off. This test is modeled in a test harness using a Test Sequence block. Run the following to open the test harness:

sltest.harness.open(Model,Harness);

Requirements Linking

The test suite contains links to the requirements document. You can view the requirements links by selecting the Fan Speed Parametric Study test suite in the Test Browser pane of the Test Manager. Then, expand the Requirements section.

The Test Sequence

Double-click the Test Sequence block to open the test sequence editor.

The T0out and T0in signals store the initial projector temperature at each test step.

PowerOnTime stores simulation time when the lamp signal activates. This facilitates subsequent data analysis.

The transition condition detects the steady-state condition. At steady-state, the system temperature change is a small fraction (Threshold) of the difference between the current projector temperature and the initial projector temperature at each step. This condition must hold for a minimum time DurationLimit, in this case 10 seconds.

You may link the steps in the test sequence blocks against prepopulated requirements in the requirements document sltestProjectorFanSpeedExampleRequirements.slreqx.

Description of the Parameter Sweep

The pre-load callbacks contain the command to set the fan speed for each test case under the Fan Speed Parametric Study test suite. The parameter overrides contain the command to recalculate fan airflow from fan speed, and then override the test harness parameter. You can view these commands in the Callbacks and Parameter Overrides section of each test case.

Run the Test

In the Test Browser, highlight Fan Speed Parametric Study and click Run. When the test suite simulation completes, open the results for each test case and select ProjectorTemp. View the results in the Test Manager.

Export the Data

You can export data from the Test Manager for post-processing. To perform post-processing in MATLAB, in the Results and Artifacts pane of the Test Manager, right-click Sim Output for each test case and select Export.

This example includes the exported data in four MAT files:

ProjectorTempFanSpeed800.mat
ProjectorTempFanSpeed1300.mat
ProjectorTempFanSpeed1800.mat
ProjectorTempFanSpeed2300.mat

Investigate Response Time and Maximum Projector Temperature

Since the test sequence transitions execute when the system reaches steady-state, and the fan speed changes the system response, the lamp activates at different simulation times for each of the four test cases. Simplify the graphical results analysis by plotting each response with the lamp activation at the same time.

Extract the lamp activation response data, and plot the system response for the four fan speeds. Evaluate the results against these criteria:

Load the results. At the command line, enter

DataAt800 = load('ProjectorTempFanSpeed800.mat');
DataAt1300 = load('ProjectorTempFanSpeed1300.mat');
DataAt1800 = load('ProjectorTempFanSpeed1800.mat');
DataAt2300 = load('ProjectorTempFanSpeed2300.mat');

The script ArrangeProjectorData.m arranges the temperature and power on data from the output for each run.

ArrangeProjectorData

The script PlotProjectorThermalResponse.m plots the thermal response of the projector after the lamp activates, for each of the fan speeds.

PlotProjectorThermalResponse

Results Interpretation

The results show that while the highest fan speed results in the lowest maximum temperature, it also takes the longest time to reach the lamp activation temperature. The lowest fan speed results in the fastest lamp activation, but the system exceeds the maximum specified temperature by a significant margin.

Fan speed = 1300 keeps the system under the maximum temperature spec, and the system also reaches lamp activation temperature approximately 3 seconds faster than with the highest fan speed.

close_system(Model,0);

clear Model;
clear Harness;
clear TestSuite;
close(figure(1));