This repository presents a high-fidelity, mathematically rigorous 6-axis robotic arm simulation. Initiated by an asynchronous and detailed GUI (based on PyQt5), the system orchestrates high-level motion using MoveIt2, leveraging TRAC-IK for robust numerical solutions to inverse kinematics and the Pilz trajectory planning engine for deterministic Cartesian paths (LIN/CIRC). At the core, Controller Chaining (based on ros2_control) interpolates these trajectories into a tight 1000Hz reference for a custom Computed Torque Control (CTC) controller. To ensure physical realism, the CTC bypasses MuJoCoโs default PID joints to run the physics engine in pure effort mode. The CTC processes the feedback states via a Kalman Filter state observation model to reconstruct noise-free joint velocities, and executes Pinocchio dynamics calculations to compute rigid-body inertia, gravity, and Coriolis matrices. Combined with friction/armature compensation (smoothly neutralizing viscous damping, Coulomb stiction, and motor rotor inertia), this architecture achieves a tracking precision of
To bridge the gap between nominal simulation parameters and real-world physical discrepancies, an automated system identification (SysID) framework is integrated to extract joint-level viscous/Coulomb friction and motor armature coefficients, which are difficult to measure directly during manufacturing. Utilizing a bounded Fourier excitation trajectory coupled with an offline Least-Squares solver, the identified parameters align closely with the ground-truth physical properties defined in the URDF and MuJoCo XML files. Furthermore, the architecture incorporates a robust software emergency stop (E-Stop) and active recovery state machine. Upon triggering, it overrides high-level commands to execute precise, limit-aware deceleration trajectories that safely halt the manipulator, while providing a seamless recovery sequence back to the active closed-loop tracking state.
To facilitate real-time diagnostics and visual verification of tracking performance, the framework incorporates a synchronized "shadow robot" virtual twin. This visualization layer is driven directly by the real-time C++ CTC controller via a lock-free RealtimePublisher. It renders the exact internal target commands (q_target_) at 1000Hz (downsampled to 200Hz for RViz), representing the manipulator's theoretical kinematic state under zero-disturbance conditions. The shadow robot depicts pure kinematic targets, while the closed-loop CTC controller forces the physical manipulator to overcome physical disturbances (gravity, Coriolis, joint friction, etc.) to track this reference. Under standard operating conditions (including software E-Stop), the physical robot tracks the shadow twin tightly (deviation is indistinguishable to the human eye). However, there will be significant visible tracking deviation during System Identification because the friction and motor armature compensations are temporarily deactivated, allowing the algorithm to isolate and extract these specific raw unmodeled dynamics.
For detailed mathematical derivations, please refer to Theory.pdf.
โ ๏ธ Notice: The robot model used in this simulation is the UFactory Lite6, sourced from the official MuJoCo Menagerie. This repository is an independent open-source project designed for research and educational purposes.
- System Architecture
- Core Features & Demonstrations
- Quick Start & Reproduction Guide
- Future Work
- Acknowledgements
- Disclaimer
- Contact
The system's architecture is illustrated in the diagram below:
graph TB
%% Style Definitions
classDef PyQt fill:#34495e,stroke:#2c3e50,stroke-width:2px,color:#fff;
classDef MoveIt fill:#2980b9,stroke:#1f3a52,stroke-width:2px,color:#fff;
classDef R2C fill:#8e44ad,stroke:#5b2c6f,stroke-width:2px,color:#fff;
classDef Custom fill:#d35400,stroke:#a04000,stroke-width:2px,color:#fff;
classDef Physics fill:#27ae60,stroke:#1e8449,stroke-width:2px,color:#fff;
classDef Util fill:#7f8c8d,stroke:#5d6d7e,stroke-width:2px,color:#fff;
%% Subgraph: PyQt5 HMI Layer (Asynchronous / Non-Real-Time)
subgraph Sub_GUI [PyQt5 HMI Layer โ Async / ~20Hz]
UI[PyQt5 GUI Window]:::PyQt
Worker[ROS2Worker Thread]:::PyQt
SysID[SysIdEngine]:::PyQt
UI <-->|Signals/Slots<br/>Bidirectional Sync| Worker
Worker <-->|Least Squares Fitting| SysID
end
%% Subgraph: MoveIt2 Planning Space (~50Hz On-Demand)
subgraph Sub_MoveIt [MoveIt2 Planning Space โ ~50Hz]
MGroup[MoveGroup Node]:::MoveIt
TracIK[TRAC-IK Kinematics Solver]:::MoveIt
Pilz[Pilz Industrial Motion Planner<br/>PTP / LIN / CIRC]:::MoveIt
MGroup <-->|Plugin API| TracIK
MGroup <-->|Plugin API| Pilz
end
%% Subgraph: Hard Real-Time ROS 2 Control Loop (1000Hz)
subgraph Sub_R2C [ROS 2 Control Manager โ Hard Real-Time 1000Hz Loop]
JTC[JointTrajectoryController<br/>Upstream Interpolator]:::R2C
CTC[Lite6CTCController<br/>Downstream CTC Controller]:::Custom
JSB[JointStateBroadcaster]:::R2C
KF[1D Kalman Filter<br/>State Observer]:::Custom
Pinocchio[Pinocchio RBD Engine<br/>RNEA Solver]:::Custom
%% Chaining Mechanism
JTC ==>|Memory Pointer Sharing<br/>q_d, dq_d, ddq_d| CTC
CTC <-->|computes M, C, G| Pinocchio
KF -->|Filtered q, dq| CTC
end
%% Subgraph: MuJoCo Simulation Environment (1000Hz)
subgraph Sub_MuJoCo [MuJoCo Simulation Environment โ 1000Hz]
HW[Lite6MujocoSystem<br/>Hardware Interface]:::Custom
Engine[MuJoCo<br/>Effort Mode]:::Physics
HW <-->|mj_step / C API| Engine
end
%% Subgraph: Shadow Robot Twin (Visualization Helper)
subgraph Sub_Display [Robot Graphical Display โ 200Hz]
ShadowRSP[Shadow Robot State Publisher]:::Util
RViz[RViz2 Visualization]:::Util
end
%% System Connections (Data Flow)
Worker -->|Action Goal:<br/>MoveGroup.Action| MGroup
MGroup ==>|Action Goal:<br/>FollowJointTrajectory| JTC
Worker -.->|Software E-Stop<br/>Injects zero-velocity pt| JTC
%% CTC Inputs & Outputs
HW -->|Hardware State Interface<br/>Raw q_measured| KF
HW -->|Hardware State Interface| JSB
CTC ==>|Hardware Command Interface<br/>Output Torques tau_cmd| HW
%% GUI State Feedback & SysID
JSB -.->|Topic /joint_states| Worker
Worker -.->|SysID Mode Commands:<br/>Int32 /system_cmd| CTC
%% Shadow Robot Flow (Direct from CTC)
CTC -.->|"RT Publisher (1000Hz -> 200Hz)"<br/>Topic /shadow/joint_states| ShadowRSP
ShadowRSP -.->|Shadow TFs| RViz
JSB -.->|Main Robot TFs| RViz
%% Frequencies & Details annotations
style Sub_GUI fill:#f4f6f9,stroke:#bdc3c7,stroke-width:2px;
style Sub_MoveIt fill:#ebf5fb,stroke:#a9cce3,stroke-width:2px;
style Sub_R2C fill:#f5eef8,stroke:#d7bde2,stroke-width:2px;
style Sub_MuJoCo fill:#eaf2f8,stroke:#a9dfbf,stroke-width:2px;
style Sub_Display fill:#f2f4f4,stroke:#ccd1d1,stroke-width:2px;
- PTP (Point-to-Point): Joint space planning to specific angles or Cartesian poses (MoveP) via TRAC-IK:
- MoveL: Deterministic linear Cartesian interpolation via the Pilz planner:
- MoveC: Circular Cartesian interpolation using an auxiliary midpoint frame:
(Note: For MoveL and MoveC, it is recommended to lower the velocity and acceleration scale, to make sure that the trajectory planning doesn't exceed the robot's physical limits; otherwise, trajectory planning may fail.)
- The automated system identification executes bounded Fourier excitation trajectories, records
$q, \dot{q}, \tau$ , and utilizes Least Squares Optimization to extract exact Armature, Viscous Friction, and Coulomb Friction matrices.
- This is a collision-free, cyan-colored digital twin that runs alongside the main robot. It reflects the controller's target commands, allowing instant visual verification of tracking error and dynamic deviations.
- The emergency stop overrides MoveIt and generates a safe deceleration trajectory based on physical kinematic limits.
You can run this project using either a containerized Docker environment (recommended for avoiding dependency conflicts) or natively on Ubuntu 22.04.
Download the repository and organize your workspace as follows:
~/lite6_ws/
โโโ mujoco-3.9.0/ # MuJoCo binaries (provided in this repository)
โโโ src/ # Source code (lite6_bringup, lite6_controllers, etc.)
โโโ Dockerfile # Docker configuration
โโโ run_docker.sh # Container boot script
โโโ README.md
This method isolates the environment and requires no local ROS 2 installation. Ensure you have a Linux host (Ubuntu 22.04 or higher), Docker and the NVIDIA Container Toolkit (recommended for GUI rendering and GPU acceleration).
1. Start the Docker Environment:
cd ~/lite6_ws
chmod +x ./run_docker.sh
./run_docker.sh(Note: The script automatically detects your GPU and chooses suitable drivers if you are using an Intel/AMD GPU or no GPU. It also automatically mounts the src folder, so you do not need to rebuild the image when modifying code).
2. Build and Launch (Inside Docker):
colcon build --symlink-install
source install/setup.bash
ros2 launch lite6_bringup system_bringup.launch.pyIf you prefer running natively, ensure you have ROS2 Humble installed on your host machine.
1. Setup MuJoCo Binaries:
Move the provided MuJoCo binaries to your home directory and link them to your .bashrc.
mkdir -p ~/.mujoco
cp -r ~/lite6_ws/mujoco-3.9.0 ~/.mujoco/
# Add MuJoCo paths to your bash profile
echo 'export MUJOCO_DIR=~/.mujoco/mujoco-3.9.0' >> ~/.bashrc
echo 'export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$MUJOCO_DIR/lib' >> ~/.bashrc
source ~/.bashrc2. Install ROS 2 Dependencies (rosdep):
Use rosdep to automatically install all required packages (Pinocchio, MoveIt2, etc.).
cd ~/lite6_ws
rosdep update
rosdep install --from-paths src --ignore-src -r -y3. Compile and Launch:
colcon build --symlink-install
source install/setup.bash
ros2 launch lite6_bringup system_bringup.launch.pyThis framework is actively evolving. Upcoming features include:
- Analytical Inverse Kinematics via Symbolic Computation: Transitioning from the current numerical IK solver to a mathematically exact, closed-form analytical IK solution. Derived using symbolic computation software, this will reduce IK computation time from milliseconds to nanoseconds, guarantee thread safety, and provide all possible joint configurations instantly.
- Deterministic Trajectory Planning: Instead of using sampling-based trajectory planners, we combine analytical IK solutions with time parameterization to achieve deterministic trajectory planning for MoveJ, MoveP, MoveL, and MoveC commands.
- Collision Detection: When calculating IK and planning trajectories, the collision mesh of the robot model is computed to prevent the robot from colliding with itself or the worktable.
(Note: These new features have been implemented; please refer to link)
I would like to sincerely thank the creators and maintainers of the following open-source projects and organizations:
This is a personal, open-source project. I am not affiliated with, sponsored by, or endorsed by any commercial entities mentioned in this repository. All trademarks and registered trademarks are the property of their respective owners. The software is provided "as is", without warranty of any kind.
If you have any questions about this repository, you can raise an issue or contact directly.





