考虑多种非线性变形的RC框架单元模型

提出一种钢筋混凝土框架结构弹塑性地震反应分析的杆件单元模型 ,着重介绍节点滑移子单元和剪切子单元的柔度矩阵的建立过程 .该单元模型把框架结构构件 ,由分布弹塑性梁子单元、节点滑移子单元和剪切子单元等 3个在节点上串联的子单元组成 .这两个子单元模型能够分别考虑纵向钢筋在节点内的粘结滑移 ,构件关键受力区域的剪切变形随复杂加载历史的变化过程 .该法从子单元柔度矩阵的计算到单元刚度矩阵的建立过程 ,力学概念明确 ,大大简化多种变形成分和复杂加载历史下 ,构件单元刚度矩阵的建立过程 ,可应用于非线性地震反应分析     

土自重及地基刚度矩阵对板土相互作用的影响

将地基取为有限分层弹性体，采用三种不同方法计算地基刚度矩阵，根据工程背景下板嵌埋于土中这一特点，考虑土自重应刀对地基形变势能的影响，借助于能量变分原理获得了不同地基刚度矩阵时不同荷载作用下的两种不同抗弯刚度的矩形板内刀、板底反力和板底位移的数值解答。着重对比分析了地基刚度矩阵的不同算法及土自重应力对板内力及位移计算结果的影响，指出用单位集中力计算地基刚度矩阵，易受Ｂｏｓｓｉｎｅｓｑ解答在集中力作用；久附近“极性”的制约，计算结果不够理想，而以单位均布力为依据计算地基刚度矩阵，所得板内力及位移能很好地符合弹性解答基本规律。     

减少单元刚度矩阵计算量的一种方法

在固体力学使用等参数单元的有限单元分析中,单元刚度矩阵的计算量是较大的。节省单元刚度矩阵的计算量也是提高程序质量的一个方面。本文给出的计算单元刚度矩阵的方法,可以减少计算量。 如果单元具有n个自由度和l个广义应变分量,在单元刚度矩阵     

对刚度矩阵元定义的一点看法

用矩阵法求解振动问题时,系统刚度矩阵的计算是首要的步骤之一,刚度矩阵的多种计算方法中,直接从刚度矩阵元k_(ij)的定义计算刚度矩阵,是最常用的方法。对于离散系统或离散化的系统,这个方法不会有什么问题,但是应用此法于刚体平面运动时,会出现不顺当的情况,且会出现错的计算结果。     

聚氨酯泡沫复合夹层板弯曲问题的有限元解法

利用“分解刚度有限元素法”建立夹层板的弯曲刚度矩阵,将一个较为复杂的问题简化为一个单层薄板的弯曲刚度矩阵和一个纯剪切刚度矩阵相加,降低了数学问题的难度,减少了计算工作量。计算表明,该方法能够在经典理论解的基础上有效地提高计算精度,可以满足工程要求。     

弹性支撑上楔形杆件的稳定计算

在文「１」、「２」建立的三维退化直梁单元刚度矩阵及切线刚度矩阵的基础上,推导出计算弹性支撑上楔形杆件结构弹性稳定问题的有限元列式。     

矩形八结点单元刚度矩阵的通用公式

本文讨论分析了现行数值计算方法在计算平面问题的矩形八结点单元的单元刚度矩阵中存在的缺陷,给出了数学计算软件MAPLE在推导矩形八结点单元的单元刚度矩阵中的应用及其相应的MAPLE通用程序、单元刚度矩阵的通用计算公式,并利用MAPLE自带的codegen库自动给出了相应的C、FORTRAN代码,得出结论,利用该通用公式来计算其单元刚度矩阵矩阵,可以节约大量的计算时间与存储空间以提高计算速度与精度,并可将其应用到C、FORTRAN中。     

一个三维有限元法拱坝应力分析程序

本程序是为了计算拱坝应力用 TQ-16BCY 语言编写的一个专用程序。本程序是在 TQ-16机上充分利用外存(磁带)的基础上,对拱坝弹性体的位移和应力进行有限单元法的分析。用四面体单元和线性位移模式,由计算机半自动形成四面体单元,对每一分块逐层地形成四面体单元同时计算单元刚度矩阵并立即叠加到总刚度矩阵中去,单元刚度矩阵的计算公式进行了改变,每块的刚度矩阵采用标记对角元的变带宽一维存贮以及选择改进的平方根法来分块求解位移方程组。计算提供:节点位移、单元应力、单元主应力及主向和约束反力。     

半空间地下结构与土相互作用的链杆法研究

根据弹性地基上的半无限平面理论,利用计算弹性地基梁的链杆法,将地基的弹性刚度矩阵与上部结构刚度矩阵相协调、综合,引进上部结构与地基土的综合刚度矩阵,解决了结构与土的相互作用问题,计算方法简单、实效。     

模具温控凝固成形的数值模拟软件系统开发

根据软件工程的理论,采用结构化分析与设计方法,开发了专用的模具温控凝固成形（MTCS）数值模拟软件系统。首先进行了软件的需求分析建模,绘制了数据流图和数据字典,直观地描述了单元信息、边界条件、刚度矩阵等数据的流动。在需求分析的基础上进行了系统设计,将程序划分成相变潜热、刚度矩阵计算及矩阵求解等10个主要模块,并编制了主模块与各子模块的程序流程图,定义了模块之间的数据接口。软件的前后处理模块的设计与实现在MSC．male软件上进行,求解器模块则采用fortran语言编程实现。最后设计了测试实例对软件进行了测试,并与MSC．marc软件的计算结果进行了对比,误差在3％以内。     

环形板元素的动刚度矩阵

本文用有限元的动刚度矩阵法推导了环形板的动刚度矩阵和传递矩阵。这样,便可用一个二阶矩阵来计算变厚度圆板或环形板的动态固有特性。     

大转动三角形平面单元有限元分析的共旋坐标法

基于能显著减小单刚计算量的场一致性原则,遵循标准的共旋坐标法,导出了三角形平面单元在大转动、小应变条件下的单元切线刚度矩阵.利用这一非对称的单元切线刚度矩阵编制程序,运用该程序对悬臂端作用有集中力的平面悬臂梁进行了计算.计算结果表明,所提出的单元切线刚度矩阵列式正确,该单元刚度矩阵虽然不对称,但计算较简单.这在非线性计算中对于减小由于计算机位数限制带来的累积舍入误差和提高迭代的收敛性具有重要意义.     

Stokes问题的Galerkin-Hermite小波方法

将圆内区域Stokes方程组的第二边值问题归化为Hadamard型强奇异自然积分方程组,然后运用Galerkin-Hermite小波方法将其变分问题的积分核离散化,获得了简单的刚度矩阵计算公式.对一个2J+3×2J+3阶的刚度矩阵,只需计算2J+3J+7个元素,而且刚度矩阵是块对角矩阵,每个块矩阵又是循环对称的,计算量大为减少,计算速度和精度也显著提高.     

基于改变结构刚度模态试验的动力学模型修改

提出了一种新的动力学有限元模型刚度与质量矩阵的修正方法.在分别测量原结构动力学系统的模态和改变结构刚度后的动力学系统模态基础上,利用矩阵计算和代数方程求解,计算出原结构有限元模型的刚度、质量矩阵的修正量.文中用实例研究了用非完整模态集进行模型修正时,实测模态阶数与模型修改后刚度、质量阵元素及模态频率、振型的计算精度之间的关系,从而证明了本文所给出方法的可用性以及当误差局限于局部小区域时该方法所具有的易用性与准确性.     

覆盖位移函数对数值流形方法刚度矩阵的影响

基于数值流形方法的覆盖位移函数理论,从覆盖位移函数对刚度矩阵形成关系的角度出发分析了其对刚度矩阵求解的影响。通过分析可知,刚度矩阵中矩阵元素数值递增速度与覆盖位移函数中绝对坐标x、y相关。并提出一种改进的覆盖位移函数来改善在采用高阶覆盖位移函数情况下刚度矩阵的局部刚度过大的问题,以利于矩阵求解计算。最后通过算例验证了该方法是有效的。     

基于类平衡矩阵的桁架构件重要性计算方法

为简化构件重要性系数的计算过程，提高计算效率，以基于刚度的构件重要性判断方法为基础进行数学推导，提出一种基于类平衡矩阵的空间结构重要构件的快速判断方法，并进行算例验证. 结果表明，当铰接杆件体系的各杆件线刚度近似相等，且不考虑结构应力刚化作用或应力刚化所贡献的刚度很小时，整体切线刚度矩阵仅与类平衡矩阵相关，构件重要性系数可直接利用类平衡矩阵的奇异值进行计算. 通过算例验证，类平衡矩阵的重要性构件判断方法与基于刚度的重要性构件判断方法计算结果一致，说明所提方法是一种切实可行的快速计算方法.     

空间过渡梁元的切线刚度矩阵

详细推导了整体拉朗日（Total Iagrange)坐标系下空间过渡梁元的切线刚度矩阵和单元抗力计算式。计算表明,过渡梁元切线刚度矩阵的, 使局部单层局部双层网壳的计算结果更趋近于实际。     

齿板滑移对轻型木桁架挠度的影响

以齿板连接轻型木桁架为研究对象,研究了齿板滑移对木桁架各结点位移的影响。采用有 限元分析方法,运用弹簧的串联原理将两个虚拟变刚度单元及一根实体杆单元整合为等效刚度杆件单元,建立单元刚度矩阵和整体刚度矩阵,再运 用普通桁架的计算理论来计算引入变刚度单元模拟齿板连接的受力情况及位移情况。根据以上原理用Fortran90编制了齿板连接木桁架MPT Design设计程序,并根据木结构试验规范要求通过木屋架试验验证该程序的正确性。     

关于轴对称有限元刚度矩阵的精确积分

叙述了轴对称问题线性位移函数的有限元算法中，计算刚度矩阵时进行精确积分的 方法．并采用精确积分刚度矩阵对一些问题进行了计算，其结果与解析解作了比较。     

沥青路面温度应力分析

把沥青路面视为多层弹性半空间轴对称体,利用热弹性力学以及Hankel和Laplace积分变换等数学方法,首先推导出任意一层沥青路面温度应力的刚度矩阵,然后按传统的有限元方法组成总体刚度矩阵.通过求解由总体刚度矩阵所构成的代数方程和Hankel及Laplace积分逆变换就可解出外荷载和温度联合作用下沥青路面温度应力的精确解.由于刚度矩阵的元素中不含有正指数项,计算时不会出现溢出现象,从而克服了传递矩阵法的缺点.在推导过程中不用人为选择应力函数,使得问题的求解更加合理,同时也为进一步研究沥青路面的湿度场、动力学等问题奠定了理论基础.计算实例证明了推导结果的正确性.