The nominal force method for truss geometry and topology optimization incorporating stability considerations

This paper presents a computationally efficient method for truss layout optimization with stability constraints. Previously proposed approaches that ensure stability of optimal frameworks are first reviewed, showing that existing studies are generally restricted to topology optimization. The present contribution aims to generalize the approach to simultaneous geometry and topology optimization. A lower-bound plastic design formulation under multiple loading will serve as basis for this purpose. The numerical difficulties associated with geometrical variations are identified and the parametrization is adapted accordingly. To avoid nodal instability, the nominal force method is adopted, which introduces artificial loading cases to simulate the effect of geometric imperfections. Hence, the truss systems with unstable nodes are eliminated from the set of optimal solutions. At the same time, the local stability of structural members is ensured via a consistent local buckling criterion. This novel formulation leads to optimal configurations that can be practically used for the preliminary design of structural frameworks. Four applications illustrate the impact of stability constraints on the solution. The importance of geometry optimization is also pointed out by comparing with results that would be unattainable by topology optimization only.     

Nonlinear buckling optimization of composite structures considering “worst” shape imperfections

Nonlinear buckling optimization is introduced as a method for doing laminate optimization on generalized composite shell structures exhibiting nonlinear behaviour where the objective is to maximize the buckling load. The method is based on geometrically nonlinear analyses and uses gradient information of the nonlinear buckling load in combination with mathematical programming to solve the problem. Thin-walled optimal laminated structures may have risk of a relatively high sensitivity to geometric imperfections. This is investigated by the concepts of “worst” imperfections and an optimization method to determine the “worst” shape imperfections is presented where the objective is to minimize the buckling load subject to imperfection amplitude constraints. The ability of the nonlinear buckling optimization formulation to solve the laminate problem and determine the “worst” shape imperfections is illustrated by several numerical examples of composite laminated structures and the application of both formulations gives useful insight into the interaction between laminate design and geometric imperfections.     

Effects of imperfections on the buckling response of compression-loaded composite shells

The results of an experimental and analytical study of the effects of initial imperfections on the buckling and postbuckling response of three unstiffened thin-walled compression-loaded graphite-epoxy cylindrical shells with different orthotropic and quasi-isotropic shell-wall laminates are presented. The results identify the effects of traditional and non-traditional initial imperfections on the non-linear response and buckling loads of the shells. The traditional imperfections include the geometric shell-wall mid-surface imperfections that are commonly discussed in the literature on thin shell buckling. The non-traditional imperfections include shell-wall thickness variations, local shell-wall ply-gaps associated with the fabrication process, shell-end geometric imperfections, non-uniform applied end loads, and variations in the boundary conditions including the effects of elastic boundary conditions. A high-fidelity non-linear shell analysis procedure that accurately accounts for the effects of these traditional and non-traditional imperfections on the non-linear responses and buckling loads of the shells is described. The analysis procedure includes a non-linear static analysis that predicts stable response characteristics of the shells and a non-linear transient analysis that predicts unstable response characteristics.     

Design of buckling constrained multiphase material structures using continuum topology optimization

This study contributes to a possibility of evaluating composite structures configuration such as steel and concrete using buckling and volume constraints based on multi-material topology optimization. A Jacobi active-phase algorithm is used to generate multiphase topology optimization. It provides a rational solution appropriated to the topology optimizer, Method of Moving Asymptotes due to the conflict in updating the design variables. A modified material interpolation scheme solving spurious buckling modes problem which occurs in the multi-material topology optimization process is given and discussed. An investigation of buckling constraint parameter is described. It allows a single-objective minimum compliance topology optimization to obtain two objectives of maximizing both structure stiffness and first buckling load factor. The optimal changing topologies of single material structure and multi-material structure corresponding to different buckling constraints are presented. Numerical examples of compression-only structures and compression-tension structures considering structural instability are performed using both single material and multiple materials to verify the efficiency and superiority of the present method.

     

Plate/shell structure topology optimization of orthotropic material for buckling problem based on independent continuous topological variables

The purpose of the present work is to study the buckling problem with plate/shell topology optimiza-tion of orthotropic material.A model of buckling topology optimization is established based on the independent,con-tinuous, and mapping method, which considers structural mass as objective and buckling critical loads as constraints. Firstly, composite exponential function (CEF) and power function(PF)as filter functions are introduced to recognize the element mass,the element stiffness matrix,and the ele-ment geometric stiffness matrix.The filter functions of the orthotropic material stiffness are deduced. Then these fil-ter functions are put into buckling topology optimization of a differential equation to analyze the design sensitiv-ity.Furthermore,the buckling constraints are approximately expressed as explicit functions with respect to the design vari-ables based on the first-order Taylor expansion.The objective function is standardized based on the second-order Taylor expansion. Therefore,the optimization model is translated into a quadratic program.Finally,the dual sequence quadratic programming(DSQP)algorithm and the global convergence method of moving asymptotes algorithm with two different filter functions(CEF and PF)are applied to solve the opti-mal model.Three numerical results show that DSQP&CEF has the best performance in the view of structural mass and discretion.     

基于激光测距仪的薄壳结构初始缺陷测量系统

薄壳结构的屈曲通常对初始缺陷十分敏感,因此对几何初始缺陷的精确测量是高质量的壳体屈曲试验的必要环节。本研究旨在开发一套基于激光测距仪的初始缺陷测量系统,以实现对壳体结构的快速、精确、非接触的测量。本文详细介绍该系统的设计思想及实现方法,包括激光测距仪的选择、转动与直线运动的实现、数据采集与计算机自动控制技术等。利用该系统对两个柱支承钢筒仓模型进行了仔细的初始缺陷三维测量,验证了系统的适用性和可靠性。文中还提出了利用二重傅立叶分解技术分析实测缺陷的方法。本文成果与方法为薄壳结构屈曲试验奠定了基础。     

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

通过对拱顶储罐罐壁承受轴向载荷、初始几何缺陷及轴压失稳状况研究,指 出在固定顶罐设计、建造和运行各阶段都应进行罐壁轴压稳定性校核. 根据圆柱薄壳稳定性 理论和轴压失稳临界应力数值分析计算结果,提出固定顶罐罐壁轴压稳定性校核方法和数学 模型,并运用回归分析方法建立罐壁轴压失稳临界应力计算公式. 对几种常用规格的拱顶罐 有初始挠度缺陷罐壁轴压稳定性分析表明：随储罐容积和罐壁初始挠度增大,罐壁轴压稳定 性呈减弱趋势.     

初始几何缺陷对网壳截面优化结果影响研究

单层网壳结构是缺陷敏感结构,初始几何缺陷不同对"网壳截面优化结果"有显著影响,该文研究了网壳杆件截面取值与缺陷分布、结构极限承载力之间的关系.结果表明:经截面优化设计的网壳结构最不利应力最早并且一直出现在初始几何缺陷较大处;优化寻优结果与初始几何缺陷分布有直接关系,增大缺陷较大处杆件截面能显著提高结构承载力;由于缺陷随机分布,经截面优化设计的网壳结构须校核其在不同缺陷下的稳定承载力以确保结构安全.     

Shape optimization and its extension to topological design based on isogeometric analysis

In most of structural optimization approaches, finite element method (FEM) has been employed for structural response analysis and sensitivity calculation. However, the approaches generally suffer certain drawbacks. In shape optimization, cumbersome parameterization of design domain is required and time consuming remeshing task is also necessary. In topology optimization, design results are generally restricted on the initial design space and additional post-processing is required for communication with CAD systems. These drawbacks are due to the use of different mathematical languages in design or geometric modeling and numerical analysis: spline basis functions are used in design and geometric modeling whereas Lagrangian and Hermitian polynomials in analysis. Isogeometric analysis is a very attractive and promising alternative to overcome the limitations resulting from the use of the conventional FEM in structural optimization. In isogeometric analysis, the same spline information such as control points and spline basis functions which represent geometries in CAD systems are also used in numerical analysis. Such unification of the mathematical languages in CAD, analysis and design optimization can resolve the issues mentioned above. In this work, structural shape optimization using isogeometric analysis is studied on 2D and shell problems. The proposed framework is extended to topology optimization using trimming techniques. New inner fronts are introduced by trimming spline curves in topology optimization. Trimmed surface analysis which was recently proposed to analyze arbitrary complex topology problems is employed for topology optimization. Some benchmarking problems in shape and topology optimization are treated using the proposed approach.     

薄壁圆筒壳初始几何缺陷不确定性量化的极大熵方法

针对薄壁圆筒壳结构轴压屈曲载荷的缺陷敏感性以及真实几何缺陷的不确定性,提出一种基于实测缺陷数据和极大熵原理的初始缺陷建模与屈曲载荷预测方法。首先,将初始几何缺陷视为二维随机场,并利用实测缺陷数据和Karhunen-Loève展开法将初始缺陷的随机场建模转化为随机向量的建模;其次,利用极大熵方法确定随机向量的概率分布;最后,基于所构建的初始缺陷随机模型,利用MCMC抽样方法和确定性屈曲分析方法,进行随机屈曲分析并给出基于可靠度的屈曲载荷折减因子。数值算例表明,与直接假设随机场相关结构的方法相比,本文方法的结果是对薄壁圆筒壳屈曲载荷的一个更无偏估计。     

Imperfection sensitivity of curved panels under combined compression and shear

The initial buckling load of curved panels under compressive loads is substantially reduced by the existence of imperfections, in particular geometric imperfections. It is therefore essential that these imperfections are considered in analysing components which incorporate such panels in order to accurately predict their buckling behaviour. Finite element analysis allows fully non-linear analysis of shells containing geometric imperfections, however, to obtain accurate results information is required on the exact size and shape of the imperfection to be modelled. In most cases this data is not available. It is therefore generally recommended that the imperfections are modelled on the first eigenmode and have an amplitude selected according to the manufacturing procedure. This paper presents the effects of varying imperfection shape and amplitude on the buckling and postbuckling behaviour of one specific case, a curved panel under combined shear and compression, to test the accuracy of such recommendations.     

基于功能度量法的桁架结构非概率可靠性拓扑优化方法研究

在结构优化中,拓扑优化相比于尺寸优化和形状优化,设计空间更加广泛,因而能够取得更大的效益.近年来,结构拓扑优化逐渐成为人们研究的热点和难点.随着科学技术的发展,工程结构越来越复杂,材料本身和外部环境的不确定性影响加剧,因此在拓扑优化中需要考虑不确定性的影响.本文研究了桁架结构的非概率可靠性拓扑优化问题,用区间模型来量化不确定性,并利用参数顶点组合法来完成不确定性的传播分析,利用基于面积比的非概率可靠性指标构建可靠性拓扑优化模型,提出了功能度量法对原可靠性约束进行等价转化,从而克服了收敛性问题.采用移动渐近方法(MMA)对优化问题进行了求解.数值算例表明,本文提出的功能度量法能够很好地适用于桁架结构的非概率可靠性拓扑优化问题.     

Repeated buckling and its influence on the geometrical imperfections of stiffened cylindrical shells under combined loading

The present experimental study aims at providing better inputs for improvement of the buckling load predictions of stiffened cylindrical shells subjected to combined loading. The work focuses on two main factors which considerably affect the combined buckling load of stiffened shells, namely geometric imperfections and boundary conditions. Six shells with nominal simple supports were tested under various combinations of axial compression and external pressure. The vibration correlation technique is employed to define the real boundary conditions. The geometric imperfections of the integrally stiffened shells are measured in the present experiments in situ and are used as inputs to a multimode analysis which yields the corresponding “knockdown” factor for various combinations of loading. Thus, when employing the repeated buckling procedure for obtaining interaction curves, each point on the curve is adjusted (using the multimode analysis) for the measured “new” surface of the shell and this results in more realistic interaction curves. The geometrical imperfections of the preloaded shells can also serve as an input to the International Imperfection Data Bank for future studies on the correlation between the manufacturing method of the shell and their geometric imperfections.     

Postbuckling of pressure-loaded shear deformable laminated cylindrical panels

IntroductionInrecentyears,fiber_reinforcedcompositelaminatedpanelshavebeenwidelyusedintheaerospace,marine ,automobileandotherengineeringindustries .Theproblemofbucklingandpostbucklingofcylindricalpanelsunderaxialcompressionortorsionhasbeenextensivelystudied .Incontrast,theliteratureoncylindricalpanelsunderpressureloadingisrelativelyspares.Thesestudiesincludealinearbucklinganalysis (Singeretal.[1]) ,anonlinearbucklinganalysi(YamadaandCroll[2 ]) ,anelastoplasticbucklinganalysisusingreducedstif…     

Nonlinear response of an imperfect microcantilever static and dynamically actuated considering uncertainties and noise

The motion of a slender, clamped-free, imperfect, electrically actuated microbeam is investigated. Special attention is given to the influence of imperfections and noise on the bifurcations and instabilities of the structure, a problem not tackled in the previous literature on the subject. To this end, a geometrically nonlinear theory is adopted for the microbeam retaining geometric nonlinear terms up to the third order and considering in a consistent way the effect of initial geometric imperfections. Also, additive white noise is considered to model forcing uncertainties, and the Galerkin discretization method, using as interpolating functions the linear vibration modes, is used to obtain a modal stochastic differential equation of Itô type, which is solved by the stochastic Runge–Kutta method. A parametric analysis clarifies the influence of geometric imperfections and noise level on the natural frequencies, resonance curves, and pull-in instability. Additionally, the global dynamics is examined through the generalized cell mapping, showing the effects of uncertainties on the attractor’s probability density functions and basins of attraction.

     

含初缺陷裂纹损伤梁的冲击动力屈曲

由Hamilton原理导出考虑初始缺陷及横向剪切变形时裂纹梁的动力屈曲控制方程;应用断裂力学中常用的线弹簧模型将裂纹引入到屈曲控制方程中;基于B-R动力屈曲判断准则,采用数值方法求解了受轴向冲击载荷作用时裂纹梁的动力屈曲;对比讨论了不同冲击速度、初始几何缺陷大小以及分布形式等因素对梁冲击动力屈曲的影响。     

加筋板的非线性相关屈曲研究

本文将线性样条有限条理论推广到大挠度弹塑性范围,建立了样条有限条非线性分析方法,成功地分析了纵向加筋板结构局部与整体稳定相关作用对其极限承载能力的影响.采用弧长法和Newton-Raphson迭代求解非线性刚度方程,可获得包括峰值点在内的一条完整的荷载—挠度曲线.分析中考虑了加筋板结构的初曲和残余应力的影响,使得计算结果更具实用价值.     

Effect of multiple localized geometric imperfections on stability of thin axisymmetric cylindrical shells under axial compression

Stability of imperfect elastic cylindrical shells which are subjected to uniform axial compression is analyzed by using the finite element method. Multiple interacting localized axisymmetric initial geometric imperfections, having either triangular or wavelet shapes, were considered. The effect of a single localized geometric imperfection was analyzed in order to assess the most adverse configuration in terms of shell aspect ratios. Then two or three geometric imperfections of a given shape and which were uniformly distributed along the shell length were introduced to quantify their global effect on the shell buckling strength. It was shown that with two or three interacting geometric imperfections further reduction of the buckling load is obtained. In the ranges of parameters that were investigated, the imperfection wavelength was found to be the major factor influencing shell stability; it is followed by the imperfection amplitude, then by the interval distance separating the localized imperfections. In a wide range of parameters this last factor was recognized to have almost no effect on buckling stresses.     

含脱层损伤复合材料层合梁的冲击动力屈曲

针对含初始缺陷和脱层损伤的复合材料层合梁的轴向冲击动力屈曲问题进行了分析。基于Hamilton原理导出了考虑初始缺陷、轴向和横向惯性、横向剪切变形以及转动惯性影响时含脱层损伤复合材料梁的非线性动力屈曲控制方程;基于B-R准则,采用有限差分方法求解了受轴向冲击载荷作用下含脱层损伤复合材料梁的动力屈曲问题;讨论了冲击速度、初始几何缺陷、铺层角度以及脱层长度等因素对复合材料层合梁动力屈曲的影响。     

Kármán-type equations for a higher-order shear deformation plate theory and its use in the thermal postbuckling analysis

Kármán-type nonlinear large deflection equations are derived occnrding to the Reddy’s higher-order shear deformation plate theory and used in the thermal postbuckling analysis The effects of initial geometric imperfections of the plate areincluded in the present study which also includes th thermal effects.Simply supported,symmetric cross-ply laminated plates subjected to uniform or nomuniform parabolictemperature distribution are considered. The analysis uses a mixed GalerkinGolerkinperlurbation technique to determine thermal buckling louds and postbucklingequilibrium paths.The effects played by transverse shear deformation plate aspeclraio, total number of plies thermal load ratio and initial geometric imperfections arealso studied.