V型缺口裂端的三维应力状态及约束分析

利用三维有限元方法，研究了有限厚度板中V型缺口根部穿透裂纹前沿的三维弹性应力场。对不同厚度、不同缺口张开角和裂纹长度对应力强度因子及裂尖附近的三维约束程度的影响进行了分析，同时还讨论了三维约束区的大小。研究结果显示：当缺口张开角小于60度时，不同缺口的应力强度因子和离面约束因子的分布基本一致，角度的影响不明显。应力强度因子是厚度的函数，板中面的应力强度因子随厚度的增加逐渐减小趋近干平面值，最大为1．08倍的平面值。当板厚超过15倍的缺口深度时，应力强度因子最大值将从中面转移至接近自由表面位置，距中面约0．4倍板厚。三维约束非常明显的区域在裂尖前沿约0.45倍厚度的范围内．二维到三维的过渡区在裂尖前沿1.5倍厚度的区域内；在中面上三维效应影响区最大，随着离中面距离的增加逐渐减小，在自由表面上降为0。     

有限厚双边U型切口板的三维弹塑性有限元分析

针对单向拉伸载荷作用下的有限厚双边 U型切口板 ,本文对其切口根部的三维效应进行了详细的三维弹塑性大变形有限元分析 ,得到了不同于平面切口问题的新结果 :( 1 )在线弹性阶段 ,三维切口根部的应力集中因子 Kt高于二维情形 ,且 Kt在板厚 ( B)约为 8倍切口根部半径 (ρ)时达到峰值 ;进入塑性后 ,切口根部的 Kt随载荷水平的增加逐渐下降 ,并介于平面应力和平面应变情形之间 ;( 2 )切口前缘等效应变εe的最大值出现在约束最高的中面 ,且可比约束最低的自由表面处的应变值高一倍 ;( 3)薄板中 ,塑性区尺寸 xp 的最大值不是出现在约束较低的自由表面 ,而是出现在约束较高的中面 ,但随 B/ρ的增大逐渐向自由表面外推 ;( 4 )面内应力比 Tx在切口根部塑性区内的分布对板厚和载荷水平不敏感 ,可以用 Hill的滑移线解很好地近似。另外 ,对三维影响区内一些重要的三维约束参量进行了详细分析 ,总结了它们的三维分布规律。     

功能梯度压电圆板自由振动问题的三维精确分析

本文对周边为广义刚性滑动和广义简支两种边界条件下的功能梯度压电材料圆板自由振动问题进行分析。根据轴对称横观各向同性压电材料基本方程，并利用有限Hankel变换得到了功能梯度压电材料圆板的状态空间方程。假设材料的机械和电学性质均沿板厚方向按统一的指数函数形式梯度分布，从而获得了周边为广义刚性滑动和广义弹性简支两种边界条件下功能梯度压电圆板自由振动问题的三维精确频率方程，该方程是一个关于自由振动频率的超越方程，通过求解该超越方程可得到在不同板厚以及不同的材料性质梯度变化情况下的圆板自由振动频率值，结果表明在相同的材料性质梯度变化情况下频率均随着板厚增加而增大，而在相同的板厚情况下频率则随材料性质梯度变化指数的增大而减小的结论。     

有限厚度板穿透裂纹前缘附近三维弹性应力场分析

通过三维有限元计算来研究有限宽度、有限厚度含有穿透裂纹板的裂纹前缘应力场，从中找出应力强度因子与板的厚度、裂纹长度之间的关系，同时还分析了裂尖的三维约束程度和三维约束区的大小。分析结果表明：应力强度因子沿厚度的分布是不均匀的，应力强度因子的最大值及其位置与厚度有关；有限厚度板中面应力强度因子(KI)m-p及最大应力强度因子(KI)max均大于平面应力或平面应变的应力强度因子。对有限厚度裂纹问题，按平面应力或平面应变来考虑是不安全的；板中面的应力强度因子(KI)m-p及最大应力强度因子(KI)max是厚度B／a的函数；板的中面离面约束系数Tx最大，自由面(z=B)Tx=0。沿厚度方向裂尖附近的离面约束系数Tx也是z／B和B／a的函数，随着厚度的增加离面约束系数Tx增大，离中面越近离面约束系数Tx越大。Tx随着x的增大急剧减小，三维约束影响区域大小大约为板厚的一半，且裂纹长度a／W对应力强度因子沿厚度变化规律及Tx影响区域大小影响较小。     

压电功能梯度板自由振动的三维解

基于三维弹性理论和压电理论，导出了有限长矩形压电功率梯度板的动力学方程及相应的边界条件，并用幂级数展开法进行了求解，得到了压电功能梯度板自由振动的三维精确解公式，求解了自由振动的固有频率，并分析了压电系数的梯度变化对不同电学边界条件下压电板的自由振动频率的影响，结果可用于校核不同的近似理论及理解压电结构的动态行为。     

板厚、温度和速度对单搭接胶接接头强度的影响

为了探讨单搭接胶接接头结构的破坏机理,采用Hopkinson拉杆实验技术,研究了板厚、温度和速度(试样端部的最大速度)对单搭接胶接接头强度的影响。实验结果表明,板越厚、温度越低、速度越大,接头的强度越大。运用有限元方法分析了胶层的应力分布,通过引入剥离应力对试样强度的影响因子,研究发现板越厚、温度越低、速度越大,剥离应力对接头强度的影响越小。     

考虑靶背自由表面和开裂影响的刚性尖头弹贯穿金属靶板模型

将靶体视为不可压缩材料,假定空腔膨胀产生塑性-弹性响应分区,构造了靶背自由表面效应的衰减函数。将衰减函数乘以可压缩幂次硬化材料的阻力方程,得到了弹体贯穿金属靶板的阻力函数。基于弹性衰减-塑性衰减-开裂3个阶段,建立了同时考虑靶体可压缩性、靶背自由表面和开裂影响的弹体贯穿有限厚金属靶板的分析模型,推导得出了弹体瞬时速度的解析方程,并采用数值方法计算得到了弹体的过载、瞬时速度和残余速度。通过与6组实验数据和已有模型的对比得到,当靶板厚度和弹体冲击速度在一定范围内时,需要考虑自由表面效应的影响。     

三维气泡与自由表面相互作用的直接数值模拟

采用VOF中的PLIC界面重构方法数值模拟了三维气泡在液流中上升并与自由表面相互作用的运动.分别考察了不同初始高度,有无来流及有无再生气泡对气泡上升高度、上升速度、压力及与自由表面相互作用等的影响.结果表明:气泡初始位置越低,顶端上升的高度越大,自由面隆起的范围更广.越靠近自由表面,底部射流横向发展越窄,而向上的压力梯度,气泡上升速度,底部射流上升高度越大,反之则反;但如果底部射流均在接近自由表面以前已横向发展充分,则差别不大.气泡外形、上升高度、破裂时间以及上升速度与来流无关.产生再生气泡后,原生气泡与再生气泡相吸,相互加速对方的上升;自由表面抬升的高度增幅较大,范围拓宽,上升速度也大大增加,且再生气泡越多,自由表面隆起的范围越大.     

拉伸载荷作用下共面表面裂纹间应力强度因子影响系数的有限元分析

采用参数化有限元方法,结合节点力法和循环迭代算法,对一有限厚矩形板表面有两个相邻共面半椭圆表面裂纹在拉伸载荷作用下进行了求解,得到了两裂纹在不同形状和相隔距离时的应力强度因子的影响系数,计算结果对含三维广布裂纹结构的剩余强度和疲劳寿命有参考意义.     

厚圆板轴对称振动的弹性力学解

本文以轴对称三维弹性力学基本方程为基础,导出厚圆板强迫振动的状态方程式。利用Ｍａｃｌａｕｒｉｎ级数和Ｓｙｌｖｅｓｔｅｒ定理,厚圆板的位移和应力可以用中面位移和应力的微分算子表示。通过载荷分解和圆板表面条件,可以得到厚圆板在对称载荷与反对称载荷作用下的振动控制方程。求解了厚圆板在周边固支和简支条件下的对称与反对称的自由振动问题。通过数值计算得到了这两类自由振动的固有频率。本文的方法适用于求解厚圆板在     

基于贝叶斯信息融合的岸桥有限元模型修正研究

为建立精确的岸桥有限元模型,研究了基于贝叶斯信息融合的模型修正方法．通过方差分析,确定待修正参数,利用中心复合试验设计获取样本点,根据有限元计算结果与实测的结果残差为目标函数获得响应样本．拟合样本点和响应样本值构建二阶多项式响应面模型,并检验响应面模型的精度．基于贝叶斯理论更新融合系数来优化响应面参数,从而获得修正模型．以宁波大榭3号岸桥为工程背景,对比修正后的模态频率和实测频率,最大频率相对误差不超过5%,进而验证了基于贝叶斯信息融合的动力学有限元模型修正方法的有效性．修正后的有限元模型可进一步应用于岸桥的健康监测和安全评估．     

The concentration of stress and strain in finite thickness elastic plate containing a circular hole

The elastic stress and strain fields of finite thickness large plate containing a hole are systematically investigated using 3D finite element method. It is found that the stress and strain concentration factors of the finite thickness plate are different even if the plate is in elasticity state except at notch root of plate surface. The maximum stress and strain do not always occur on the mid plane of plate. They occur on the mid plane only in thin plate. The maximum stress and strain concentration factors are not on mid plane and the locations of maximum stress and strain concentration factors are different in thick plate. The maximum stress and strain concentration factors of notch root increase from their plane stress value to their peak values, then decrease gradually with increasing thickness and tend to each constant related to Poisson’s ratio of plate, respectively. The stress and strain concentration factors at notch root of plate surface are the same and are the monotonic descent functions of thickness. Their values decrease rapidly and tend to each constant related to Poisson’s ratio with plate thickness increasing. The difference between maximum and surface value of stress concentration factor is a monotonic ascent function of thickness. The thicker the plate is or the larger the Poisson’s ratio is, the larger the difference is. The corresponding difference of strain concentration factor is similar to the one of stress concentration factor. But the difference magnitude of stress concentration factor is larger than that of strain concentration factor in same plate.     

On the Measurement and Prediction of the Out-of-Plane Displacement Surrounding Cold-Expanded Holes

Experimental measurements of the out-of-plane displacement surrounding cold-expanded holes in a 6082-T6 aluminum alloy were made with a 3D optical scanner using the technique of the encoded light-pattern projection in white light. The measured surface profiles have shown the thickness discontinuity along the hole edge due to the effect of the split in the sleeve. An analytical–numerical solution of the out-of-plane displacement is presented based on existing analytical models. Thus, the results given by the analytical model were then compared with the experimental data and with a finite element (FE) model that simulates the cold-expansion process. The location of the elastic–plastic boundary was estimated as the point at which no change in thickness was observed; a good agreement was found in the comparison of measured, FE and analytical results. The measured surface profiles agreed with those predicted by the FE model and analytical solution. The proposed experimental approach can be used together with FE analysis for predicting the radial and circumferential residual stresses in cold-expanded hole. It is quite versatile and can also be used as quality-control technique in the manufacturing processes of cold-expanded holes.     

FRP-混凝土界面粘结行为的参数影响研究

FRP-混凝土界面的粘结性能对FRP加固混凝土结构力学行为和破坏模式有着重要影响。本文对表征FRP-混凝土界面粘结性能的三个重要参数(界面初始刚度、最大剪应力、界面破坏能)开展研究,通过13个单剪试件的试验考察了混凝土强度、胶层厚度和粘结长度等因素对界面粘结行为的影响,根据试验结果拟合了界面破坏能、最大剪切应力与胶层剪切刚度、混凝土强度之间的函数关系。在试验研究基础上,构建了外贴FRP-混凝土界面粘结的有限元模型。通过有限元分析考察了界面破坏能等三个参数不变的前提下,不同的局部粘结滑移本构关系对界面粘结行为的影响;进而研究了其中一个参数变化时引起的界面粘结性能改变。研究结果表明:界面粘结承载力随着胶层厚度增加而逐渐提高;胶层厚度与界面破坏能成正比,与峰值剪应力成反比;当界面破坏能等三个参数保持不变时,局部粘结滑移本构关系对FRP-混凝土界面粘结性能的影响较小;三个参数中的一个增大时将延缓界面破坏的过程。     

An Application of the Kane and Mindlin Theory to Crack Problems in Plates of Arbitrary Thickness

Classical plane solutions of the theory of elasticity, which are sometimes more than 100 years old, are still used today and provide a framework for the analysis of many practical problems. But, strictly speaking, these analytical solutions are only applicable to plates with vanishing thickness or infinite thickness, where the stress state could be classified as plane stress or plane strain, respectively. However, the through-the-thickness stresses that exist in a plate of given thickness have a significant impact in a number of practical applications; and these stresses are often inevitably ignored due to the lack of analytical tools. This paper presents new analytical results for crack tip opening displacement (CTOD) for the through-the-thickness crack in infinite plates with various thicknesses. These results are based on the solution for an edge dislocation in infinite plate of arbitrary thickness and an application of the distributed dislocation technique. The analytical predictions of the CTOD and the constraint factor are compared with the three-dimensional elasto-plastic finite element (FE) results. It is shown that both analytical and numerical results are in good agreement when the numerical calculations are not affected by the size of the FE mesh and by the boundaries of the FE model.     

金属蜂窝夹层结构抗水下爆炸特性

金属蜂窝夹层结构是一种新型的舰船防护结构,在舰船防护领域具有广阔的应用前景,但目前缺乏对其在实际水下爆炸载荷作用下动态响应的研究。为研究金属蜂窝夹层结构在水下爆炸载荷作用下的动态响应及防护性能,设计并制备了背板加筋蜂窝夹层结构样件以及相应的浮箱,在大型露天水池中进行了水下实爆 实验;通过声固耦合算法对结构响应进行模拟,实验结果与模拟结果吻合良好,随后分析了蜂窝夹层板的变形过程及能量吸收特性,量化了载荷参数（冲击因子）及结构参数（前后面板厚度比和芯体相对密度）对结构动态响应的影响;最后,以蜂窝夹层板的面密度和后面板中心点最大变形的无量纲量为目标函数,使用NSGA-Ⅱ遗传算法对结构进行了多目标优化,得到对应的Pareto前沿。结果表明,随着冲击因子的增大,蜂窝夹层板整体变形显著增大,蜂窝芯体始终是主要的吸能构件,但其吸能占比逐渐降低;随着前后面板厚度比或芯体相对密度的增加,蜂窝夹层结构的最大变形呈现先降低后升高的趋势,同时呈现不同的变形模式,芯体相对密度对结构变形的影响更为显著;对蜂窝夹层结构开展多目标优化可有效降低结构的面密度及最大变形,优化结果可为蜂窝夹层结构的设计选型提供参考。     

Mechanical Stimuli for Left Ventricular Growth During Pressure Overload

Background

The mechanical stimulus (i.e., stress or stretch) for growth occurring in the pressure-overloaded left ventricle (LV) is not exactly known.

Objective

To address this issue, we investigate the correlation between local ventricular growth (indexed by local wall thickness) and the local acute changes in mechanical stimuli after aortic banding.

Methods

LV geometric data were extracted from 3D echo measurements at baseline and 2 weeks in the aortic banding swine model (n?=?4). We developed and calibrated animal-specific finite element (FE) model of LV mechanics against pressure and volume waveforms measured at baseline. After simulation of the acute effects of pressure-overload, the local changes of maximum, mean and minimum myocardial stretches and stresses in three orthogonal material directions (i.e., fiber, sheet and sheet-normal) over a cardiac cycle were quantified. Correlation between mechanical quantities and the corresponding measured local changes in wall thickness was quantified using the Pearson correlation number (PCN) and Spearman rank correlation number (SCN).

Results

At 2 weeks after banding, the average septum thickness decreased from 10.6?±?2.92 mm to 9.49?±?2.02 mm, whereas the LV free-wall thickness increased from 8.69?±?1.64 mm to 9.4?±?1.22 mm. The FE results show strong correlation of growth with the changes in maximum fiber stress (PCN?=?0.5471, SCN?=?0.5111) and changes in the mean sheet-normal stress (PCN?=?0.5266, SCN?=?0.5256). Myocardial stretches, however, do not have good correlation with growth.

Conclusion

These results suggest that fiber stress is the mechanical stimuli for LV growth in pressure-overload.

     

Non-linear finite element analysis of cone penetration in layered sandy loam soil – Considering precompression stress state

Axisymmetric finite element (FE) method was developed to simulate cone penetration process in layered granular soil. The FE was modeled using ABAQUS/Explicit, a commercially available package. Soil was considered as a non-linear elastic plastic material which was modeled using variable elastic parameters of Young’s Modulus and Poisson’s ratio and Drucker–Prager criterion with yield stress dependent material hardening property. The material hardening parameters of the model were estimated from the USDA-ARS National Soil Dynamics Laboratory – Auburn University (NSDL-AU) soil compaction model. The stress–strain relationship in the NSDLAU compaction model was modified to account for the different soil moisture conditions and the influence of precompression stress states of the soil layers. A surface contact pair (‘slave-master’) algorithm in ABAQUS/Explicit was used to simulate the insertion of a rigid cone (RAX2 ABAQUS element) into deformable and layered soil medium (CAX4R ABAQUS element). The FE formulation was verified using cone penetration data collected on a soil chamber of Norfolk sandy loam soil which was prepared in two compaction treatments that varied in bulk density in the hardpan layer of (1) 1.64 Mg m⁻³ and (2) 1.71 Mg m⁻³. The FE model successfully simulated the trend of cone penetration in layered soils indicating the location of the sub-soil compacted (hardpan) layer and peak cone penetration resistance. Modification of the NSDL-AU model to account for the actual soil moisture content and inclusion of the influence of precompression stress into the strain behavior of the NSDL-AU model improved the performance of FE in predicting the peak cone penetration resistance. Modification of the NSDL-AU model resulted in an improvement of about 42% in the finite element-predicted soil cone penetration forces compared with the FE results that used the NSDL-AU ‘virgin’ model.     

A coupled IEM/FEM approach for solving elastic problems with multiple cracks

In this paper, the detailed two-dimensional infinite element method (IEM) formulation with infinite element (IE)–finite element (FE) coupling scheme for investigating mode I stress intensity factor in elastic problems with imbedded geometric singularities (e.g. cracks) is presented. The IE–FE coupling algorithm is also successfully extended to solve multiple crack problems. In this method, the domain of the primary problem is subdivided into two sub-domains modeled separately using the IEM for the multiple crack sub-domain, and the FEM for the uncracked sub-domain. In the IE sub-domain, the similarity partition concept together with the IEM formulation are employed to automatically generate a large number of infinitesimal elements, layer by layer, around the tip of each crack. All degrees of freedom related to the IE sub-domain, except for those associated with the coupling interface, are condensed and transformed to form a finite master IE for each crack with master nodes on sub-domain boundary only. All of the stiffness matrices constructed in the IE sub-domains are assembled into the system stiffness matrix for the FE sub-domain. The resultant FE solution with a symmetrical stiffness matrix, having the singularity effect of imbedded cracks in IEs, is required only for solving multiple crack problems.Using these efficient numerical techniques a very fine mesh pattern can be established around each crack tip without increasing the degree of freedom of the global FEM solution. One is easily allowed to conduct parametric analyses for various crack sizes without changing the FE mesh. Numerical examples are presented to show the performance of the proposed method and compared with the corresponding known results where available.     

Agricultural tire deformation in the 2D case by finite element methods

The mechanical characteristics of the rubber tire and the interaction between a tire and a rigid surface were investigated by a two-dimensional (2D) finite element (FE) model. The FE model consists of a rigid rim and a rigid contact surface which interact with the elastic tire. Four distinct sets of elastic parameters are used to represent beads, sidewall, tread and lugs. Several sets of tire loads and inflation pressures were applied to the FE model as boundary conditions, together with various displacements and friction conditions. The deformation of the tire profile, the tire displacements in the vertical and lateral directions, the normal contact pressures, the frictional forces and the stress distribution of the tire components were investigated by the 2D FE model under the above boundary conditions. The calculated tire deflections were compared with the measured data. The results show a good fit between calculated and measured data, especially at high load and inflation pressure. The comparison shows that the FE analysis is suitable to predict aspects of the tire performance like its deflection and interactions with the contact surface. Compared with the experimental methods, the FE methods show many advantages in the prediction of tire deformation, contact pressure and stress distribution.