CN106709264B - 双起重机系统变幅角响应建模算法及随机响应域预测方法 - Google Patents
双起重机系统变幅角响应建模算法及随机响应域预测方法 Download PDFInfo
- Publication number
- CN106709264B CN106709264B CN201710019654.2A CN201710019654A CN106709264B CN 106709264 B CN106709264 B CN 106709264B CN 201710019654 A CN201710019654 A CN 201710019654A CN 106709264 B CN106709264 B CN 106709264B
- Authority
- CN
- China
- Prior art keywords
- random
- autocrane
- phase angle
- variable phase
- luffing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Control And Safety Of Cranes (AREA)
- Jib Cranes (AREA)
Abstract
本发明公开了一种双起重机系统变幅角响应建模算法及随机响应域预测方法,变幅角响应建模算法步骤如下:1、建立系统几何模型,并给出各点的坐标;2、建立系统约束方程;3、结合系统约束方程,建立变幅角响应模型。随机响应域预测方法步骤如下:1、根据随机参数模型来描述载荷的不确定性;2、基于变幅角响应模型和随机参数模型,建立随机变幅角响应等效方程;3、结合复合函数特性和随机摄动法,提出了摄动随机复合函数法,并求解出随机变幅角响应表达式;4、根据随机变量函数矩法,进一步求解出变幅角响应的期望和方差。本发明可解决带有随机参数下变幅角响应域的预测问题,具有快速性和精度高的特点,保证了系统作业的可靠性。
Description
技术领域
本发明涉及可靠性技术领域,具体是一种双起重机系统变幅角响应建模算法 及随机响应域预测方法。
背景技术
汽车起重机是实现物料搬运机械化的重要工具之一,广泛应用于大型重载起 重吊装作业中。近年来,随着单台起重机在结构和控制等方面的快速发展,使得 双台汽车起重机系统联合吊装作业成为可能。然而,这也增加了起吊作业的复杂 程度和危险性。根据GB6067.1-2010《起重机械安全规程》,“特殊情况需要两台 或两台以上的起重机共同起吊重物时,各台起重机的升降、运行应保持同步”。 而这类问题的前提和关键是首先要确定多台起重机系统在特定运动下的响应模 型。
现有关于多台起重机系统在某一运动下的响应模型,针对的都是确定性载荷 参数来分析,而没有考虑载荷参数的不确定性。根据GB6067.1-2010《起重机械 安全规程》中多台起重机械联合起升总则,“在多台起重机械的联合起升操作中, 由于起重机械之间的相互运动可能产生作用于起重机械、物品和吊索具上的附加 载荷,而这些附加载荷的监控是困难的”。根据GB3811-2008《起重机设计规范》 中起重机整体抗倾覆稳定性,“对于固定的起重机,在具体使用现场或地区如有 地震或其他的基础外部激励效应,则在相应的工作状态或非工作状态抗倾覆稳定 性的核算中,将其作为附加的载荷情况加以考虑”。另外,对于一个复杂的系统 而言,输入参数即使发生很小的变化,也会给系统响应造成明显地波动,甚至导 致严重的事故。因此,将不确定载荷参数建立成随机模型,并研究随机模型下的响应分布情况具有重要的意义。随机理论在其他领域,如结构学、热学以及声学 等已经取得了一定的成果,然而在多起重机工程应用领域才刚刚起步。另外,以 摄动理论为基础的随机摄动法在多层复合函数的应用也是一片空白。
因此,如何计算双台汽车起重机系统变幅运动下变幅角响应以及快速预测随 机载荷下变幅角响应域,是目前起重机可靠性学术领域的一个重点,对于弥补现 有双台汽车起重机系统变幅运动中确定性参数下响应模型和随机参数下响应域
发明内容
本发明的目的是提供一种双起重机系统变幅角响应建模算法及随机响应域 预测方法,以解决现有技术中双台汽车起重机系统变幅运动中确定性参数下响应 模型难以建立及随机参数下响应域难以预测的问题。
为了达到上述目的,本发明所采用的技术方案为:
双起重机系统变幅角响应建模算法,其特征在于:包括以下步骤:
(1.1)、建立双台汽车起重机系统的几何模型及设定坐标系,并给出各点的 位置向量如下:
第一台汽车起重机系统中吊臂A1B1与第一台汽车起重机系统中转台的铰接 点A1在基坐标系{B}下的位置向量为:
第二台汽车起重机系统中吊臂A2B2与第二台汽车起重机系统中转台的铰接 点A2在基坐标系{B}下的位置向量为:
吊臂A1B1与第一台汽车起重机系统中吊绳B1C1的铰接点B1在基坐标系{B} 下的位置向量为:
吊臂A2B2与第二台汽车起重机系统中吊绳B2C2的铰接点B2在基坐标系{B} 下的位置向量为:
吊绳B1C1与负载C1C2的铰接点C1在动坐标系{P}下的位置向量为:
吊绳B2C2与负载C1C2的铰接点C2在动坐标系{P}下的位置向量为:
负载C1C2的重心Op在基坐标系{B}下的位置向量为:
其中,D和d分别为汽车起重机系统间距A1A2和负载C1C2的长度;基坐标系 {B}:O-YZ坐落于A1A2连接点的中心;动坐标系{P}:Op-YpZp坐落于C1C2连接 点的中心;L1和L2分别是吊臂A1B1和吊臂A2B2的长度;γ1和γ2分别是吊臂A1B1和 吊臂A2B2的变幅角;y和z分别是负载C1C2中心Op的沿Y轴和Z轴的笛卡尔坐标值;
吊绳B1C1与负载C1C2的铰接点C1在基坐标系{B}下的位置向量为:
吊绳B2C2与负载C1C2的铰接点C2在基坐标系{B}下的位置向量为:
(1.2)、根据步骤(1.1)各点的位置向量,建立系统约束方程,其中:
吊绳B1C1的约束方程为:
吊绳B2C2的约束方程为:
根据上述所有方程,整理可得:
其中θ表示动坐标系{P}相对于基坐标系{B}的旋转的角度。S1和S2分别为吊 绳B1C1和吊绳B2C2的长度;
上式可重新写成系统约束方程为:
K1isinγi+K2icosγi+Kai=0,i=1,2;
(1.3)、求解步骤(1.2)所得到的系统约束方程,得到双台汽车起重机系统 变幅运动下的变幅角响应方程:
针对上式系统约束方程求解,可得双台汽车起重机系统变幅运动下的变幅角 方程为:
其中:
根据变幅角的非负性,双台汽车起重机系统变幅运动下的变幅角响应方程可 重新写成为:
(1.4)、根据步骤(1.3)所得到的双台汽车起重机系统变幅运动下的变幅角 响应方程,进一步建立双台汽车起重机系统变幅运动下的变幅角响应模型:
根据上述分析,进一步构建双台汽车起重机系统变幅运动下的变幅角响应模 型:
M=Tγ
其中系统随机向量M={M1 M2}T,
系统随机矩阵T1=K31-K21,T2=K32-K22;
系统变幅角响应向量γ={γ1 γ2}T,
一种基于变幅角响应模型的双起重机系统变幅角随机响应域预测方法,其特 征在于:包括以下步骤:
(2.1)、在吊运负载过程中,复杂的环境因素及绳索的振动往往导致载荷参 数具有随机性,因此,建立随机参数模型如下:用随机向量 X={x1,x2,...,xr,...,xn}T表示双台汽车起重机系统变幅运动中所有随机负载参 数,其中n为随机参数的个数,双台汽车起重机系统随机参数包括:负载沿Y轴 的位置坐标y,负载沿Z轴的位置坐标z,动坐标系{P}相对于基坐标系{B}的旋转 的角度θ;
(2.2)、基于步骤(1.4)得到的双台汽车起重机系统变幅运动下的变幅角响 应模型,结合步骤(2.1)中的随机参数模型,建立带有随机参数模型的双台汽 车起重机系统变幅运动下的变幅角响应等效方程:
Mi(Ki(X))=Ti(Ki(X))γi,i=1,2,
其中,复合函数向量Mi(Ki(X))=Mi,复合函数矩阵Ti(Ki(X))-Ti。 Ki(X)={K1i(X),K2i(X),Kai(X)}T为引入的关系函数矩阵;
(2.3)、根据摄动随机复合函数法对步骤(2.2)中的带有随机参数模型的双 台汽车起重机系统变幅运动下的变幅角响应等效方程进行求解,得到带有随机载 荷参数的双台汽车起重机系统变幅运动下的变幅角表达式:
首先,根据一阶泰勒展开式和复合函数微分法,可以将复合函数向量 Mi(Ki(X))和复合函数矩阵Ti(Ki(X))分别展开表示为:
其中:
表示复合函数向量Mi(Ki(X))的期望,Δ1Mi表示复合函数向量 Mi(Ki(X))的增量,表示复合函数向量Ti(Ki(X))的期望,Δ1Ti表示复合函数 向量Ti(Ki(X))的增量,Xo表示随机向量X的期望,表示随机参数xr的期望;
因此,双台汽车起重机系统变幅运动下的随机变幅角响应方程可表示为:
用纽曼级数展开式对展开可表示为:
忽略高阶项,根据随机摄动法,将上式代入双台汽车起重机系统变幅运动下 的随机变幅角响应方程可得:
进一步整理可得:
(2.4)、根据随机变量函数矩法对步骤S3中的双台汽车起重机系统变幅运 动下的随机变幅角响应方程进行求解,可得随机载荷下变幅角响应域的期望的方 程为:
E(γi)=(Ti(Ki(Xe)))-1Mi(Ki(Xe))r i=1,2,
随机载荷下变幅角响应域的方差的方程为:
本发明提供的双台汽车起重机系统变幅角响应模型计算方法中,操作人员能 够根据确定性结构参数,结合负载运动轨迹曲线即确定性载荷参数,然后根据双 台汽车起重机系统变幅运动下的变幅角响应方程求得负载在任意位置下的变幅 角。针对随机载荷下的变幅角响应域数学特征的预测问题,本发明在双台汽车起 重机系统变幅角响应模型获得的基础上,基于摄动随机复合函数法和随机函数矩 法,根据双台汽车起重机系统变幅运动下随机载荷下的变幅角响应域的预测方 法,以获得变幅角响应域的数学特征公式。基于此,本发明还给出了双台汽车起 重机系统变幅运动下随机载荷下的变幅角响应域的预测方法在计算机中实施的 两套方案。本发明在充分考虑随机载荷参数的情况下设计出变幅角响应域的预测 方法,以提高双台汽车起重机系统随机响应域预测能力,具有快速性和精度高的 特点,以提高双台汽车起重机系统作业的可靠性。具体有益技术效果如下:
1)本发明提供的双台汽车起重机系统变幅运动下的变幅角响应模型的计算 方法,具体包括各结构参数、载荷参数和变幅角的运动关系,可以快速给出双台 汽车起重机系统变幅运动下的各输入端、输出端的运动状态,对大型起吊作业工 程具有重要的参考意义。
2)与传统的多起重机系统响应分析方法相比,本发明提供的双台汽车起重 机系统变幅运动下随机载荷下的变幅角响应域的预测方法充分考虑了载荷参数 的随机性,计算得到的变幅角响应域的数学特征对系统的可靠性设计及参数优化 具有重要的指导意义。
3)本发明提供的双台汽车起重机系统变幅运动下随机载荷下的变幅角响应 域的预测方法在计算机的实施方案一和实施方案二,可根据工程问题的复杂程度 合理选取合适的方案。具体而言,实施方案一中针对的是已知各类确定性结构参 数和详细的随机性载荷参数分布的场合,求解结果精度更高,但是计算时间较长; 而实施方案二中所提出的摄动随机复合函数法充分考虑变幅运动下的变幅角响 应方程的多层复合函数关系,并结合随机函数矩法推导出变幅角响应域的数学特 征,同时,针对载荷参数的小不确定性的场合,采用一阶泰勒级数展开式和一阶 纽曼级数展开式进行求解,既充分考虑了工程问题的复杂性,又保证了计算结果 的相对精确性,特别地,计算时间较短且计算效率显著提高。
附图说明
图1是双台汽车起重机系统三维模型示意图;图中示出了第一台汽车起重机 系统中转台1、第二台汽车起重机系统中转台2、第一台汽车起重机系统中吊臂 A1B1、第二台汽车起重机系统中吊臂A2B2、第一台汽车起重机系统中吊绳B1C1、 第二台汽车起重机系统中吊绳B2C2、负载C1C2、负载重心Op、铰接点A1、A2、B1、 B2、C1、C2及其位置关系。
图2是双台汽车起重机系统坐标结构示意图,图中示出了基坐标系{B}、动 坐标系{P}、吊臂A1B1、吊臂A2B2、吊绳B1C1、吊绳B2C2、负载C1C2、负载重心 Op、铰接点A1、A2、B1、B2、C1、C2及其几何关系。
图3是本发明提供的双台汽车起重机系统变幅运动下随机载荷下的变幅角 响应域的预测方法的流程图。
图4是本发明提供的随机载荷参数的方差系数范围为[0,0.1]时,采用双台汽 车起重机系统变幅运动下随机载荷下的变幅角响应域的预测方法在计算机中实 施方案一和实施方案二所计算的第一台汽车起重机系统变幅角响应域的方差曲 线图。
图5是本发明提供的随机载荷参数的方差系数范围为[0,0.1]时,采用双台汽 车起重机系统变幅运动下随机载荷下的变幅角响应域的预测方法在计算机中实 施方案一和实施方案二所计算的第一台汽车起重机系统变幅角响应域的方差曲 线图。
具体实施方式
参见图2,双台汽车起重机系统变幅运动下的变幅角响应模型的计算方法, 按如下步骤进行:
步骤一:建立几何模型及设定坐标系,并给出各点的位置向量。
吊臂A1B1与第一台汽车起重机系统的转台1的铰接点A1在基坐标系{B}下的 位置向量为:
吊臂A2B2与第二台汽车起重机系统的转台2的铰接点A2在基坐标系{B}下的 位置向量为:
吊臂A1B1与吊绳B1C1的铰接点B1在基坐标系{B}下的位置向量为:
吊臂A2B2与吊绳B2C2的铰接点B2在基坐标系{B}下的位置向量为:
吊绳B1C1与负载C1C2的铰接点C1在动坐标系{P}下的位置向量为:
吊绳B2C2与负载C1C2的铰接点C2在动坐标系{P}下的位置向量为:
负载C1C2的重心Op在基坐标系{B}下的位置向量为:
其中D和d分别为起重机间距A1A2和负载C1C2的长度。基坐标系{B}:O-YZ 坐落于A1A2连接点的中心。动坐标系{P}:Op-YpZp坐落于C1C2连接点的中心。L1和L2分别是吊臂A1B1和吊臂A2B2的长度。γ1和γ2分别是吊臂A1B1和吊臂A2B2的 变幅角。y和z分别是负载C1C2中心Op的沿Y轴和Z轴的笛卡尔坐标值。
吊绳B1C1与负载C1C2的铰接点C1在基坐标系{B}下的位置向量为:
吊绳B2C2与负载C1C2的铰接点C2在基坐标系{B}下的位置向量为:
步骤二:根据步骤一各点的位置向量,建立系统约束方程。
吊绳B1C1的约束方程为:
吊绳B2C2的约束方程为:
根据上述所有方程,整理可得:
其中θ表示动坐标系{P}相对于基坐标系{B}的旋转的角度。S1和S2分别为吊 绳B1C1和吊绳B2C2的长度。
上式可重新写成系统约束方程为:
K1isinγi+K2icosγi+Kat=0,i=1,2
步骤三:求解步骤二所得到的系统约束方程,得到双台汽车起重机系统变幅 运动下的变幅角响应方程。
针对上式系统约束方程求解,可得双台汽车起重机系统变幅运动下的变幅角 方程为:
其中,
根据变幅角的非负性,双台汽车起重机系统变幅运动下的变幅角响应方程可 重新写成为:
步骤四:根据步骤三所得到的双台汽车起重机系统变幅运动下的变幅角响应 方程,进一步建立双台汽车起重机系统变幅运动下的变幅角响应模型。
根据上述分析,进一步构建双台汽车起重机系统变幅运动下的变幅角响应模 型:
M=Tγ
其中系统随机向量M={M1 M2}T, 系统随机矩 阵T1=K31-K21,T2=K32-K22。系统变幅角响应向量 γ={γ1 γ2}T,
图1是双台汽车起重机系统三维模型示意图,包括第一台汽车起重机系统的 转台1、第二台汽车起重机系统的转台2、第一台汽车起重机系统的吊臂A1B1、 第二台汽车起重机系统的吊臂A2B2、第一台汽车起重机系统的吊绳B1C1、第二 台汽车起重机系统的吊绳B2C2、负载C1C2、负载重心Op、铰接点A1、A2、B1、B2、 C1、C2。变幅运动中,转台1(转台2)保持静止状态,即并不通过各自的回转 机构实现负载C1C2绕起重机回转中心轴线转动的运动;吊臂A1B1(吊臂A2B2) 保持静止状态,包括多节相互套接的伸缩臂,即伸缩臂并不通过伸缩驱动机构的 伸缩作用产生相对运动,即并不改变吊臂A1B1(吊臂A2B2)的长度以调整汽车 起重机的作业半径;吊绳B1C1(吊绳B2C2)保持静止状态,即并不通过变幅机构 中起升机构中吊绳B1C1(吊绳B2C2)的伸缩动作来实现负载C1C2在竖直平面内的 升降运动。变幅油缸D1E1(变幅油缸D2E2)一端与转台1(转台2)铰接,另一 端与吊臂A1B1(吊臂A2B2)铰接,通过调节变幅机构中变幅油缸D1E1(变幅油 缸D2E2)的长度,进一步实现吊臂A1B1(吊臂A2B2)在竖直平面内绕变幅油缸 D1E1(变幅油缸D2E2)与转台1(转台2)铰接点处做旋转运动以改变吊臂 A1B1(吊臂A2B2)仰角的变化,从而改变汽车起重机的变幅角度。对于上述双 台汽车起重机系统,以下对本发明提供的双台汽车起重机系统变幅运动下随机载 荷下的变幅角响应域的预测方法进行描述。
参见图3,该图是本发明提供的双台汽车起重机系统变幅运动下随机载荷下 的变幅角响应域的预测方法的流程图,按如下步骤进行:
S1:在吊运负载过程中,复杂的环境因素及绳索的振动往往导致载荷参数具 有随机性。因此,建立随机参数模型如下:用随机向量X={x1,x2,...,xr,...,xn}T表示双台汽车起重机系统变幅运动中所有随机负载参数,其中n为随机参数的个 数。双台汽车起重机系统随机参数包括:负载沿Y轴的位置坐标y,负载沿Z轴 的位置坐标z,动坐标系{P}相对于基坐标系{B}的旋转的角度θ。
S2:基于双台汽车起重机系统变幅运动下的变幅角响应模型的计算方法中步 骤四的双台汽车起重机系统变幅运动下的变幅角响应方程,结合步骤S1中的随 机参数模型,建立带有随机参数模型的双台汽车起重机系统变幅运动下的变幅角 响应等效方程:
Mi(Ki))=Ti(Ki(X))γi,i=1,2
其中,复合函数向量Mi(Ki(x))=Mi,复合函数矩阵Ti(Ki(X))=Ti。 Ki(X)={K1i(X),K2i(X),K3i(X)}T为引入的关系函数矩阵。
S3:根据摄动随机复合函数法(包括随机摄动法和复合函数微分法)对步骤 S2中的带有随机参数模型的双台汽车起重机系统变幅运动下的变幅角响应等效 方程进行求解,得到带有随机载荷参数的双台汽车起重机系统变幅运动下的变幅 角表达式。
首先,根据一阶泰勒展开式和复合函数微分法,可以将复合函数向量 Mi(Ki(X))和复合函数矩阵Ti(Ki(X))分别展开表示为:
其中:
其中,表示复合函数向量Mi(Ki(X))的期望,Δ1Mi表示复合函数向量 Mi(Ki(X))的增量,表示复合函数向量Ti(Ki(X))的期望,Δ1Ti表示复合函数 向量Ti(Ki(X))的增量,Xo表示随机向量X的期望,表示随机参数xr的期望。
因此,双台汽车起重机系统变幅运动下的随机变幅角响应方程可表示为:
用纽曼级数展开式对展开可表示为:
忽略高阶项,根据随机摄动法,将上式代入双台汽车起重机系统变幅运动下 的随机变幅角响应方程可得:
进一步整理可得:
S4:根据随机变量函数矩法对步骤S3中的双台汽车起重机系统变幅运动下 的随机变幅角响应方程进行求解,可得随机载荷下变幅角响应域的期望的方程 为:
E(γi)=(Ti(Ki(Xe)))-1Mi(Ki(Xe)),i=1,2
随机载荷下变幅角响应域的方差的方程为:
对于双台汽车起重机系统变幅运动下确定性响应模型的随机性响应域预测 方法在计算机中实施的步骤,进一步说明如下:
根据起重机设计参数和工况要求,确定各结构参数和载荷参数的确定值;
在上述各结构参数和载荷参数获得的前提下,利用MATLAB编程将确定性 结构参数和确定性载荷参数依次带入双台汽车起重机系统变幅运动下的变幅角 响应方程。
因此,得到确定性参数下的双台汽车起重机系统变幅运动下的变幅角响应。
确定性结构参数包括起重机间距A1A2的长度D,负载C1C2的长度d,吊臂A1B1的长度L1,吊臂A2B2的长度L2,吊绳B1C1的长度S1,吊绳B2C2的长度S2。
确定性载荷参数包括负载C1C2中心Op的位姿坐标,包括沿Y轴的笛卡尔坐标 值y、沿Z轴的笛卡尔坐标值z以及动坐标系{P}相对于基坐标系{B}的旋转的角度 θ。
对于双台汽车起重机系统变幅运动下随机载荷下的变幅角响应域的预测方 法在计算机中实施方案一的步骤,进一步说明如下:
根据起重机设计参数和工况要求,确定各结构参数的确定值和载荷参数随机 分布值的数学特征;
利用MATLAB编程将各结构参数的确定值和载荷参数随机分布值的数学特 征依次带入随机载荷下变幅角响应域的期望的公式和随机载荷下变幅角响应域 的方差的公式。
因此,得到随机载荷参数下的双台汽车起重机系统变幅运动下的变幅角响应 域的数学特征。
载荷参数随机分布值的数学特征包括载荷参数随机分布值的期望和方差。
变幅角响应域的数学特征包括变幅角响应域的期望和方差。
对于双台汽车起重机系统变幅运动下随机载荷下的变幅角响应域的预测方 法在计算机中实施方案二的步骤,进一步说明如下:
根据起重机设计参数和工况要求,确定各结构参数的确定值和载荷参数随机 分布值的数学特征;
在上述各结构参数的确定值和载荷参数随机分布的数学特征获得的前提下, 从每个载荷参数的随机分布值中任选一个随机值,并输入到MATLAB程序;
利用MATLAB编程将各结构参数的确定值和载荷参数的随机值依次带入双 台汽车起重机系统变幅运动下的变幅角响应方程。
因此,得到随机载荷参数下的双台汽车起重机系统变幅运动下的变幅角响 应。
重复上述过程至次数i=10000次,并输出随机载荷参数下的双台汽车起重机 系统变幅运动下的变幅角响应域分布曲线,并根据计算机指令输出随机载荷参数 下的双台汽车起重机系统变幅运动下的变幅角响应域的数学特征。
参见图4,该图是本发明提供的随机载荷参数的方差系数范围为[0,0.1]时, 采用双台汽车起重机系统变幅运动下随机载荷下的变幅角响应域的预测方法在 计算机中实施方案一和实施方案二所计算的第一台汽车起重机系统变幅角响应 域的方差曲线图。
引入随机参数模型来表示不同类型的随机载荷参数,并采用实施方案一和实 施方案二分别计算双台汽车起重机系统中第一台汽车起重机系统的变幅角响应 域的期望和方差,如表1~2所示。
表1第一台汽车起重机系统变幅角响应域的期望.
表2 第一台汽车起重机系统变幅角响应域的方差.
以第一台汽车起重机系统为研究对象,从图1~2和表1~2所示可知,在小范 围的随机载荷参数的方差系数下,双台汽车起重机系统变幅运动下随机载荷下的 变幅角响应域的预测方法在计算机中实施方案一和实施方案二计算的结果基本 保持一致;随着随机载荷参数的方差系数的增加,变幅角响应域的期望和方差的 误差增加。实施方案一的特点是:计算精度较高,计算效率高,计算耗时短,此 功能是一般计算软件不能实现的。实施方案二的特点是:计算效率低,耗时长。
参见图5,该图是本发明提供的随机载荷参数的方差系数范围为[0,0.1]时, 采用双台汽车起重机系统变幅运动下随机载荷下的变幅角响应域的预测方法在 计算机中实施方案一和实施方案二所计算的第二台汽车起重机系统变幅角响应 域的方差曲线图。
引入随机参数模型来表示不同类型的随机载荷参数,并采用实施方案一和实 施方案二分别计算双台汽车起重机系统中第二台汽车起重机系统的变幅角响应 域的期望和方差,如表3~4所示。
表3第二台汽车起重机系统变幅角响应域的期望.
表4 第二台汽车起重机系统变幅角响应域的方差.
以第二台汽车起重机系统为研究对象,从图3~4和表3~4所示可知,在小范 围的随机载荷参数的方差系数下,双台汽车起重机系统变幅运动下随机载荷下的 变幅角响应域的预测方法在计算机中实施方案一和实施方案二计算的结果基本 保持一致;随着随机载荷参数的方差系数的增加,变幅角响应域的期望和方差的 误差增加。实施方案一的特点是:计算精度较高,计算效率高,计算耗时短,此 功能是一般计算软件不能实现的。实施方案二的特点是:计算效率低,耗时长。
因此,本发明可以解决双台汽车起重机系统变幅运动中随机载荷下的变幅角 响应域数学特征的预测问题。两套实施方案在计算精度和计算效率方面各有特 点。
上述实施计算例仅仅为本发明的典型实施例,本发明不仅仅局限于上述实施例,凡在本发明的原理和内容之内所作的改动均应包含在本发明的保护范围之内。
Claims (2)
1.双起重机系统变幅角响应建模算法,其特征在于:包括以下步骤:
(1.1)、建立双台汽车起重机系统的几何模型及设定坐标系,并给出各点的位置向量如下:
第一台汽车起重机系统中吊臂A1B1与第一台汽车起重机系统中转台的铰接点A1在基坐标系{B}下的位置向量为:
第二台汽车起重机系统中吊臂A2B2与第二台汽车起重机系统中转台的铰接点A2在基坐标系{B}下的位置向量为:
吊臂A1B1与第一台汽车起重机系统中吊绳B1C1的铰接点B1在基坐标系{B}下的位置向量为:
吊臂A2B2与第二台汽车起重机系统中吊绳B2C2的铰接点B2在基坐标系{B}下的位置向量为:
吊绳B1C1与负载C1C2的铰接点C1在动坐标系{P}下的位置向量为:
吊绳B2C2与负载C1C2的铰接点C2在动坐标系{P}下的位置向量为:
负载C1C2的重心Op在基坐标系{B}下的位置向量为:
其中,D和d分别为汽车起重机系统间距A1A2和负载C1C2的长度;基坐标系{B}:O-YZ坐落于A1A2连接点的中心;动坐标系{P}:Op-YpZp坐落于C1C2连接点的中心;L1和L2分别是吊臂A1B1和吊臂A2B2的长度;γ1和γ2分别是吊臂A1B1和吊臂A2B2的变幅角;y和z分别是负载C1C2中心Op的沿Y轴和Z轴的笛卡尔坐标值;
吊绳B1C1与负载C1C2的铰接点C1在基坐标系{B}下的位置向量为:
吊绳B2C2与负载C1C2的铰接点C2在基坐标系{B}下的位置向量为:
(1.2)、根据步骤(1.1)各点的位置向量,建立系统约束方程,其中:
吊绳B1C1的约束方程为:
吊绳B2C2的约束方程为:
根据上述所有方程,整理可得:
其中θ表示动坐标系{P}相对于基坐标系{B}的旋转的角度,S1和S2分别为吊绳B1C1和吊绳B2C2的长度;
上式可重新写成系统约束方程为:
K1isinγi+K2icosγi+K3i=0,i=1,2
(1.3)、求解步骤(1.2)所得到的系统约束方程,得到双台汽车起重机系统变幅运动下的变幅角响应方程:
针对上式系统约束方程求解,可得双台汽车起重机系统变幅运动下的变幅角方程为:
其中:
根据变幅角的非负性,双台汽车起重机系统变幅运动下的变幅角响应方程可重新写成为:
(1.4)、根据步骤(1.3)所得到的双台汽车起重机系统变幅运动下的变幅角响应方程,进一步建立双台汽车起重机系统变幅运动下的变幅角响应模型;
根据上述分析,进一步构建双台汽车起重机系统变幅运动下的变幅角响应模型:
M=Tγ
其中系统随机向量:
M={M1 M2}T
其中系统随机矩阵:
T1=K31-K21
T2=K32-K22
其中系统变幅角响应向量:
γ={γ1 γ2}T
2.一种基于权利要求1建立的变幅角响应模型的双起重机系统变幅角随机响应域预测方法,其特征在于:包括以下步骤:
(2.1)、在吊运负载过程中,复杂的环境因素及绳索的振动往往导致载荷参数具有随机性,因此,建立随机参数模型如下:用随机向量X={x1,x2,...,xr,...,xn}T表示双台汽车起重机系统变幅运动中所有随机负载参数,其中n为随机参数的个数,双台汽车起重机系统随机参数包括:负载沿Y轴的位置坐标y,负载沿Z轴的位置坐标z,动坐标系{P}相对于基坐标系{B}的旋转的角度θ;
(2.2)、基于步骤(1.4)得到的双台汽车起重机系统变幅运动下的变幅角响应模型,结合步骤(2.1)中的随机参数模型,建立带有随机参数模型的双台汽车起重机系统变幅运动下的变幅角响应等效方程:
Mi(Ki(X))=Ti(Ki(X))γi i=1,2
其中,复合函数向量Mi(Ki(X))=Mi,复合函数矩阵Ti(Ki(X))=Ti,Ki(X)={K1i(X),K2i(X),K3i(X)}T为引入的关系函数矩阵;
(2.3)、根据摄动随机复合函数法对步骤(2.2)中的带有随机参数模型的双台汽车起重机系统变幅运动下的变幅角响应等效方程进行求解,得到带有随机载荷参数的双台汽车起重机系统变幅运动下的变幅角表达式:
首先,根据一阶泰勒展开式和复合函数微分法,可以将复合函数向量Mi(Ki(X))和复合函数矩阵Ti(Ki(X))分别展开表示为:
其中,
表示复合函数向量Mi(Ki(X))的期望,Δ1Mi表示复合函数向量Mi(Ki(X))的增量,表示复合函数向量Ti(Ki(X))的期望,Δ1Ti表示复合函数向量Ti(Ki(X))的增量,Xe表示随机向量X的期望,表示随机参数xr的期望;
因此,双台汽车起重机系统变幅运动下的随机变幅角响应方程可表示为:
用纽曼级数展开式对展开可表示为:
忽略高阶项,根据随机摄动法,将上式代入双台汽车起重机系统变幅运动下的随机变幅角响应方程可得:
进一步整理可得:
(2.4)、根据随机变量函数矩法对步骤(2.3)中的双台汽车起重机系统变幅运动下的随机变幅角响应方程进行求解,可得随机载荷下变幅角响应域的期望的方程为:
E(γi)=(Ti(Ki(Xe)))-1Mi(Ki(Xe))i=1,2
随机载荷下变幅角响应域的方差的方程为:
,
其中var(xr)指的是随机负载参数xr的方差值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710019654.2A CN106709264B (zh) | 2017-01-11 | 2017-01-11 | 双起重机系统变幅角响应建模算法及随机响应域预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710019654.2A CN106709264B (zh) | 2017-01-11 | 2017-01-11 | 双起重机系统变幅角响应建模算法及随机响应域预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106709264A CN106709264A (zh) | 2017-05-24 |
CN106709264B true CN106709264B (zh) | 2019-01-01 |
Family
ID=58907270
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710019654.2A Active CN106709264B (zh) | 2017-01-11 | 2017-01-11 | 双起重机系统变幅角响应建模算法及随机响应域预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106709264B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107741997A (zh) * | 2017-08-31 | 2018-02-27 | 合肥工业大学 | 适用于小区间参数下起重机系统变幅角响应域获取方法 |
CN107679277A (zh) * | 2017-08-31 | 2018-02-09 | 合肥工业大学 | 双起重机系统动力学建模算法及非奇异区间参数设计方法 |
CN108229045B (zh) * | 2018-01-16 | 2021-05-04 | 西南交通大学 | 一种基于灵敏度分析的高速受电弓关键参数辨识方法 |
CN111460591A (zh) * | 2020-03-31 | 2020-07-28 | 合肥工业大学 | 基于子区间理论的大不确定性起重机系统变幅角预测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105731279A (zh) * | 2016-03-23 | 2016-07-06 | 合肥工业大学 | 起重机参数类型的识别方法、响应误差控制方法和系统 |
-
2017
- 2017-01-11 CN CN201710019654.2A patent/CN106709264B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105731279A (zh) * | 2016-03-23 | 2016-07-06 | 合肥工业大学 | 起重机参数类型的识别方法、响应误差控制方法和系统 |
Non-Patent Citations (4)
Title |
---|
A modified hybrid uncertain analysis method for dynamic response field of the LSOAAC with random and interval parameters;Bin Zi等;《Journal of Sound and Vibration》;20160721;第374卷;第111-137页 |
Dynamics and trajectory tracking control of cooperative multiple mobile cranes;Sen Qian等;《Nonlinear Dynamics》;20150819;第83卷(第1-2期);第89-108页 |
Inverse kinematic control of a dual crane system experiencing base motion;Frank A. Leban等;《IEEE Transactions on Control Systems Technology》;20140425;第23卷(第1期);第331-339页 |
Kinematics and error analysis of cooperative cable parallel manipulators for multiple mobile cranes;Sen Qian等;《International Journal of Mechanics and Materials in Design》;20140318;第10卷(第4期);第395-409页 |
Also Published As
Publication number | Publication date |
---|---|
CN106709264A (zh) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106709264B (zh) | 双起重机系统变幅角响应建模算法及随机响应域预测方法 | |
Chang et al. | A fast path planning method for single and dual crane erections | |
Qian et al. | Dynamics and trajectory tracking control of cooperative multiple mobile cranes | |
CN109911773B (zh) | 一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法 | |
CN103466458A (zh) | 基于行为的多塔机三维空间防碰撞方法 | |
Zhang et al. | Research on trajectory planning and autodig of hydraulic excavator | |
CN110451408B (zh) | 一种基于能耗优化的塔机运动规划方法 | |
CN107741997A (zh) | 适用于小区间参数下起重机系统变幅角响应域获取方法 | |
CN110228754B (zh) | 一种自适应速度规划的工业吊车控制方法及系统 | |
Wang et al. | Collision-free path planning for mobile cranes based on ant colony algorithm | |
An et al. | Re-optimization strategy for truck crane lift-path planning | |
Pambudi et al. | Simulation design of trajectory planning robot manipulator | |
CN110342400B (zh) | 一种基于负载能量耦合的桥式起重机定位消摆控制方法 | |
Antic et al. | Anti-swing fuzzy controller applied in a 3D crane system | |
CN110407095A (zh) | 一种基于在线轨迹规划的桥式起重机定位消摆控制方法 | |
Hung et al. | Strategies to accelerate the computation of erection paths for construction cranes | |
FI20215750A1 (en) | DYNAMIC FLEXIBILITY COMPENSATION, COORDINATED LIFTING DEVICE CONTROL, AND ANTI-SWING IN LOAD HANDLING DEVICES | |
Li et al. | Study and implementation of a cooperative hoisting for two crawler cranes | |
Tubaileh | Working time optimal planning of construction site served by a single tower crane | |
CN107679277A (zh) | 双起重机系统动力学建模算法及非奇异区间参数设计方法 | |
CN102662331B (zh) | 一种基于虚拟现实的汽车吊臂挠度的仿真方法 | |
Wu et al. | Optimizing GPS-guidance transit route for cable crane collision avoidance using artificial immune algorithm | |
Huang et al. | Point cloud based path planning for tower crane lifting | |
Lin et al. | Lift path planning for tower cranes based on environmental point clouds | |
CN102530725B (zh) | 汽车起重机防摆控制技术 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |