地基土—多自由度非线性结构相互作用体系的非平稳随机动力响应分析

本文应用一个改进的等效线性化方法探讨了地基土——多自由度非线性结构相互作用体系在随机地震荷载作用下的非平稳动力响应问题。对非线性恢复力模型采用一个具有三线性滞回特性的非线性系统。最后将在平稳Gauss过滤白噪声激励下的非平稳随机响应与Monte-carlo法的统计结果进行了比较,得出了较为满意的结果。     

求解结构动应力的振型函数差分法

本文利用常规有限元方法的计算结果,结合数值计算方法对振型函数进行［Ｌ］算子的微分计算,从而可方便迅速获得到复杂结构动应力响应,并对梁和板进行了计算,计算结果表明该方法具有较高的精度,较一般的动态有限元具有通用性强,计算简单等特点。     

Some fast 3D finite element solvers for the generalized Stokes problem

This paper is devoted to a comparison of various iterative solvers for the Stokes problem, based on the preconditioned Uzawa approach. In the first section the basic equations and general results of gradient-like methods are recalled. Then a new class of preconditioners, whose optimality will be shown, is introduced. In the last section numerical experiments and comparisons with multigrid methods prove the quality of these schemes, whose discretization is detailed.     

梁中复合应力波的传播

采用有限差分法讨论了梁中复合应力波的传播．给出了粘塑性悬臂梁当自由端受突加弯矩载荷作用时梁内复合应力波传播的基本图象。指出,在冲击早期响应阶段．截面横向转动惯性效应起着重要作用,是不可忽视的。标志弹塑性边界的塑性效,一开始由自由端向固定端运动,但在反射卸载波的迎面作用下,会出现回退现象。在波动早期阶段,固定端主要处于弹性变形状态。此外,还对弹塑性梁中复合应力波的控制方程进行了必要的讨论。     

粘弹性大变形动力响应的有限元分析

本文基于ＴｏｔａｌＬａｇｒａｎｇｉａｎ增量叠加方法,采用Ｋｉｒｃｈｈｏｆｆ应力增量和Ｇｒｅｅｎ应变增量表示的动力虚功方程和Ｋｉｒｃｈｈｏｆｆ应力－Ｇｒｅｅｎ应变的单积分型本构关系,导出粘弹性大变形的动力变分方程。依此采用Ｎｅｗｍａｒｋ法和八节点轴对称等参数元与二十节点三维等参数元编制了轴对称及三维问题的动力响应计算程序,典型例题的计算结果表明分析符合结构的物理性质。     

新型p型板单元及其在结构振动分析中的应用

利用Legendre正交多项式作为形函数基底函数,开发了两种新型的通用p型板单元.单元矩阵的解析积分保证了p型有限单元解的精确性及单调收敛性,计算实例表明所开发的p型有限单元计算结果随基底函数中附加项数量的增加而快速收敛,且它们的计算精度远高于一般线性单元.另外,p型板单元不使用缩减积分也能分析薄板的振动问题,利用它们收敛率高的特点,分析了结构破坏的时频特性.p型有限单元仿真结果与实测结果良好的吻合证明了它们用于结构振动响应分析的有效性.     

辛解析奇异元在动荷载作用下V型切口问题中的应用

利用辛解析奇异单元,结合时域精细算法,研究了动荷载作用下的含平面V型切口问题。时域上,采用时域精细算法,并结合自适应算法控制展开项数,保证了计算精度。空间域上,切口尖端附近采用辛解析奇异单元,其余区域采用常规有限单元,避免了局部网格加密。本文使用的辛解析奇异单元不需要过渡单元和局部网格加密,且能够通过奇异单元内部的参数关系直接给出切口尖端的应力强度因子,不需要复杂的后处理过程。数值结果表明,本文方法具有良好的精度和稳定性,可以准确地计算动态应力强度因子。     

三维五向编织复合材料渐进损伤分析及强度预测

基于五向编织复合材料连续体细观结构单胞,提出了材料的三维渐进损伤分析模型,采用非线性有限元方法并结合均匀化平均思想,建立了三维五向编织复合材料的强度预测胞元模型.经研究典型编织角材料在拉伸载荷作用下细观损伤的发生及演化过程,分析了材料的细观失效机理,获得了材料的宏观拉伸应力应变曲线和极限破坏强度,并详细探讨了工艺参数编织角对材料宏观力学性能的影响规律.     

Three-dimensional elliptic crack under impact loading

The dynamic stress intensity factor of a three-dimensional elliptic crack under impact loading is determined with the finite element method. The computation results can take into account the influence of time and the ratio of the wave speeds on the stress intensity factor. The present method is suitable not only for three-dimensional dynamic crack, but also for three-dimensional dynamic contact. Project supported by the National Natural Science Foundation of China (No. K19672007).