开口和闭口薄壁杆件基本微分方程的统一和它的有限元分析

本文将开口和闭口薄壁理论统一在乌曼斯基假设之下,建立开口和闭口薄壁杆件统一的基本微分方程式。通过基本微分方程式建立它们的刚度矩阵,进行有限元分析、对开口和闭口联接处的位移谐调,应用最小二乘法处理。 本文还指出文[3]中的错误。     

薄壁曲杆的有限元法

本文提出了开口和闭口截面薄壁曲杆的刚度矩阵.位移模式采用多项式逼近,并考虑了温度影响.对于闭口截面的薄壁曲杆,利用了乌曼斯基假设.文中还利用截面线元素,按线性分布规律逼近,讨论了圣维南翘曲函数的离散化计算.作了不同于藤谷义信的分析.     

考虑约束扭转的薄壁梁单元刚度矩阵

推导了薄壁空间梁单元刚度矩阵 ,考虑了双向弯曲及截面约束扭转对杆件轴向变形的影响 ;计算了截面的翘曲变形 ,以及二次剪应力对翘曲变形的影响 ,可适用于任意截面 (包括开口、闭口和混合剖面 )的薄壁杆件。计算结果表明 ,考虑约束扭转的薄壁梁单元刚度矩阵有相当好的精确度 ,可以用于薄壁杆件的静动力分析。     

一种新的薄壁杆件单元扭转刚度矩阵

本文提出一种新的薄壁杆件单元扭转刚度矩阵,它能够计及二次剪应力对翘曲变形的影响,并适用于任意剖面(包括开口,闭口和混合剖面)的薄壁杆件。计算表明,这个新的扭转刚度矩阵有相当好的精确度,可以代替Kawai或Gunnlaugsson-Pedersen的刚度矩阵,用于薄壁杆件的有限元静动力分析。     

弹塑性有限元分析的拟线性法

本文在理论分析的基础上,建立了弹塑性有限元分析的拟线性法.这种方法可用于解算理想弹塑性材料、强化和软化材料的弹塑性问题.计算时可只采用一个弹性刚度矩阵.选用了三个算例以说明本法的计算精度.     

剪力滞效应对箱形梁挠度影响的研究

从剪力滞翘曲应力的轴向平衡条件出发,选取双室箱梁的合理翘曲位移函数,引入相应于剪力滞翘曲变形的惯性矩和惯性积等几何特性,用能量变分法建立薄壁箱梁剪力滞效应分析的控制微分方程。通过求解控制微分方程,导出集中荷载和均布荷载作用下简支箱梁和悬臂箱梁的挠度公式及有限梁段单元刚度矩阵,模型试验和ANSYS壳单元计算结果证实了其正确性。结合简支、悬臂和连续箱梁数值算例,具体分析剪力滞效应对箱梁挠度的提高程度。结果表明,无论在集中荷载还是均布荷载作用下,剪力滞效应对简支箱梁的挠度均有显著提高。在集中荷载作用下,剪力滞效应对连续箱梁挠度的提高可达14%;对于跨宽比约为4.0～6.0的简支箱梁,可将按初等梁计算的跨中挠度乘以提高系数1.05～1.11;计算悬臂箱梁的挠度时,一般可以忽略剪力滞效应的影响。     

开口薄壁杆件结构稳定分析的精确单元和两步求解算法

从控制微分方程的通解出发,构造受偏心压力作用开口薄壁杆件的精确形函数,建立用于开口薄壁杆件结构稳定性分析的精确有限元,得到了单元刚度矩阵和几何刚度矩阵的显式表达,提出了计算给定区间内各阶临界荷载以及相应失稳模态的两步计算方法。计算结果表明,与常规单元相比,采用精确单元无需进行网格细分就可以获得精确的数值结果,结合本文的两步求解算法,可以准确获得给定区间内全部临界荷载和失稳模态。     

预应力混凝土箱梁横向框架效应有限元分析

以预应力混凝土箱梁桥顶板纵向开裂为工程背景,以箱梁顶板与腹板交接处两个角点的转角为基本未知量,基于箱梁纵向一维单元,利用最小势能原理推导出箱梁横向框架效应有限单元法的单元刚度矩阵,及相应的单元等效节点力,对比计算表明,该有限单元法计算精度与理论分析基本一致。实桥分析表明：现行规范计算箱梁角点弯矩相对安全,对于中点弯矩相对偏于不安全,设计过程不充分考虑混凝土箱梁横向框架效应可能导致顶板的纵向开裂。本文提出的有限单元法具有概念简单、运算速度快等优点,弥补了目前商业程序的不足,对工程设计具有一定的参考意义和推广价值。     

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

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

单支开口薄壁梁屈曲分析的Riccati有限条传递矩阵法

针对单支开口薄壁梁屈曲问题,通过采用半解析有限条法推导薄壁梁的条元控制方程、传递矩阵基础形成单元传递方程、Riccati变换改善数值计算稳定性,提出了开口单支薄壁梁屈曲分析的Riccati有限条传递矩阵法。为了验证该方法的正确性与高效性,通过提出的Riccati有限条传递矩阵法和有限元法分析了一般支撑条件下两种不同截面形式的薄壁梁屈曲临界载荷与屈曲模态,计算发现两种方法的结果具有很好的一致性,且与传递矩阵法相比,Riccati有限条传递矩阵法数值稳定性好。因此,Riccati有限条传递矩阵法可计算广义边界条件的薄壁结构屈曲问题。     

Improved torsional analysis of laminated box beams

The improved torsional analysis of the laminated box beams with single- and double-celled sections subjected to a torsional moment is performed by introducing 14 displacement parameters. For this, a thin-walled laminated box beam theory considering the effects of shear and elastic couplings is presented. The governing equations and the force-displacement relations are derived from the variation of the strain energy. The system of linear algebraic equations with non-symmetric matrix is constructed by introducing the displacement parameters and by transforming the higher order simultaneous differential equations into first order ones. This numerical technique determines eigenmodes corresponding to 12 zero and 2 non-zero eigenvalues and derives displacement functions for displacement parameters based on the undetermined parameter method. Finally, the element stiffness matrix is determined using the member force-displacement relations. The theory developed by this study is validated by comparing several torsional responses from the present approach with those from the finite element beam model using the Lagrangian interpolation polynomials and three-dimensional analysis results using the shell elements of ABAQUS for coupled laminated beams with single- and double-celled sections.     

基于新型广义位移的箱形梁扭转单元及应用

为了改进变截面连续箱梁桥的扭转分析理论,将截面总扭转角分解为自由翘曲扭转角和约束剪切扭转角,选取自由翘曲转角扭率作为广义位移,提出一个2节点8自由度的扭转梁段单元。从约束扭转控制微分方程出发,推导单元刚度矩阵及等效节点荷载列阵。引入应力增大系数,以反映约束扭转对初等梁应力的增大效应。数值算例验证了本文梁段单元的可靠性。最后对一个三跨变截面连续箱梁桥进行分析,结果表明,双力矩影响线与弯矩影响线较为类似,按双力矩影响线进行最不利荷载加载时最大应力值偏小;应力增大系数在集中荷载作用截面出现极值,均发生在腹板与顶板交点处;利用偏载放大系数来考虑扭转附加效应时,不宜考虑弯曲正应力较小及翘曲正应力出现极值的梁段区域。     

采用有限单元法计算梁任意形状截面特性

采用将梁截面离散化的方式,用数值积分计算截面的几何特性,并根据梁剪切变形和扭转理论,利用变分原理建立截面的有限元法方程,求解任意形状截面的扭转常数、剪切中心以及剪切面积修正系数等特性.本方法适用于各种形式的截面,具有计算精度高及适应性强的特点.根据上述理论编制了相应程序,按照不同的单元划分方式,分别计算出矩形截面截面特性,与理论解进行比较;又对舟山市定海长峙至岙山预应力混凝土连续箱梁截面进行了计算,并与Ansys结果进行比较,均证明采用本文的计算方法能得到满意的结果,且该方法适用于各种形状的截面形式.     

Out-of-plane responses of a circular curved Timoshenko beam due to a moving load

For the cases of using the finite curved beam elements and taking the effects of both the shear deformation and rotary inertias into consideration, the literature regarding either free or forced vibration analysis of the curved beams is rare. Thus, this paper tries to determine the dynamic responses of a circular curved Timoshenko beam due to a moving load using the curved beam elements. By taking account of the effect of shear deformation and that of rotary inertias due to bending and torsional vibrations, the stiffness matrix and the mass matrix of the curved beam element were obtained from the force–displacement relations and the kinetic energy equations, respectively. Since all the element property matrices for the curved beam element are derived based on the local polar coordinate system (rather than the local Cartesian one), their coefficients are invariant for any curved beam element with constant radius of curvature and subtended angle and one does not need to transform the property matrices of each curved beam element from the local coordinate system to the global one to achieve the overall property matrices for the entire curved beam structure before they are assembled. The availability of the presented approach has been verified by both the existing analytical solutions for the entire continuum curved beam and the numerical solutions for the entire discretized curved beam composed of the conventional straight beam elements based on either the consistent-mass model or the lumped-mass model. In addition to the typical circular curved beams, a hybrid curved beam composed of one curved-beam segment and two identical straight-beam segments subjected to a moving load was also studied. Influence on the dynamic responses of the curved beams of the slenderness ratio, moving-load speed, shear deformation and rotary inertias was investigated.     

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.     

钢-混凝土组合箱梁梁段有限元法研究

组合梁界面滑移将减小组合梁刚度,增大变形,影响构件性能;剪应力沿截面横向分布不均匀,造成其弯曲正应力的横向分布呈曲线形状,即剪力滞效应;同时组合梁往往重载,具有较小的跨高比,剪切变形不可忽略.根据虚功原理,建立了同时考虑滑移效应、剪力滞及剪切变形效应的组合梁单元刚度矩阵及等效节点力向量,并在此基础上编制了组合梁梁段有限元程序.利用本文程序对现有组合梁试件的混凝土顶板应力、钢梁底板应力、跨中挠度和梁端滑移进行了计算,并将本文计算结果与解析法计算结果及试验结果进行了比较.结果表明,本文计算结果与解析法计算结果及试验结果吻合良好;同时,本文计算结果具有较好的稳定性,验证了本文计算方法的正确性.本文所建立的梁单元刚度矩阵同时考虑了剪切变形、剪力滞及滑移效应的影响,符合工程实际,为有限梁段法分析组合箱梁提供了理论基础.     

Measurement of Torsionally Induced Shear Stresses in Shrink-Fit Assemblies

Shear stresses along the shaft/hub interface in shrink-fit components, generated by torsional loads, can drive premature failure through fretting mechanisms. It is difficult to numerically predict these shear stresses, and the associated circumferential slip along the shaft/hub interface, due to uncertainties in frictional behaviour and the presence of steep stress gradients which can cause meshing problems. This paper attempts to provide validation of a numerical modelling methodology, based on finite element analysis, so the procedure may be used with confidence in fitness-for-purpose cases. Very few experimental techniques offer the potential to make measurements of stress and residual stress interior to metallic components. Even fewer techniques provide the possibility of measuring shear stresses. This paper reports the results of neutron diffraction measurements of shear stress and residual shear stress in a bespoke test specimen containing a shrink-fit. One set of measurements was made with a torsional load ‘locked-in’. A second set of measurements was made to determine the residual shear stress when the torsional load had been applied and removed. Overall, measurement results were consistent with numerical models, but the necessity for a small test specimen to allow penetration of the neutron beam to the measurement locations meant the magnitude of shear stresses was at the limits of what could be measured experimentally.     

钢－混凝土组合梁非线性有限元分析

提出了部分剪力连接钢-混凝土组合梁非线性有限元分析的简化模型，该模型同时考虑了钢、混凝土材料的非线性本构关系和剪力连接件的滑移非线性对组合梁刚度和强度的影响，推导了考虑滑移的剪力连接件单元刚度矩阵。利用该模型计算了连续组合梁的极限荷载、挠度、应力以及滑移，与已有的理论计算结果和实验结果吻合，证明本方法分析钢-混凝土组合梁非线性的可靠性和实用性。     

Shear deformable finite beam elements for composite box beams

The shear deformable thin-walled composite beams with closed cross-sections have been developed for coupled flexural, torsional, and buckling analyses. A theoretical model applicable to the thin-walled laminated composite box beams is presented by taking into account all the structural couplings coming from the material anisotropy and the shear deformation effects. The current composite beam includes the transverse shear and the restrained warping induced shear deformation by using the first-order shear deformation beam theory. Seven governing equations are derived for the coupled axial-flexural-torsional-shearing buckling based on the principle of minimum total potential energy. Based on the present analytical model, three different types of finite composite beam elements, namely, linear, quadratic and cubic elements are developed to analyze the flexural, torsional, and buckling problems. In order to demonstrate the accuracy and superiority of the beam theory and the finite beam elements developed by this study,numerical solutions are presented and compared with the results obtained by other researchers and the detailed threedimensional analysis results using the shell elements of ABAQUS. Especially, the influences of the modulus ratio and the simplified assumptions in stress–strain relations on the deflection, twisting angle, and critical buckling loads of composite box beams are investigated.     

Computational analysis of multi-stepped beams and beams with linearly-varying heights implementing closed-form finite element formulation for multi-cracked beam elements

The model where the cracks are represented by means of internal hinges endowed with rotational springs has been shown to enable simple and effective representation of transversely-cracked slender Euler–Bernoulli beams subjected to small deflections. It, namely, provides reliable results when compared to detailed 2D and 3D models even if the basic linear moment–rotation constitutive law is adopted.This paper extends the utilisation of this model as it presents the derivation of a closed-form stiffness matrix and a load vector for slender multi-stepped beams and beams with linearly-varying heights. The principle of virtual work allows for the simple inclusion of an arbitrary number of transverse cracks. The derived at matrix and vector define an ‘exact’ finite element for the utilised simplified computational model. The presented element can be implemented for analysing multi-cracked beams by using just one finite element per structural beam member. The presented expressions for a stepped-beam are not exclusively limited to this kind of height variation, as by proper discretisation an arbitrary variation of a cross-section’s height can be adequately modelled.The accurate displacement functions presented for both types of considered beams complete the derivations. All the presented expressions can be easily utilised for achieving computationally-efficient and truthful analyses.