一种考虑横向剪切变形的三角形板弯曲单元

本文从部分协调的三角形薄板弯曲单元出发，并假设横向剪切应变在单元内线性变化，提出了一种考虑横向剪切变形有具有１５个自由度的三角平板弯曲单元。该单元应用于薄板和中等厚度板分析均有较高的精度，计算效率高，可用于工程中具有复杂形状的薄板和中等厚度板结构的分析。     

热载荷作用下梯度多孔材料梁弯曲及过屈曲力学行为的研究

基于Bernoulli-Euler梁理论,引入物理中面解耦了复合材料结构的面内变形与横向弯曲特性,研究了梯度多孔材料矩形截面梁在热载荷作用下的弯曲及过屈曲力学行为．假设沿梁厚度方向材料的性质是连续变化的,利用能量法推导了矩形截面梁的控制微分方程和边界条件,并用打靶法对无量纲化的控制方程进行数值求解．利用计算得到的结果分析了材料的性质、热载荷、边界条件对矩形截面梁非线性力学行为的影响．结果表明,对称材料模型下,固支梁与简支梁均显示出了典型的分支屈曲行为特征,而其临界屈曲热载荷值均会随着孔隙率系数的增加而单调增加．非对称材料模型下,固支梁仍显示出分支屈曲行为特征,但其临界屈曲热载荷不再随着孔隙率系数的变化而单调变化;而对于两端简支梁,发生了弯曲变形,弯曲挠度随载荷的增大而增大．     

不同拉压特性的双层厚壁圆筒极限承载力解答

采用统一强度理论并考虑材料拉伸与压缩弹性模量的差异性，建立均匀内压作用下双层厚壁圆筒的应力表达式，获得了其内压相应的弹性极限解答、塑性极限解答，并分析拉压强度比、拉压模量系数、统一强度理论参数、半径比及分层半径对弹性、塑性极限内压的影响规律．研究结果表明：弹性、塑性极限内压随拉压强度比的增加而减小，但随统一强度理论参数、半径比的增加而增大；弹性极限内压随分层半径的增加呈现先增大后减小变化，随拉压模量系数的增加而一直减小；塑性极限内压与拉压模量系数、分层半径无关．应用于实际工程时，可根据所得结果选择合理的壁厚及分层半径，再根据材料特性确定其他参数，以便更加准确地计算结构的受力状况．     

三维疲劳弹塑性弯曲裂纹塑性区研究

主要研究疲劳载荷作用下的三维弹塑性弯曲裂纹尖端的塑性区问题.用数值解法计算出三维弹塑性弯曲裂纹尖端交变塑性区于裂纹直线部分延长线上的投影长度的最大值和变化幅值,作图分析了三维弹塑性弯曲裂纹尖端交变塑性区尺寸的最大值和变化幅值与三维裂纹体几何尺寸之间的关系.三维弹塑性弯曲裂纹尖端交变塑性区的最大值和变化幅值随着三维裂纹体厚度的增大而减小,随着三维裂纹体厚度的均匀增大,三维弹塑性弯曲裂纹尖端塑性区的最大值和变化幅值不断减小,减小的幅度越来越小,最终趋于平面应变状态下的弹塑性弯曲裂纹尖端塑性区尺寸最大值和变化幅值.当三维裂纹体几何尺寸相同时,三维弯曲裂纹尖端塑性区的最大值和变化幅值随外载荷的不断增大而逐渐增大.建立了一个计算三维弹塑性弯曲裂纹交变塑性区的最大值和变化幅值的崭新理论模型.     

一种考虑横向剪切变形的三角形板弯曲单元

本文从部分协调的三角形薄板弯曲单元出发,并假设横向剪切应变在单元内线性变化,提出了一种考虑横向剪切变形的具有15个自由度的三角平板弯曲单元。该单元应用于薄板和中等厚度板分析均有较高的精度,计算效率高,可用于工程中具有复杂形状的薄板和中等厚度板结构的分析。     

功能梯度矩形板的三维弹性分析

将功能梯度三维矩形板的位移变量按双三角级数展开，以弹性力学的平衡方程为基础．导出位移形式的平衡方程。引入状态空间方法，以三个位移分量及位移分量的一阶导数为状态变量，建立状态方程。考虑四边简支的边界条件，由状态方程得到了功能梯度三维矩形板的静力弯曲问题和自由振动问题的精确解。由给出的均匀矩形板自由振动问题的计算结果表明．与已有的理论解以及有限元方法的计算结果相吻合。假设功能梯度三维矩形板的材料常数沿板的厚度方向按照指数函数的规律变化．进一步给出了功能梯度三维矩形板的自由振动问题和静力弯曲问题的算例分析，并讨论了材料性质的梯度变化对板的动力响应和静力响应的影响。     

轴对称件超塑约束胀形中的摩擦效应

尽管超塑胀形作为一种金属成形方法正日渐受到世界各国的普遍关注，但对超塑约束胀形理论的研究报道却还很少，尤其对胀形过程的有限元模拟研究就更为罕见。针对这种情况，采用大变形刚粘塑性有限元法模拟了轴对称零件向圆筒形凹模内超塑约束胀形的变形过程，着重研究了工具工件之界面摩擦对胀形件厚度分布不均匀性和胀形板料向凹模角部充填性的影响．结果表明，随着摩擦的降低，胀形件的侧向较厚部分能有所减薄，可以改善整个胀形件的厚度均匀性，但当摩擦因子Ａｍ≤０．２（相当于摩擦系数μ≤０．１２）时，胀形件极顶部分的减薄过大；摩擦较小时，胀形板料向凹模角部的充填性较好；在考虑到极点附近厚度适度减薄和胀形板料对凹模角部充填性好的前提下，工艺上应当适当减小摩擦，其最佳状态是μ值约为０．３．为了检验所用刚粘塑性有限元法模拟的可靠性，将计算结果与试验结果作了对比，发现两者相当吻合。     

固支及简支圆板在均布侧压作用下的弹脆塑性承载能力分析

圆板是工程中常见的一种结构形式。本文采用弹脆塑性木构模型，对弹塑性载荷作用下固支圆板的弯曲作了详细的分析研究，得出了对应材料残余强度系数不同取值范围的图板的承载能力解析表达式。文中对加载过程中屈服面的变化情况作了探讨。文未还给出了相应的简支圆板的弹脆塑性承载能力表达式。     

三维冲击荷载弹塑性弯曲裂纹张开位移

主要研究冲击载荷作用下的三维弹塑性弯曲裂纹尖端的张开位移问题.综合考虑了冲击作用应力,三维塑性区域边界上正应力与剪应力,利用二阶摄动方法计算了三维弹塑性弯曲裂纹尖端的张开位移.用数值解法计算出三维弹塑性弯曲裂纹尖端张开位移,作图分析了三维弹塑性弯曲裂纹尖端张开位移与三维裂纹体几何尺寸之间的变化关系.三维弹塑性弯曲裂纹尖端张开位移随着三维裂纹体厚度的增大而减小,随着三维裂纹体厚度的均匀增大,三维弹塑性弯曲裂纹尖端张开位移尺寸不断减小,减小的幅度越来越小,最终趋于平面应变状态下的弹塑性弯曲裂纹尖端张开位移尺寸.当三维裂纹体几何尺寸相同时,三维弯曲裂纹尖端动态张开位移随外部冲击载荷的不断增大而逐渐增大,三维弯曲裂纹尖端动态张开位移随动荷系数的增大而迅速增大,建立了一个计算三维弹塑性弯曲裂纹尖端动态张开位移的崭新理论模型.     

功能梯度与均匀圆板弯曲解的线性转换关系

基于一阶剪切理论, 研究了功能梯度材料圆板与均匀圆板轴对称弯曲解之间 的线性转换关系. 通过理论分析和比较 功能梯度材料圆板和均匀圆板在一阶剪切理论下的位 移形式的轴对称弯曲控制方程, 发现了功能梯度材料圆板的转角与均匀圆板的转角之间 的相似转换关系. 在假设材料性质沿板厚连续变化的情况下, 给出了相似转换系数的解析表 达式. 在此基础上, 进一步导出了一阶剪切理论下功能梯度圆板的挠度与经典理论 下, 均匀圆板的挠度之间的线性关系. 从而, 可将功能梯度材料圆板在一阶剪切理 论下的弯曲问题求解, 转化为相应均匀薄圆板在经典理论下的弯曲问题求解, 以及 转换系数的计算问题. 这一方法为功能梯度非均匀中厚度圆板的求解提供了简捷途 径, 而且更便于工程应用. 作为例子, 采用上述方法分别求得了周边简支和夹紧条 件下, 梯度圆板在均布横向载荷作用下的弯曲解析解, 该解答与Reddy得到的结果 完全吻合.     

Study of an MHD Flow of the Carreau Fluid Flow Over a Stretching Sheet with a Variable Thickness by Using an Implicit Finite Difference Scheme

The present analysis deals with a two-dimensional MHD flow of the Carreau fluid over a stretching sheet with a variable thickness. The governing partial differential equations are converted into an ordinary differential equation by using the similarity approach. The solution of the differential equation is calculated by using the Keller box method. The solution is studied for different values of the Hartmann number, Weissenberg number, wall thickness parameter, and power-law index. The skin friction coefficient is calculated. The present results are compared with available relevant data.     

The effects of sheet thickness on the oscillation of an air-blasted liquid sheet

In the present work, experiments have been performed to study the influence of the sheet thickness on the oscillatory movement of an air-blasted liquid sheet, based on measurements of oscillation frequencies and wavelengths. To measure the wavelengths, a novel double-beam diffraction technique has been implemented. It has been determined that the sheet oscillation frequency depends on the inverse of the square root of the sheet thickness, while the wavelength increases linearly as sheet thickness is increased. Based on these dependencies, suitable non-dimensional groups are proposed to describe the problem. Wave propagation velocities have also been calculated and discussed.     

Strain tensor for large three-dimensional surface deformation of sheet metal from an object grating

Grating techniques are used to determine the three-dimensional deformation and the tangential strain of sheet metal. A grating is fixed on the surface and taken by stereo CCD cameras in different deformation states. By suitable line-following software, the grating coordinates in the images are determined with subpixel accuracy. Using photogrammetric methods, the three-dimensional coordinates are calculated from the image coordinates. The strain usually is determined by means of a deformation gradient, which is calculated from every deformed triangle. In this paper, the gradient is determined in the center of four neighboring meshes using a polynomial approximation of the displacement function in a reference position. The influence of the nontangential deformation is considered. By simulation, a flat sheet metal is deformed to a rotational symmetric surface. The difference of the known exact strain is compared with the numerically derived strain with respect to different grating pitches. The proposed method yields good results even in the case of large spatial deformation. It is applied to the deformation of a hatlike test specimen.     

Analytical springback model for lightweight hexagonal close-packed sheet metal

Experiments have shown that magnesium alloy sheet a common hexagonal close-packed metal, exhibits mechanical behavior unlike that of sheets made of cubic metals (X.Y. Lou et al., 2007, Int. J. Plasticity, 24, 44). The unique stress–strain response includes a strong asymmetry in the initial yield and subsequent plastic hardening. In other words, the stress–strain curves in tension and compression are significantly different. A proper representation of the constitutive relationships is crucial for the accurate evaluation of springback, which occurs due to the residual moment distribution through the sheet thickness after bending. In this paper, we propose an analytical model for asymmetric elasto-plastic bending under tension followed by elastic unloading in order to evaluate the bending moment, which is equivalent to the springback amount. To simplify the calculations, the experimentally measured stress–strain curve of the magnesium alloy sheet was approximated with discrete linear hardening in each deformation region, and the material properties were characterized according to several simplifying assumptions. The bending moment was calculated analytically using the approximate asymmetric stress–strain relationship up to the prescribed curvature corresponding to the radius of the tool in sheet metal forming operations. A numerical example showed an unusual springback increase, even with an increase in the applied force; this is an unexpected result for conventional symmetric materials. We also compared the calculated springback amounts with the results of physical measurements. This showed that the proposed model predicts the main trends of the springback in magnesium alloy sheets reasonably well considering the simplicity of the analytical approach.     

The propagation of vorticity in a viscoelastic fluid

A layer of constant vorticity exists in an infinite space of incompressible isotropic viscoelastic fluid for which the shear stress for rectilinear shearing flows depends linearly on the history of the velocity gradient. At some instant of time the forces maintaining this flow are removed. The subsequent time-dependent vorticity field is calculated explicitly in the particular cases when the fluid is Newtonian and when it is Maxwellian. The limiting case in which the vorticity layer becomes a vortex sheet is also calculated.     

Forming simulation of aluminum sheets using an anisotropic yield function coupled with crystal plasticity theory

This paper describes the application of a coupled crystal plasticity based microstructural model with an anisotropic yield criterion to compute a 3D yield surface of a textured aluminum sheet (continuous cast AA5754 aluminum sheet). Both the in-plane and out-of-plane deformation characteristics of the sheet material have been generated from the measured initial texture and the uniaxial tensile curve along the rolling direction of the sheet by employing a rate-dependent crystal plasticity model. It is shown that the stress–strain curves and R-value distribution in all orientations of the sheet surface can be modeled accurately by crystal plasticity if a “finite element per grain” unit cell model is used that accounts for non-uniform deformation as well as grain interactions. In particular, the polycrystal calculation using the Bassani and Wu (1991) single crystal hardening law and experimental electron backscatter data as input has been shown to be accurate enough to substitute experimental data by crystal plasticity data for calibration of macroscopic yield functions. The macroscopic anisotropic yield criterion CPB06ex2 (Plunkett et al., 2008) has been calibrated using the results of the polycrystal calculations and the experimental data from mechanical tests. The coupled model is validated by comparing its predictions with the anisotropy in the experimental yield stress ratio and strain ratios at 15% tensile deformation. The biaxial section of the 3D yield surface calculated directly by crystal plasticity model and that predicted by the phenomenological model calibrated with experimental and crystal plasticity data are also compared. The good agreement shows the strength of the approach. Although in this paper, the Plunkett et al. (2008) yield function is used, the proposed methodology is general and can be applied to any yield function. The results presented here represent a robust demonstration of implementing microscale crystal plasticity simulation with measured texture data and hardening laws in macroscale yield criterion simulations in an accurate manner.     

Fracture prediction in stretch forming using finite element simulation combined with ductile fracture criterion

Summary  A criterion for ductile fracture is introduced in the finite element simulation of sheet metal forming. From the calculated histories of stress and strain in each element, the fracture initiation site and the critical stroke are predicted by means of the ductile fracture criterion. The calculations are carried out for axisymmetric stretch forming of various aluminium alloy sheets and their laminates clad by mild steel sheets. The predictions so obtained are compared with experimental observations. The results show that the combination of the finite element simulation and the ductile fracture criterion enables the prediction of forming limit in a wide range of sheet metals. Accepted for publication 11 August 1996     

冰排与斜面结构作用时的破断长度分析

研究了冰排与斜面结构的相互作用过程．文中将冰排分为在斜面附近脱离水面部分和较远处漂浮于水面部分两个区域，建立了冰排因弯曲而发生破断的数学模型，根据弯矩分布分析了冰排可能的破断位置，从而得到冰排与斜面结构相互作用时的破断长度解析表达式．结果表明，冰排破断长度既不是冰厚的固定倍数，也不与冰排特征长度直接相关．计算结果与前人实验结果吻合较好，可以为近海冰区抗冰结构的设计提供参考．