层合圆柱厚壳热应力分析的状态方程及其半解析法

通过对Ｈｅｌｌｉｎｇｅｒ－Ｒｅｉｓｓｎｅｒ变分原理的修正，导出了变温作用下层合圆柱厚壳的状态方程及其半解析法，该方法在ｚ－θ曲面内采用通常的有限元离散，而沿壳厚（ｒ）方向采用状态空间法给出解析解答，且通过采用状态转移矩阵，建立了层合圆柱壳内外表面应力和位移之间的关系式，然后利用打靶法进行求解，从而大大降低了计算中的未知量数目。     

球形厚壳的应力分析

文中将球形厚壳视为三维弹性体，按三维弹性理论问题求解其应力及位移。     

厚壳结构弹塑性分析的半解析有限元解

本文提出了对厚壳弹塑性分析的一种有效的、简便的分析方法,该法适用于圆柱形、任意截面柱形、旋转形等各种中厚壳,强厚壳的物理非线性分析。文中给出了厚壳弹塑性阶段的变形、应力状态和塑性区发展与破坏规律,为建立厚壳弹塑性理论提供了必要的依据,为厚壳非线性分析的实用计算探索了一条新的途径。     

三维弹塑性问题的半解析解

1.前言半解析数值方法是一类在各科学计算与工程分析领域中有着广泛应用前景的新的计算方法,在各种工程力学分析中已有广泛应用。但这种方法目前只限于线性问题。很少有物理非线性问题结果,特别是三维问题。半解析元法综合了解析法与有限元法的特点,是比较适合于求解弹塑性问题的一种方法。由于它在空间离散中引入了解析函数,使其在空间上由三维降为一维数值问题(大约省二个数量级),使弹塑性分析也得以在一般微机上实现,为固体力学中三维弹塑性问题求解探索一条较为简便而实用的途径。本文探     

半圆柱壳挠曲问题解析解精确性的研究

壳体力学已于上世纪由多位专家发展成熟,其中简支柱壳挠曲问题采用改进莱维解法的三角级数法解出,但是其解法复杂,手算难以完成．为讨论其结果的精确性,通过编写运行基于MATLAB的运算程序导出实例化解析解,与基于力学基本理论的推想假设对比,再引入有限元计算结果进行比较研究．最终发现,理论解析解应力和位移具有分布形式大致准确性,但仍存在不容忽视的细节与局部性问题．研究表明,理论解法工程意义有限,结果尚需改进．     

厚圆柱壳的非轴对称振动分析

本文对厚圆柱壳的非轴对称振动进行了分析，其中除包含通常的薄膜和弯曲效应外，还反映了转动惯性，横向剪切变形和横向挤压的影响，数值结果表明：对于厚圆柱壳来说存在着频率密集区，频率位置发生移动，横向挤压的影响必须要考虑。     

大型纵向加筋截断锥壳稳定性的半离散分析

如何计算大型纵向加筋截断锥壳的稳定问题，尚未见到有关报导。本文在沿纵向和环向分别采用Ｈｅｒｍｉｔｅ插值和三角级数插值的基础上，建立了几何非线带筋环单元，进而导出了半离散分析法，并经过实验检验了它的正确与可靠。这样，本文为建立各类工程结构的半离散单元提出了一种新的途径。     

两端固支叠层圆柱厚壳轴对称问题的精确解析解

基于柱坐标系下的三维弹性力学基本方程,采用状态空间法得到两端固支单层与叠层圆柱厚壳轴对称问题的精确解析解。为严格满足固支端的边界条件,将固支端的边界位移函数作为状态变量引入状态方程,采用增维方法把非齐次状态方程变为齐次状态方程,并通过层合渐近技术将变系数状态矩阵转为常系数矩阵进行求解。所得到的解不仅严格满足三维弹性力学基本方程,而且严格满足固支边界条件,是真正意义上的三维精确解。算例表明,本研究解与有限元解吻合,具有很高的精度,且关于级数项数和分层数具有很好的收敛性。另外,通过圆柱厚壳各力学量沿径向和轴向的精确分布规律分析了厚径比和跨径比变化对位移和应力分布的影响。     

拉压异性线性等向强化材料厚壁球壳极限分析

运用Mohr屈服准则对承受内压的拉压屈服强度不同的线性等向强化材料的厚壁球壳进行了极限载荷分析，得到了依赖于拉压比和强化模量的厚壁球壳极限载荷解析式和依赖于拉压比和强化模量的厚壁球壳极限载荷分析式。结果表明，材料拉压屈服强度的不同和强化特性对厚壁球壳极限载荷均有一定的影响。     

轴压圆柱壳后屈曲的有限元分析

采用流动的曲线坐标建立了大位移下的样条曲壳元，利用它分析了轴压圆柱壳的后屈曲问题。一个跨越和求取分支点的技巧和对弧长法的一个改进在文中提出。计算所得的载荷－轴向缩短曲线较解析解更接近于实验结果，所得的后屈曲波形也和实验相当符合。     

NUMERICAL ANALYSIS OF MINDLIN SHELL BY MESHLESS LOCAL PETROV-GALERKIN METHOD

The objectives of this study are to employ the meshless local Petrov-Galerkin method （MLPGM） to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors＇ sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.     

浸没的球面各向同性球壳自由振动的中厚壳理论

采用六模态的中厚壳理论研究了可压缩流体中球面各向同性球壳的自由振动，即在分析中考虑了剪切变形，旋转惯量和横向正应变的影响，引入５个辅助变量可以理到两类振动及其频率方程，对频率方程作了简化并算例进行了相应的探讨。     

???????????????????????

对具有不同拉压模量的厚壁球壳,采用双剪统一强度理论推导了其扩张问题的应力及位移的 统一解. 分析了不同模量、不同模型控制参数对厚壁球壳扩张时的扩张压力和应力场的影响. 结果表明：厚壁球壳弹性极限压力、应力场、位移场等均随着模量控制参数、模型参数的变 化而变化,在$\alpha<1$的情况下(即$E^ + < E^ - $),可以明显提高球壳的弹 性极限压力$p_e $; 厚壁球壳塑性极限压力与材料的拉压模量无关,与模型参数$\eta $有关,且随$\eta$的增加,先增大后减小. 因此若采用经典的弹性理论和单一的 模型参数对厚壁球壳进行设计计算,会带来较大的误差.     

不同拉压模量厚壁球壳弹塑性分析

对具有不同拉压模量的厚壁球壳,采用双剪统一强度理论推导了其扩张问题的应力及位移的统一解. 分析了不同模量、不同模型控制参数对厚壁球壳扩张时的扩张压力和应力场的影响.结果表明：厚壁球壳弹性极限压力、应力场、位移场等均随着模量控制参数、模型参数的变化而变化,在$\alpha<1$的情况下(即$E^ + < E^ - $),可以明显提高球壳的弹性极限压力$p_e $; 厚壁球壳塑性极限压力与材料的拉压模量无关,与模型参数$\eta$有关,且随$\eta$的增加,先增大后减小. 因此若采用经典的弹性理论和单一的模型参数对厚壁球壳进行设计计算,会带来较大的误差.     

厚板弹塑性混合状态Hamiltonian元的半解析法

采用初应力法并通过对增量形式Ｈｅｌｌｉｎｇｅｒ－Ｒｅｉｓｓｎｅｒ变分原理的修正，给出了厚板弹塑性分析的Ｈａｍｉｌｔｏｎ正则方程及其相应的半解析法。此法在厚板平面内采用通常的有限元离散，而沿板厚方向采用状态空间法给出半解析解答。     

层合厚板混合状态Hamiltonian动力元及其半解析解

文中在时间方向采用Ｌａｐｌａｃｅ变换，给出了层合厚板动力学问题混合状态Ｈａｍｉｌｔｏｎ正则方程及其半解析法．该方法在层板平面内采用通常的有限元离散，而沿板厚方向采用状态控制方程给出解析解答．在层与层之间采用迁移矩阵法，给出层合板上下表面力学量之间的关系式．利用打靶法得到响应在象空间的一般解，然后再利用拉氏逆变换的数值解求出层合板的瞬时位移场和应力场．     

四边简支各向同性矩形厚板的插值矩阵法解

针对强厚度矩形板四边简支情况,论文根据状态变量法思想,基于三维弹性理论基本方程,以3个位移分量及3个应力分量按双三角级数展开,将三维弹性力学控制方程转化为常微分方程边值问题.尽管一些各向异性弹性矩形厚板早已由状态空间法获得分析解,可是各向同性厚板的分析解至今难以获得,因为状态空间解法中特征方程有重根问题而不易于收敛.论文提出采用插值矩阵法直接对常微分方程进行求解,获得各向同性矩形厚板在四边简支边界条件下三维理论的位移和应力解,并与有限元精细结果进行比较,证明了本文解的准确性.     

电磁波导的半解析辛分析

根据电磁波导的Hamilton体系，辛几何可用于任意各向异性材料，而且便于处理不同区段的界面条件，横向的电场和磁场构成了对偶向量．基于Hamilton变分原理用半解析法进行横向离散应当保持体系的辛结构．离散后可以运用应用力学的有效算法，求解其辛本征值问题．每段波导可以引入两端Riccati矩阵，用精细积分法求解其方程组．