基于-Matlab-手动编程了解有限元的基本原理.docVIP

基于 Matlab 手动编程了解有限元的基本原理 卢鸣飞 复旦大学力学与工程科学系 1 背景阐述 在有限元方法飞速发展的今天，各类有限元软件也越来越普及，功能越来越强大。并且，针对各个不同行业 的专业有限元软件也犹如雨后春笋般地涌现出来，大大降低了有限元软件的使用门槛。但其背后的理论基础 却是大同小异。 本文试图通过梳理有限元最基本的原理，让大家了解如何通过数值计算来求解一个简单的实际问题。通过原 理来了解数值计算本质，能够更有利于大家理解、掌握有限元技术。 2 问题描述 一电线塔如图 1 所示，所有杆件间的连接均为铰接。已知弹性模量 2.0x1011 Pa，密度为 7.8x103 kg/m3。杆 件的截面积 A=2.827433 x10-3 m2。求解：系统的前十阶固有频率。 F (t) 1 1 1 F (t) 1 1 1 1111输出点 1 1 1 1 0 0.08 0.1 1图 1 电线塔模型 1 3 理论分析 由于是杆件系统，所以只需要用到杆件系统的质量矩阵与刚度矩阵即可： ?c2 EA ? cs ?c2 s2 ?cs ?cs? ?s2 ? ?2 0 1 0? ? Al ?0 2 0 1? [k (e) ] ? ? ? [m(e) ] ? ? ? l ? 对称 c2 cs ? 6 ?1 0 2 0? ? s2 ? ?0 1 0 2?