Mixed-mode I/II fields around a crack with a cohesive zone ahead of the crack tip

Considering a cohesive zone ahead of the crack tip, the linear elastic crack problem under Mode I, Mode II or mixed-mode conditions is formulated in an elliptic coordinate system, so that the cohesive surfaces are conveniently represented by straight line segments. It is shown that the displacement and stress fields around the crack tip and the cohesive zone, expressed in terms of elliptic coordinates, have a simple mathematical form, which does not contain a stress singularity at the crack tip due to the existence of the cohesive zone.     

Mode-3 spontaneous crack propagation along functionally graded bimaterial interfaces

The effects of combining functionally graded materials (FGMs) of different inhomogeneous property gradients on the mode-3 propagation characteristics of an interfacial crack are numerically investigated. Spontaneous interfacial crack propagation simulations were performed using the newly developed spectral scheme. The numerical scheme derived and implemented in the present work can efficiently simulate planar crack propagation along functionally graded bimaterial interfaces. The material property inhomogeneity was assumed to be in the direction normal to the interface. Various bimaterial combinations were simulated by varying the material property inhomogeneity length scale. Our parametric study showed that the inclusion of a softening type FGM in the bimaterial system leads to a reduction in the fracture resistance indicated by the increase in crack propagation velocity and power absorbed. An opposite trend of increased fracture resistance was predicted when a hardening material was included in the bimaterial system. The cohesive tractions and crack opening displacements were altered due to the material property inhomogeneity, but the stresses ahead of the cohesive zone remained unaffected.     

数学网格和物理网格分离的有限单元法(II):粘聚裂纹扩展问题中的应用

强化有限单元法将物理网格与数学网格分离开来,可以方便地描述非连续变形;粘聚区域模型是模拟断裂过程区作用最简单有效的方法,且可以避免裂纹尖端的应力奇异性.本文以平面问题为例,将强化有限单元法与粘聚区域模型相结合,利用富集数学节点描述任意粘聚裂纹扩展过程中的非连续变形问题,提出了裂纹扩展过程中数学节点富集和数学单元定义的方法.本文还导出了与平面4～8节点平面等参单元对应的8～16节点粘聚裂纹单元列武.最后,通过三点弯梁的裂纹扩展过程模拟验证了本文提出的粘聚裂纹扩展模拟方法的有效性.     

用内聚力模型研究PBGA封装焊点/IMC界面的分层开裂

塑料球栅阵列封装PBGA的可靠性分析中,考虑封装过程中SnAgCu焊料与铜焊盘界面间产生的金属间化合物(Intermetallic compound,IMC)的影响,并引入内聚力模型(Cohesive zone model,CZM),利用ANSYS对热循环作用下焊点/IMC界面的脱层开裂情况进行研究.结果表明:热循环作用下,在封装器件中焊点承受较大的应力应变,且远离中心的外侧焊点具有比内侧焊点更大的应力应变.IMC的存在极大的降低了焊点的可靠性.界面分层最先发生在最外侧的IMC/焊点界面的两端,随着热循环次数的增加,分层逐渐沿着界面两端向里扩展.在热循环的前几个阶段,各个界面的最大损伤值增大较快,随着热循环的继续加载,界面最大损伤值逐渐趋于稳定.整个过程中四号焊点界面的损伤值始终最大.     

A mode split,Godunov‐type model for nonhydrostatic,free surface flow

A previously developed model for nonhydrostatic, free surface flow is redesigned to improve computational efficiency without sacrificing accuracy. Both models solve the Reynolds averaged Navier–Stokes equations in a fractional step manner with the pressure split into hydrostatic and nonhydrostatic components. The hydrostatic equations are first solved with an approximate Riemann solver. The hydrostatic solution is then corrected by including the nonhydrostatic pressure and requiring the velocity field to obey the incompressibility constraint. The original model requires the solution of a Riemann problem at every cell face for each vertical layer of cells, which is computationally expensive. The redesigned model instead solves the shallow water (long wave) equations for the hydrostatic solution. Vertical shear is computed by subtracting the shallow water equations from the full three dimensional equations, which removes the hydrostatic thrust terms. Therefore, the required fluxes may be more efficiently computed with velocity based upwind differencing rather than solving a Riemann problem in each vertical layer of cells. This approach is termed mode splitting and has been used in hydrostatic coastal and ocean circulation models, but not surf zone models. Numerical predictions are compared with analytical solutions and experimental data to show that the mode split model is as accurate as the original model, but requires significantly less computational effort especially for large numbers of cell layers. Published 2014. This article is a U.S. Government work and is in the public domain in the USA.     

Two and three dimensional modeling of fluidized bed with multiple jets in a DEM–CFD framework

Fluidized beds with multiple jets have widespread industrial applications. The objective of this paper is to investigate the jet interactions and hydrodynamics of a fluidized bed with multiple jets. Discrete element modeling coupled with in-house CFD code GenIDLEST has been used to simulate a bed with nine jets. The results are compared with published experiments. Mono dispersed particles of size 550 μm are used with 1.4 times the minimum fluidization velocity of the particles. Both two and three dimensional computations have been performed. To the best of our knowledge, the results presented in this paper are the first full 3D simulations of a fluidized bed performed with multiple jets. Discrepancies between the experiment and simulations are discussed in the context of the dimensionality of the simulations. The 2D solid fraction profile compares well with the experiment close to the distributor plate. At higher heights, the 2D simulation over-predicts the solid fraction profiles near the walls. The 3D simulation on the other hand is better able to capture the solid fraction profile higher up in the bed compared to that near the distributor plate. Similarly, the normalized particle velocities and the particle fluxes compare well with the experiment closer to the distributor plate for the 2D simulation and the freeboard for the 3D simulation, respectively. A lower expanded bed height is predicted in the 2D simulation compared to the 3D simulation and the experiment. The results obtained from DEM computations show that a 2D simulation can be used to capture essential jetting trends near the distributor plate regions, whereas a full scale 3D simulation is needed to capture the bubbles near the freeboard regions. These serve as validations for the experiment and help us understand the complex jet interaction and solid circulation patterns in a multiple jet fluidized bed system.     

贝雷砂岩射孔孔道流动效应正交计算分析

对射孔孔道流动效应进行了量纲分析,建立了轴向和径向流量与影响流量因素间的工程表达式。综合考虑影响流动效应的主要因素和取值范围,设计了数值模拟计算工况的正交计算方案,并利用ANSYS/FLOTRAN流体程序进行了带损伤区贝雷砂岩流量数值模拟研究,获得了各个计算方案的径向流量和轴向流量,用以确定出流量公式的系数。最终确定的流量公式表明,各要素对流量影响的重要顺序依次为固有渗透率、损害区厚度、孔深、渗透率剩余百分数和孔径。     

东海区海岸带景观格局变化对生态系统服务价值的影响

对大尺度、长时序的海岸带景观格局和生态系统服务价值(ecosystem services value,ESV)的研究,有利于更好地促进区域生态系统保护工作的开展。以1990,1995,2000,2005,2010,2015年6期遥感影像解译结果为基础,对东海区海岸带景观格局变化特征和ESV时空演变进行了分析,并探讨了ESV与景观格局变化二者间的关系。结果表明,1990—2015年的25 a间,受人类活动的影响,东海区海岸带农田和海域面积大量减少,建设用地面积大幅增加;在所研究的5类地貌类型中,除中起伏山地外,其余地貌类型景观面积转移都较为剧烈。25 a间东海区海岸带斑块数量、斑块密度、边界密度、shannon多样性指数等不断增加,景观破碎化、复杂化、多样化增强。农田景观的最大斑块指数、形态指数、边界密度等变化均较大。东海区海岸带ESV呈先增加后减小的趋势,其中森林景观对ESV总量贡献最大;水文调节、生物多样性维持、气候调节、土壤保持、废物处理和气体调节是东海区海岸带主要的生态系统服务功能;在空间上表现出以单位面积中、低ESV为主导的空间分布,高、极高和极低ESV分布较少。研究区ESV总量与农田和海域以及水域和建设用地面积显著相关;与斑块密度、边界密度和形态指数呈显著负相关。同时,食物生产、气体调节、气候调节、土壤保持和生物多样性维持的服务价值与景观指数间均显著相关。     

Effect of hot ion temperature on obliquely propagating ion-acoustic solitons and double layers in an auroral plasma

Properties of obliquely propagating ion-acoustic solitons and double layers in a magnetized auroral plasma composed of hot adiabatic ions and two types of, cool and hot Maxwellian electrons are studied using Sagdeev pseudo-potential technique and assuming the quasi-neutrality condition. The new and surprising result which emerges from the model is that in contrast to the case of cold ions where ion-acoustic solitons and double layers are found for subsonic Mach numbers only, the hot ions case allows these nonlinear structures to exist for both subsonic and supersonic Mach number regimes. The double layers exist at lower angle of propagation as hot ion temperature is increased. The soliton electric field amplitudes are increased but their width and pulse duration are decreased with the increase in hot ion temperature. For the auroral zone parameters, the maximum electric field amplitude, width, pulse duration and speed for the solitons come out to be in the range ∼ (0.3–15) mV/m, ∼ (195–455) m, (7–20) ms and (22–26) km/s, respectively. The results seem to be in agreement with the Viking satellite observations in the auroral zone.     

Thermoelastoplastic deformation of a thick-walled cylinder with a radial crack

The thermoelastoplastic fracture mechanics problem of a thick-walled cylinder subjected to internal pressure and a nonuniform temperature field is solved by the method of elastic solutions combined with the finite-element method. The correctness of the solution is provided by using the Barenblatt crack model, in which the stress and strain fields are regular. The elastoplastic problem of a cracked cylinder subjected to internal pressure and a nonuniform temperature field are solved. The calculation results are compared with available data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 173–183, May–June, 2008.