如何在Modelica中建模带y-z平面约束弹簧质量块的刚体？

我是Modelica新手，希望建模一个长圆柱刚体，通过弹簧连接一个点质量块。设定刚体固连坐标系的x轴为圆柱长轴，要求该点质量块仅能在y-z平面内平移。请问如何使用Modelica标准库实现该系统建模？以下是我的尝试代码（无法编译），望指点：

model RigidBody6DOFWithSpring
  // Parameters
  parameter Real mass_slosh = 90; // Mass of a single slosh mass
  parameter Real springStiffness = 1e4; // Spring stiffness
  parameter Real springRestLength = 1; // Spring rest length
  
  // Rigid body properties
  parameter Real mass_rigid = 10; // Mass of the rigid body
  parameter Modelica.Units.SI.Position r_CM[3](start={4,0,0});
  parameter Modelica.Units.SI.Inertia Ixx = 100;
  parameter Modelica.Units.SI.Inertia Iyy = 100;
  parameter Modelica.Units.SI.Inertia Izz = Iyy;

  // Model a world with no gravity
  inner Modelica.Mechanics.MultiBody.World world(g=0);

  // Rocket structure
  Modelica.Mechanics.MultiBody.Parts.Body rigidBody(m=mass_rigid,
                                                    r_CM = r_CM,
                                                    cylinderDiameter = 0.0, // Disable vis
                                                    sphereDiameter=0.5, // Disable vis
                                                    I_11 = Ixx,
                                                    I_22 = Iyy,
                                                    I_33 = Izz);
  
  // Joints to constrain point mass motion to body frame y-z motion only
  import Modelica.Mechanics.MultiBody.Joints;
  Joints.Constraints.Prismatic slosh1_constraint(
    x_locked=true,
    y_locked=false,
    z_locked=false);
  
  
  // Point masses
  parameter Real slosh1_x = 8;
  Modelica.Mechanics.MultiBody.Parts.Body sloshMass1(m=mass_slosh,
                                                     sphereDiameter=1);
  Modelica.Mechanics.MultiBody.Forces.Spring spring1(c=springStiffness,
                                                     s_unstretched=springRestLength);
  Modelica.Mechanics.MultiBody.Parts.FixedTranslation attachmentPoint1(
    r={slosh1_x,0,0}
  );
  
initial equation
  // Initial rigid body state
  rigidBody.w_0_start = {0, 0, 0};

  // Initial slosh positions in world space
  sloshMass1.r_0 = {slosh1_x, 0, 0};

  // Spring lengths
  spring1.length = 2;

equation
  // Define translation of spring mounting frame w.r.t rigid body frame
  connect(rigidBody.frame_a, attachmentPoint1.frame_a);
  connect(spring1.frame_a, attachmentPoint1.frame_b);
  connect(spring1.frame_b, sloshMass1.frame_a);
  
  connect(slosh1_constraint.frame_a, rigidBody.frame_a);
  connect(slosh1_constraint.frame_b, sloshMass1.frame_a);
  connect(slosh1_constraint.frame_b, attachmentPoint1.frame_b);
  connect(slosh1_constraint.frame_b, spring1.frame_b);
end RigidBody6DOFWithSpring;

问题诊断 #

• 连接冲突：同一frame端口被多个组件重复连接（slosh1_constraint.frame_b同时连了3个端口），违反MultiBody组件的连接规则。
• 初始条件错误：rigidBody.w_0_start是内部隐藏变量，无法直接赋值；spring1.length是输出变量，不能直接设置初始值。
• 可视化参数矛盾：同时禁用圆柱显示、启用球体显示，不符合长圆柱刚体的建模需求。

修正后的代码 #

model RigidBody6DOFWithSpring
  // 参数定义
  parameter Real mass_slosh = 90; // 点质量块质量
  parameter Real springStiffness = 1e4; // 弹簧刚度
  parameter Real springRestLength = 1; // 弹簧原长
  
  // 刚体属性
  parameter Real mass_rigid = 10; // 刚体质量
  parameter Modelica.Units.SI.Position r_CM[3](start={4,0,0});
  parameter Modelica.Units.SI.Inertia Ixx = 100; // x轴转动惯量（长轴）
  parameter Modelica.Units.SI.Inertia Iyy = 100; // y轴转动惯量
  parameter Modelica.Units.SI.Inertia Izz = Iyy; // z轴转动惯量

  // 无重力世界坐标系
  inner Modelica.Mechanics.MultiBody.World world(g=0);

  // 长圆柱刚体（x轴为长轴）
  Modelica.Mechanics.MultiBody.Parts.Body rigidBody(
    m=mass_rigid,
    r_CM = r_CM,
    cylinderDiameter = 0.5, // 圆柱直径，启用可视化
    cylinderLength = 10, // 圆柱长度，体现长轴特性
    sphereDiameter=0, // 禁用球体可视化
    I_11 = Ixx,
    I_22 = Iyy,
    I_33 = Izz);
  
  // 约束点质量仅在刚体y-z平面平移
  Modelica.Mechanics.MultiBody.Joints.Constraints.Prismatic sloshConstraint(
    x_locked=true, // 锁定x轴平移（刚体固连坐标系）
    y_locked=false,
    z_locked=false);
  
  // 点质量块
  Modelica.Mechanics.MultiBody.Parts.Body sloshMass(
    m=mass_slosh,
    sphereDiameter=1);
  
  // 弹簧
  Modelica.Mechanics.MultiBody.Forces.Spring spring(
    c=springStiffness,
    s_unstretched=springRestLength);
  
  // 刚体上的弹簧连接点（沿x轴偏移slosh_x）
  parameter Real slosh_x = 8;
  Modelica.Mechanics.MultiBody.Parts.FixedTranslation attachmentPoint(
    r={slosh_x,0,0});
  
initial equation
  // 初始化刚体角速度
  rigidBody.frame_a.w_0 = {0, 0, 0};
  
  // 初始化点质量在约束坐标系中的位置（y=0,z=0，弹簧处于原长）
  sloshConstraint.s = {0, 0, 0};

equation
  // 连接刚体与弹簧挂载点
  connect(rigidBody.frame_a, attachmentPoint.frame_a);
  
  // 弹簧一端连挂载点，另一端连约束的活动端
  connect(spring.frame_a, attachmentPoint.frame_b);
  connect(spring.frame_b, sloshConstraint.frame_b);
  
  // 约束固定端连刚体，活动端连点质量
  connect(sloshConstraint.frame_a, rigidBody.frame_a);
  connect(sloshConstraint.frame_b, sloshMass.frame_a);
end RigidBody6DOFWithSpring;

关键修正说明 #

• 解决连接冲突：每个frame端口仅连接一个组件，通过约束的活动端串联弹簧和点质量，避免重复约束。
• 修正初始条件：通过sloshConstraint.s（约束坐标系内的位移）初始化点质量位置，通过rigidBody.frame_a.w_0设置刚体初始角速度，均为合法可赋值变量。
• 优化可视化：启用圆柱显示并设置长度，明确体现长圆柱刚体的结构，禁用球体显示。
• 约束逻辑验证：Prismatic约束的frame_a与刚体固连，因此x_locked=true直接锁定刚体固连坐标系的x轴平移，实现点质量仅在y-z平面运动的要求。

内容的提问来源于stack exchange，提问作者Tombombadilly