广义几何规划的完全二阶原算法

对于广义几何规划问题,对目标函数和约束函数皆采用二阶缩并,化为对数变量空间中的约束二次的二次规划,采用文[7]的广义二次规划解法求解之。算法收敛快而稳定。     

用几何规划求解梁截面尺寸优化设计

梁的优化设计是一个非常复杂的问题,本文直接取截面尺寸为设计变量,用广义多项式处理方法化梁优化模型为符号几何规划,再用几何规划的切割平面法进行了求解,实现了多变量梁截面尺寸的优化设计,并用几个算例进行了验证。     

工程结构优化的序列二次规划

用序列二次规划可以高效率地解决许多工程结构的优化问题。一类问题将原来的非线性规划用力学处理和Taylor展开化为一系列二次规划问题,在目标函数可分离的情况下,通过Kuhn-Tucker条件进一步化为以Lagrange乘子为变量的准无约束的二次规划;在目标函数不可分离情况下,则直接求解原变量的近似二次规划,经过序列迭代,可以求解工程结构的静、动力优化问题。另一类问题是将原问题化为正多项式的几何规划,再将其对偶问题化为二次规划,序列迭代求解之,可以求解目标函数、约束函数非线性程度高的工程问题,对按设计规范公式进行设计的一类问题尤为有效。     

粘弹性阻尼结构系统振动特性的迭代分析

含粘弹性阻尼材料结构系统的振动特性分析是相当复杂的,求解比较困难,本文给出了一种有效的迭代解法。文中首先定义了广义Rayleigh商,指出求广义Rayleigh商的驻值问题等价于原特征问题的求解;由近似估计得到的特征向量,广义Rayleigh商给出一个较好的特征值近似估计。据此构造了迭代方法,并成功地计算了一些例子。     

多困难度几何规划的一个简单解法

推导出方向因子的计算式和一维搜索的迭代式,为求解多困难度几何规划问题提供了一个简单的计算方法。     

用广义变分原理求解多节点双模量静不定桁架

在外载荷作用下的多节点双模量静不定桁架平衡问题,是任意有限多个自变量的多元函数在任意有限多个约束条件下的极值问题,采用广义变分原理可以方便求解多节点双模量静不定桁架内力．通过求解多节点双模量静不定桁架内力的几个算例,阐述广义变分原理在计算多节点双模量静不定桁架内力中的应用．研究结果表明：采用广义变分原理求解多节点双模量静不定桁架内力的通用性较强,所求的结果是精确解析解．采用广义变分原理求解多节点双模量静不定桁架内力的方法不但克服了常规方法需利用几何关系建立协调方程的缺陷,且具有力学概念清晰直观、计算过程简便、便于工程设计人员在实际中掌握和计算等优点．     

随机桁架结构几何非线性问题的混合摄动-伽辽金法求解

提出应用混合摄动-伽辽金法求解随机桁架结构的几何非线性问题.将含位移项的随机割线弹性模量以及随机响应表示为幂多项式展开,利用高阶摄动方法确定随机结构几何非线性响应的幂多项式展开的各项系数.将随机响应的各阶摄动项假定为伽辽金试函数,运用伽辽金投影对试函数系数进行求解,从而得到随机桁架结构几何非线性响应的显式表达式.同已有的随机伽辽金法相比,本文所给的试函数由摄动解的线性组合而成,在求解非线性问题时,试函数的获取具有自适应性.数值算例结果表明,对于具有不同概率分布的多随机变量问题,本文方法无需对随机变量的概率分布形式进行转换,避免了转换误差,因而比同阶的广义正交多项式方法 (generalized polynomial chaos, GPC)计算精度高.同时,在结果精度相当时,和GPC方法相比,本文方法得到的试函数系数的非线性方程维度不大,方程的求解工作量小且更易求解.当随机量涨落较大时,混合摄动-伽辽金法计算所得的结构响应的各阶统计矩比高阶摄动法所得结果更逼近于蒙特卡洛模拟结果,显示了该方法对几何非线性随机问题求解的有效性.     

无额外自由度广义有限元非线性分析

将无额外自由度的广义有限元法由线弹性分析扩展到弹塑性大变形分析.局部强化函数的构建依赖于已有节点,不引入额外自由度,避免了线性相关性问题.在更新拉格朗日框架下,通过控制方程弱形式的线性化推导得到了节点内力的率形式,并分为材料和几何两部分.考虑超弹性和亚弹-塑性两种材料模型,采用Newton-Raphson迭代求解,给出了相关的一致切线刚度阵.三个典型算例的数值结果表明,本文发展的非线性无额外自由度广义有限元方法不仅能够准确求解超弹性和弹塑性大变形问题,同时相比于传统的线性有限元方法具有更高的精度.本文工作进一步拓宽了无额外自由度广义有限元方法的应用领域.     

非线性有限元分析的非协调模式及存在的问题

利用非协调模式提高非线性有限元分析广泛采用的低阶单元的精度和性能,是国际计算力学界研究的热点和难点.阐述了国际上在非线性有限元分析中已广泛采用的增广假设应变法方法(the enhanced assumed strain, EAS)的基本原理,详细讨论了非协调模式用于非线性有限元分析保证收敛、稳定的条件及增广假设应变场插值函数的构造方法.介绍了国内学者关于几何非线性非协调模式的研究方法和研究成果: (1)从Hellinger-Reissner广义变分原理出发,提出了几何非线性非协调模式的收敛条件,并采用非线性计算的若干简化措施建立几何非线性非协调元的简化模型;(2)一类放松单元间协调要求的非线性广义变分原理,对几何非线性问题可以选择事先无协调约束的非协调函数建立非协调元,收敛性可以保证,并根据此非线性广义变分原理可建立C$^1$或C$^0$类几何非线性广义杂交元,C$^1$或C$^0$类精化杂交元和精化直接刚度法.指出了EAS方法用于非线性有限元分析存在的问题,即本构关系和求解方法的限制,并对非协调元应用于非线性有限元分析提出了展望.      

杆系结构几何优化中的广义变量近似方法

提出一种杆系结构几何优化的广义中间变量近似方法，首先引入一套广义中间变量，包括各杆件局部柔度特性、方向余弦及内力与位移关系矩阵各元素，结构响应与这些广义中间变量的关系比与设计变量本身成更好的线性关系。因此以广义中间变量做一介泰勒展开近似位移和内力。应力约束和屈曲约束由近似内力计算，近似问题由优化器在设计变量空间内求解。最后给出了几个算例，结果表明本文的方法。是十分有效的。     

最优目标敏度紧系的桁架几何两步优化方法

对于桁架几何优化问题，采用由最优敏度紧系的序列两步优化解法：在完成截面优化之后，求出其最优目标对节点坐标变量的一、二阶导数，以此构造形状优化模型，对较难建模和求解的第二步问题建立了近似显式并能用二次规划有效地求解，取得了满意的结果．     

经典悬链线理论精确解与近似解的非线性数值计算

针对经典悬链线数学解中存在两个未知参数,即水平张力h和广义倾角α迄今尚未妥善解决的问题,进行了深入细致的分析。利用悬链线两点边值约束条件和不可拉伸假设,推导出求解隐含独立未知量水平张力的超越方程。引进互逆的无量纲参数求解超越方程中的水平张力,使得水平张力形式上具有最简单的参数依赖关系。探讨了广义倾角β,α和θ与几何参数的相互关系,得出广义倾角α不是独立未知参数的结论。提出了水平距离趋于0和趋于极限距离的各种近似解、在真小数全局计算范围内的近似解以及这些近似解关于精确解的误差程度,其结果在工程上具有应用价值。     

随机桁架结构几何非线性问题的混合摄动-伽辽金法求解

提出应用混合摄动$\!$-$\!$-$\!$伽辽金法求解随机桁架结构的几何非线性问题.将含位移项的随机割线弹性模量以及随机响应表示为幂多项式展开,利用高阶摄动方法确定随机结构几何非线性响应的幂多项式展开的各项系数.将随机响应的各阶摄动项假定为伽辽金试函数,运用伽辽金投影对试函数系数进行求解,从而得到随机桁架结构几何非线性响应的显式表达式.同已有的随机伽辽金法相比,本文所给的试函数由摄动解的线性组合而成,在求解非线性问题时,试函数的获取具有自适应性.数值算例结果表明,对于具有不同概率分布的多随机变量问题,本文方法无需对随机变量的概率分布形式进行转换,避免了转换误差,因而比同阶的广义正交多项式方法(generalizedpolynomial chaos, GPC)计算精度高.同时,在结果精度相当时,和GPC方法相比,本文方法得到的试函数系数的非线性方程维度不大,方程的求解工作量小且更易求解.当随机量涨落较大时,混合摄动$\!$-$\!$-$\!$伽辽金法计算所得的结构响应的各阶统计矩比高阶摄动法所得结果更逼近于蒙特卡洛模拟结果,显示了该方法对几何非线性随机问题求解的有效性.      

杆系结构几何优化中的广义变量近似方法

提出一种杆系结构几何优化的广义中间变量近似方法。首先引入一套广义中间变量，包括各杆件局部柔度特性、方向余弦及内力与位移关系矩阵各元素，结构响应与这些广义中间变量的关系比与设计变量本身（即节点坐标和截面积）成更好的线性关系。因此，以广义中间变量做一阶泰勒展开近似位移和内力。应力约束和屈曲约束由近似内力计算。近似问题由优化器在设计变量空间内求解。最后给出了几个算例，结果表明了本文的方法是十分有效的。     

非线性有限元的若干基本问题

本文介绍了非线性有限元中的若干基本问题。其中包括有关应变、应力和非线性平衡方程的一些基本概念,基于不同非线性广义变分原理的位移模式、杂交模式和拟协调模式几何非线性有限元及其在壳体屈曲问题中的应用等。      

广义变分原理在求解杆系装配应力中的应用

利用拉格朗日乘数法建立广义变分原理以求解有误差杆件结构装配应力.引入拉格朗日乘数并结合静力平衡方程,构造无条件广义变分原理的新泛函,求解新泛函的极值问题,获得超静定的变形协调方程,从而计算有误差杆件结构的装配应力.结果表明:该方法求解装配应力的通用性较强,不但克服传统几何方法建立变形协调方程的缺陷,而且具有计算过程简洁以及便于掌握等优点.     

Application of generalized variational principle for solving assembly stress of rod system

利用拉格朗日乘数法建立广义变分原理以求解有误差杆件结构装配应力. 引入拉格朗日乘数并结合静力平衡方程,构造无条件广义变分原理的新泛函,求解新泛函的极值问题,获得超静定的变形协调方程,从而计算有误差杆件结构的装配应力. 结果表明：该方法求解装配应力的通用性较强,不但克服传统几何方法建立变形协调方程的缺陷,而且具有计算过程简洁以及便于掌握等优点.     

直接有限元法求解广义磁热弹二维旋转问题

为了验证直接有限元法求解广义磁热弹耦合旋转问题的有效性及准确性,该文基于Lord和Shulman(L-S)广义热弹性理论,采用直接有限元方法,求解了置于磁场中的旋转半无限大体受热冲击作用的动态响应问题.文中给出了L-S型广义磁热弹耦合旋转问题的控制方程,建立了L-S型广义磁热弹旋转问题的虚位移原理,推导得到了相应的有限...     

无限元方法及其应用

限元是几何上趋于无穷的单元，它是一种特殊的有限元,也是对有限元在求解无界域 问题上的有效补充, 并可实现与有限元间的无缝连接.无限元分为映射无限元和非映射 无限元:映射无限元需要引入几何映射,在局部坐标系中构造插值形状函数,如Bettess 元和Astley元;非映射无限元则直接在整体坐标系中构造插值形状函数,如Burnett元. 本文评述求解无界域问题的无限元方法的研究现状和最新发展.首先介绍无限单元的概念 和无限元方法的特点;围绕求解以Helmholtz方程控制的波动问题,评述几种常规无限单 元的优劣,这些单元包括Bettess元、Astley元和Burnett元.然后介绍新近提出的广义 无限元方法,以及与常规无限元方法的区别与联系.最后对无限元方法在各种问题中的 应用做了总结.     

简支梁大变形的优化算法

提出一种求解几何非线性问题的优化算法,并研究了简支梁的几何非线性大变形问题。首先取简支梁大变形后的平衡状态为研究对象,分别创建它的变分模型和微分模型;然后基于微分模型,通过动坐标的迭代关系式求得微段端点坐标,构建微段端点未知坐标的目标函数;最后确定简支梁几何非线性大变形的最优化问题,并编制相应优化程序进行求解。通过分析典型算例,并同有限元方法的计算结果相比较,表明提出的优化算法在求解强几何非线性大变形问题中的正确性,为处理几何非线性大变形问题提供了一种新方法。