薄壁杆系结构的梁元分析模型

本文导出了用于薄壁杆系结构弹性分析的薄壁梁元分析模型，在空间梁元分析模型^[3]的基础上，采用了一种改进的位移模式，考察了薄壁杆件可能发生的拉压，剪切，弯曲，扭转和翘曲等各变形形式以及它们的耦合效应，得出了相应的单元形函数，同时从工程应变的定义出发，采用Taylor级数展开的方法，建立了单元的五阶近似正交变表达式，并建立了相应的薄壁单元刚度方程，从而得出了局部坐标系下单元刚度矩阵的显式，根据本文所导出的薄壁梁元分析模型，编制了相应的结构计算程序，通过算例验证了本文所推导的单元刚度矩阵，同时通过与传统空间梁元计算模型计算结果的比较，阐述了截面翘曲对薄壁杆系结构的影响。     

An asymptotically correct model for initially curved and twisted thin-walled composite beams

The variational–asymptotic method has been applied to develop an asymptotically correct model for initially curved and twisted, thin-walled, composite beams of arbitrary cross-sectional shapes and arbitrary anisotropic materials. In a two-step asymptotic reduction procedure, the three-dimensional strain energy is asymptotically reduced first to a two-dimensional shell strain energy and then to a one-dimensional beam strain energy. This is a new attempt where initially curved and twisted, thin-walled, composite beams, with open or closed sections, have been modeled in an asymptotically correct unified framework.     

Asymptotic theory for static behavior of elastic anisotropic I-beams

End effects for prismatic anisotropic beams with thin-walled, open cross-sections are analyzed by the variational-asymptotic method. The decay rates for disturbances at the ends of prismatic beams are evaluated, and the most influential end disturbances are incorporated into a refined beam theory. Thus, the foundations of Vlasovs theory, as well as restrictions on its applicability, are obtained from the variational-asymptotic point of view. Vlasovs theory is proved to be asymptotically correct for isotropic I-beams. The asymptotically correct generalization of Vlasovs theory for static behavior of anisotropic beams is presented. In light of this development, various published generalizations of Vlasovs theory for thin-walled anisotropic beams are discussed. Comparisons with a numerical 3-D analysis are provided, showing that the present approach gives the closest agreement of all published theories. The procedure can be applied to any thin-walled beam with open cross-sections.     

A non-linear viscoelastic model for ceramics at high temperatures

High-temperature mechanical behavior of ceramics is characterized by non-linear rate dependent responses, asymmetric behavior in tension and compression, and nucleation and coalescence of voids leading to rupture. Moreover, rupture experiments show considerable scatter or randomness in fatigue lives of nominally equal specimens. To capture the non-linear, asymmetric time-dependent behavior, a new non-linear viscoelastic model is proposed. Non-linearity and asymmetry are introduced in the volumetric component. To model the random formation and coalescence of voids, each element is assigned a failure strain sampled from a lognormal distribution. An element is deleted when its volumetric strain exceeds its failure strain. Temporal increases in strains produce a sequential loss of elements (a model for void nucleation and growth), which in turn leads to failure. Non-linear viscoelastic model parameters are determined from uniaxial tensile and compressive creep experiments on silicon nitride. The model is then used to predict the deformation of four-point bending and ball-on-ring specimens. Simulation is used to predict statistical moments of rupture lives. Numerical simulation results compare well with results of four-point bending experiments.     

Geometrically non-linear dynamics of composite four-bar mechanisms

This work intends to demonstrate the importance of a geometrically nonlinear cross-sectional analysis of certain composite beam-based four-bar mechanisms in predicting system dynamic characteristics. All component bars of the mechanism are made of fiber reinforced laminates and have thin rectangular cross-sections. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. They are linked to each other by means of revolute joints. We restrict ourselves to linear materials with small strains within each elastic body (beam). Each component of the mechanism is modeled as a beam based on geometrically non-linear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and non-linear 1-D analyses along the three beam reference curves. For the thin rectangular cross-sections considered here, the 2-D cross-sectional non-linearity is also overwhelming. This can be perceived from the fact that such sections constitute a limiting case between thin-walled open and closed sections, thus inviting the non-linear phenomena observed in both. The strong elastic couplings of anisotropic composite laminates complicate the model further. However, a powerful mathematical tool called the Variational Asymptotic Method (VAM) not only enables such a dimensional reduction, but also provides asymptotically correct analytical solutions to the non-linear cross-sectional analysis. Such closed-form solutions are used here in conjunction with numerical techniques for the rest of the problem to predict multi-body dynamic responses more quickly and accurately than would otherwise be possible. The analysis methodology can be viewed as a three-step procedure: First, the cross-sectional properties of each bar of the mechanism is determined analytically based on an asymptotic procedure, starting from Classical Laminated Shell Theory (CLST) and taking advantage of its thin strip geometry. Second, the dynamic response of the non-linear, flexible four-bar mechanism is simulated by treating each bar as a 1-D beam, discretized using finite elements, and employing energy-preserving and -decaying time integration schemes for unconditional stability. Finally, local 3-D deformations and stresses in the entire system are recovered, based on the 1-D responses predicted in the previous step. With the model, tools and procedure in place, we identify and investigate a few four-bar mechanism problems where the cross-sectional non-linearities are significant in predicting better and critical system dynamic characteristics. This is carried out by varying stacking sequences (i.e. the arrangement of ply orientations within a laminate) and material properties, and speculating on the dominating diagonal and coupling terms in the closed-form non-linear beam stiffness matrix. A numerical example is presented which illustrates the importance of 2-D cross-sectional non-linearities and the behavior of the system is also observed by using commercial software (I-DEAS + NASTRAN + ADAMS).     

A new beam element for analyzing geometrical and physical nonlinearity

Based on Timoshenko＇s beam theory and Vlasov＇s thin-walled member theory, a new model of spatial thin-walled beam element is developed for analyzing geometrical and physical nonlinearity, which incorporates an interior node and independent interpolations of bending angles and warp and takes diversified factors into consideration, such as traverse shear deformation, torsional shear deformation and their coupling, coupling of flexure and torsion, and the second shear stress. The geometrical nonlinear strain is formulated in updated Lagarange （UL） and the corresponding stiffness matrix is derived. The perfectly plastic model is used to account for physical nonlinearity, and the yield rule of von Mises and incremental relationship of Prandtle-Reuss are adopted. Elastoplastic stiffness matrix is obtained by numerical integration based on the finite segment method, and a finite element program is compiled. Numerical examples manifest that the proposed model is accurate and feasible in the analysis of thin-walled structures.     

A simple meso-macro model based on SQP for the non-linear analysis of masonry double curvature structures

A 3D model for the evaluation of the non-linear behavior of masonry double curvature structures is presented. In the model, the heterogeneous assemblage of blocks is substituted with a macroscopically equivalent homogeneous non-linear material. At the meso-scale, a curved running bond representative element of volume (REV) constituted by a central block interconnected with its six neighbors is discretized through of a few six-noded rigid wedge elements and rectangular interfaces. Non linearity is concentrated exclusively on joints reduced to interface, exhibiting a frictional behavior with limited tensile and compressive strength with softening. The macroscopic homogenous masonry behavior is then evaluated on the REV imposing separately increasing internal actions (in-plane membrane actions, meridian and parallel bending, torsion and out-of-plane shear). This simplified approach allows to estimate heuristically the macroscopic stress–strain behavior of masonry at the meso-scale. The non-linear behavior so obtained is then implemented at a structural level in a novel FE non-linear code, relying on an assemblage of rigid infinitely resistant six-noded wedge elements and non-linear interfaces, exhibiting deterioration of the mechanical properties.Several numerical examples are analyzed, consisting of two different typologies of masonry arches (a parabolic vault and an arch in a so-called “skew” disposition), a ribbed cross vault, a hemispherical dome and a cloister vault. To fully assess numerical results, additional non-linear FE analyses are presented. In particular, a simplified model is proposed, which relies in performing at a structural level a preliminary limit analysis – which allows to identify the failure mechanism – and subsequently in modeling masonry through elastic elements and non-linear interfaces placed only in correspondence or near the failure mechanism provided by limit analysis. Simulations performed through an equivalent macroscopic material with orthotropic behavior and possible softening are also presented, along with existing experimental evidences (where available), in order to have a full insight into the capabilities and limitations of the approach proposed.     

泡沫铝填充薄壁方形铝管的静态弯曲崩毁行为

研究了泡沫铝填充的方形铝管的准静态三点弯曲行为.实验表明:泡沫铝填充有效地改变了铝管的局部崩毁变形模式;界面粘接提高了填充结构的抗弯刚度,但使结构易在较小转角下发生破坏.最后,基于实验提出了一个分析填充结构弯曲崩塌行为的理论方法,在小转角下给出了与实验相当吻合的结果.     

Length dependence of critical measures in single-walled carbon nanotubes

This paper presents an investigation on the buckling behaviour of single-walled carbon nanotubes under various loading conditions (compression, bending and torsion) and unveils several aspects concerning the dependence of critical measures (axial strain, bending curvature and twisting angle) on the nanotube length. The buckling results are obtained by means of an atomistic-scale generalized beam theory (GBT) that incorporates local deformation of the nanotube cross-section by means of independent and orthogonal deformation modes. Moreover, some estimates are also obtained by means of non-linear shell finite element analyses using Abaqus code. After classifying the buckling modes of thin-walled tubes (global, local and distortional), the paper addresses the importance of the two-wave distortional mode (flattening or ovalization mode) in their structural behaviour. Then, the well known expression to determine the critical strain of compressed nanotubes, which is based on Donnell theory for shallow shells, is shown to be inadequate for moderately long tubes due to warping displacements appearing in the distortional buckling modes. After that, an in-depth study on the buckling behaviour of nanotubes under compression, bending and torsion is presented. The variation of the critical kinematic measures (axial strain, bending curvature and twisting angle) with the tube length is thoroughly investigated. Concerning this dependence, some uncertainties that exist in the specific literature are meticulously explained, a few useful expressions to determine critical measures of nanotubes are proposed and the results are compared with available data collected from several published works (most of them, obtained from molecular dynamics simulations).     

钢-混凝土组合梁负弯矩区畸变屈曲弯矩计算公式

侧向弯曲屈曲及侧向弯扭屈曲均为钢-混凝土组合梁负弯矩区钢梁的重要屈曲模式,而现有计算方法通常只考虑到侧向弯曲屈曲,未考虑到侧向弯扭屈曲,因此有一定局限性。本文在钢梁腹板对钢梁下翼缘的转动约束刚度及侧向约束刚度的计算公式上,结合弹性介质中薄壁杆件的屈曲理论推导了工形钢-混凝土组合梁负弯矩区的侧向弯曲屈曲弯矩及侧向弯扭屈曲弯矩计算公式。实例分析表明,现有计算方法均存在一定理论缺陷,本文计算方法更为合理;同时,本文计算式比现有文献中同类型计算式更为简洁,便于手算,实用性强并适于工程应用。     

GEOMETRICALLY NONLINEAR FINITE ELEMENT MODEL OF SPATIAL THIN-WALLED BEAMS WITH GENERAL OPEN CROSS SECTION

Based on the theory of Timoshenko and thin-walled beams, a new finite element model of spatial thin-walled beams with general open cross sections is presented in the paper, in which several factors are included such as lateral shear deformation, warp generated by nonuni- form torsion and second-order shear stress, coupling of flexure and torsion, and large displacement with small strain. With an additional internal node in the element, the element stiffness matrix is deduced by incremental virtual work in updated Lagrangian （UL） formulation. Numerical examples demonstrate that the presented model well describes the geometrically nonlinear property of spatial thin-walled beams.     

缩径焊接空心球压弯下节点刚度研究

以缩径焊接空心球节点为研究对象,运用正交设计法设计了 9个试件,试验研究了压弯受力状态下节点的名义抗弯刚度.利用ABAQUS软件建立缩径焊接空心球节点分析模型,采用线性强化应力-应变关系,考虑几何和材料非线性的影响,对压弯受力状态下的缩径焊接空心球节点的抗弯刚度进行了有限元分析,绘制了节点弯矩-转角曲线,得到了节点名义...     

Theoretical prediction of stiffness and strength of three-dimensional and four-directional braided composites

Based on unit cell model, the 3D 4-directional braided composites can be simplified as unidirectional composites with different local axial coordinate system and the compliance matrix of unidirectional composites can be defined utilizing the bridge model. The total stiffness matrix of braided composites can be obtained by the volume average stiffness of unidirectional composites with different local axial coordinate system and the engineering elastic constants of braided composites were computed further. Based on the iso-strain assumption and the bridge model, the stress distribution of fiber bundle and matrix of different unidirectional composites can be determined and the tensile strength of 3D 4-directional braided composites was predicted by means of the Hoffman＇s failure criterion for the fiber bundle and Mises＇ failure criterion for the matrix.     

Non-linear buckling and postbuckling behavior of thin-walled beams considering shear deformation

The static stability of thin-walled composite beams, considering shear deformation and geometrical non-linear coupling, subjected to transverse external force has been investigated in this paper. The theory is formulated in the context of large displacements and rotations, through the adoption of a shear deformable displacement field (accounting for bending and warping shear) considering moderate bending rotations and large twist. This non-linear formulation is used for analyzing the prebuckling and postbuckling behavior of simply supported, cantilever and fixed-end beams subjected to different load condition. Ritz's method is applied in order to discretize the non-linear differential system and the resultant algebraic equations are solved by means of an incremental Newton-Rapshon method. The numerical results show that the beam loses its stability through a stable symmetric bifurcation point and the postbuckling strength is in relation with the buckling load value. Classical predictions of lateral buckling are conservative when the prebuckling displacements are not negligible and the non-linear buckling analysis is required for reliable solutions. The analysis is supplemented by investigating the effects of the variation of load height parameter. In addition, the critical load values and postbuckling response obtained with the present beam model are compared with the results obtained with a shell finite element model (Abaqus).     

ANCF索梁单元应变耦合问题与模型解耦

索粱结构在土木工程、航空航天等领域有着广泛的应用.在各类索梁动力学建模方法中,由于绝对节点坐标方法(absolute nodal coordinate formulation,ANCF)能够描述柔性体的大变形和大转动问题,因此非常适合大变形索梁结构的动力学建模.对绝对节点坐标索梁单元的应变进行分析可知,弯曲变形会引起单元内部轴向应变的不均匀分布,即单元轴向应变与弯曲应变相互耦合.这种应变耦合效应使单元产生伪应变能,导致单元刚度增大,造成单元失真.分析不同弯曲角下的单元应变及应变能可知,弯曲变形越大,单元失真越严重.通过构造等效一维杆单元重新描述轴向应变,实现了轴向应变与弯曲应变解耦.在此基础上推导广义弹性力,得到了绝对节点坐标索梁单元的应变解耦模型.对解耦前后的两种梁模型进行静力学和动力学仿真,结果表明;解耦模型消除了单元伪应变,相比原模型表现出更好的收敛性和曲率连续性,在相同单元数目下具有更高的精度.同时由于解耦模型降低了单元刚度,因此相比原模型,速度曲线中不再有高频振动.     

A new finite element of spatial thin-walled beams

Based on the theories of Timoshenko＇s beams and Vlasov＇s thin-walled members, a new spatial thin-walled beam element with an interior node is developed. By independently interpolating bending angles and warp, factors such as transverse shear deformation, torsional shear deformation and their Coupling, coupling of flexure and torsion, and second shear stress are considered. According to the generalized variational theory of Hellinger-Reissner, the element stiffness matrix is derived. Examples show that the developed model is accurate and can be applied in the finite element analysis of thinwalled structures.     

Buckling and post-buckling of long pressurized elastic thin-walled tubes under in-plane bending

The present paper focuses on the structural stability of long uniformly pressurized thin elastic tubular shells subjected to in-plane bending. Using a special-purpose non-linear finite element technique, bifurcation on the pre-buckling ovalization equilibrium path is detected, and the post-buckling path is traced. Furthermore, the influence of pressure (internal and/or external) as well as the effects of radius-to-thickness ratio, initial curvature and initial ovality on the bifurcation moment, curvature and the corresponding wavelength, are examined. The local character of buckling in the circumferential direction is also demonstrated, especially for thin-walled tubes. This observation motivates the development of a simplified analytical formulation for tube bifurcation, which considers the presence of pressure, initial curvature and ovality, and results in closed-form expressions of very good accuracy, for tubes with relatively small initial curvature. Finally, aspects of tube bifurcation are illustrated using a simple mechanical model, which considers the ovalized pre-buckling state and the effects of pressure.     

Post-buckling behaviour of slender structures with a bi-linear bending moment-curvature relationship

Certain classes of slender structures of complex cross-section or fabricated from specialised materials can exhibit a bi-linear bending moment-curvature relationship that has a strong influence on their global structural behaviour. This condition may be encountered, for instance, in (a) non-linear elastic or inelastic post-buckling problems if the cross-section stiffness may be well approximated by a bi-linear model; (b) multi-layered structures such as stranded cables, power transmission lines, umbilical cables and flexible pipes where the drop in the bending stiffness is associated with an internal friction mechanism. This paper presents a mathematical formulation and an analytical solution for such slender structures with a bi-linear bending moment versus curvature constitutive behaviour and subjected to axial terminal forces. A set of five first-order non-linear ordinary differential equations are derived from considering geometrical compatibility, equilibrium of forces and moments and constitutive equations, with hinged boundary conditions prescribed at both ends, resulting a complex two-point boundary value problem. The variables are non-dimensionalised and solutions are developed for monotonic and unloading conditions. The results are presented in non-dimensional graphs for a range of critical curvatures and reductions in bending stiffness, and it is shown how these parameters affect the structure's post-buckling behaviour.     

Theoretical statistical solution and numerical simulation of heterogeneous brittle materials

The analytical stress-strain relation with heterogeneous parameters is derived for theheterogeneous brittle materials under a uniaxial extensional load, in which the distributions of theelastic modulus and the failure strength are assumed to be statistically independent. This theoreticalsolution gives an approximate estimate of the equivalent stress-strain relations for 3-D heterogeneousmaterials. In one-dimensional cases it may provide comparatively accurate results. The theoreticalsolution can help us to explain how the heterogeneity influences the mechanical behaviors, Further, anumerical approach is developed to model the non-linear behavior of three-dimensional heterogeneousbrittle materials. The lattice approach and statistical techniques are applied to simulate the initialheterogeneity of heterogeneous materials. The load increment in each loading stage is adaptivelydetermined so that the better approximation of the failure process can be realized. When the maximumtensile principal strain exceeds the failure strain, the elements are considered to be broken, which canbe carried out by replacing its Young‘s modulus with a very small value. A 3-D heterogeneous brittlematerial specimen is simulated during a full failure process. The numerical results are in good agreementwith the analytical solutions and experimental data.     

低速冲击下泡沫金属填充薄壁圆管的弯曲行为

采用实验方法研究了低速冲击下泡沫金属填充薄壁圆管的弯曲行为,详细说明了实验方法和原 理。通过与准静态实验结果的比较发现,冲击加载使泡沫金属填充圆管跨中截面的局部压入变形增大,跨中 截面高度变小,结构下缘拉裂破坏延迟。由于结构的惯性效应,锤头总冲击力高于准静态加载时的对应值。