连续梁的多结点力矩分配法

 计算了连续梁在结点力偶荷载作用下柔度矩阵的表达式，在此基础上，提出了 多结点力矩分配法. 对于一个承受任意荷载的$N$个结点的连续梁，该方法同时松开多个 结点，通过迭代获得各结点转角值，然后利用杆件的转角位移方程就可获得杆端弯矩的值. 该方法综合了位移法和力矩分配法的优点，较传统力矩分配法可显著减少计算工作量.     

多功能力矩分配法

提出了新的力矩分配思想，将角位移和线位移综合考虑，使力矩分配法能直接应用于有侧移结构，且在计算内力的同时还能计算结构的位移，使力矩分配法具有多功能. 实例计算表明，该方法正确、简单，有利用价值.     

回传波矩阵法在杆系结构静力分析中的应用

回传波矩阵法最初是由Pao等人分析二维框架结构动力响应时提出的。对于三维杆系结构的静力分析，为了确定结构的位移和内力，先要建立传递分配矩阵和载荷源向量，这可通过列出所有节点的静力平衡方程和位移协调方程来实现。同时，通过分析每根杆近端位移和远端位移的关系，建立结构的回传波矩阵（重分配矩阵）。在此基础上求解线性方程组，就可以得到结构的位移和内力。本文推导了空间杆系结构的有关矩阵方程式，并给出了一固定梁的两端弯矩求解算例。     

带有全铰节点的三维杆系结构静力分析的回传波射矩阵法研究

推导出含有全铰节点的三维杆系结构的回传波射矩阵表达式,完善了具有任意连接和约束的空间杆系结构静力分析的回传波矩阵法.基于节点平衡方程和协调方程,推导出表达杆件近端位移和远端位移关系的传递分配矩阵及载荷源向量,并通过由对偶坐标系下近端位移和远端位移的关系获得结构的总体相位矩阵,再引入转列矩阵,进而推导出结构的回传波射矩阵,在此基础上求解以杆端位移为基本未知量的线性方程组,最终得到精确确定所有杆件的杆端位移及杆端内力的矩阵列式.给出了空间杆系结构算例分析,与有限元结果比较,验证了回传波射矩阵法的计算精度.     

斜交连续梁桥数值分析研究

在分析了矩阵位移法和传递矩阵法的理论基础上，将二者组合而成矩阵混合法。矩阵混合法的实质就是将结构物的一端未知向量通过矩阵的简单乘法传递到另一端，再将已知条件代入求解的过程。在矩阵传递的过程中引入位移法中的前进代入的概念，即将未知向量不断地以新的未知向量替换，可以防止误差的传递和积累。将矩阵混合法运用于斜交连续梁桥内力计算，完成了矩阵混合法在斜交连续梁桥内力计算中的理论推导和程序的设计，给出了实际算例并与有限元法进行了比较。     

含任意椭圆核各向异性板杂交应力有限元

基于含椭圆核有限大各向异性板弹性问题的复变函数级数解,应用杂交变分原理建立了一种与常规有限元相协调的含任意椭圆核各向异性板杂交应力有限元.单元内的应力场和位移场采用满足平衡方程、几何方程与物理方程的复变函数级数解,假设的复变函数级数解精确满足椭圆核边界处的位移协调条件和应力连续条件,单元外边界上的位移场按常规有限元位移场假设,单元内椭圆核的长轴可以与材料主轴不重合.单元刚度矩阵采用Gauss积分求得,并给出了建立刚度矩阵的主要公式和推倒过程.数值计算结果表明该单元具有计算精度高、计算工作量小等优点.     

不可压缩平面问题的位移-压力混合重心插值配点法

引入人工压力变量,将弹性本构方程以应力、应变和压力表达,建立求解不可压缩平面弹性问题的位移-压力方程和不可压缩条件方程的耦合偏微分方程组。利用张量积型重心Lagrange插值近似二元函数,得到计算插值节点处偏导数的偏微分矩阵。采用配点法离散不可压缩弹性控制方程,利用偏微分矩阵直接离散弹性力学控制方程为矩阵形式方程组。利用插值公式离散位移和应力边界条件,将离散边界条件与离散控制方程组合为新的方程组,得到求解弹性问题的过约束线性代数方程组;利用最小二乘法求解线性方程组,得到弹性力学问题位移数值解。数值算例验证了所提方法的数值计算精度为10-14~10-10。     

一种结构动力时程分析的积分求微方法

传统采用微分求积(differential quadrature,DQ)法求解动力问题时都是以位移响应作为基本未知量,而将速度响应和加速度响应表示为位移响应的加权和的形式.如此做法需要处理线性方程组或者矩阵方程(Sylvester方程)才能求得动力响应,导出的算法一般为有条件稳定算法.本文利用动力响应的Duhamel积分解,逆用DQ原理,提出了一种计算卷积的高精度显式算法.该算法可以逐时段地求解出动力时程响应,当各时段内DQ节点分布完全一致时,仅须进行一次Vandermonde矩阵求逆计算即可应用于各个时段,一次性获得时段内多个时刻的位移响应值,因而具有计算效率高的优点.通过分析动力方程积分格式,证明本文动力算法传递矩阵的谱半径恒等于1,因而该算法具有无条件稳定特性,且计算过程中不会产生数值耗散. 本文算法的数值精度取决于分析时段内布置的DQ节点数量$N$,具有$N-1$阶代数精度.实际操作时可以取10个甚至更多的DQ节点数,从而获得比较高的数值精度.      

一种结构动力时程分析的积分求微方法

传统采用微分求积(differential quadrature, DQ)法求解动力问题时都是以位移响应作为基本未知量,而将速度响应和加速度响应表示为位移响应的加权和的形式.如此做法需要处理线性方程组或者矩阵方程(Sylvester方程)才能求得动力响应,导出的算法一般为有条件稳定算法.本文利用动力响应的Duhamel积分解,逆用DQ原理,提出了一种计算卷积的高精度显式算法.该算法可以逐时段地求解出动力时程响应,当各时段内DQ节点分布完全一致时,仅须进行一次Vandermonde矩阵求逆计算即可应用于各个时段,一次性获得时段内多个时刻的位移响应值,因而具有计算效率高的优点.通过分析动力方程积分格式,证明本文动力算法传递矩阵的谱半径恒等于1,因而该算法具有无条件稳定特性,且计算过程中不会产生数值耗散.本文算法的数值精度取决于分析时段内布置的DQ节点数量N,具有N-1阶代数精度.实际操作时可以取10个甚至更多的DQ节点数,从而获得比较高的数值精度.     

考虑在地基固结和流变的上部结构和地基基础的共同作用问题

本文提出一种能同时考虑地基的瞬时变形,团结变形和流变的地基基础与上部结构共同作用的一种新的计算方法,在这种方法中,土的计算考虑了Biot固结理论和流变,上部结构利用子结构法进行凝聚可得边界刚度矩阵和边界力,利用位移协讯条件则可以得到共同作用的控制方程,对上述方程进行Laplace变换,可以得到Laplace变换域内的控制方程,在Laplace变换域内对上述方程进行数值求解并进行相应的Laplace逆变换则可以得到时间境内任意时刻的解。     

拱桥悬臂施工过程中面内特征值的传递矩阵法

采用传递矩阵法对拱桥悬臂施工过程中的面内特征值问题进行求解,建立了该类桥型施工过程中面内竖弯刚度的评估方法.首先,将索和拱分别视为无垂度的张紧弦和欧拉伯努利梁,基于传递矩阵法基本理论推导了系统的总传递矩阵,考虑拱和索的边界条件以及索拱节点的位移连续性条件得到系统的特征值方程,进而计算出系统的频率和模态.同时,采用有限元...     

力矩分配法在对称结构中的应用

以力矩分配法为基础，探索了对称结构的简化计算方法．为了克服对称结构计算时通常方法所产生的不便，提出了新的力矩分配概念，找出了新的分配系数和传递系数．应用改进的力矩分配法，对对称结构进行了计算，算例表明该方法简化了计算，加快了收敛速度．     

具有任意连接和约束的空间杆系结构静力分析的回传波矩阵法

本文发展了具有任意连接和约束的空间杆系结构静力分析的回传波矩阵法。以杆端位移和转角为基本未知量,通过结构所有节点的平衡方程和位移协调条件,推导出传递分配矩阵和载荷源向量,并进一步利用设定的同一杆件两个局部坐标系下杆端位移之间的关系,最终得到结构的回传矩阵。据此可求出结构所有杆件的杆端位移及杆端内力。对不同的杆件连接形式,如刚接、铰接、半刚接,以及不同的约束情况,如固定支座、铰支座、定向支座等,本文推导出了空间杆系结构的回传波矩阵表达式,可直接用于相应空间杆系结构内力的计算。同时,针对一个具体刚架结构进行了算例分析,并通过与弯矩分配法和有限元结果的比较,验证了本文方法的精确度。     

集中载荷作用下薄壁开口曲杆的稳定性

本文研究了薄壁开口曲杆的稳定性问题.该曲杆的边界条件比较复杂,一端固支,另一端既有对位移的限制,又有集中力和集中力矩作用.本文从薄壁开口曲杆的平衡方程出发,应用假想分布载荷法,导出一组关于曲杆稳定性的变系数常微分方程.为了解特征值问题,应用有限差分法,将特征方程用矩阵形式表达.对某一产品结构中的实际构件进行了数值计算与稳定性试验,构件为一Ⅰ字形截面的薄壁开口曲杆,对特征方程应用矩阵迭代法求得临界载荷的理论计算值.应用Southwell法结合最小二乘法整理了重复试验所获得的数据,得到临界载荷的实验值,理论计算值与实验值接近.     

地震作用下空间框架结构层间相对位移响应梯度和海赛矩阵的计算

提出一种在地震作用下空间框架结构层间相对位移响应梯度和海赛矩阵的精确计算方法。在有限单元法和纽马克-β法的基础上推导出在地震作用下空间框架结构层间相对位移响应梯度和海赛矩阵的计算公式,给出在地震作用下空间框架结构层间相对位移响应梯度和海赛矩阵的计算步骤,用Matlab语言编制了空间框架结构层间相对位移响应梯度和海赛矩阵的计算程序,实现了空间框架结构层间相对位移响应梯度和海赛矩阵的精确计算。最后通过一个二层空间框架结构的计算实例表明本文所提出的计算方法是有效的。     

阶梯梁的位移计算方法对比

为了加深学生对位移法及阶梯梁位移计算的理解，以阶梯式简支梁为例，分别利用积分法、常规图乘法、刚度叠加法及位移法求解了阶梯式简支梁的位移。计算结果表明：利用位移法求解阶梯式简支梁的位移具有独特的优势，可显著减少计算工作量并提高计算准确率，值得在教学过程中推广。     

平面弹性问题的位移-应力混合重心插值配点法

提出数值分析平面弹性问题的位移-应力混合重心插值配点法。将弹性力学控制方程表达为位移和应力的耦合偏微分方程组,采用重心插值近似未知量,利用重心插值微分矩阵得到平面问题控制方程的矩阵形式离散表达式。使用重心插值离散位移和应力边界条件,采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,应用最小二乘法求解过约束方程组,得到平面弹性问题位移和应力数值解。数值算例结果表明,重心Lagrange插值方法的计算精度可达到10~(-10)量级。位移-应力混合重心插值配点法的计算公式简单、程序实施方便,是一种高精度的无网格数值分析方法。     

不规则区域平面弹性问题的正则区域重心插值配点法

将不规则区域嵌入到规则的矩形区域,在矩形区域上将弹性平面问题的控制方程采用重心Lagrange插值离散,得到控制方程矩阵形式的离散表达式。在边界节点上利用重心插值离散边界条件,规则区域采用置换法施加边界条件,不规则区域采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法进行求解,得到整个规则区域上的位移数值解。利用重心插值计算得到不规则区域内任意节点的位移值,计算精度可到10-14以上。数值算例验证了所建立方法的有效性和计算精度。     

基于等效刚度的矩形孔蜂窝梁钢框架结构内力计算方法

蜂窝梁钢框架结构因梁截面沿长度周期性变化，不能直接采用普通钢框架结构矩阵位移法计算框架内力和位移．本文基于等效刚度法推导了矩形孔蜂窝梁的等效抗弯刚度、抗剪刚度和轴向刚度，建立了矩形孔蜂窝梁的单元刚度方程，提出了矩形孔蜂窝梁钢框架内力和位移计算方法．算例理论计算结果与有限元分析结果表明，两种方法计算结果非常接近．本文提出的等效刚度法概念清晰，准确性好，适用于计算蜂窝梁钢框架结构的内力和位移．     

不规则区域平面弹性问题的正则区域重心插值配点法

将不规则区域嵌入到规则的矩形区域,在矩形区域上将弹性平面问题的控制方程采用重心Lagrange插值离散,得到控制方程矩阵形式的离散表达式。在边界节点上利用重心插值离散边界条件,规则区域采用置换法施加边界条件,不规则区域采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法进行求解,得到整个规则区域上的位移数值解。利用重心插值计算得到不规则区域内任意节点的位移值,计算精度可到10^-14以上。数值算例验证了所建立方法的有效性和计算精度。