CN102750261A - 一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法 - Google Patents

一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法 Download PDF

Info

Publication number
CN102750261A
CN102750261A CN2012101747855A CN201210174785A CN102750261A CN 102750261 A CN102750261 A CN 102750261A CN 2012101747855 A CN2012101747855 A CN 2012101747855A CN 201210174785 A CN201210174785 A CN 201210174785A CN 102750261 A CN102750261 A CN 102750261A
Authority
CN
China
Prior art keywords
matrix
node
singularity
driving
nodes
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.)
Pending
Application number
CN2012101747855A
Other languages
English (en)
Inventor
袁行飞
周练
李阿龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2012101747855A priority Critical patent/CN102750261A/zh
Publication of CN102750261A publication Critical patent/CN102750261A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Structure Of Telephone Exchanges (AREA)

Abstract

本发明提供了一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法。它包括以下步骤：1）确定铰接杆系机构的自由节点中的驱动节点和从动节点；2）针对铰接杆系在运动过程中产生的构型建立铰接杆系机构平衡矩阵，再由奇异值分解得到整体机构位移模态矩阵：3）将整体机构位移模态矩阵分解为关于驱动节点对应的机构位移模态矩阵 U a (2l×m)和从动节点对应的机构位移模态矩阵 U p (2k×m)；4）计算机构位移模态子矩阵的奇异值，如其最小奇异值趋于零，则可判定此铰接杆系机构发生奇异和运动分岔。本发明所提供的方法具有较强的操作性，为新型空间结构—机构提供了结构设计和施工设计的依据，具有广阔的应用前景。

Description

一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法
技术领域
本发明涉及新型空间结构—机构的奇异和运动分岔的确定方法，属于新型空间结构工程分析与设计技术领域。
背景技术
随着结构形式的不断发展,现代结构中出现了许多可动体系，如张拉整体结构、索穹顶、攀达穹顶以及可展结构等。这些新体系具有可变的初始形态，施工成型过程伴随着区别于弹性位移的几何位移，因此常与“机构”联系在一起。已有研究发现机构运动存在奇异和分岔现象，机构在奇异和分岔构型处可能出现运动“死点”，也可能按照非设计运动路径运动失去可控性。当机构按照非设计路径运动时，由于控制与实际运动之间的矛盾，将在系统中产生较大的附加荷载，导致整个结构系统失稳，造成严重后果。因此，对机构奇异和运动分岔的研究具有非常重要的实际意义。
发明内容
本发明所要解决的技术问题是提供一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法，具有较强的操作性。
为此，本发明采用以下技术方案： 
1）、确定铰接杆系机构的自由节点中的驱动节点和从动节点，所述驱动节点为与铰接杆系机构驱动杆件直接相连的自由节点，所述从动节点为除驱动节点外的其余自由节点；
2）、针对铰接杆系在运动过程中产生的构型建立铰接杆系机构平衡矩阵，所述铰接杆系机构平衡矩阵为2n×b阶矩阵，再由奇异值分解得到整体机构位移模态矩阵：
                       (1)
其中n为自由节点数，b为平面铰接机构杆件数，m为铰接杆系机构位移模态数，u为机构位移模态向量；
3）、将步骤2）获得的整体机构位移模态矩阵分解为关于驱动节点对应的机构位移模态矩阵U a (2l×m)和从动节点对应的机构位移模态矩阵U p (2k×m)，其中l为驱动节点数，k为从动节点数，并定义U a 的转置矩阵U a T和U p 的转置矩阵U p T分别为关于驱动节点的机构位移模态子矩阵和关于从动节点的机构位移模态子矩阵；
4）、计算机构位移模态子矩阵的奇异值，如其最小奇异值趋于零，则可判定此铰接杆系机构发生奇异和运动分岔。
进一步地，设驱动节点对应约束力矢量为F a ，为阶矩阵；从动节点对应约束力矢量为F p ，为
Figure 2012101747855100002DEST_PATH_IMAGE006
阶矩阵，由于约束力属于内力，不做功，可得：
Figure 2012101747855100002DEST_PATH_IMAGE008
                               (2)
当U a T非行满秩而U p T行满秩时，对应U a T的零空间，F p 等于零，即从动节点不受约束，铰接杆系机构发生输入奇异；
当U p T非行满秩而U a T行满秩时，对应U p T的零空间，F a 等于零，即驱动节点不受约束，铰接杆系机构发生输出奇异；
当U a T，U p T同时非行满秩时，铰接杆系机构发生结构奇异，驱动节点和从动节点不受约束，铰接杆系机构同时具有输入奇异和输出奇异的特点。
本发明通过考察铰接杆系机构的机构位移模态子矩阵是否行满秩来判定铰接杆系机构的奇异和运动分岔。
当驱动节点或从动节点对应的机构位移模态转置矩阵（机构位移模态子矩阵）出现非行满秩现象时，驱动节点或从动节点失去约束，获得自由度，机构运动发生奇异。本发明讨论的运动分岔为广义的分岔，即当机构在某个构型的运动状态发生奇异，则认为此构型为分岔构型。因此，可通过考察机构驱动节点和从动节点的机构位移模态子矩阵是否行满秩来判定机构奇异与运动分岔。
本发明所提供的方法具有较强的操作性，为新型空间结构—机构提供了结构设计和施工设计的依据，具有广阔的应用前景。
附图说明
图1为一平面铰接四杆体系的初始构型示意图。
图2为平面铰接四杆体系的顶升运动路径示意图。
图3为对应驱动节点的机构位移模态子矩阵最小奇异值变化曲线图。
图4为对应从动节点的机构位移模态子矩阵最小奇异值变化曲线图。
图5a为平面铰接四杆机构两种输入奇异分岔构型中第一种分岔构型示意图，
图5b为平面铰接四杆机构两种输入奇异分岔构型中第二种分岔构型示意图。
图6为输出奇异分岔构型示意图，图中，Eamin和Epmin分别为驱动节点和从动节点的机构位移模态子矩阵的最小奇异值，n为迭代步。
具体实施方式
以图1所示一平面铰接四杆体系为例，各点坐标为A (0, 0)，B (3, -4)，C (5, -20/3)，D (7, -4)，E (10, 0)， 杆件长度AB=5，BC =10/3，CD =10/3，DE =5，AE =10。
选择AB、DE杆件为驱动杆，则节点B、D为驱动节点，C节点为从动节点。在B、D节点上分别作用一竖直向上的单位力，通过MATLAB编程模拟其顶升运动过程，其运动路径如图2所示。
针对运动路径中的构型建立铰接杆系机构平衡矩阵A，所述铰接杆系机构平衡矩阵A为2n×b阶矩阵：
Figure 2012101747855100002DEST_PATH_IMAGE010
                 (3)
再由奇异值分解得到整体机构位移模态矩阵： 
Figure 755505DEST_PATH_IMAGE002
                        (4)
其中n为平面杆系自由节点数，b为杆件数，m为机构位移模态数，u为机构位移模态向量。
将整体机构位移模态矩阵分解为驱动节点对应的机构位移模态矩阵U a (2l×m)和从动节点对应的机构位移模态矩阵U p (2k×m)，其中l为驱动节点数，k为从动节点数，并定义U a 的转置矩阵U a T和U p 的转置矩阵U p T分别为关于驱动节点的机构位移模态子矩阵和关于从动节点的机构位移模态子矩阵；
计算机构位移模态子矩阵的特征值（或奇异值），如其最小特征值（或奇异值）趋于零，则可判定此铰接杆系机构发生奇异和运动分岔。
四杆体系顶升过程中，机构位移模态子矩阵的最小奇异值变化曲线如图3、4所示。
图3中，当迭代步n=50时，机构位移模态矩阵为：
， 
Figure 2012101747855100002DEST_PATH_IMAGE014
驱动节点的机构位移模态子矩阵U a T的最小奇异值Eami趋近于零，此时U a T非行满秩，U p T行满秩，机构发生输入奇异，奇异构型如图5(a)中实线所示。
图3中，当迭代步n=170时，机构位移模态子矩阵为:
Figure 2012101747855100002DEST_PATH_IMAGE016
，
Figure 2012101747855100002DEST_PATH_IMAGE018
驱动节点的机构位移模态子矩阵U a T的最小奇异值Eami趋于零。此时U a T非行满秩，U p T行满秩，机构发生输入奇异，对应奇异构型如图5(b)中实线所示。
图5(a)和5(b)中两奇异构型处，对应U a T的零空间F p 均为零，从动节点C均失去约束，获得自由度。当驱动节点B、D固定时，从动节点C仍具有如图5(a)和5(b)中虚线所示的可动性,因此此两奇异构型均为分岔构型。
图4中，当迭代步n=0时，机构位移模态子矩阵为：
Figure 2012101747855100002DEST_PATH_IMAGE020
，
从动节点的机构位移模态子矩阵最小奇异值Epmi趋于零，此时U a T行满秩，U p T非行满秩，机构发生输出奇异，奇异构型如图6中实线所示。对应U p T 的零空间F a 为零，驱动节点B，D失去约束。当从动节点C固定时，驱动节点B，D仍具有如图6中虚线所示的可动性。因此，此奇异构型也为分岔构型。

Claims (2)

1.一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法，其特征在于，它通过考察机构位移模态子矩阵是否非行满秩来判定铰接杆系机构的奇异和运动分岔，并包括以下步骤：
1）、确定铰接杆系机构的自由节点中的驱动节点和从动节点，所述驱动节点为与铰接杆系机构驱动杆件直接相连的自由节点，所述从动节点为除驱动节点外的其余自由节点；
2）、建立铰接杆系机构平衡矩阵，该矩阵为2n×b阶矩阵，对平衡矩阵进行奇异值分解得到整体机构位移模态矩阵： 
Figure 2012101747855100001DEST_PATH_IMAGE002
                          (1)
其中n为平面杆系自由节点数，b为杆件数，m为机构位移模态数，u为机构位移模态向量；
3）、将步骤2）获得的整体机构位移模态矩阵分解为驱动节点对应的机构位移模态矩阵U a (2l×m)和从动节点对应的机构位移模态矩阵U p (2k×m)，其中l为驱动节点数，k为从动节点数，并定义U a 的转置矩阵U a T和U p 的转置矩阵U p T分别为关于驱动节点的机构位移模态子矩阵和关于从动节点的机构位移模态子矩阵；
4）、计算机构位移模态子矩阵的奇异值，如其最小奇异值趋于零，则判定此铰接杆系机构发生奇异和运动分岔。
2.如权利要求1所述的一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法，其特征在于： 
设驱动节点对应约束力矢量为F a ，为
Figure 2012101747855100001DEST_PATH_IMAGE004
阶矩阵；从动节点对应约束力矢量为F p ，为
Figure 2012101747855100001DEST_PATH_IMAGE006
阶矩阵，由于约束力属于内力，不做功，可得：
Figure 2012101747855100001DEST_PATH_IMAGE008
                               (2)
当U a T非行满秩而U p T行满秩时，对应U a T的零空间，F p 等于零，即从动节点不受约束，铰接杆系机构发生输入奇异；
当U p T非行满秩而U a T行满秩时，对应U p T的零空间，F a 等于零，即驱动节点不受约束，铰接杆系机构发生输出奇异；
当U a T，U p T同时非行满秩时，铰接杆系机构发生结构奇异，驱动节点和从动节点不受约束，铰接杆系机构同时具有输入奇异和输出奇异的特点。
CN2012101747855A 2012-05-28 2012-05-28 一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法 Pending CN102750261A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012101747855A CN102750261A (zh) 2012-05-28 2012-05-28 一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012101747855A CN102750261A (zh) 2012-05-28 2012-05-28 一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法

Publications (1)

Publication Number Publication Date
CN102750261A true CN102750261A (zh) 2012-10-24

Family

ID=47030457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012101747855A Pending CN102750261A (zh) 2012-05-28 2012-05-28 一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法

Country Status (1)

Country Link
CN (1) CN102750261A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106021705A (zh) * 2016-05-17 2016-10-12 浙江大学 一种定量评价空间铰接杆系结构构件重要性的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1710723A2 (en) * 2005-04-08 2006-10-11 Dassault Systèmes Solver for a restrained deformable system with released degrees of freedom
CN101582095A (zh) * 2009-06-19 2009-11-18 东南大学 确定索杆系静力平衡状态的非线性动力有限元法
CN102393876A (zh) * 2011-07-13 2012-03-28 燕山大学 用杆组自由度计算并联机构自由度的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1710723A2 (en) * 2005-04-08 2006-10-11 Dassault Systèmes Solver for a restrained deformable system with released degrees of freedom
CN101582095A (zh) * 2009-06-19 2009-11-18 东南大学 确定索杆系静力平衡状态的非线性动力有限元法
CN102393876A (zh) * 2011-07-13 2012-03-28 燕山大学 用杆组自由度计算并联机构自由度的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐静: "铰接杆系机构的运动形态研究", 《中国优秀硕士学位论文全文数据库（工程科技Ⅱ辑）2007年》 *
马刚 等: "平面多自由度机构奇异位形分析新方法与仿真验证", 《机械科学与技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106021705A (zh) * 2016-05-17 2016-10-12 浙江大学 一种定量评价空间铰接杆系结构构件重要性的方法
CN106021705B (zh) * 2016-05-17 2019-12-06 浙江大学 一种定量评价空间铰接杆系结构构件重要性的方法

Similar Documents

Publication Publication Date Title
Kan et al. Nonlinear dynamic and deployment analysis of clustered tensegrity structures using a positional formulation FEM
CN110414165B (zh) 一种基于全局应力约束的多相材料柔顺机构拓扑优化方法
Ibrahimbegovic et al. Energy conserving/decaying implicit time-stepping scheme for nonlinear dynamics of three-dimensional beams undergoing finite rotations
CN104504284B (zh) 一种基于悬链线单元的松弛索网找形方法
Zhou et al. A valid inhomogeneous cell-based smoothed finite element model for the transient characteristics of functionally graded magneto-electro-elastic structures
Thai et al. Large deflection inelastic analysis of space trusses using generalized displacement control method
Luo et al. PD control for vibration attenuation in Hoop truss structure based on a novel piezoelectric bending actuator
Li et al. Earthquake input for finite element analysis of soil-structure interaction on rigid bedrock
CN108038324B (zh) 一种面向各向异性材料的连续体结构双材料拓扑优化方法
Zhu et al. Design of fully decoupled compliant mechanisms with multiple degrees of freedom using topology optimization
Yun et al. Self-learning simulation method for inverse nonlinear modeling of cyclic behavior of connections
CN107563005B (zh) 一种拉压不同刚度结构振动瞬时最优控制方法
Yoon et al. Nonlinear performance of continuum mechanics based beam elements focusing on large twisting behaviors
Chen et al. Topology optimization of bistable mechanisms with maximized differences between switching forces in forward and backward direction
Kim et al. Hysteretic mechanical–informational modeling of bolted steel frame connections
CN111368466A (zh) 一种基于频响函数参数修正的机械振动预测方法
Gao et al. An exact block-based reanalysis method for local modifications
CN112784489A (zh) 一种针对连续体结构的高效的动力学鲁棒性拓扑优化方法
Yang et al. Recent developments in geometrically nonlinear and postbuckling analysis of framed structures
Rong et al. Structural modal reanalysis for topological modifications with extended Kirsch method
Yang et al. Hybrid simulation of a zipper‐braced steel frame under earthquake excitation
CN103942381B (zh) 用于飞机铝合金结构性能预测的状态近场动力学方法
Ansola et al. Evolutionary optimization of compliant mechanisms subjected to non-uniform thermal effects
CN102750261A (zh) 一种采用机构位移模态子矩阵法检测铰接杆系机构奇异和运动分岔的方法
Vaezzadeh et al. Three-dimensional nonlinear dynamic analysis of slack cable structures using node Relaxation method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20121024