因为公式是变量的矩阵表示方式，这个公式可以通过总结所有的方程来得到，可以表示如下这里分别是质量阻尼和刚度矩阵，是负载向量。.稳定状态基础上的数值模拟曲柄的转速是,该灵活曲柄滑块机构的各项数值表示如下.，.,.,.,.。这里和分别是曲柄和耦合器的长度，是滑块的质量。通过曲柄和耦合器的个运动周期，可以看出稳态横向位移和中点弯曲应力的变化情况，以及分析本课题的结果。可以通过增加物理阻尼矩阵提高稳定性，被称作瑞利阻尼这里和是两个常数，可以从中对应于两个不同频率的振动的阻尼比得到。本文中和的值取决于自然频率。通过在运动方程中增加物理阻尼，也可以通过时间步骤观测超过个周期的运动，从而得到分析结果。当采用数值时间积分是出示条件从零开始。误差可以表示为这里和分别表示以有限元方法和参考方法为基础的两个值，总的来说，可以建立时间方程，而且很容易被接受，比如能量位移弯矩等等。和指的是时间积分的间隔，通常指的是稳态条件下的以个周期。因为没有个合适的准确的方法，在本文中可以通过个五次多项式表示个元件链接为基础的位移有限元方法得到参考值。图总能量的时间响应，耦合器的量纲中点挠度，耦合器在稳态条件下的中点应变。.数值模拟在这节中，我们讨论刚性曲柄机构。耦合器是唯的个灵活的连杆。在第六节中以以梁单元为基础，该梁单元可以做刚性轴运动，但是存在横向挠度。在第三节中讨论以有限元为基础的方法时，很有必要考虑模型的边界条件和形函数的相近程度，我们粗略的建立了耦合器应变线性分布方程，而且在弯矩不为零的条件下考虑耦合器的边界条件。在下面这个例子中，我们认为耦合器是由两个三个四个或者五个元件组成的，同时它的曲率分布可以表示为线性方程于是，时间响应和总能量误差，耦合器的中点挠度和应变都可以通过以应力为基础的有限元方法得到。同时，也评估了第自然频率。曲柄的转速为，该灵活的曲柄滑块机构中各个部件的值可以表示如下.，.,.,.,.。这里和分别是曲柄和耦合器的长度，是滑块的质量。为了通过以位移为基础的有限元方法比较误差，我们同样要用它建立个机构，结果可以参考文献。表两种有限元方法的第自然频率误差元件数目第自然频率以位移为基础的有限元方法以应力为基础的有限元方法.自由度数目表两种有限元方法的总能量误差元件数目第自然频率以位移为基础的有限元方法以应力为基础的有限元方法.自由度数目图显示了总能量的时间响应，耦合器的量纲中点挠度，耦合器在稳态条件下的中点应变。表分别比较了以位移为基础和以应力为基础的有限元方法的第自然频率误差总能量耦合器的中点挠度量纲以及耦合器的中点应变。误差可以由公式得到。结果表明，当两种方法中的元件数目相同时，以应力为基础的方法误差较以位移为基础的误差大。但是，当自由度的数目相同时，以应力为基础的有限元方法的误差比以位移为基础的有限元方法的误差小很多。同时，我们注意到当元件相同，除去第自然频率误差时，以应力为基础的有限元方法的误差也比以位移为基础的有限元方法的小很多。这说明以应力为基础的有限元方法可以提供大量精确的解决动态弹塑性问题的方法。表两种有限元方法的耦合器中点挠度误差元件数目第自然频率以位移为基础的有限元方法以应力为基础的有限元方法.自由度数目表两种有限元方法的耦合器中点应变误差元件数目第自然频率以位移为基础的有限元方法以应力为基础的有限元方法.自由度数目.结论本文提出了种新的以应力为基础的有限元方法来解决欧拉拉格朗日梁问题。该方法尤其适用于解决动态弹塑性问题。并且提出了梁的近似曲率。然后我们可以通过整合近似曲率得到横向挠度和应力分布。在整合过程中，有必要使梁单元的边界条件得到满足，从而可以得到整合常数。本文中，我们提出了在高速运作下解决灵活曲柄滑块机构问题。结果表明，在同样的自由度下，以应力为基础的有限元方法的误差小于常规方法的误差，常规方法也就是以位移为基础的有限元方法。同样，在元件数目相同的条件下，以应力为基础的有限元方法可以提供更多准确的解决方法。参考文献.

（其他） 回转体质心位置检测设备的设计论文.doc

（图纸） 机架.dwg

（图纸） 脚轮.dwg

（其他） 任务书.doc

（其他） 外文翻译--应力为基础的有限元方法应用于灵活的曲柄滑块机构.doc

（图纸） 摇杆.dwg

（图纸） 装配图.dwg