共查询到20条相似文献,搜索用时 15 毫秒
1.
移动荷载作用下地基动力分析的有限元方法 总被引:3,自引:0,他引:3
通过对地基动力问题的基本方程进行变换,把基本方程变换到随荷载移动的运动坐标系中,通过加权残数法推导了相应的单元刚度矩阵,从而建立了移动问题的有限元格式,并发现移动荷载问题的单元刚度矩阵是对相应静力问题单元刚度矩阵的修正,在静力单元刚度矩阵的主对角元素上增加与移动速度有关的项,即可得到移动问题有限元的单元刚度矩阵,这样就将动力学问题转化为“拟静力”问题处理。文中用移动问题有限元方法计算了地基的动力响应,并与解析解进行了对比,以说明本方法具有较好的精度。 相似文献
2.
一种考虑初始垂度影响的非线性索单元 总被引:2,自引:0,他引:2
本文首先运用Mathematic的符号运算功能,通过解偏微分方程给出了不等高单索的显式解析解,然后把该解析解应用于索微元的应变计算中,在此基础上用虚功原理推导了考虑初始垂度影响的两节点非线性索单元,并给出了索单元的单刚矩阵的具体形式,通过与两节点直线索单元的单铡矩阵的形式的比较,明确了该曲线非线性索单元的修正项,并给出了垂度影响因子的变化曲线。比较直观的给出了垂度对索单元刚度各项的影响程度。该非线性索单元既有多节点索单元精度高的特点,又有节点少,刚度元素较易求解以及有限元列式简洁等特点。本文通过两个数值算例表明,本文的非线性索单元是正确的,也表明了所编制的非线性计算程序是正确和可靠的。 相似文献
3.
构造了一种3节点三角形膜单元,以适用于平面薄膜横向振动的有限元分析.在给出单元形函数的基础上,根据最小势能原理建立了薄膜自由振动方程,并推导了单元刚度矩阵和单元质量矩阵.研究结果表明,单元刚度矩阵和单元质量矩阵形式简单,且自由度少;通过两个典型算例,证明3节点三角形膜单元的计算结果非常接近理论解,同时可以达到很高的精度... 相似文献
4.
构建金属桁架结构航天器陨落再入气动热环境有限元传热模型,是准确预测在轨服役期满大型航天器陨落再入解体过程温度分布的关键。本文采用四节点四面体单元对空间进行离散,依据泛函理论,将传热控制方程离散为代数方程组;利用有限单元法总体合成得到具有对称正定、高度稀疏和非0元素分布的规则性刚度矩阵,发展一维变带宽压缩存贮技术,有效解决大型稀疏矩阵的数据存贮问题;为有效抑制求解过程出现的温度在时间和空间上的振荡问题,发展集中热容矩阵系数处理方法,将热容矩阵的同行或同列元素相加代替对角线元素,使非对角线元素化为0,构造求解三维瞬态温度场的两点向后差分格式、Crank-Nicolson格式和Galerkin格式。通过对正方体瞬态传热计算验证分析,在相同条件下,采用以上三种格式均可获得一致稳定的温度解,并得到与现有ANSYS有限元软件较为吻合的计算结果,验证了所建立三维瞬态传热有限元计算模型的准确可靠性。在此基础上,对铝合金低轨航天器薄壳结构进行了传热计算,给出了类天宫飞行器两舱体陨落飞行107.5 km~90 km不同高度的瞬态温度分布,为这类寿命末期航天器陨落再入解体预报提供理论支撑与可计算模型。 相似文献
5.
构建金属桁架结构航天器陨落再入气动热环境有限元传热模型,是准确预测在轨服役期满大型航天器陨落再入解体过程温度分布的关键。本文采用四节点四面体单元对空间进行离散,依据泛函理论,将传热控制方程离散为代数方程组;利用有限单元法总体合成得到具有对称正定、高度稀疏和非0元素分布的规则性刚度矩阵,发展一维变带宽压缩存贮技术,有效解决大型稀疏矩阵的数据存贮问题;为有效抑制求解过程出现的温度在时间和空间上的振荡问题,发展集中热容矩阵系数处理方法,将热容矩阵的同行或同列元素相加代替对角线元素,使非对角线元素化为0,构造求解三维瞬态温度场的两点向后差分格式、Crank-Nicolson格式和Galerkin格式。通过对正方体瞬态传热计算验证分析,在相同条件下,采用以上三种格式均可获得一致稳定的温度解,并得到与现有ANSYS有限元软件较为吻合的计算结果,验证了所建立三维瞬态传热有限元计算模型的准确可靠性。在此基础上,对铝合金低轨航天器薄壳结构进行了传热计算,给出了类天宫飞行器两舱体陨落飞行107.5 km~90 km不同高度的瞬态温度分布,为这类寿命末期航天器陨落再入解体预报提供理论支撑与可计算模型。 相似文献
6.
在有限元位移法中,单元矩阵(刚度矩阵、几何矩阵、一致质量矩阵、载荷列阵等)往往是在单元局部坐标下计算的(以下简称局部方位);结构整体矩阵的形成,单元应力和变形的计算,都要求确定总体节点参数(自由度)和单 ... 相似文献
7.
本文基于拟协调元的基本原理,引入薄壳分层子单元的概念,提出了用于圆柱薄壳结构弹-塑性分析的三角形和矩形拟协调分层圆柱薄壳单元的模型,给出了第k层薄壳分层子单元的弹-塑性矩阵和物理矩阵[D~(k)]的表达式,并最后给出了整个薄壳单元的物理矩阵和刚度矩阵的表达式。 根据上述原理和方法用Fortran语言编制了圆柱薄壳结构弹-塑性有限元分析程序,实例计算与实验结果较为吻合,从而说明用本方法进行圆柱薄壳结构的弹-塑性分析是可靠而行之有效的。 相似文献
8.
基于含椭圆核有限大各向异性板弹性问题的复变函数级数解,应用杂交变分原理建立了一种与常规有限元相协调的含任意椭圆核各向异性板杂交应力有限元.单元内的应力场和位移场采用满足平衡方程、几何方程与物理方程的复变函数级数解,假设的复变函数级数解精确满足椭圆核边界处的位移协调条件和应力连续条件,单元外边界上的位移场按常规有限元位移场假设,单元内椭圆核的长轴可以与材料主轴不重合.单元刚度矩阵采用Gauss积分求得,并给出了建立刚度矩阵的主要公式和推倒过程.数值计算结果表明该单元具有计算精度高、计算工作量小等优点. 相似文献
9.
10.
11.
面向对象的有限元程序设计 总被引:23,自引:2,他引:21
本文按照面向对象的程序设计方法,遵循有限元分析的本质,采用C++语言,建立了有关描述有限元模型的类,给出了类的描述和它的实现方法。相关的类和方法包括处理矩阵的类、节点类、单元类、材料类、形函数类等。据此编制了有限元分析的数值计算程序,并给出了一个实例。 相似文献
12.
为明确ANSYS中基于铁摩辛柯梁理论的Beam188单元的单刚矩阵,以形函数为基础,根据最小势能原理系统给出了铁摩辛柯梁单元的刚度矩阵推导过程,尤其将常规2节点单元的刚度矩阵扩展至3节点单元,包括等截面梁和变截面梁,并给出了各自的单刚矩阵理论表达式。以矩形、圆形、箱形和圆环四种截面为例,分别给出单刚矩阵理论值与ANSYS中Beam188单元刚度矩阵进行对比。结果表明,3节点等/变截面梁的单刚矩阵与理论推导所得结果一致,但2节点梁的单刚矩阵与理论推导所得结果存在偏差,且此偏差仅表现在与弯曲相关的元素中,并且随着单元长度的减小,两者的偏差也越来越小。另外,对比时应注意,ANSYS所得截面参数如剪切系数和极惯性矩与理论值存在较大偏差。 相似文献
13.
无单元法在薄板稳定问题中的应用 总被引:8,自引:1,他引:8
用无单元法研究了薄板的弹性稳定问题,从滑动最小二乘法和变分原理出发导出了薄板的无单元法几何刚度矩阵,编制了相应的计算程序,并给出了算例,结果表明,方法合理可行,且精度高于有限元。 相似文献
14.
本文在文献[1]给出的放松应力增量平衡约束的修正余能广义变分原理基础上,提出一种高效率的弹塑性有限元分析的新方法——增量杂交/混合修正弦线模量法。该法保持了文[1]方法的全部优点,而在迭代过程中,依据材料的单向拉伸应力—应变关系,不断改变过渡区和塑性区单元柔度矩阵和塑性矩阵中的弹性模量;并在体积压缩模量不变假设下,相应地改变过渡区单元矩阵中的泊松系数。从而大大降低了迭代收敛次数和单刚计算量,提高了多类变量弹塑性有限元分析的计算效率和收敛精度。 相似文献
15.
四、装配成结构刚度矩阵 在结构矩阵分析的过程中,把复杂结构的计算问题转化为简单的单元分析和整体装配问题。在第一讲中已经介绍了单元分析,现在就来讨论如何把各个单元刚度矩阵装配成结构刚度矩阵(或简称为总刚)。关键的问题是要知道单刚的各个元素在总刚中的下标,即单刚各元素在总刚中的行号和列号。对应于每一个单刚来讲,应有一个下标矩阵。 相似文献
16.
模拟火炮后坐运动的非线性弹簧-阻尼单元 总被引:1,自引:0,他引:1
根据非线性有限元法的基本理论,建立了火炮发射动力学的有限元运动方程,推导出单元耗散为向量和单元阻尼矩阵的计算公式。提出采用非线性弹簧-阻尼单元来模拟反后坐装置对炮身后坐的阻尼作用,并给出了该单元的刚度矩阵和阻尼矩阵。根据炮身的后坐过程的三个阶段,提出采用分段拟合的方法确定非线性单元的参数以模拟火炮后坐阻力-位移曲线。采用非线性有限元单元模拟后坐装置,能够较准确地描述炮身部分的后坐运动,同时在大简化了全炮模型的建立,为全炮的动力学仿真提出了一条新的思路。 相似文献
17.
18.
本文构造了一族圆筒壳的部分正交升阶谱协调元素,具体给出了位移函数、均匀圆筒壳刚度与质量矩阵中对应于高阶插值函数的自由度部份的统一显表达式,以保证刚度与质量矩阵的计算精度和推迟矩阵方程出现病态情况,从而能充分发挥该元素的p-收敛特性的优点.然后,用该族元素计算了等厚度、轴向与环向变厚度圆筒壳的固有振动频率和振型,并给出了它们随厚度变化而改变的趋势.结果表明,与其它元素相比,本文元素具有带宽较小、收敛速度较快、计算效率较高等优点. 相似文献
19.