OpenModelica仿真中关节位置限制及部件碰撞属性实现方法问询

一、转动副关节角度限位实现 #

针对你需要将转动副角度限制在-π/2到π/2区间的需求，有两种成熟实现方式：

• 直接调用标准库现成组件
  替换原有Revolute转动副为Modelica.Mechanics.MultiBody.Joints.Assemblies.RevoluteWithLimits组件，在参数面板直接配置如下参数即可：
  • phi_min：设置为-Modelica.Constants.pi/2
  • phi_max：设置为Modelica.Constants.pi/2
  • useLimits：设置为true开启限位功能
    若需要硬限位碰撞效果，可同时开启stopAtLimit参数，角度达到阈值时会直接锁死关节；如果需要更平滑的数值表现，保持默认的软限位模式即可，组件会自动施加弹簧阻尼力限制角度超出范围。
• 自定义限位逻辑实现
  如果不想修改原有关节结构，可以在转动副的两个旋转法兰之间添加自定义力矩逻辑，示例代码如下：

  model CustomRevoluteLimit
    import Modelica.Constants.pi;
    // 转动副两侧法兰接口
    Modelica.Mechanics.Rotational.Interfaces.Flange_a flange_a annotation(Placement(transformation(extent={{-110,-10},{-90,10}})));
    Modelica.Mechanics.Rotational.Interfaces.Flange_b flange_b annotation(Placement(transformation(extent={{90,-10},{110,10}})));
    // 限位参数配置
    parameter Real phi_min = -pi/2 "最小转角限制";
    parameter Real phi_max = pi/2 "最大转角限制";
    parameter Real limit_stiffness = 1e6 "限位刚度，根据系统惯量调整";
    parameter Real limit_damping = 1e3 "限位阻尼，根据系统惯量调整";
    // 内部变量
    Real rel_phi = flange_a.phi - flange_b.phi "相对转角";
    Real limit_torque "限位力矩";
  equation
    // 超出限位范围时施加反向抑制力矩
    limit_torque = if rel_phi < phi_min then limit_stiffness*(phi_min - rel_phi) + limit_damping*der(rel_phi)
                   else if rel_phi > phi_max then limit_stiffness*(phi_max - rel_phi) + limit_damping*der(rel_phi)
                   else 0;
    flange_a.tau + flange_b.tau = limit_torque;
  end CustomRevoluteLimit;
  

二、多体部件碰撞属性配置 #

MSL默认的多体示例没有开启碰撞检测，所以会出现实体穿透的情况，OpenModelica环境下有两种实现碰撞拦截的方案：

• 基础包围盒碰撞实现
  OpenModelica 1.18及以上版本自带基础碰撞检测支持，直接找到对应实体的Shape可视化组件，开启参数面板中的enableCollision选项，系统会自动根据实体的几何形状生成AABB轴向包围盒，检测到碰撞时自动施加接触力阻止穿透。需要注意根据系统实际情况调整contactStiffness和contactDamping参数，参数过大会导致仿真数值刚性过高、步长过小甚至收敛失败，参数过小会出现明显穿透。
• 自定义复杂碰撞逻辑
  如果需要实现非规则形状的精确碰撞，可以自己编写碰撞检测方程：首先计算两个实体表面的最小距离，当距离小于0时根据Hertz接触模型施加法向接触力和切向摩擦力，法向力计算公式参考：fn = contact_stiffness * delta^1.5 + contact_damping * der(delta)，其中delta为穿透深度。

内容的提问来源于stack exchange，提问作者Roberto A.