CN105260499A - 一种三维柔性梁系统力学的运动仿真方法及系统 - Google Patents

一种三维柔性梁系统力学的运动仿真方法及系统 Download PDF

Info

Publication number
CN105260499A
CN105260499A CN201510515412.3A CN201510515412A CN105260499A CN 105260499 A CN105260499 A CN 105260499A CN 201510515412 A CN201510515412 A CN 201510515412A CN 105260499 A CN105260499 A CN 105260499A
Authority
CN
China
Prior art keywords
flexible beam
beam system
dimensional flexible
equation
dimensional
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
CN201510515412.3A
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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201510515412.3A priority Critical patent/CN105260499A/zh
Publication of CN105260499A publication Critical patent/CN105260499A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出一种三维柔性梁系统力学的运动仿真方法及系统。首先设定建立动力学方程的条件,然后根据设定的条件建动力学方程,将每个物体的位置和梯度作为动力学方程待求解的未知数;最后求解动力学方程,获得三维柔性梁系统中每个物体在每个时刻的位置,根据每个物体在每个时刻的位置信息构建三维柔性梁系统动画图像并显示出来。本发明可以对三维柔性梁进行动力学仿真并方便、快捷地设计三维柔性梁系统。

Description

一种三维柔性梁系统力学的运动仿真方法及系统
技术领域
本发明属于力学与机械技术领域,具体涉及一种三维柔性梁系统力学的运动仿真方法。
背景技术
从20世纪70-80年代起,传统意义上的CAD/CAE/CAM技术开始进入实用阶段,它们主要关注产品零部件质量和性能,通过采用结构设计、工程分析和制造过程控制的软件或工具,以达到设计和制造高质量零部件的目的。具体地说,传统的CAD技术基于三维实体几何造型技术,支持产品零部件的详细结构设计和形态分析。传统的CAE技术主要指应用有限元软件,完成产品零部件的结构分析、热分析、振动特性等功能分析问题。传统的CAM技术旨在提高产品零部件的可制造性,提供对机床、机器人、铸造过程、冲压过程、锻造加工等方面更好的控制。在过去的几十年里,传统的CAD/CAE/CAM技术在主要的工业领域(汽车、航空、通用机械、机械电子等)得到了广泛的应用,并且取得了巨大的成效。以汽车工业来说,在1995-1999的五年里,零部件故障率降低了40%,与之相伴的,是产品开发和制造成本的相应降低。
机械系统力学的运动仿真技术作为一种虚拟样机技术与CAD/CAE/CAM技术,向系统的设计/分析/制造、以提高产品整体质量和性能并降低开发与制造成本,提高零部件的质量和性能产生很好的作用。
现阶段国内的产品主要还是依靠于实验以及经验,对于使用仿真方法进行设计还是比较少,当前,安全、稳定、节能成为机械制造业发展的主题,为了更好的制造,对于系统进行动力学仿真是非常重要的,这有助于对于机械产品进行设计以及修改,进而进行生产销售。国内对于动力学仿真方法急需自主发展。如何进行自有仿真方法的开发,并把它投入到工业界进行学习以及发展,是当前的一个研究重点,以及急需发展的一个重要环节。
发明内容
本发明的目的在于提出一种三维柔性梁系统的动力学仿真方法,可以对三维柔性梁系统的运动进行仿真。
为了解决上述技术问题,本发明提供一种三维柔性梁系统力学的运动仿真方法,包括以下步骤:
第一步,设定建立动力学方程的条件,包括
(1)设定三维柔性梁系统中各物体的材料参数(弹性模量,密度等)、尺寸参数(长度,宽度,高度)以及初始位置参数;
(2)选择三维柔性梁系统中各个物体间的约束方式建立三维的柔性梁系统模型,约束方式包括固定约束,移动约束等;
(3)给定三维柔性梁系统模型所受的外力;
(4)根据需要设定三维柔性梁系统动力学仿真的时间,以及计算的时间步长。
第二步,构建动力学方程,将每个物体的位置和梯度作为动力学方程待求解的未知数;
第三步,求解动力学方程,获得三维柔性梁系统中每个物体在每个时刻(即每个步长)的位置;根据每个物体在每个时刻的位置信息构建三维柔性梁系统动画图像并显示出来。这样让用户有一个清晰的认识,判断该设计是否满足自己的需要,用户可以根据需要修改设计。
进一步,体三维柔性梁系统动力学方程如式(1)所示:
M q ·· Φ q T Φ q 0 q ·· λ = Q γ - - - ( 1 )
式(1)中,Φq为Φ关于q的导数阵,q为系统的广义坐标,Φ=0为三维柔性梁系统的约束方程,λ和γ分别为拉格朗日乘子矩阵和加速度约束方程的右项;为单体三维柔性梁系统的动力学方程;M=diag(M1M2……MN)为总的质量阵;Q=diag(Q1Q2……QN)为总的外力阵。
进一步,使用龙格库塔法对动力学方程进行求解。
本发明还提供一种三维柔性梁系统力学的运动仿真系统,包含以下模块:
模型参数输入模块,用于输入三维柔性梁系统各物体的尺寸参数,材料参数以及各物体的初始位置参数;
约束方式选定模块,用于在预先存储的多种约束方式中选定各物体间需要的约束方式;
外力输入模块,用于输入模型所受的外力;
计算时间以及计算步长输入模块,用于输入三维柔性梁系统力学仿真的时间,以及计算的时间步长;
动力学计算模块,用于进行动力学分析,其包括动力学方程构建子模块和计算子模块;动力学方程构建子模块用于根据模型参数输入、连接方式选定、外力模块输入,建立动力学方程,将每个物体的位置和梯度作为待求解的未知数;计算子模块用于求解三维的柔性梁系统运动的力学方程,获得三维的柔性梁系统中每个物体的在每个时刻(即时间步长)的位置;
显示模块,用于根据每个时刻系统中每个物体的位置,显示出三维柔性梁系统的图像画面。
本发明与现有技术相比,其显著优点在于,本发明提供了输入环节,计算环节,以及结果显示功能,能全面地对三维柔性梁系统进行动力学模拟仿真;本发明的功能在国内属于先进水平,很好的填补国内空白。
附图说明
图1是三维梁单元广义坐标表示示意图。
图2是柔性单摆示意图。
图3是本发明仿真试验一中对柔性梁自由下落的动力学问题进行仿真的结果图,其中(a)为使用本发明方法获得的仿真结果图,(b)为使用现有技术获得仿真结果图。
图4是三维柔性梁系统自由下落示意图。
图5是本发明仿真试验二中末端点A和B的位置示意图,在梁的刚性比较大的情况,梁的弹性模量为E=2.0e10N/m2,其中,(a)为第一根梁末端A点的位置示意图,(b)为第二根梁末端B点的位置示意图。
图6是本发明仿真试验二中对三维柔性梁系统自由下落的动力学问题进行仿真的结果示意图,其中E=2.0e10N/m2。
图7是本发明仿真试验二中第一根梁末端A点的位置仿真结果图,其中,E=2.0e6N/m2。
图8是本发明仿真试验二中对三维柔性梁系统自由下落的动力学问题进行仿真的结果图,其中E=2.0e6N/m2。
图9是本发明仿真试验二中第一根梁末端A的位置仿真结果图,其中,E=2.0e5N/m2。
图10是本发明仿真试验二中对三维柔性梁系统自由下落的动力学问题进行仿真的结果图,其中E=2.0e5N/m2。
具体实施方式
一、基于绝对节点坐标法的三维梁单元
1.1梁单元的表达式
结合图1,在绝对节点坐标法中,三维柔性梁单元的广义坐标是由梁单元的两个节点的全局位置坐标和梯度组成的,梁单元和两个节点的表达式如式(1)和(2)所示:
e=[eAeB]T(1)
e j = r T r x T r y T r z T j T j = A , B - - - ( 2 )
式(1)和(2)中,e为梁单元的广义坐标,ej分别为两节点的广义坐标。
其中(2)式中的r为节点A和B的全局位置坐标,表示梁单元的节点位置坐标,rxryrz为节点A和B的全局位置坐标对x,y,z方向的梯度,A和B表示梁单元的A点和B点,这样ej为12×1的向量。
而梁单元中任意点P的位置坐标可以通过单元坐标以及式(3)所示形函数来表示:
r=Se(3)
式(3)中S为形函数,S的表达式为:
S=[S1I3S2I3S3I3…S7I3S8I3]
S1=1-3ξ2+2ξ3,S2=l(1-2ξ23),
S3=l(η-ξη),S4=l(ζ-ξζ),
S5=3ξ2-2ξ3,S6=l(-ξ23),
S7=lξη,S8=lξζ
其中,I3为3阶单位矩阵,ξ=x/l,η=y/l,ζ=z/l,l表示梁单元的长度,x,y,z表示P点在梁单元坐标中的位置。
1.2梁单元的质量阵
将式(3)对时间求导,可以得到梁单元中任一点的绝对速度公式,如式(4)所示,
r · = S e · - - - ( 4 )
式(4)中,表示梁单元中任一点的绝对速度,表示广义坐标的导数。利用(4)式可以求得梁单元的动能,动能如式(5)所示,
T e = 1 2 ∫ V ρ r · T r · d V = 1 2 e · T ∫ V ρS T S d V e · = 1 2 e · T m e e · - - - ( 5 )
式(5)中,Te为梁单元的动能,me=∫VρSTSdV为梁单元的质量阵,ρ为材料的质量密度,V为梁单元的体积。由(5)可以看出,质量矩阵是一个常量矩阵。
1.3梁单元的刚度阵
梁单元对x,y,z方向的位移梯度J如式(6)所示,
J = ∂ r ∂ x = S 1 x e S 1 y e S 1 z e S 2 x e S 2 y e S 2 z e S 3 x e S 3 y e S 3 z e - - - ( 6 )
式(6)中,x=[xyz]T S i x = ∂ S i / ∂ x , S i y = ∂ S i / ∂ y , S i z = ∂ S i / ∂ z , Si为单元形函数的第i行。
利用位移梯度,可以获得梁单元的格林-柯西应变张量εG如式(7)所示,
ϵ G = 1 2 ( J T J - I ) = 1 2 ( e T a 1 e - 1 ) e T a 4 e e T a 5 e e T a 4 e ( e T a 2 e - 1 ) e T a 6 e e T a 5 e e T a 6 e ( e T a 3 e - 1 ) - - - ( 7 )
式(7)中,
a 1 = S 1 x T S 1 x + S 2 x T S 2 x + S 3 x T S 3 x , a 2 = S 1 y T S 1 y + S 2 y T S 2 y + S 3 y T S 3 y , a 3 = S 1 z T S 1 z + S 2 z T S 2 z + S 3 z T S 3 z , a 4 = S 1 x T S 1 y + S 2 x T S 2 y + S 3 x T S 3 y , a 5 = S 1 x T S 1 z + S 2 x T S 2 z + S 3 x T S 3 z , a 6 = S 1 y T S 1 z + S 2 y T S 2 z + S 3 y T S 3 z
其中,a1,a2,..,a6均为24×24的矩阵。
应变张量的列矢量ε的形式如式(8)所示,
ε=[εG11εG22εG33G12G13G23]T(8)
式(8)中,εG11,...,εG23表示式子(7)中εG的第1行第1列,...,第2行第三列等。
根据材料的本构关系,单元应力σ如式(9)所示,
σ=Dε(9)
式(9)中,D为弹性模量矩阵,其表达式为:
D = λ + 2 μ λ λ + 2 μ s y m m e t r i c λ λ λ + 2 μ 0 0 0 μ 0 0 0 0 μ 0 0 0 0 0 μ
其中,γ为梁的泊松比,其表达式如式(10)所示,
λ = γ E ( 1 + γ ) ( 1 - 2 γ ) , μ = E 2 ( 1 + γ ) - - - ( 10 )
式(10)中,E为梁的弹性模量。
梁单元的弹性力所作的虚功δUe的表达式如式(11)所示:
δU e = δ ( - 1 2 ∫ V σ T ϵ d V ) = - 1 2 δ ∫ V ϵ T D ϵ d V = - ∫ V ( D ϵ ) T δ ϵ d V = ( Q ϵ e ) T δ e - - - ( 11 )
其中,为梁的广义弹性力,其表达式如式(12)所示:
Q ϵ e = - ∫ V a 1 e a 2 e a 3 e 2 a 4 e 2 a 5 e 2 a 6 e ( D ϵ ) d V - - - ( 12 )
1.4单元的广义外力阵
设FP是作用于梁单元上P点的外力,则外力所做的虚功如式(13)所示:
δP F e = δr T F P = δe T Q F e - - - ( 13 )
其中,为梁单元的集中力的广义阵。
重力的虚功如式(14)所示:
δP g e = δr T ∫ V ρ g d V = δe T Q g e - - - ( 14 )
其中为梁单元的重力的广义力阵,g为重力加速度矢量,令总的广义外力 Q e = Q F e + Q g e + Q ϵ e .
2梁系统的动力学方程
设Be为单元节点坐标和总体节点坐标的转换阵,则有qe=Beqi,其中,qi为梁Bi总体节点坐标列,qe为单元的节点坐标列阵,则Bi的质量阵Mi和外力阵Qi如式(15)所示:
M i = Σ j = 1 n B i T m j e B i Q i = Σ j = 1 n B i T Q j e - - - ( 15 )
设柔性梁系统由梁Bi(i=1,…,N)组成,则单体系统的动力学方程如式(16)所示,
M q ·· = Q - - - ( 16 )
式(16)中:
M=diag(M1M2……MN)为总的质量阵;
Q=diag(Q1Q2……QN)为总的外力阵;
q为系统的广义坐标。
设Φ=0为系统的约束方程,则多体系统的动力学方程如式(17)所示:
M q ·· Φ q T Φ q 0 q ·· λ = Q γ - - - ( 17 )
式(17)中,Φq为Φ关于q的导数阵;λ和γ分别为拉格朗日乘子矩阵和加速度约束方程的右项。
3动力学求解与图像表示
式(17)是一组微分方程,使用龙格库塔法对式(17)进行求解,这样就可以计算得到系统中各个物体的每个时刻的位置q。
通过所得到的位置q的结果,对于某一时刻,画出图形。根据梁节点的位置,求得梁中每个点的位置,得到每个梁的图形。从而得到每个时刻系统的图形。
为了说明本发明,做了一下仿真实验。
仿真实验一
对有代表意义的柔性单摆的自由下落问题进行了计算,如图2,柔性单摆的参数为:L=1m,宽度和高度b=h=0.01m,泊松比v=0.3。弹性模量为E=2.0e6N/m2。
本仿真实验把仿真的结果画在x-y平面上并与现有技术]进行了比较,如图3。两者的结果一致,这样验证了本发明方法的有效性。
仿真实验二
对三维柔性梁系统的自由下落的例子进行了计算,如图4,其中,第一个梁与地面,第一个梁与第二个梁采用移动约束,两个梁刚开始成直角。其中梁的参数为L=1m,宽度和高度b=h=0.02m,泊松比v=0.3。
通过改变梁的弹性模量E可以计算不同结果,仿真中E分别取为2.0e10N/m2,2.0e6N/m2直到2.0e5N/m2。计算结果如图5至图10中所示,图中ADAMS结果取自多体动力学软件MSC.ADAMS与ANSYS联合使用得到的计算结果。从图5中可以发现,在刚性较大时,两者结果一致。
而当梁的柔性比较大时,弹性模量达到2.0e6N/m2时(图7),ADAMS计算结果出现异常。随着弹性模量的继续减小,当弹性模量达到2.0e5N/m2时(图9),ADAMS计算失败。这个结果提醒我们这样一个事实,在使用ADAMS软件进行动力学计算时,一定要谨慎看待所得的结果,或许这个结果有较大误差,甚至是错误的结论。
图8与图10可以明显的表明本发明方法可以很好的计算三维柔性梁系统的动力学问题。

Claims (4)

1.一种三维柔性梁系统的动力学仿真方法,其特征在于,包括以下步骤:
第一步,设定建立动力学方程的条件,包括:
设定三维的柔性梁系统中各物体的材料参数、尺寸参数以及初始位置参数、选择三维柔性梁系统中各个物体间的约束方式;
给定三维柔性梁系统模型所受的外力;
根据需要设定三维的柔性梁系统动力学分析的时间,以及计算的时间步长;
第二步,构建动力学方程,将每个物体的位置和梯度作为动力学方程待求解的未知数;
第三步,求解动力学方程,获得三维柔性梁系统中每个物体在每个时刻的位置;根据每个物体在每个时刻的位置信息构建三维柔性梁系统动画图像并显示出来。
2.如权利要求1所述三维柔性梁系统力学的运动仿真方法,其特征在于,多体三维柔性梁系统动力学方程如式(1)所示:
M q ·· Φ q T Φ q 0 q ·· λ = Q γ - - - ( 1 )
式(1)中,Φq为Φ关于q的导数阵,q为系统的广义坐标,Φ=0为三维柔性梁系统的约束方程,λ和γ分别为拉格朗日乘子矩阵和加速度约束方程的右项;为单体三维柔性梁系统的动力学方程;M=diag(M1M2……MN)为总的质量阵;Q=diag(Q1Q2……QN)为总的外力阵。
3.如权利要求2所述三维柔性梁系统力学的运动仿真方法,其特征在于,使用龙格库塔法对动力学方程进行求解。
4.一种三维柔性梁系统力学的运动仿真系统,其特征在于,包含以下模块:
模型参数输入模块,用于输入三维柔性梁系统各物体的尺寸参数,材料参数以及各物体的初始位置参数;
约束方式选定模块,用于在预先存储的多种约束方式中选定各物体间需要的约束方式;
外力输入模块,用于输入模型所受的外力值;
计算时间以及计算步长输入模块;
动力学计算模块,用于进行动力学分析,其包括动力学方程构建子模块和计算子模块;动力学方程构建子模块用于构建三维柔性梁系统运动的力学方程,将每个物体的位置和梯度作为待求解的未知数;计算子模块用于求解三维柔性梁系统运动的动力学方程,获得三维柔性梁系统运动中每个物体的在每个时刻的位置;
显示模块,用于根据每个时刻系统中每个物体的位置,显示出三维柔性梁系统的图像画面。
CN201510515412.3A 2015-08-20 2015-08-20 一种三维柔性梁系统力学的运动仿真方法及系统 Pending CN105260499A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510515412.3A CN105260499A (zh) 2015-08-20 2015-08-20 一种三维柔性梁系统力学的运动仿真方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510515412.3A CN105260499A (zh) 2015-08-20 2015-08-20 一种三维柔性梁系统力学的运动仿真方法及系统

Publications (1)

Publication Number Publication Date
CN105260499A true CN105260499A (zh) 2016-01-20

Family

ID=55100188

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510515412.3A Pending CN105260499A (zh) 2015-08-20 2015-08-20 一种三维柔性梁系统力学的运动仿真方法及系统

Country Status (1)

Country Link
CN (1) CN105260499A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107122515A (zh) * 2017-03-17 2017-09-01 北京航空航天大学 基于绝对节点坐标法的绳系运输系统的动力学分析方法
CN107220421A (zh) * 2017-05-18 2017-09-29 北京理工大学 一种空间复杂柔性结构多体系统动力学建模与计算方法
CN108182330A (zh) * 2018-01-15 2018-06-19 南京理工大学 一种基于b样条计算柔性矩形薄板刚柔耦合动力学响应的方法
CN108427322A (zh) * 2018-03-16 2018-08-21 南京航空航天大学 一种大柔性飞行器基于在线辨识的建模方法
CN109514556A (zh) * 2018-12-10 2019-03-26 南京理工大学 柔性仿人机械手手指摩擦碰撞瞬态响应的计算方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090228256A1 (en) * 2008-03-10 2009-09-10 Min Gyu Choi Method and system of graphical simulation of thin shell objects
US8222581B2 (en) * 2005-04-08 2012-07-17 Dassault Systemes Solver for a restrained deformable system with released degrees of freedom
CN103558769A (zh) * 2013-09-23 2014-02-05 广东工业大学 带末端质量体的柔性梁系统动力学建模方法及其控制方法
CN104091003A (zh) * 2014-06-30 2014-10-08 南京航空航天大学 一种基础运动时柔性壳结构大变形响应的有限元建模方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8222581B2 (en) * 2005-04-08 2012-07-17 Dassault Systemes Solver for a restrained deformable system with released degrees of freedom
US20090228256A1 (en) * 2008-03-10 2009-09-10 Min Gyu Choi Method and system of graphical simulation of thin shell objects
CN103558769A (zh) * 2013-09-23 2014-02-05 广东工业大学 带末端质量体的柔性梁系统动力学建模方法及其控制方法
CN104091003A (zh) * 2014-06-30 2014-10-08 南京航空航天大学 一种基础运动时柔性壳结构大变形响应的有限元建模方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
JINYANG LIU等: "Nonlinear formulation for flexible multibody system with large deformation", 《ACTA MECHANICA SINICA》 *
REFAAT Y. YAKOUB等: "Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements:Implementation and Applications", 《JOURNAL OF MECHANICAL DESIGN》 *
张斌: "基于绝对坐标法的柔性旋转梁的动力学研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
李永亮: "基于绝对节点坐标法的柔性梁动力学建模与分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
蒋依宁: "大变形柔性多体系统建模理论与实验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
郑彤等: "航空发动机叶片刚柔耦合动力学分析", 《机械工程学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107122515A (zh) * 2017-03-17 2017-09-01 北京航空航天大学 基于绝对节点坐标法的绳系运输系统的动力学分析方法
CN107122515B (zh) * 2017-03-17 2020-06-23 北京航空航天大学 基于绝对节点坐标法的绳系运输系统的动力学分析方法
CN107220421A (zh) * 2017-05-18 2017-09-29 北京理工大学 一种空间复杂柔性结构多体系统动力学建模与计算方法
CN108182330A (zh) * 2018-01-15 2018-06-19 南京理工大学 一种基于b样条计算柔性矩形薄板刚柔耦合动力学响应的方法
CN108427322A (zh) * 2018-03-16 2018-08-21 南京航空航天大学 一种大柔性飞行器基于在线辨识的建模方法
CN108427322B (zh) * 2018-03-16 2020-04-17 南京航空航天大学 一种大柔性飞行器基于在线辨识的建模方法
CN109514556A (zh) * 2018-12-10 2019-03-26 南京理工大学 柔性仿人机械手手指摩擦碰撞瞬态响应的计算方法
CN109514556B (zh) * 2018-12-10 2022-04-08 南京理工大学 柔性仿人机械手手指摩擦碰撞瞬态响应的计算方法

Similar Documents

Publication Publication Date Title
CN105260499A (zh) 一种三维柔性梁系统力学的运动仿真方法及系统
Olshevskiy et al. Three-dimensional solid brick element using slopes in the absolute nodal coordinate formulation
CN109902404B (zh) 不同阻尼形式的结构时程响应积分的统一递推计算方法
Zhang et al. Extremum response surface method of reliability analysis on two-link flexible robot manipulator
Pappalardo et al. A comparative study of the principal methods for the analytical formulation and the numerical solution of the equations of motion of rigid multibody systems
Lum et al. Structural optimization for flexure-based parallel mechanisms–Towards achieving optimal dynamic and stiffness properties
Aginaga et al. Improving static stiffness of the 6-RUS parallel manipulator using inverse singularities
US20090012751A1 (en) Method and program for structure analysis by finite element method
CN109634111B (zh) 一种高速重载机器人动态变形计算方法
Liu et al. Simple formulations of imposing moments and evaluating joint reaction forces for rigid-flexible multibody systems
Bascetta et al. Closed form Newton–Euler dynamic model of flexible manipulators
Moberg et al. Modeling and parameter estimation of robot manipulators using extended flexible joint models
CN104281730A (zh) 一种大转动变形的板壳结构动响应的有限元分析方法
Angel et al. Adams/matlab co-simulation: dynamic systems analysis and control tool
Stravopodis et al. Rectilinear tasks optimization of a modular serial metamorphic manipulator
CN117171864B (zh) 一种梁结构线性振动预测方法
Scaglioni et al. Closed-form control oriented model of highly flexible manipulators
Hurskainen et al. A planar beam finite-element formulation with individually interpolated shear deformation
Guo et al. A transformation method to generate the workspace of an n (3rrs) serial–parallel manipulator
Angel et al. Modeling, simulation and control of a differential steering type mobile robot
CN103217924B (zh) 面向实时控制应用的过约束重型并联机床动力学建模方法
Mathew et al. Dynamics of suspended cable driven parallel robots using the geometric variable strain approach
Balasubramanian Plate analysis with different geometries and arbitrary boundary conditions
Malik et al. Static and dynamic analysis of absolute nodal coordinate formulation planar elements
Nandihal et al. Dynamics of Rigid-Flexible Robots and Multibody Systems

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160120