Asymptotic field of mode Ⅱ dynamic growing crack in visco-elastic material

A mechanical model of visco-elastic material is established in order to investigate viscous effect in dynamic growing crack-tip field of mode Ⅱ. It is shown that in stable creep growing phase, elastic deformation and viscous deformation are equally dominant in the near-tip field, the stress and strain have the same singularity, namely, (σε) ∝r^-1/(n-1). The asymptotic solutions of separatied variables of stress, stain and displacement in crack-tip field are obtained by asymptotic analysis, and the results of numerical value of stress and strain in crack-tip field are obtained by shooting method. Through numerical calculation, it is shown that the near-tip fields are mainly governed by the creep exponent n and Mach number M. By the asymptotic analysis to the crack-tip field, the fracture criterion of mode Ⅱ dynamic growing crack of visco-elastic materials is put forward from the point of view of strain.     

刚性-压力敏感性材料界面I型准静态扩展裂纹尖端场

为了研究船用工程复合材料的界面裂纹特性,建立了刚性-压力敏感粘弹塑性材料Ⅰ型准静态扩展裂纹的力学模型.在稳态扩展阶段,应力和应变具有相同的奇异量级,即(σ,ε)∝γ<'-1/(n-1)>.引入Airy应力函数,通过渐近分析得出了裂纹尖端应力和应变的分离变量形式的渐近解,并采用打靶法求得了裂纹尖端应力和应变的数值结果.数值计算结果表明,界而裂尖场主要受材料的泊松比和幂硬化指数的控制.通过对裂纹尖端场的渐近分析,从应变角度出发,提出了刚性一压力敏感性材料界面Ⅰ型准静态扩展裂纹的断裂判据.     

The crack tip fields of Mode Ⅱ stationary growth crack on bimaterial interface

A mechanical model is established for mode Ⅱ interfacial crack static growing along an elastic-elastic power law creeping bimaterial interface. For frictional contact of boundary conditions on crack faces, asymptotic solutions of the stresses and strains of near tip-crack are got. It was shown that in stable creep growing phase, elastic deformation and viscous deformation are equally dominant in the near-tip field, the stress and strain have the same singularity and there is not the oscillatory singularity the field. Through numerical calculation, it is shown that the frictional coefficient η notably influence the crack-tip field.     

Self-propulsion computations using a speed controller and a discretized propeller with dynamic overset grids

A method that can be used to perform self-propulsion computations of surface ships is presented. The propeller is gridded as an overset object with a rotational velocity that is imposed by a speed controller, which finds the self-propulsion point when the ship reaches the target Froude number in a single transient computation. Dynamic overset grids are used to allow different dynamic groups to move independently, including the hull and appendages, the propeller, and the background (where the far-field boundary conditions are imposed). Predicted integral quantities include propeller rotational speed, propeller forces, and ship’s attitude, along with the complete flow field. The fluid flow is solved by employing a single-phase level set approach to model the free surface, along with a blended k−ω/k−ɛ based DES model for turbulence. Three ship hulls are evaluated: the single-propeller KVLCC1 tanker appended with a rudder, the twin propeller fully appended surface combatant model DTMB 5613, and the KCS container ship without a rudder, and the results are compared with experimental data obtained at the model scale. In the case of KCS, a more complete comparison with propulsion data is performed. It is shown that direct computation of self-propelled ships is feasible, and though very resource intensive, it provides a tool for obtaining vast flow detail.     

片体间距对高速三体船水动力性能影响的比较分析(英文)

This paper numerically investigates the influence of separation variation of the outriggers on the hydrodynamic performance of a high speed trimaran (HST) aiming at improving its applicability in diverse realistic disciplines. The present investigation was performed within the framework of the 2-D slender body method (SBM) by calculating the resistance of three symmetric trimaran series moving in a calm free surface of deep water. Each trimaran series comprises of 4681 configurations generated by considering 151 staggers (50%≤α≤+100%), and 31 separations (100%≤β≤400%) for 81 Froude numbers (0.20≤Fn≤1.0). In developing the three trimaran series, Wigley-st, AMECRC-09, and NPL-4a models were used separately for both the main and side hulls of each individual series models. A computer macro named Tri-PL was created using the Visual Basic for Applications. Tri-PL sequentially interfaced Maxsurf then Hullspeed to generate the models of the three trimaran series together with their detailed hydrostatic particulars, followed by their resistance components. The numerical results were partially validated against the available published numerical calculations and experimental results, to benchmark the Tri-PL macro and hence to rely on the analysis outcomes. A graph template was created within the framework of SigmaPlot to visualize the significant results of the Tri-PL properly.     

Influence of turbulence closure models on the vortical flow field around a submarine body undergoing steady drift

Manoeuvring underwater vehicles experience complex three-dimensional flow. Features include stagnation and boundary layer separation along a convex surface. The resulting free vortex sheet rolls up to form a pair of streamwise body vortices. The track and strength of the body vortex pair results in a nonlinear increase in lift as body incidence increases. Consequently, accurate capture of the body vortex pair is essential if the flow field around a manoeuvring submarine and the resulting hydrodynamic loading is to be correctly found. This work highlights the importance of both grid convergence and turbulence closure models (TCMs) to the strength and path of the crossflow-induced body vortices experienced by the DOR submarine model at an incidence angle of 15°. Five TCMs are considered; Spalart–Allmaras, k-ε, k-ω, shear stress transport, and the SSG Reynolds stress model. The SSG Reynolds stress model shows potential improvements in predicting both the path and strength of the body vortex over standard one- and two-equation TCMs based on an eddy viscosity approach.     

Asymptotic analysis of mode I propagating crack-tip field in a creeping material

Adopting an elastic-viscoplastic, the asymptotic problem of mode I propagating crack-tip field is investigated. Various asymptotic solutions resulting from the analysis of crack growing programs are presented. The analysis results show that the quasistatically growing crack solutions are the special case of the dynamic propagating solutions. Therefore these two asymptotic solutions can be unified.     

Asymptotic analysis of mode Ⅰ propagating crack-tip field in a creeping material

0　INTRODUCTIONWeknowthatthemetalmaterialwillbegincreepwhenthetemperatureishigherthan 30 per centoftheirsmeltingpoint.Foraturbine,adieselengine,anaerospacestructureandothermotivepowermechanisms,thematerialsstillworkwhentheirtemperaturesovershootthecreepthresholdval ue.Thecreepdeformationwilloccuratnormaltem peratureforthepolymerstructures(e.g .plasticspipingelements) .Inthecaseofcreepingthefrac tureofthestructureoccursduetocreepfractureandcreepcrackgrowth.Thereforetheresearchonthecrackgro…     

Strip method for a laterally drifting ship in waves

Lateral drift occurs due to the effects of wind forces, wave drifting forces, or both on ships sailing in actual seas. It is important therefore to investigate the influence of lateral drift on seakeeping performance for improved ship operation. The velocity potential was expanded as an asymptotic power series in terms of the lateral speed parameter, τ, defined as ω _e V ₀/g, where ω _e is the frequency of wave encounter; V ₀ denotes the lateral velocity, which is assumed to be sufficiently small; and g is the acceleration due to gravity. By combining this technique with the strip method, two sets of motion equations of all the hydrodynamic force coefficients for ship seakeeping were derived. The first set is for ships without lateral drift and is the same as the equations in the new strip method, and the second set is for the additional motions induced by lateral drift. It was found that all ship motion modes except surge are coupled when a ship drifts laterally in waves.     

舰船气泡尾流数值模拟(英文)

Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship (KCS) was first numerically simulated. Then with a rigid-lid-free-surface method, the underwater flow field was computed based on the mixture multiphase model to simulate the bubbly wake around the KCS hull. The realizable k-ε two-equation turbulence model and Reynolds stress model were used to analyze the effects of turbulence model on the ship bubbly wake. The air entrainment model, which is relative to the normal velocity gradient of the free surface, and the solving method were verified by the qualitatively reasonable computed results.     

深水环境下粘弹性复合材料夹层结构蠕变特性研究

文章结合单轴压缩蠕变试验和粘弹性材料广义Maxwell建模理论得到吸声和浮体两种粘弹性填充材料的松弛模量Prony级数系数,开展了深水环境下粘弹性复合材料夹层结构蠕变特性试验,并将之与仿真研究结果进行对比,得到如下结论：复合材料夹层结构长时蠕变初期,芯材蠕变较表层复合材料占主导地位;随着时间增加,表层复合材料基体的粘弹性特性所表现出的蠕变和松弛现象交织出现,表层应变出现波动;芯材蠕变对复合材料夹层结构长时蠕变变形的贡献约为60%。     

Numerical simulations of fatigue crack initiation and propagation based on re-tensile plastic zone generating load criterion for in-plane gusset welded joints

Many accidents are caused by fatigue in welded built-up steel structures, and so it is important to estimate the fatigue lives of such structures quantitatively for safety reasons. By assuming that fatigue cracks cannot grow without an accumulation of alternating tensile/compressional plastic strain, one of the authors identified an improved effective stress intensity factor range ΔK _RPG based on the re-tensile plastic zone generating (RPG) load, which represents the driving force for fatigue cracks, and suggested that ΔK _RPG should be used as the parameter to describe fatigue crack growth behavior. The “FLARP” numerical simulation code in which ΔK _RPG is implemented as the fatigue crack growth parameter, was developed in order to predict fatigue crack initiation and propagation behavior. In this paper, it is demonstrated that FLARP gives accurate estimates for fatigue life by comparing the estimated fatigue crack growth curves and S–N curves with the experimental results for in-plane gusset welded joints, which are used in many welded steel structures. Moreover, the effect of induced bending moment due to the linear misalignment in the out of plane direction on the fatigue strength of in-plane gusset welded joints is investigated through numerical simulations.     

Depth estimation of surface cracks based on the measurement of crack opening deformation and numerical–experimental iteration

We propose a new method to estimate the depth of a surface crack based on the measurement of crack opening deformation (COD) by using strain gauges. Through finite-element (FE) analysis of several surface cracks with different crack depth, it was found that the distribution of crack depth along the crack line can be approximated by multiplying a certain proportional α to the distribution of COD per unit nominal strain (normalized COD). The strain gauges are cemented just on the crack line and at a reference position, and the normalized CODs are measured under the impact load excited by hammer punching. The surface crack depth is estimated from the normalized COD measurements by a numerical–experimental iteration method based on FE analysis. The estimated distribution of the surface crack depth along the crack line shows good agreement with the shape of a real crack depth.     

Study of fatigue crack nucleation processes by crystalline FE analysis

 In order to accumulate knowledge about how material compositions and manufacturing methods affect fatigue strength, this paper investigates the relationship between the swiftness of the changes in the macroscopic indices of the crack initiation process and the hardening properties of a material. This is done by calculating the cyclic deformation behavior of a f.c.c. single crystal. The relationship between the swiftness of the changes and the crystal geometries is also examined by calculating the deformation behavior of a f.c.c. crystal with a high Schmid factor buried in a large crystal with a low Schmid factor. In addition, a multiscale hardening rule based on forest theory is developed in order to examine the microscopic mechanisms of fatigue slip band (persistent slip bands, PSBs) formation. The validity of the theories presented is examined by comparing the changing nature of the measured and calculated hysteresis loop shapes of f.c.c. single crystals. Inhomogenous slip deformation through the crystal is also investigated, and inferences are then drawn about the microscopic mechanisms of cyclic hardening and PSB formation. Received: August 5, 2002 / Accepted: December 18, 2002 Address correspondence to: N. Osawa (osawa@naoe.eng.osaka-u.ac.jp) Updated from the Japanese original, which won the 2002 SNAJ prize (J Soc Nav Archit Jpn 1998;184:351–363, 1999;185:283–292 and 186:535–544, 2001;190:539–551)     

URANS and DES analysis for a Wigley hull at extreme drift angles

Vortical structures and associated instabilities for flows around the Wigley hull for a wide range of drift angles (10° ≤ α ≤ 60°) with free surface are identified and analyzed. Quantitative verification and validation are conducted on three systematically refined grids with comparison to the experimental data for α = 10°. Analysis of the flow pattern shows a strong correlation between the vortical structures and free-surface wave elevation. For α = 10° and 30°, the flows remain steady and vortices are generated at the keel and fore and aft perpendiculars of the hull. The strength and complexity of these vortices increase with increasing α. At α = 45°, flow becomes unsteady without any significant change in the main flow pattern. At α = 60°, a complex and unsteady flow field on the leeward side of the hull is formed with a large recirculation region from the aft to the fore end, which prevents the flow coming from below the keel from moving up and generating the keel vortices observed at lower drift angles. Karman-like and helical instabilities are analyzed. The effect of Froude number is more apparent for large than for small drift angles.     

A review of the effect of a／W ratio on fracture toughness （Ⅰ）—experimental investigation in EPFM

Many experimental investigations have previously been performed and recently done on different shipbuilding structural steels where the specimens size and crack depth/specimen width (a/W) were varied. A series of interesting results have been gained. It is worthwhile to have a review on the effect of a/W ratio on fracture toughness, and further theoretical analysis is necessary. In this paper, experimental work in elasticplastic fracture mechanics (EPFM) was discussed. Tests had been carried out on 10 kinds of strength steels. Results showed that J _i and δ₁ values increased with decreasing a/W when a/W<0.3 for three-point bend specimens and that shallow crack specimens which have less constrained flow field give markedly higher values of toughness than deeply notched specimens. However, for a/W>0.3, the toughness was found to be independent of a/W. Slip line field analysis shows that for shallow cracks, the hydrostatic stress is lower than that from standard deeply cracked bend specimen which develops a high level of crack tip constraint, provides a lower bound estimate of toughness, and will ensure an unduly conservative approach when applied to structure defects especially if initiation values of COD/J-integral are used.     

Numerical simulation of interference effects for a high-speed catamaran

The simulations of the flow around a high-speed vessel in both catamaran and monohull configurations are carried out by the numerical solution of the Reynold averaged Navier–Stokes (RANS) equations. The goal of the analysis is the investigation of the interference phenomena between the two hulls, with focus on its dependence on the Reynolds number (Re). To this aim, numerical simulations are carried out for values of Re ranging from 10⁶ to 10⁸ for two different values of the Froude number (Fr = 0.30, 0.45). Wave patterns, wave profiles, limiting streamlines, surface pressure and velocity fields are analyzed; comparison is made between the catamaran and the monohull configurations. Dependence of the pressure and viscous resistance coefficients, as well as of the interference factor, on the Reynolds number is investigated. Verification and validation for both resistance coefficients and wave cuts is also performed.     

不同焊接顺序下加筋板变形及残余应力的三维有限元研究

针对焊接过程的二维有限元计算与实际情况存在一定差别的问题，使用三维热弹-塑性有限元法对不同焊接顺序下加筋板焊接过程进行了仿真，获得了加筋板焊接引起的温度场、位移场和应力场。结果表明：在4点约束条件下，加筋板横截面的变形为中垂变形，纵筋的变形为中拱变形，方案1横截面变形更小，方案2纵筋变形更小。焊接引起的加筋板残余应力主要表现为横向应力，其在近焊缝区为拉应力，达到材料屈服强度，远离焊缝区表现为压应力，达到0.2倍材料屈服强度。加强筋横向应力峰值出现在起弧端和收弧端，约为0.85倍材料屈服强度，纵向应力峰值出现在焊接起弧端，约为0.3倍材料屈服强度。在加筋板横截面位置，焊接顺序主要影响加强筋处的残余应力；在加强筋位置，焊接顺序主要影响纵向应力。每组焊缝同时同向焊接，且每根纵筋从左向右依次焊接的焊接方案产生更小的残余应力。     

450MPa级船用钢冲击实验研究及Cowper-Symonds本构模型

应用SHPB实验技术,对450MPa级船用钢进行冲击实验研究,建立了材料的Cowper-Symonds本构模型.冲击实验发现,450MPa级船用钢的动态流动应力在变形过程中,具有二次非线性强化规律、存在强化转折点.当材料的变形小于转折点的应变时,动态流动应力随应变的增大而增大;当材料的变形大于转折点的应变时,动态流动应力随应变的增大而减小.文中还对Cowper-Symonds本构模型的适用范围进行了讨论.     

花岗岩残积砂质黏性土的蠕变特性与本构模型

花岗岩残积砂质黏性土越来越多地被选作建(构)筑物天然地基持力层,但目前关于花岗岩残积砂质黏性土的研究主要停留在短期的强度和变形阶段。为探讨花岗岩残积砂质黏性土在长期荷载作用下的变形特性,对潮州市花岗岩残积砂质黏性土进行蠕变试验。结果表明,花岗岩残积砂质黏性土蠕变变形随时间的增加不断增大,但变形速率不断减小,直至变形稳定;随着剪应力及固结压力的增加,蠕变特性更加显著。引入邓肯-张模型对整数阶Poynting-Thomson流变模型进行修正,修正后的P-T模型可以较好地描述花岗岩残积砂质黏性土蠕变特性。研究结果有助于进一步认识花岗岩残积砂质黏性土的蠕变特性,为花岗岩残积砂质黏性土地基处理提供参考。