Thread Subject: Simscape - Wrong solution to implicit quadratic equation

Subject: Simscape - Wrong solution to implicit quadratic equation

From: John Kreuder

Date: 11 Jun, 2009 19:56:02

Message: 1 of 1

I am using Simscape to solve for the steady state of a thermal system (Brayton cycle). This is being done as proof of concept so that I can do some detailed thermal analysis in the future. I have created a domain in which mass flow is a through variable and temperature and pressure are across variables. Performance curves for the compressor and turbine are available in the form of quadratic curve fit relationships. In simscape, I have the equation:

 input power [kw] == 1020-.383*outlet_pressure[Pa] + .00513*(outlet pressure)[Pa])^2
for an input pressure and temperature of 101kPa and 25 C.

The steady state solution for the system is known (solved outside of Simscape) and the input power is 1530 kw, which yeilds outlet pressure solutions of -280kPa and 355kPa. The compressor performance curve is only valid for discharge pressures of 175 to 490kPa, but I haven't found a way to work this into the simscape component equation. as a result, the simulation runs and solves for the wrong root of the quadratic equation above, causing the simulation to fail to converge.

Do you know of any way that I can set an interval over which to search for solutions to a quadratic equation in a Simscape component file? Do you have any other suggestions for how to ensure that my simulation always picks the solution that I need when solving equations with multiple solutions?

Tags for this Thread

Everyone's Tags:

Add a New Tag:

Separated by commas
Ex.: root locus, bode

What are tags?

A tag is like a keyword or category label associated with each thread. Tags make it easier for you to find threads of interest.

Anyone can tag a thread. Tags are public and visible to everyone.

Tag Activity for This Thread
Tag Applied By Date/Time
quadratic John Kreuder 11 Jun, 2009 15:59:09
wrong root John Kreuder 11 Jun, 2009 15:59:09
simscape John Kreuder 11 Jun, 2009 15:59:09
rssFeed for this Thread

Contact us at files@mathworks.com