Add rotary mechanical elements to PyPHS about pyphs

[afalaize]

Hi @payam-rigi;
This warning indicates that the graph has been reduced by merging two components.

The default in PyPHS is to not merge except if a realizability issue is detected (i.e. the core PHS structure can not be derived from the graph due to components controlled by the same quantity).

To verify that this is your case, you can turn off the solving of anti-realizable components (abbreviated as arc) as follows (where netlist_filename should be replaced by the path to your netlist):

# read in Netlist object
netlist = Netlist(netlist_filename)

# Build Graph object
graph = netlist.to_graph()

# Build Core object
core = graph.to_core(solve_arc=False)

Can you provide your netlist and the result of the above code?
Also, you could make a PR for us to check your implementation in context.

Best regards;
Antoine

[payam-rigi]

netlist:
mechanics.mass M_chassis ('chassis',): M_chassis=20.0;
mechanics.inertia J_chassis ('chassis',): J_chassis=2.00;
mechanics.mass M_arm ('arm',): M_arm=0.5;
mechanics.inertia J_arm ('arm',): J_arm=2.0;
mechanics.mass M_wheel ('wheel',): M_wheel=0.5;
mechanics.stiffnessAngular Kr ('p1','arm'): Kr=10000;
mechanics.damperAngular Ar ('p1','arm'): Ar=1.0;
mechanics.damperAngular Conector ('arm','wheel'): Ar=0.0;
mechanics.source Sf_I ('#','wheel'): type=force;
mechanics.sourceAngular Svr_O ('chassis', 'p1'): type=angularvelocity;

Result of the above code:
The storage components [xM_arm, xJ_arm] have been merged into xM_armJ_arm
warnings.warn(merge_warning("storage", self.x_arc, X), Warning)
The storage components [xM_chassis, xJ_chassis] have been merged into xM_chassisJ_chassis
warnings.warn(merge_warning("storage", self.x_arc, X), Warning)

[payam-rigi]

rotaryDic.zip