输电杆塔螺栓连接区域局部结构有限元分析

利用ABAQUS有限元软件建立了500k V同塔双回输电线路塔线耦合体系模型,在不同脱冰率情况下对塔线体系的动力响应进行了数值模拟。分析了杆塔在各脱冰率下的受力情况并结合现场倒塔情况,确定了杆塔危险区域,建立了该危险区域局部结构三维有限元模型。将由塔线体系模型计算得到的危险区域边界的力施加于三维局部模型上,计算了该区域在导线脱冰过程中的应力。模拟研究连接螺栓预紧力不同的情况下危险区域的应力和变形,发现螺栓预紧力减小时杆塔容易破坏,这一结论对输电杆塔的设计和安装具有指导意义。     

三维自适应FE-SPH耦合算法在多层间隔金属靶侵彻问题中的应用

鉴于有限元算法不能有效地模拟侵彻过程所产生的金属碎片, 本文中基于三维自适应FE-SPH耦合算法的基本理论, 自主开发了模拟多层间隔金属靶侵彻问题的三维FE-SPH耦合计算程序。该程序采用四面体单元对多层间隔金属靶侵彻模型进行初始离散, 计算过程中, 当四面体单元等效塑性应变超过某一设定值时, 单元自动转化为SPH粒子, 并引入有限单元-粒子接触算法和耦合算法, 实现大变形和破碎区域采用SPH方法计算, 克服有限元法单元畸变存在的问题。多层间隔靶侵彻算例分析表明, 三维FE-SPH耦合计算程序采用等效塑性应变作为转化判据计算结果较稳定, 并且能够有效地再现侵彻过程中所产生的碎片, 能够模拟侵彻碎片对后层靶的毁伤效应。     

横向撞击下短梁的剪切破坏研究

应用实验和数值模拟手段研究了 2 0 # 钢梁在横向撞击下的双剪破坏问题。实验中通过测量加速度 ,分析了梁破坏需要吸收的能量 ;数值分析中采用Cauchy应力对数应变Cowper Symonds过应力本构方程模拟材料率相关行为 ,累积塑性应变破坏准则与单元死活技术结合模拟材料弱化行为 ,数值分析给出的结构破坏吸收能量与实验结果基本吻合 ,由计算得到的剪切塑性铰长度与S .B .Menkes等的实验结果基本一致。     

岩石试件端面摩擦效应数值模拟研究

试件端面摩擦效应直接影响试件内的塑性等效应变、侧向位移的分布和单元应力应变曲线。本文运用ANSYS中的接触单元模拟了平面应变状态下端面摩擦效应对塑性等效应变、侧向位移和单元应力应变曲线的影响,得到了不同摩擦系数时塑性等效应变及侧向位移的渐进变化形式。当接触面摩擦较小时,塑性等效应变图案为上下两个X形网络,侧向位移上下分布均匀;当接触面摩擦增大时,塑性等效应变网络向中部靠拢并且明显增大,侧向位移上下分布不均匀,中部较上下端面位移大;当试件端面侧向位移被限制,即摩擦力很大时,塑性等效应变网络变为一个X形局部化带,侧向位移分布更加不均匀,中部明显隆起。     

基于应变能的复杂应力场低周疲劳分析

塑性应变能使材料微观组织结构发生不可逆变化,从而引起等效宏观应力,该应力随循环加载而增大．假定材料疲劳源处破坏是由最大拉应力引起的,最大等效宏观应力与外加应力叠加达到材料本征断裂应力时形成微裂纹．微裂纹引起上述两部分应力变化,继续加载直至宏观裂纹出现,从而得到材料的疲劳寿命．本文所建立的多轴疲劳寿命公式包含材料参数、拉应力以及塑性应变能等,以上数据可通过单轴疲劳数据和有限元方法获得．通过对SM45C材料的计算验证,表明该模型对多轴随机应变加载低周疲劳寿命,具有良好的预测结果．     

叠加型A-F类随动强化模型塑性应变的数值计算法

以Chaboche随动强化模型为例,在M isses屈服准则及正交流动准则的前提下,推导了叠加型A rm-strong-F rederick(A-F)类随动强化模型塑性应变的数值计算法,联合利用四阶龙格-库塔法与径向返回法实现数值计算中的内部平衡迭代。同时推导了统一切向矩阵以便确定每一平衡迭代后的试算应变。利用AN SY S提供的U PF s将算法嵌入到AN SY S有限元程序,实现了叠加型A-F类随动强化模型塑性应变的数值计算,并利用四边形单元模拟了单轴循环加载时的棘轮应变,计算结果能够很好地与实验值吻合。     

饱和多孔介质中的混合有限元法和有限应变下应变局部化分析

对基于Biot理论的饱和多孔介质中动力-渗流耦合分析提出了一个耦合场混合元．固相位移、应变和有效应力以及流相压力、压力梯度和Darcy速度在单元内均处理为独立变量分别插值．基于胡海昌-Washizu三变量广义变分原理给出的饱和多孔介质动力-渗流耦合问题控制方程的单元弱形式，导出了单元公式．进一步导出了考虑压力相关非关联塑性的非线性单元公式和发展了相应的一致性算法．对几何非线性分析，采用了共旋公式途径．数值结果例题显示所发展耦合场混合元模拟大应变下由应变软化引起以应变局部化为特征的渐进破坏现象的性能．     

耦合蠕变损伤的Chaboche粘塑性模型的研究

Chaboche粘塑性统一本构模型由于没有包含损伤变量,在应力恒定时只能得到不变的塑性应变速率,因此不能描述蠕变第三阶段的加速过程.该文按照Lemaitre有效应力的概念,采用Kachanov损伤演化方程,推导了耦合损伤的Chaboche粘塑性流动方程和硬化方程,并链接到有限元软件ANSYS中,将之用于高温合金钢P91的蠕变寿命计算,结果表明耦合损伤模型的模拟值与实验数据基本吻合;为该模型描述蠕变损伤和疲劳交互作用奠定了基础.     

蠕变硬化材料中Ⅱ型扩展裂纹尖端场特性研究

采用弹牯塑性力学模型,对蠕变硬化材料中平面应变扩展裂纹尖端场进行了渐近分析.假设人工粘性系数与等效塑性应变率的幂次成反比,通过量级匹配表明应力和应变均具有幂奇异性,奇异性指数由粘性系数中等效塑性应变率的幂指数唯一确定.通过数值计算讨论了Ⅱ型准静态扩展裂纹尖端场的分区构造以及裂纹尖端应力和应变场的特性随各材料参数的变化规律,结果表明裂尖场由材料的粘性和塑性共同主导.当硬化系数为零时裂尖场可退化为相应的HR场.     

一种新的弹塑性杂交应力元法

本文基于一个改进的弹塑性的Hellinger/Reis■ner 混合变分原理构造了一种用于解弹塑性问题的四节点等参杂交应力元.新的模型中,在单元内增加了等效应力增量、塑性等效应变增量及不协调位移变量,从而使单元内的屈服准则及流动法则平均得到满足,不协调位移改进了单元应力精度.计算表明,新的模型可以提高弹塑性杂交法的精度和计算效率.     

梁柱简化模型在自立塔塔线体系动态响应中的适用性研究

对输电塔进行合理简化可以提高塔线体系动力学仿真的效率。本文给出自立塔梁柱简化模型的计算方法,并提出利用梁柱简化模型计算方法建立自立塔塔线体系整体模型,同时采用桁梁混合模型建立精细化塔线体系整体模型,对两种模型塔线体系静力特性及振型和固有频率等动力特性进行对比分析。以脱冰工况为例,采用生死单元技术将施加在输电线节点上的集中质量单元杀死来模拟脱冰,实现对塔线体系动力学响应的有限元模拟,研究塔线体系简化模型在动态响应中的适用性。结果表明,两种模型弯曲变形误差小,低阶的振型相同,固有频率值误差小,动力特性基本相同;脱冰工况下,自立塔节点位移和塔材内力时程曲线一致,在提高计算效率的情况下,能有效保证计算精度。     

山地输电线动刚度计算及模态分析

为了揭示地形变化对连续档输电线动态特性的影响，本文提出了基于任意高差角的动刚度理论．首先将单档导线视为连续档导线的子结构，构建了两悬挂点不等高的单档导线振动方程，并在该理论方程中引入了二次项，能够考虑重力沿弦向分量的影响，推导出绝缘子串绕其悬挂点转动的动刚度计算方法，通过动刚度理论研究连续档导线横向振动模态和频率，并应用有限元软件ABAQUS 验证理论模态和频率公式的准确性．所得结果表明，动刚度的计算方法具有较高的精度，理论模态与有限元输出结果吻合较好，同时获得的频率和模态可为山地高压输电线路防舞技术的研究提供依据．     

激光刻痕聚丙烯塑料失效破坏的参数识别

采用等效截面积法建立了激光加工微孔的有限元模型及带激光刻痕聚丙烯塑料板的冲击模型,并 通过参数识别方法确定了材料的失效应变。通过极差分析得出该塑料材料的冲击断裂强度随撕裂线直径的 增大而增大,随冲击速度的增大而减小。这与实验结果一致。     

An analysis of the shear failure of rigid-linear hardening beams under impulsive loading

A theoretical rigid-plastic analysis for the dynamic shear failure of beams under impulsive loading is presented when using a travelling plastic shear hinge model which takes into account material strain hardening. The maximum dynamic shear strain and shear strain-rate can be predicted in addition to the permanent transverse deflections and other parameters. The conditions for the three modes of shear failure, i.e., excess deflection failure, excess shear strain failure and adiabatic shear failure are analyzed. The special case of an infinitesimally small plastic zone is discussed and compared with Nonaka's solution for a rigid, perfectly plastic material. The results can also be generalized to examine the dynamic response of fibre-reinforced beams.     

Experimental and FEM study of windshield subjected to high speed bird impact

Both experimental and finite element model (FEM) results are presented for the dynamic strength behavior of windshield subjected to bird impact. The experimental data taken from a series of high speed photographs are compared with the numerical results predicted by using FEM in which the windshield was modeled entirely with solid elements and the bird body was approximately simulated by an elastic-plastic material with failure element behavior. Effective plastic strain and element pressure were adopted as the failure criteria and once the pressure or the effective plastic strain of an element reached the critical value, the element would lose the tensile resistance capability completely. The deflection and stress distribution in the windshield were obtained. It is shown that the result from the finite element analysis agrees with those from the full-scale bird impact test. The project supported by the National Natural Science Foundation of China (10272011)     

不同加载状态下TA2钛合金绝热剪切破坏响应特性

一般认为绝热剪切现象在宏观上表现为材料动态本构失稳,即热软化大于应变硬化.本文采用帽型受迫剪切试样研究TA2钛合金的动态力学特性和本构失稳过程.首先对剪切区加载应力状态进行理论和数值分析,通过合理设计帽型试样,剪切区变形可近似按剪切状态处理;结合二维数字图像相关法(two-dimensional digital image correlation,DIC-2D)直接测试试样剪切区应变演化,给出帽型受迫剪切实验的等效应力-应变响应曲线.进一步,利用Hopkinson压杆对TA2钛合金开展动态压缩及帽型剪切对比试验研究,比较压缩、剪切试验得到的等效应力-应变曲线,采用"冻结"试样方法分析试样中绝热剪切局域化演化过程,探讨不同加载状态下TA2钛合金的绝热剪切破坏现象及其动态力学响应特性.实验结果表明,在塑性变形初始阶段,动态压缩及剪切加载下的等效应力-应变曲线符合较好,但随塑性损伤发展及绝热剪切带形成,两者出现分离,表明损伤及绝热剪切演化过程与应力状态相关.剪切试样实验得到的本构"软化"特性能够反映绝热剪切带起始、破坏演化过程的力学响应特性,而在动态压缩实验中,即使试样中已出现双锥形的绝热剪切带及局部裂纹分布,其表观等效应力-应变曲线并不出现软化特征,动态压缩实验无法得到关于绝热剪切起始、发展以及破坏的本构软化响应特性.     

Equivalent nonlinear beam model for the 3-D analysis of shear-type buildings: Application to aeroelastic instability

Tower buildings can be very sensitive to dynamic actions and their dynamic analysis is usually carried out numerically through sophisticated finite element models. In this paper, an equivalent nonlinear one-dimensional shear–shear–torsional beam model immersed in a three-dimensional space is introduced to reproduce, in an approximate way, the dynamic behavior of tower buildings. It represents an extension of a linear beam model recently introduced by the authors, accounting for nonlinearities generated by the stretching of the columns. The constitutive law of the beam is identified from a discrete model of a 3D-frame, via a homogenization process, which accounts for the rotation of the floors around the tower axis. The macroscopic shear strain in the equivalent beam is produced by the bending of columns, accompanied by negligible rotation of the floors. A coupled nonlinear shear–torsional mechanical model is thus obtained. The coupling between shear and torsion is related to a non-symmetric layout of the columns, while mechanical nonlinearities are proportional to the slenderness of the columns. The model can be used for the analysis of the response of tower buildings to any kind of dynamic and static excitation. A first application is here presented to investigate the effect of mechanical and aerodynamic coupling on the critical galloping conditions and on the postcritical behavior of tower buildings, based on a quasi-steady model of aerodynamic forces.     

Numerical simulation of penetration of aluminum targets by spherical-nose steel rods

Numerical simulations of the penetration processes in aluminium blocks by spherical-nose steel rods were performed in this study. The specific impact configuration of this study involves 152-mm diameter 6061-T651 aluminum bars impacted by spherical-nose projectiles machined from T-200 maraging steel rods at nominal impact velocities between 300 and 1000 m/s. The transient dynamic finite element code LS-DYNA2D was used for the numerical analysis. The erosion capability in LS-DYNA2D was exercised in conjunction with the maximum equivalent plastic strain criterion to carry out failure simulations in the target. Calculated results were compared to the experimental data. Good correlation was obtained.     

等径通道挤压过程三维有限元模拟

利用三维有限元模型对等径通道挤压过程中变形分布的不均匀性进行了模拟,通过对不同摩擦系数下截面等效塑性应变分布比较发现:沿三个方向截面上的变形分布都不是均匀的,这说明二维有限元模型不能真实模拟等径通道挤压过程中试样中的变形分布。此外,摩擦对等效塑性应变分布及挤压力的影响较大,截面变形不均匀参数和最大等效塑性应变出现的位置随摩擦系数的增大而变化。压力-位移曲线的变化可以结合试样的变形过程来解释。