Difference between revisions of "MGED CMD simulate"

From BRL-CAD
(mged_simulation tutorial)
 
(set units to meters)
(26 intermediate revisions by 4 users not shown)
Line 1: Line 1:
 
[[category:tutorials]]
 
[[category:tutorials]]
  
BRL-CAD recently integrated a new command in its mged tool, the simulate command. Those are the steps one should follow in order to run a simulation in the mged tool.
+
BRL-CAD recently integrated a new command into '''mged''', the '''simulate''' command, which allows for rigid-body physics simulations. This tutorial demonstrates the usage of the '''simulate''' command with a simple scene.
  
First open the mged tool, if you have BRL-CAD installed on your Linux system it should be easy, just type mged in a terminal and a mged terminal should open.
+
{|align="right"
 +
|[[Image:Simulation_1.png|thumb|256px|Cube falling to the ground plane modeled and rendered in BRL-CAD]]
 +
|}
  
 +
First, open the '''mged''' interface and create a new database. Next we want to add some geometry to the database. In the command window we will now type the following series of commands:
 +
<br>
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
brlcad@brlcad: mged
+
mged> units m
 +
mged> in cube.s rpp -1 1 -1 1 -1 1
 +
mged> r cube.r u cube.s
 +
mged> attr set cube.r simulate::type region
 
</pre>
 
</pre>
 
|}
 
|}
 +
<br style="clear: both" />
 +
This will create a cube primitive and a containing region. We will use the region later to apply material appearances. We also set the '''simulate::type''' attribute so that the '''simulate''' command regards all primitives within this region as one object.
  
In the mged terminal select the File menu, and then on the New... (file dialog button). A window will appear and a name for your new database will be required. For the purpose of this tutorial we will name the new database, sim1. Introduce the chosen name (sim1) in the "File name:" field and hit save.
+
<br style="clear: both" />
Now you have opened a new database. Next we want to add some geometry to the database. Something that we will have to simulate on.  
+
Now we want to add a fixed "ground" object to our simulation. BRL-CAD's simulated gravity acts in the ''x''-''y'' axis, so our ground plane should be parallel to this axis.
In the mged command window we will now type the following series of commands :
+
<br>
 
 
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
mged> in
+
mged> in ground.s rpp -15 15 -15 15 -1 1
Enter name of solid: box
+
mged> r ground.r u ground.s
Enter solid type: rpp
+
mged> attr set ground.r simulate::type region
Enter XMIN, XMAX, YMIN, YMAX, ZMIN, ZMAX: -1 1 0 5 -1 1
 
 
</pre>
 
</pre>
 
|}
 
|}
  
The "in" command creates new geometry, it allows the user to type in the arguments to create a shape with the name provided by the second argument, in the example above is the name provided when "Enter name of solid:" appears on the screen. The arguments of the in command can be provided one by one, as in the example above or in one single line, like in the example below:
+
<br style="clear: both" />
  
 +
We will also add an attribute specifying that this object's position is fixed:
 +
<br>
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
mged> in box rpp -1 1 0 5 -1 1
+
mged> attr set ground.r simulate::mass 0.0
 
</pre>
 
</pre>
 
|}
 
|}
 +
<br style="clear: both" />
  
We have so far a box, now we want to add a ground plane to our simulation.
+
We can also give the cube an initial angular velocity (linear velocity can be applied in a similar way):
 
+
<br>
 
 
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
mged> in ground rpp -15 15 -1 1 -15 15
+
mged> attr set cube.r simulate::angular_velocity <2.0,-1.0,3.0>
 
</pre>
 
</pre>
 
|}
 
|}
 +
<br style="clear: both" />
  
You should see now 2 boxes on your screen, one perpendicular on the other one, with a little overlap between them.
+
Now we will add some material appearances to the objects:
Using the sed command we will select the box and translate it 100 units on the Y axis.
+
<br>
 
 
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
mged> sed box
+
mged> mater cube.r "plastic {tr 0.2 re 0.2}" 255 0 0 y
mged> tra 0 100 0
+
mged> mater ground.r "plastic {tr 0.2 re 0.2}" 0 0 255 y
mged> accept
 
 
</pre>
 
</pre>
 
|}
 
|}
 +
<br style="clear: both" />
  
To be able to run the simulate command our geometry shapes should be grouped in regions to do this we will use the "r" command and create 2 regions for the box and the ground shapes.
+
To specify our simulated system we will define a combination containing both objects. This combination will be modified in-place during the simulation.
 
+
We will also change the units of gravity for the simulation to be
 +
acceleration in mm/s/s due to the small scale of our model.
 +
<br>
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
mged> r region1.r u box
+
mged> comb system.c u cube.r u ground.r
mged> r region2.r u ground
+
mged> attr set system.c simulate::gravity <0,0,-9.80665>
 
</pre>
 
</pre>
 
|}
 
|}
 +
<br style="clear: both" />
  
Now we are ready to run the simulate command:
+
We will now select the cube and translate it 50 units "up" along the ''y''-axis.
 
+
<br>
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
mged> simulate 100
+
mged> B system.c
 +
mged> oed /system.c/cube.r cube.s
 +
mged> tra 0 0 50
 +
mged> accept
 
</pre>
 
</pre>
 
|}
 
|}
 +
<br style="clear: both" />
  
After the simulation runs its course, we will use the "who" command to find out what objects are being displayed.
+
In order to see what we did so far we should set the view for the geometry we have so far. We will use an "ae 35 15", that means the azimuth angle is at 35 degrees and the elevation at 15 degrees. The angle of elevation of an object as seen by an observer is the angle between the horizontal and the line from the object to the observer's eye (the line of sight).The azimuth of an object is the angular distance along the horizon to the location of the object.We will use "autoview" command to reset the view size and the view center such that all displayed objects are within the view. And then we will do a translation on the Y axis.
 
+
<br>
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
mged> who
+
mged> ae 35 15
 +
mged> autoview
 +
mged> tra 0 10 0
 
</pre>
 
</pre>
 
|}
 
|}
 
+
<br style="clear: both" />
You should get the following answer:
+
Now we can save this view for later:
 
+
<br>
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
box ground sim.c
+
mged> saveview ae3515.rt
 
</pre>
 
</pre>
 
|}
 
|}
 +
<br style="clear: both" />
  
Now, we don't really want to display the box and the ground shapes, we already have them in the sim.c group so we will unload all the geometry(using the "Z" command) and will load only the simulate geometry(using the "draw sim.c" command).
+
Now we are ready to run the simulate command. To preserve our configuration we will first create a copy of our combination.
 
+
<br>
 
{|align="left"
 
{|align="left"
 
  |<pre>
 
  |<pre>
mged> Z
+
mged> cp system.c output.c
mged> draw sim.c
+
mged> B output.c
 +
mged> simulate output.c 10.0
 
</pre>
 
</pre>
 
|}
 
|}
 +
<br style="clear: both" />
 +
This will evaluate 10 seconds of simulated physics.
  
 
+
In order to create an animation from your simulation, you will have to run the script below (note that it requires ImageMagick): (in progress)
Now to actually see the simulation we will raytrace it, using the "rt" command.
 
 
 
{|align="left"
 
|<pre>
 
mged> rt
 
</pre>
 
|}
 

Revision as of 16:34, 24 March 2017


BRL-CAD recently integrated a new command into mged, the simulate command, which allows for rigid-body physics simulations. This tutorial demonstrates the usage of the simulate command with a simple scene.

Cube falling to the ground plane modeled and rendered in BRL-CAD

First, open the mged interface and create a new database. Next we want to add some geometry to the database. In the command window we will now type the following series of commands:

mged> units m
mged> in cube.s rpp -1 1 -1 1 -1 1
mged> r cube.r u cube.s
mged> attr set cube.r simulate::type region


This will create a cube primitive and a containing region. We will use the region later to apply material appearances. We also set the simulate::type attribute so that the simulate command regards all primitives within this region as one object.


Now we want to add a fixed "ground" object to our simulation. BRL-CAD's simulated gravity acts in the x-y axis, so our ground plane should be parallel to this axis.

mged> in ground.s rpp -15 15 -15 15 -1 1
mged> r ground.r u ground.s
mged> attr set ground.r simulate::type region


We will also add an attribute specifying that this object's position is fixed:

mged> attr set ground.r simulate::mass 0.0


We can also give the cube an initial angular velocity (linear velocity can be applied in a similar way):

mged> attr set cube.r simulate::angular_velocity <2.0,-1.0,3.0>


Now we will add some material appearances to the objects:

mged> mater cube.r "plastic {tr 0.2 re 0.2}" 255 0 0 y
mged> mater ground.r "plastic {tr 0.2 re 0.2}" 0 0 255 y


To specify our simulated system we will define a combination containing both objects. This combination will be modified in-place during the simulation. We will also change the units of gravity for the simulation to be acceleration in mm/s/s due to the small scale of our model.

mged> comb system.c u cube.r u ground.r
mged> attr set system.c simulate::gravity <0,0,-9.80665>


We will now select the cube and translate it 50 units "up" along the y-axis.

mged> B system.c
mged> oed /system.c/cube.r cube.s
mged> tra 0 0 50
mged> accept


In order to see what we did so far we should set the view for the geometry we have so far. We will use an "ae 35 15", that means the azimuth angle is at 35 degrees and the elevation at 15 degrees. The angle of elevation of an object as seen by an observer is the angle between the horizontal and the line from the object to the observer's eye (the line of sight).The azimuth of an object is the angular distance along the horizon to the location of the object.We will use "autoview" command to reset the view size and the view center such that all displayed objects are within the view. And then we will do a translation on the Y axis.

mged> ae 35 15
mged> autoview
mged> tra 0 10 0


Now we can save this view for later:

mged> saveview ae3515.rt


Now we are ready to run the simulate command. To preserve our configuration we will first create a copy of our combination.

mged> cp system.c output.c
mged> B output.c
mged> simulate output.c 10.0


This will evaluate 10 seconds of simulated physics.

In order to create an animation from your simulation, you will have to run the script below (note that it requires ImageMagick): (in progress)