Sensitivity analysis for 3-D elastoplastic structural reliability

By taking the elastoplastic effect of structural material into account and based on 3-D elastoplastic stochastic finite element method, methods for sensitivity analysis with respect to both the distribution parameters of random variables and parameters in the limit state function are suggested. In the incremental iterative calculation, the sub-increment changingK, method and the corresponding formulas for accelerating convergence are used. The sensitivity of 3-D structural system reliability with respect to random variables is also studied.     

一种基于应力的双方向结构拓扑优化算法

基于传统的渐进结构优化方法,提出了一种基于应力的双方向渐进拓扑优化算法。该方法是对传统方法和目前的双方向法的改进和完善。算例表明该方法能避免目前的双方向法具有的解的振荡问题,具有更高的可靠性,能获得更佳的拓扑结构。     

一类有限元结构分析的EBE计算方法

本文利用ＥＢＥ策略￣［１］的基本思想，给出一类有限元结构分析的ＥＢＥ计算方法，即ＥＢＥ共轭梯度法ＥＢＥ－ＣＧ和ＥＢＥ预处理共轭梯度法ＥＢＥ－ＰＣＧ，这类方法避免了传统有限元结构分析中总刚度阵的组集而可大大降低存储量要求。同时它们还特别适合在各种粒度下的多处理机系统上实现。初步数值试验结果表明：这类ＥＢＥ有限元结构分析方法对串行和并行计算都是很有效的。     

基于虚荷载变量的形状优化和灵敏度分析

基于选择施加在结构“控制点”上的虚荷载作为优化设计变量,针对一种新的承受约束的形状优化数值方法进行了研究。借助于节点位移与虚荷载之间的线性关系,提出了一种新的计算灵敏度系数的解析方法。利用节点移动速度域概念构造了优化新形状产生的计算公式,以结构中节点的最大应力最小化作为优化目标,通过控制网格结点的最大位移量,较好地解决了单元网格在形状优化中的扭曲问题。对三个不同的实例成功地完成了形状优化。     

基于Monte-Carlo方法的结构系统可靠度计算及敏度分析

基于可靠性分析理论,将结构失效概率对随机变量均值的敏度表示成失效概率与正则化随机变量在失效域上期望的乘积,并利用敏度分析的结果给出了结构线性等效安全余量的表达式.通过等概率转换,使得该方法可以应用于服从任意分布的随机变量.该方法在给出失效概率的同时,能够给出失效概率对随机变量均值的敏度,而无需重新对结构进行计算,提高了...     

基于随机有限元的非线性结构稳健性优化设计

结合结构优化技术和摄动随机有限元方法研究了非线性结构稳健设计问题。将结构稳健性优化设计问题构造为双目标优化问题。优化目标包含结构性能函数的期望值和标准差。约束函数的变异也给予考虑,并采用基于函数梯度的算法进行求解。为对具有路径相关特征的非线性结构性能及结构响应的平均值及标准差进行分析。本文采用缩减的随机变量,提出了基于增量法的摄动随机有限元计算格式。在此框架下,进一步提出以一般泛函形式表达的结构性能的平均值和方差及其灵敏度的计算格式。为显示方法的有效性。文中给出几个数值算例。     

一种精确结构静力重分析法

本文提出了一种结构静力重分析方法。通过引入结构刚体位移特征向量，可以导出结构广义柔度矩阵，原阶数较高的刚度方程被转化成一阶数较小的线性系统，位移一般解可以在边界条件尚未引入结构刚度矩阵之前导出，对于有局部变化的结构，新的结构广义柔度矩阵可以迅速进行修改。这种静力重分析可以用在载荷条件、边界条件、结构单元同时或分别改变时的静力分析之中，文中提供了两个算例，以证明此方法的有效性     

超级有限元法及其在结构工程中的应用

本文探讨一种基于半连续半离散思想，适用于复杂结构（如高层框架、剪力墙、桁架、网架等结构系统）工程分析的超级有限元，其结构数值分析是按连续体进行，但又按单个构件进行有限元计算。这种按整体系统进行离散所获得的单元内部包含众多构件，有别于一般常见的实体有限元，称为“超级有限元”。这种方法自由度数比一般有限元法少很多，又与单元内部所含构件数多少无关，并可求取结构内每个构件的内力值。     

非确定性结构静动态特性稳健优化设计

本文研究了考虑参数随机性的结构静动态特性稳健性优化设计问题的数学模型和数值求解。在考虑结构设计变量和其研究了考虑参数随机性的结构静动态特性稳健性优化设计问题的数学模型和数值求 解. 在考虑结构设计变量和其他参数随机分布的前二阶矩的条件下,采用基于二阶摄动法的 随机有限元方法对结构响应的平均值和方差进行近似求解. 在摄动法有限元分析的框架下, 提出以一般函数形式表达的结构性能的平均值和标准差及其灵敏度的计算格式. 将结构 稳健性优化设计问题构造为双目标优化问题,优化目标包含结构性能函数的期望值和标准 差,约束函数的变异也给予考虑. 优化问题采用基于函数梯度的算法进行求解. 文中给出的数值算例显示了方法的有效性.     

Sensitivity analysis of soil parameters based on interval

Interval analysis is a new uncertainty analysis method for engineering structures. In this paper, a new sensitivity analysis method is presented by introducing interval analysis which can expand applications of the interval analysis method. The interval analysis process of sensitivity factor matrix of soil parameters is given. A method of parameter intervals and decision-making target intervals is given according to the interval analysis method. With FEM, secondary developments are done for Marc and the Duncan-Chang nonlinear elastic model. Mutual transfer between FORTRAN and Marc is implemented. With practial examples, rationality and feasibility are validated. Comparison is made with some published results.     

Modification of an implicit approach based on nonstandard rules for structural dynamics analysis

The present study aims to modify a recently suggested implicit approach consisted of the approximate Euler method and closed-form exponential mapping (herein referred to as the Liu scheme) for the dynamic analysis of structures. Such modification has been developed based upon nonstandard rules. The equation of motion is formulated in the augmented dynamic space to apply the exponential mapping as a group preserving scheme. The formulation of the proposed method involves the hyperbolic sine and cosine functions. The method is therefore prone to divergence due to the behavior of the hyperbolic functions in structures with a high ratio of stiffness to mass. In the present study, to consider the properties of the structural equation into the formulation of the time step size and thereby avoid the divergence, a parameter, known as stability parameter, is thus derived from the exact solution of the equation of motion based on nonstandard rules. Embedding this parameter into the proposed method improves its stability. Afterward, for evaluating the performance of the proposed method, it is applied to several structures with different loading patterns while implemented in programing environment of the Matlab software. The results are compared to those of several commonly used numerical methods in structural applications. It is found that the proposed method has acceptable convergence and accuracy, and low time consumption compared to several commonly used methods. Furthermore, its stability is guaranteed by embedding the stability parameter into the proposed method.     

非确定结构系统区间分析的泛灰求解方法

工程中的不确定性问题可以用区间分析、概率理论或模糊理论来求解。采用泛灰区间分析法来处理结构静力分析和设计中的不确定性问题。将结构系统中的不确定性参数用区间数来表示，用有限元法建立系统的控制方程。该控制方程是线性区间方程组。然后，在概述泛灰数的概念及其运算规则的基础上，介绍了泛灰数与区间数的转化，利用泛灰数的可扩展性对区间进行分析，研究了泛灰线性方程求解，然后将它应用于结构静力分析和设计中的不确定性问题，泛灰数不仅具有区间分析的功能，而且能解决区间分析所不能解决的问题。文中给出了两个算例，列出了本文算法与其他算法的结果比较。     

Homotopy solution of the inverse generalized eigenvalue problems in structural dynamics

The structural dynamics problems, such as structural design, parameter identification and model correction, are considered as a kind of the inverse generalized eigenvalue problems mathematically. The inverse eigenvalue problems are nonlinear. In general, they could be transformed into nonlinear equations to solve. The structural dynamics inverse problems were treated as quasi multiplicative inverse eigenalue problems which were solved by homotopy method for nonlinear equations. This method had no requirements for initial value essentially because of the homotopy path to solution. Numerical examples were presented to illustrate the homotopy method.     

结构有限元分析中阶跃型弹性约束的一种有效处理方法

本文通过引入弹性约束刚度矩阵和结构位移约束列阵，提出了结构有限元分析中处理阶跃型弹性约束的一种有效方法。该法通过改变弹性约束系数及位移非约束量大小，可方便有效地处理结点常弹性约束，刚性约束，阶跃型弹性的约束及阶跃型刚性约束等问题。     

A mathematical programming algorithm for limit analysis

This paper deals with the limit analyses of perfect rigid-plastic continua. Based on the kinematic theorem of the limit analysis theory, a mathematical programming finite element formula for determining the upper bound load multiplier has been established, and an iteration algorithm proposed accordingly. In this algorithm the plastic and rigid zones are distinguished for every iteration step, and the goal function is modified gradually. The difficulties caused by the nonsmoothness of the goal function are overcome. Some examples solved by this algorithm are presented. The project supported by National Natural Science Foundation of China.     

基于等几何分析的比例边界有限元方法

提出了一种具有比例边界有限元的半解析特性和等几何分析的几何特性的新方法。该新方法是在比例边界有限元框架中用NURBS曲线或曲面精确描述域边界几何形状,同时域边界位移场采用描述几何形状的NURBS形函数等参构造。这种新方法具有比例边界有限元固有的径向解析特性和NURBS的高阶连续性的优点。数值算例显示,与传统的比例边界有限元相比,基于等几何分析的比例边界有限元方法提高了域边界单元和域内应力场的连续性,减少了计算自由度。应用此方法可以用较少的计算自由度获得更高连续阶和更高精度的位移、应力和应变场。     

基于整体-局部技术的结构形状优化策略

针对大尺度结构局部形状优化设计中的求解效率与精度之间的矛盾,提出了一种基于整体-局部技术的结构形状优化方法.首先,分析并讨论了是否考虑对切割边界的边界条件进行修正的两种优化方案的优缺点.然后,提出了一种将整体-局部技术与无梯度优化法相结合的双循环优化程序.最后,通过一系列实例对双循环优化程序的实用性、效率和结果精度进行了验证分析.结果 表明,提出的基于整体-局部技术的结构形状优化策略对解决大尺度结构局部形状优化问题,无论是在计算效率还是在结果精度方面,都具有很好的效果.     

结构动力学方程的二次加速度积分法

本文提出了结构动力学方程求解的一类二次加速度逐步积分法，推导了计算公式，分析了积分稳定性和精度。通过理论分析和具体算例表明，这种方法具有相当高的积分精度，但积分是条件稳定的。     

间隙元在钻柱瞬态动力学分析中的应用

为了描述钻柱与井壁的随机碰撞接触状态，本文构造了动力间隙元，并推导了动力间隙元的相对压缩量和迭代求解格式。把动力间隙元与梁单元相结合，采用Newmark直接积分法和冲量定理，进行了钻柱碰撞接触非线性瞬态动力学分析。经工程应用表明，构造的动力间隙元能够描述钻柱的碰撞接触状态，所设计的偏心防斜钻具已在大庆油田得到应用，取得了明显的经济和社会效益。     

预锻模具形状优化设计与有限元灵敏度分析

采用刚－粘塑性有限元灵敏度分析方法研究锻造过程的预锻模具形状优化设计问题。表示预锻模具形状的三次Ｂ样条曲线的控制系数用作设计变量。以实际终锻件与理想终锻件的形状差异为目标函数,给出了与设计变量有关的目标函数,节点坐标和节点速度等方面的灵敏度及其数学关系,对于典型的轴对称锻造过程,优化设计的预锻模具形状可获得理想的的终锻件形状,为实现净成形锻造提供了一种有效方法。