二次型问题

2020年09月25日 星期五

在应用变分原理后有限元问题最后会转化对泛函的求导。本文将讲解对于泛函中被积函数具有二次项形式时所得到方程的形式。……

二次型问题

​ 在应用变分原理后有限元问题最后会转化对泛函的求导。本文将讲解对于泛函中被积函数具有二次项形式时所得到方程的形式。

1 基本公式

为了后续推导,需要有下述数学公式。(相关公式可以网上搜索得到)

二次型关于向量求导:

二次型

标量和对向量求导的公式: 向量积对标量求导公式:

2 单元刚度矩阵推导

采用变分原理时,如果泛函积分项中被积函数为关于整体位移 及其导数的二次函数,即

求导可应用标量和对向量求导的公式: 故这里首先推导上述二次型中第一项的导数。

1)整体坐标系下

上述可以表示为 在整体坐标系下,的关系式为: 故进一步可表示为: 这里为节点转换矩阵。

​ 此时可表示为: 相对于的导数按照二次项导数公式为: 相应的可记上述积分项为整体坐标系下单元刚度矩阵为: 整体刚度矩阵: 可以看出,若是对称的,则整体刚度也应是对称的。

2)局部坐标系

​ 一般积分都是在局部坐标系进行的,此时上述均应是在局部坐标系下,后续推导过程一致,不过的关系式应为: 对于可参看相关对于坐标转换的讲解。

​ 此时相应可得到相对于的导数按照二次项导数公式为: 记局部坐标系下单元刚度矩阵为: 相应的可以看到此时在集成整体坐标系过程中,需要做两步:

①坐标转换,②节点转换

3)导数部分

对于二次型的一次导数部分,利用向量积对标量求导公式: 有: 相应的可知局部坐标系下单元刚度矩阵为: 对于n阶导数类似可以总结为: 4)节点荷载

一次项类似的可表示为在局部坐标系下: 相应可以得到单元等效节点荷载为:

整体等效节点荷载: 综上,得到有限元方程: 对于原本荷载就施加在节点上的情形,其在泛函中可表达为: 其对于导数为

3 平面杆系结构应用

​ 对于平面杆系结构而言,可建立平面直角坐标系。记局部坐标系下杆件沿轴向位移为,沿法向为,势能在可表示为: 注意到实际上上述势能项是局部坐标和整体坐标的混搭,而积分项是在基于杆件的局部坐标系下推导得到的。

若记: 则第一项内力势能被积函数可以表示成二次项为: 其中: 形函数取: 则局部坐标系下单元刚度矩阵为:


更新于2020年9月25日 星期五 天气晴

精选博客

关于日本钢结构相关信息的摘取

应公司领导要求,需要给出日本有关钢结构建筑、桥梁及钢结构用量的统计数据。参考国内的情形,在日本国土交通省找到了相关数据。而且相对于国内的,数据更加全面。……

继 续 阅 读

板件宽厚比和加劲肋

这里主要想讲的是钢标中宽厚比取值的来源及其应用,以及局部稳定加劲肋的原理。对于四边简支单向受压板件,类似于轴压杆欧拉荷载,有一个临界应力,其表达式可以通过板的弹性平衡微分方程求解得到,这里直接给出(出处可见于各类钢结构书籍讲到板的局部稳定部分,如《稳定》8.3节或是《基理》5.6节)……

继 续 阅 读

五帝本纪

黄帝者,少典之子,姓公孙,名曰轩辕。生而神灵,弱而能言,幼而徇齐,长而敦敏,成而聪明。轩辕之时,神农氏世衰。诸侯相侵伐,暴虐百姓,而神农氏弗能征。於是轩辕乃习用干戈,以征不享,诸侯咸来宾从。而蚩尤最为暴,莫能伐。炎帝欲侵陵诸侯,诸侯咸归轩辕。轩辕乃修德振兵,治五气,艺五种,抚万民,度四方,教熊罴貔貅貙虎,以与炎帝战于阪泉之野。三战,然后得其志。蚩尤作乱,不用帝命。于是黄帝乃征师诸侯,与蚩尤战于涿……

继 续 阅 读