oestersoetebier / IdealizedContact

A Modelica library for idealized simulation of contact phenomena
6 stars 6 forks source link

The initial conditions are not fully specified. #23

Open tbeu opened 9 years ago

tbeu commented 9 years ago

The initial conditions are not fully specified. This is the Dymola 2015 FD01 output of Example5.

Translation of IdealizedContact.Examples.Example5:

The DAE has 39119 scalar unknowns and 39119 scalar equations.

The initial conditions for variables of type Real are not fully specified.
Dymola has selected default initial conditions.

Assuming fixed start value for the continuous states:
contact.contactDefinition.geometry[1].contact_point_cylinder.v_0[1](start = 0)
contact.contactDefinition.geometry[1].contact_point_cylinder.v_0[2](start = 0)
contact.contactDefinition.geometry[1].pos_offset.v(start = 0)
contact.contactDefinition.geometry[2].contact_point_cylinder.v_0[1](start = 0)
contact.contactDefinition.geometry[2].pos_offset.v(start = 0)
cylinder.body.w_a[1](start = cylinder.body.R_start.T[1, 1]*cylinder.body.w_0_start[1] +cylinder.body.R_start.T[1, 2]*cylinder.body.w_0_start[2]+cylinder.body.R_start.T[1, 3] *cylinder.body.w_0_start[3])
cylinder.body.w_a[2](start = cylinder.body.R_start.T[2, 1]*cylinder.body.w_0_start[1] +cylinder.body.R_start.T[2, 2]*cylinder.body.w_0_start[2]+cylinder.body.R_start.T[2, 3] *cylinder.body.w_0_start[3])
cylinder.body.w_a[3](start = cylinder.body.R_start.T[3, 1]*cylinder.body.w_0_start[1] +cylinder.body.R_start.T[3, 2]*cylinder.body.w_0_start[2]+cylinder.body.R_start.T[3, 3] *cylinder.body.w_0_start[3])

Statistics

Original Model
Number of components: 526
Variables: 3269
Constants: 17 (17 scalars)
Parameters: 593 (847 scalars)
Unknowns: 2659 (39119 scalars)
Differentiated variables: 93 scalars
Equations: 2057
Nontrivial: 1604

Translated Model
Constants: 3904 scalars
Free parameters: 157 scalars
Parameter depending: 30933 scalars
Continuous time states: 20 scalars
Time-varying variables: 912 scalars
Alias variables: 4077 scalars
Assumed default initial conditions: 8
Number of mixed real/discrete systems of equations: 0
Sizes of linear systems of equations: {197, 159}
Sizes after manipulation of the linear systems: {26, 7}
Sizes of nonlinear systems of equations: {191}
Sizes after manipulation of the nonlinear systems: {4}
Number of numerical Jacobians: 0

Initialization problem
Sizes of linear systems of equations: {2, 2, 4, 49}
Sizes after manipulation of the linear systems: {2, 2, 4, 3}
Sizes of nonlinear systems of equations: {13}
Sizes after manipulation of the nonlinear systems: {4}
Number of numerical Jacobians: 0

Selected continuous time states

Statically selected continuous time states
contact.contactDefinition.geometry[1].filter.x[2]
contact.contactDefinition.geometry[1].pos_offset.s
contact.contactDefinition.geometry[2].filter.x[2]
contact.contactDefinition.geometry[2].pos_offset.s
freeMotion.r_rel_a[1]
freeMotion.r_rel_a[2]
freeMotion.r_rel_a[3]

Dynamically selected continuous time states
There are 2 sets of dynamic state selection.
From set 1 there are 3 states to be selected from:
cylinder.body.Q[1]
cylinder.body.Q[2]
cylinder.body.Q[3]
cylinder.body.Q[4]

From set 2 there are 10 states to be selected from:
contact.contactDefinition.geometry[1].contact_point_cylinder.v_0[1]
contact.contactDefinition.geometry[1].contact_point_cylinder.v_0[2]
contact.contactDefinition.geometry[1].contact_point_cylinder.v_0[3]
contact.contactDefinition.geometry[1].contact_point_rectangle.v_0[1]
contact.contactDefinition.geometry[1].contact_point_rectangle.v_0[2]
contact.contactDefinition.geometry[1].contact_point_rectangle.v_0[3]
contact.contactDefinition.geometry[1].filter.x[1]
contact.contactDefinition.geometry[1].pos_offset.v
contact.contactDefinition.geometry[2].contact_point_cylinder.v_0[1]
contact.contactDefinition.geometry[2].contact_point_cylinder.v_0[2]
contact.contactDefinition.geometry[2].contact_point_cylinder.v_0[3]
contact.contactDefinition.geometry[2].contact_point_rectangle.v_0[1]
contact.contactDefinition.geometry[2].contact_point_rectangle.v_0[2]
contact.contactDefinition.geometry[2].contact_point_rectangle.v_0[3]
contact.contactDefinition.geometry[2].filter.x[1]
contact.contactDefinition.geometry[2].pos_offset.v
cylinder.body.w_a[1]
cylinder.body.w_a[2]
cylinder.body.w_a[3]
der(cylinder.body.Q[1])
der(cylinder.body.Q[2])
der(cylinder.body.Q[3])
der(cylinder.body.Q[4])
freeMotion.v_rel_a[1]
freeMotion.v_rel_a[2]
freeMotion.v_rel_a[3]

Finished

// experiment StartTime=0 StopTime=1

WARNING: 1 warning was issued
oestersoetebier commented 9 years ago

Thanks. These warnings occur quite often. But since the default conditions are normally ok I didnt take care of it yet. It'll be done in the next release at the latest.

tbeu commented 9 years ago

Yes, that's fine for me.