计算力学第2课.pptVIP

计算力学 Computational Mechanics (2) 哈尔滨工程大学航天与建筑工程学院 王滨生 E-mail: wang_bs@126.com，heu.tulixue@163.com 基本内容 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 有限元法引论 总结步骤 利用变分法直接进行近似计算的几种方法： 里兹法：该方法是假设一位移函数v(x)，只令其先满足位移边界条件，然后再通过dP=0(最小势能原理)去近似满足力边界条件和平衡方程(控制方程)式。 伽辽金法：该方法是在里兹法的基础上发展起来的，其特点是在设位移函数时除了使v(x)满足位移边界条件外，还要求它满足力边界条件，然后通过dP=0使其近似满足控制方程式。 例：基于最小势能原理的变分法(以直梁为例) 求单元位移函数ve(x) 现设单元位移函数 (插值函数)为 由边界条件可求出ai。 边界条件为： 解得 所以 Ni叫做形状函数。在梁单元中，它表示一个两端固定梁只产生一个单位位移时梁弯曲成的形状。 求单元总势能及其变分 梁的弯曲变形能 对第一项进行积分 对第二项进行积分 小结 《小论文大作业》 数值计算或某一数值计算方法综述方面； 两种以上算法的耦合计算方面：如有限元与边界元耦合； 随机数值计算方面：如随机有限元、随机边界元； 颗粒物质(散体材料)力学性能的研究现状及方法方面； 数值计算应用领域方面； 其他数值计算方面，如新方法、新应用领域。 公共信箱：heu.tulixue@163.com 密码：045182519201 1.有限元的发展 50年代有限元首次用于飞机设计中，1956年M.J.Turner和R.W.Clough等人用矩阵对飞机结构进行了受力和变形分析。 1960年R.W.Clough首次提出“有限元”这个名词，有限元发作为一种数值分析方法正式出现于工程技术领域。 Clough（Berkeley土木系）1952，Clough被西雅图波音飞机公司聘请为他们summer faculty项目的成员。在结构动力学小组，小组的领导人是Turner。Turner让他从事一种delta翼结构的振动分析工作。由于它是三角平面形状，这个问题不能用基于标准梁理论的方法解决；Turner建议通过组装三角或矩形形状的平面应力小块来建立机翼的刚度特性列式，……。波音小组称这种方法为直接刚度法。 【插曲】航天与建筑工程学院：林同炎、兄林同骅、弟林同骥。兄林同骅：交通大学唐山学院土木工程学士（1933）、美国麻省理工航空工程硕士（1936）、美国米歇根大学工程力学博士（1953）。1944年8月18日，以林同骅为总设计师完成我国自行设计制造试飞成功的首架运输机，被定名为“C1010运输机”。 六、七十年代中国在设计水坝中(1964，冯康[U1]?) 创立了数值求解偏微分方程的有限元方法，独立的发展了有特色的数学理论基础。 陈省身、华罗庚、冯康。 法国科学院院士J. I. Lions教授在70年代末就已对此给出了很高的评价说：“有限元方法意义重大，中国学者在对外隔绝的环境下，独立创造了有限元方法” 上角标：工况 下角标：第一个下标表示节点编号，第二个下标表示单元编号 刚度 yx：沿x方向产生单位位移需要在y方向施加的力 这样把求解桁架节点位移的问题转变成求解代数方程组，已知{P}和[K]，求得{u}以后，桁架各杆的伸长、应变、内力和应力等都很容易得到解决。 可理解为 非变分的←→真实的运动或位移。 变分原理←→满足边界条件的可能的运动或位移。 微分形式←→某一瞬时。 积分形式←→某一有限的时间间隔或过程 驻值：极值、拐点 举例：稳定平衡、不稳定平衡 有势力又称保守力。如果作用在物体上的力所做之功仅与力作用点的起始位置和终了位置有关，而与其作用点经过的路径无关，这种力称为有势力或保守力。重力、弹簧力、静电库仑力、万有引力等都具有这一特征，因而都是有势力。 自由度 自由度 外力对系统做正功=系统势能减少，如重物下落。