Republic of the Philippines
BATANGAS STATE UNIVERSITY
The National Engineering University
Alangilan Campus
Golden Country Homes, Alangilan Batangas City, Batangas, Philippines 4200
Tel Nos.: (+63 43) 425-0139 local 2121 / 2221
E-mail Address:[Link]@[Link] | Website Address: [Link]
CHAPTER 3 - APPLICATION OF STATE-SPACE REPRESENTATION
PRE-LABORATORY SIMULATION (MATLAB OFFLINE)
Procedures:
1. Open the script in MATLAB and familiarize yourself with the initial system parameters
defined for the RLC circuit, Pneumatic Wastegate, and Fuel Solenoid Valve.
2. Navigate to Part 1 to complete the State-Space matrices for all three systems.
3. For the RLC Circuit, construct the A_rlc and B_rlc matrices using A_rlc = [0, 1/C_init;
-1/L_init, -R_init/L_init]; and B_rlc = [0; 1/L_init];.
4. For the Pneumatic Wastegate, construct the A_pneu and B_pneu matrices using A_pneu
= [0, 1; -k_init/m_init, -b_init/m_init]; and B_pneu = [0; 1/m_init];.
5. For the Fuel Solenoid Valve, construct the A_sol and B_sol matrices using A_sol = [0, 1;
-ks_init/ms_init, -(1 + 10/Rs_init)/ms_init]; and B_sol = [0; 2/(Rs_init*ms_init)];.
6. Run the script to launch the "ICE Multiphysics Digital Twin" dashboard. Use the
dropdown menu in the "System Selector & Controls" panel to switch between the three
different multiphysics domains.
7. Use the "Dynamic Physical Parameters" sliders to adjust the physical properties of the
currently selected system and observe the real-time changes in the State-Space matrices,
the 3D Render, the 2D Schematic, and the Time Domain Step Response.
Guide Questions:
1. Using the System Selector to switch between the three different systems, what
mathematical similarities do you observe in the structure of the A matrix across the
electrical, mechanical, and electromechanical domains?
2. While running the "3. Fuel Solenoid Valve" simulation, how does increasing the Spring
Stiffness (ks) parameter affect the maximum displacement (Plunger Pos) shown on the
Time Domain Step Response plot?
Republic of the Philippines
BATANGAS STATE UNIVERSITY
The National Engineering University
Alangilan Campus
Golden Country Homes, Alangilan Batangas City, Batangas, Philippines 4200
Tel Nos.: (+63 43) 425-0139 local 2121 / 2221
E-mail Address:[Link]@[Link] | Website Address: [Link]
LABORATORY SIMULATION (MATLAB ONLINE)
Procedures:
1. Open the script within your MATLAB Online workspace.
2. Review STEP 1 to identify the physical parameters defining the electrical (R,L,Ke) and
mechanical (Kt, J,b ) domains of the DC motor.
3. Navigate to STEP 2 to observe how these distinct physical domains are mathematically
coupled into a single continuous-time system matrix, A_baseline. Pay attention to the
mathematical formulas used to construct each row.
4. Check STEP 3 to see how simulating physical wear and tear (doubling the mechanical
friction b) translates into modifying specific elements of the mathematical matrix to
create A_modified.
5. Run the script. Review the numerical matrices printed in your Command Window, then
analyze the generated Figure window to visually compare the shift in the system's
eigenvalues (poles) due to the increased friction.
Guide Questions:
1. Based on the mathematical equations defining A_baseline in STEP 2, which specific
matrix element (Row, Column) represents the direct physical coupling where the
𝑑ω
electrical Armature Current (i) influences the mechanical angular acceleration ( 𝑑𝑡
), and
mathematically, what would happen to this matrix element's value if you significantly
increased the motor's Rotor Inertia (J)?
2. Observing the plotted complex plane in STEP 4, along which specific axis do the red
modified poles shift compared to the blue baseline poles when the mechanical friction is
doubled, and what does this geometric movement further away from the imaginary axis (
𝑗ω) imply about how quickly the motor's transient response will settle?