Simulation of a system a mechanical subsystem and a hydraulic subsystem

• US 20040020356A1
• Filed: 08/20/2003
• Published: 02/05/2004
• Est. Priority Date: 10/17/2000
• Status: Abandoned Application

• Current Assignee: Lumeo Software Oy
• Original Assignee: Lumeo Software Oy
• Inventors: Ojanen, Harri
• Application Number:
  US10/398,174
• Publication Number:
  US 20040020356A1
• Time in Patent Office: Days
• US Class Current: 91/361
• CPC Class Codes: G05B 17/02 electric

Abstract

A method (90) for simulating a system having a mechanical subsystem and a hydraulic subsystem is presented. The method comprises the steps of: describing (92) the mechanical subsystem and the effect of the hydraulic subsystem thereupon by integrating a first differential equation system of an n-dimensional first variable using a first numerical method, resulting in a sequence of values for the first variable; and describing (94) the hydraulic subsystem and the effect of the mechanical subsystem thereupon by integrating a second differential equation system of an m-dimensional second variable using a second numerical method, resulting in a sequence of values for the second variable. A corresponding arrangement having a first integrator (1) and a second integrator (4) is also presented.

Claims

1. A method (90) for simulating a system having a mechanical subsystem and a hydraulic subsystem, said mechanical subsystem having a number of bodies and said hydraulic subsystem having a number of hydraulic actuators with hydraulic cylinders, said number of bodies being at least one and said number of hydraulic actuators being at least one, characterized in that

the mechanical subsystem and the effect of the hydraulic subsystem thereupon is described by a first differential equation system of an n-dimensional first variable, the first variable comprising coordinates defining the configuration of said mechanical subsystem and the effect of the hydraulic subsystem on the mechanical subsystem being taken into account using actuator forces,

the hydraulic subsystem and the effect of the mechanical subsystem thereupon is described by a second differential equation system of an nm-dimensional second variable, the second variable relating to pressures in said hydraulic cylinders,

the first differential equation system is discretized using a first numerical integration method,

the second differential equation system is discretized using a second numerical integration method,

and in that for each time step of the simulation, after an initialization step, the method comprises the following steps, performed by means of a computer:

determining a current value for the second variable using the discretization of the second differential equation system, a previous value for the first variable and a previous value for the second variable, said previous values having been determined in a previous time step,

computing current actuator forces using at least the current value of the second variable, and

determining a current value for the first variable using the discretization of the first differential equation system, the current actuator forces and the previous value for the first variable.

2. A method according to claim 1, wherein

the hydraulic cylinders are coupled to a number of tanks and to a number of pumps with a number of control valves and flows relating to the control valves are calculated using a predefined flow formula,

the second differential equation system is uncoupled (93) by taking into account the interactions within the hydraulic subsystem using discrete event simulation, and

analytical expression for the discretization of the second differential equation system is determined (94a) using a predefined differential equation for the pressure in a chamber of a hydraulic cylinder and the flow corresponding to said hydraulic cylinder,

and wherein the step of determining a current value for the second variable comprises the substeps of:

calculating (23) current flows relating to the control valves by evaluating the predefined flow formula, and

evaluating (94b, 24) the current value for the second variable using the analytical expression and the current flows,

and wherein the step of computing current actuator forces comprises the substep of;

computing current cylinder forces, which form part of the actuator forces, by evaluating a predefined cylinder force formula and the current value for the second variable.

3. A method according to claim 2, wherein the predefined flow formula takes as inputs for a control valve in a hydraulic cylinder at least previously determined pressure in said hydraulic cylinder and information about previous value for the first variable,

4. A method according to claim 2 or 3, wherein a set of the hydraulic cylinders of said hydraulic subsystem are connected to a single pump and the step of determining a current value for the second variable comprises the substep of:

modifying (23) said flows relating to said set of hydraulic cylinders to take into account the couplings within the hydraulic subsystem and the capacity of the pump.

5. A method according to claim 4, wherein said flows are modified so that the ratio of the modified flows remains the same as the ratio of the unmodified flows and the sum of the modified flows is equal to a source flow produced by said pump.

6. A method according to claim 4, wherein said flows are modified so that a number of largest flows, the sum of said largest flows being less than or equal to a source flow produced by said pump and the number of largest flows being as large as possible, remains unmodified and modified flows corresponding to unmodified flows not belonging to said number of largest flows are zero flows.

7. A method according to claim 4, wherein said flows are modified so that according to a preference order said flows are retained as long as the sum of said retained flows is less than or equal to a source flow produced by said pump.

8. A method according to claim 4, wherein said flows are modified by using predetermined information about the flows in said set of hydraulic cylinders given certain positions for said control valves.

9. A method according to any of the claims 2 to 8, wherein the step of computing current actuator forces comprises the substep of:

calculating motion limiting forces and adding the calculated motion limiting forces to the cylinder forces.

10. A method according to any of the preceding claims, wherein the first differential equation system is a second order differential equation system and the step of determining current value for the first variable comprises a substep of:

calculating current numerical value for the acceleration of the first variable using said current actuator forces.

11. A method according to claim 10, wherein

the second order first differential equation system is modified into a first order differential equation system by replacing the n-dimensional first variable with a 2n-dimensional third variable containing the first variable and the first derivative of the first variable,

and the step of determining current value for the first variable comprises a substep of:

evaluating Jacobian relating to said first order differential equation system, and

calculating a current value for the third variable using said Jacobian and said current numerical value for the acceleration of the first variable.

12. A method according to claim 11, wherein the Jacobian is evaluated using divided differences.

13. A method according to claim 11 or claim 12, wherein the numerical value for the acceleration of the first variable is calculated using the Featherstone algorithm.

14. A method according to claim 11 or claim 12, wherein the numerical value for the acceleration of the first variable is derived from equations describing explicitly said mechanical subsystem.

15. A method according to any of the preceding claims, wherein the second numerical method is different from the first numerical method.

16. A method according to claim 15, wherein the first numerical method is a semi-implicit method for numerical integration and the second numerical method is an implicit method for numerical integration.

17. A method according to claim 15, wherein the first numerical method is a semi-implicit numerical method and the second numerical method is a semi-implicit method for numerical integration.

18. A method according to claim 15, wherein the first numerical method is an implicit method for numerical integration and the second numerical method is an implicit method for numerical integration.

19. A method according to claim 15, wherein the first numerical method is an implicit numerical method and the second numerical method is a semi-implicit method for numerical integration.

20. A method according to any of the preceding claims, wherein a second time step, which is shorter than the time step for the simulation, is used for determining intermediate values for the second variable in the step of determining the current value for the second variable, and wherein the time step for simulation is used in determining the current value for the first variable.

21. A method according to any of the preceding claims, further comprising a step of:

visualizing (18, 96) the behavior of said system during simulation.

22. A method according to any of the preceding claims, further comprising a step of:

controlling (91) said hydraulic subsystem using control signals.

23. A method according to claim 22, wherein the control signals are obtained from a separate control unit and the method further comprising a step of:

sending (95) information relating to current value of at least one of said first and second variables to said separate control unit for at least part of the simulation time steps.

24. A method according to claim 22 or 23, wherein the hydraulic cylinders of the hydraulic subsystem are connected to a number of tanks and to a number of pumps with a number of control valves and the control signals are controlling (17, 20, 21) the relative openings of said control valves.

25. A computer program for simulating a system having a mechanical subsystem and a hydraulic subsystem, said mechanical subsystem having a number of bodies and said hydraulic subsystem having a number of hydraulic actuators with hydraulic cylinders, said number of bodies being at least one and said number of hydraulic actuators being at least one, characterized in that said computer program comprises computer program code means adapted to perform, for each time step of the simulation after an initialization step, the following steps when the program is run on a computer:

determining a current value for a second variable using a discretization of a second differential equation system, a previous value for the second variable and a previous value for a first variable, where the second differential equation system of the m-dimensional second variable, which relates to pressures in said hydraulic cylinders, describes the hydraulic subsystem and the effect of the mechanical subsystem thereupon, the second differential equation system is being discretized using a second numerical integration method and said previous values of the first and second variables have been determined in a previous time step,

computing current actuator forces using at least the current value of the second variable, and

determining a current value for the first variable using a discretization of a first differential equation system, the current actuator forces and the previous value for the first variable, where the first differential equation system of the n-dimensional first variable, which comprises coordinates defining the configuration of said mechanical subsystem, describes the mechanical subsystem and the effect of the hydraulic subsystem thereupon, the effect of the hydraulic subsystem on the mechanical subsystem is being taken into account using the actuator forces, and the first differential equation system is being discretized using a first numerical integration method.

26. A computer program as claimed in claim 25 embodied on a computer readable medium.

27. A simulation arrangement for simulating a system having a mechanical subsystem and a hydraulic subsystem, said mechanical subsystem having a number of bodies and said hydraulic subsystem having a number of hydraulic actuators with hydraulic cylinders, said number of bodies being at least one and said number of hydraulic actuators being at least one, characterized in that it comprises

means (4) for determining a current value for a second variable using a discretization of a second differential equation system, a previous value for the second variable and a previous value for a first variable, where the second differential equation system of the nm-dimensional second variable, which relates to pressures in said hydraulic cylinders, describes the hydraulic subsystem and the effect of the mechanical subsystem thereupon, the second differential equation system is being discretized using a second numerical integration method and said previous values of the first and second variables have been determined in a previous time step,

means for computing current actuator forces using at least the current value of the second variable, and

means (1, 3) for determining a current value for the first variable using a discretization of a first differential equation system, the current actuator forces and the previous value for the first variable, where the first differential equation system of the n-dimensional first variable, which comprises coordinates defining the configuration of said mechanical subsystem, describes the mechanical subsystem and the effect of the hydraulic subsystem thereupon, the effect of the hydraulic subsystem on the mechanical subsystem is being taken into account using actuator forces, and the first differential equation system is being discretized using a first numerical integration method,

and in that said means for determining the current value for the second variable, said means for computing the current actuator forces the and said means for determining the current value for the first variable are adapted to be active for each time step of the simulation after an initialization step.

28. An arrangement according to claim 27, further comprising

means (2) for visualizing behavior of said system.

29. An arrangement according to claim 27 or 28, further comprising

means (6) for controlling said system using control signals, and

means (5) for communicating said control signals to said means (4) for determining the current value for the second variable,

and wherein said means (4) for determining the current value for the second variable comprise means for mapping said control signals to openings of control valves belonging to said hydraulic part.

30. An arrangement according to claim 29, wherein said means (6) for controlling said system is a computer program.

31. An arrangement according to claim 29, wherein said means (6) for controlling said system is a control unit for controlling a physical device, the control unit being connected to other means forming said arrangement.

32. An arrangement according to claim 31, wherein said other means forming said arrangement are realized as a set of computer programs.

33. An arrangement according to any one of claims 27 to 32, wherein at least one of said means (1, 3) for determining the current value for the first variable and said means (4) for determining the current value for the second variable is arranged to send information relating to the current value of at least one of said first and second variable for at least part of the simulation time steps to said means (6) for controlling said system.

34. An arrangement according to any one of claims 27 to 33, wherein said arrangement is realized as a set of computer programs.