受载时船舶屈曲临界载荷数值的数学模型设计

船舶航行中船体会受到来自于水体的巨大冲击力,一旦超过屈曲载荷临界值就会增加船体断裂的风险。针对传统载荷控制模型对于环境参数变化不敏感、反应时滞过长的弊端,设计一种用户船舶屈曲临界值计算分析的数学模型。分析船舶行进中轴向速度和径向速度,并从速度与应力的比例关系中识别出船体受力综合载荷值,以此建立屈曲临界载荷数学模型。基于临界载荷数学模型分析不同载荷值条件下径向、轴向力的变化与临界误差之间的关系,进而准确地识别出屈曲临界载荷数据的变化规律。仿真结果表明,提出数学模型具有更高的载荷值传递效率,能够及时做出主动调整,并有效降低各种作用力对于船体的冲击。

“边坡-地基”的边坡失稳极限荷载与地基极限承载力研究

为研究"边坡-地基"所存在的地基承载力和边坡失稳问题,考虑不同基础宽度B与坡顶距De,建立了12个"边坡-地基"模型,采用2种解析法和2种有限元法,分析了"边坡-地基"的变形破坏模式,讨论了边坡失稳极限荷载和地基极限承载力。结果表明:"边坡-地基"的破坏模式随荷载的增加发生3种模式的变化;对于"边坡-地基"问题,传统水平地基承载力算法并不适用,且随着坡顶距增大,Saran解与其他解的误差也逐渐增大,基于塑性理论的解析法无法较好地全面考虑边坡对地基的影响问题;基于弹塑性理论的有限元强度折减法和有限元载荷数值试验方法得到的极限荷载基本一致,两者算得的极限荷载可以认为是"边坡-地基"严格意义上的极限承载力。

FINITE ELEMENT MODELING OF MODERATELY THICK COMPOSITE PLATE WITH PIEZOELECTRIC SENSORS AND ACTUATORS

A rectangular finite element for laminated plate with bonded and/or embedded piezoelectric sensors and actuators is developed based on the variational principle and the first order shear deformation theory. The element has four-node, 20-degrees-of-freedom with one potential degree of freedom for each piezoelectric layer to represent the piezoelectric behavior. The higher order derivation of deflection is obtained by using the normal rotation expressions to take the effects of transverse shear deformation into considerations. The finite element can accurately simulate the deformation of both thin and moderately thick plates. A Fortran program is written and a number of benchmark tests are exercised to verify its effectiveness. Results are compared well with the existing data. The unbalanced composite with piezoelectric layers is then analyzed by using the model. Results show that the changes of the ratio between the thickness of positive angle layers and the negative angle layers have an effect on the deformation of the structure under the same electric loading.

Analysis of granular assembly deformation using discrete element method

The discrete element method is used to simulate specimens under three different loading conditions（conventional triaxial compression,plane strain,and direct shear）with different initial conditions to explore the underlying mechanics of the specimen deformation from a microscale perspective.Deformations of specimens with different initial void ratios at different confining stresses under different loading conditions are studied.Results show that the discrete element models successfully capture the specimen deformation and the strain localization.Particle behaviors including particle rotation and displacement and the mesoscale void ratio distributions are used to explain the strain localization and specimen deformation.It is found that the loading condition is one of the most important factors controlling the specimen deformation mode.Microscale behavior of the granular soil is the driving mechanics of the macroscale deformation of the granular assembly.

Numerical simulation of dynamic fracture properties of rocks under different static stress conditions

When underground cavities are subjected to explosive stress waves,a uniquely damaged zone may appear due to the combined effect of dynamic loading and static pre-load stress.In this study,a rate-dependent two-dimensional rock dynamic constitutive model was established to investigate the dynamic fractures of rocks under different static stress conditions.The effects of the loading rate and peak amplitude of the blasting wave under different confining pressures and the vertical compressive coefficient(K_(0))were considered.The numerical simulated results reproduced the initiation and further propagation of primary radial crack fractures,which were in agreement with the experimental results.The dynamic loading rate,peak amplitude,static vertical compressive coefficient(K_(0))and confining pressure affected the evolution of fractures around the borehole.The heterogeneity parameter(m)plays an important role in the evolution of fractures around the borehole.The crack propagation path became more discontinuous and rougher in a smallerheterogeneity parameter case.

Numerical Simulation of Dynamic Response of an Existing Subway Station Subjected to Internal Blast Loading

In order to design and retrofit a subway station to resist an internal blast, the distribution of blast loading and its effects on structures should be investigated firstly. In this paper, the behavior of a typical subway station subjected to different internal blast Ioadings was analyzed. It briefly introduced the geometric characteristics and material constitutive model of an existing two-layer and three-span frame subway station. Then three cases of different explosive charges were consid- ered to analyze the dynamic responses of the structure. Finally, the maximum principal stress, dis- placement and velocity of the columns in the three cases were obtained and discussed. It con- cluded that the responses of the columns are sensitive to the charge of explosive and the distance from the detonation. It＇s also found that the stairs between the two layers have significant effects on the distribution of the maximum principal stress of the columns in the upper layer. The explicit dynamic nonlinear finite element software ANSYS/LS-DYNA was used in this study.

Cracks coalescence mechanism and cracks propagation paths in rock-like specimens containing pre-existing random cracks under compression

The mechanism of cracks propagation and cracks coalescence due to compressive loading of the brittle substances containing pre-existing cracks （flaws） was modeled experimentally using specially made rock-like specimens from Portland Pozzolana Cement （PPC）. The breakage process of the specimens was studied by inserting single and double flaws with different inclination angles at the center and applying uniaxial compressive stress at both ends of the specimen. The first crack was oriented at 50＆#176; from the horizontal direction and kept constant throughout the analysis while the orientation of the second crack was changed. It is experimentally observed that the wing cracks are produced at the first stage of loading and start their propagation toward the direction of uniaxial compressive loading. The secondary cracks may also be produced in form of quasi-coplanar and/or oblique cracks in a stable manner. The secondary cracks may eventually continue their propagation in the direction of maximum principle stress. These experimental works were also simulated numerically by a modified higher order displacement discontinuity method and the cracks propagation and cracks coalescence were studied based on Mode I and Mode II stress intensity factors （SIFs）. It is concluded that the wing cracks initiation stresses for the specimens change from 11.3 to 14.1 MPain the case of numerical simulations and from 7.3 to 13.8 MPa in the case of experimental works. It is observed that cracks coalescence stresses change from 21.8 to 25.3 MPa and from 19.5 to 21.8 MPa in the numerical and experimental analyses, respectively. Comparing some of the numerical and experimental results with those recently cited in the literature validates the results obtained by the proposed study. Finally, a numerical simulation was accomplished to study the effect of confining pressure on the crack propagation process, showing that the SIFs increase and the crack initiation angles change in this case.

Numerical Simulation of Dynamic Response and Collapse for Steel Frame Structures Subjected to Blast Load

The progressive collapse of steel frame structures under the blast load was investigated using LS-DYNA. The multi-material Eulerian and Lagrangian coupling algorithm was adopted. A flu-id-structure coupling finite element model was established which consists of Lagrange element for simulating steel frame structures and concrete ground, multiple ALE element for simulating air and TNT explosive material. Numerical simulations of the blast pressure wave propagation, struc-tural dynamic responses and deformation, and progressive collapse of a five-story steel frame structure in the event of an explosion near above ground were performed. The numerical analysis showed that the Lagrangian and Eulerian coupling algorithm gave good simulations of the shock wave propagation in the mediums and blast load effects on the structure. The columns subjected to blast load may collapse by shear yielding rather than by flexural deformation. The columns and joints of steel beam to column in the front steel frame structure generated enormous plastic defor-mation subjected to intensive blast waves, and columns lost carrying capacity, subsequently lead-ing to the collapse of the whole structure. The approach coupling influence between struc-tural deformation and fluid load well simulated the progressive collapse process of structures, and provided an effective tool for analyzing the collapse mechanism of the steel frame structure under blast load.

A stability study of goaf based on mechanical properties degradation of rock caused by rheological and disturbing loads

Based on the classical static theory and static numerical simulation,the static method could not accurately reflect the stability of goaf where the rocks on the pillar and roof are influenced by Theological and blasting disturbance for a long time.According to the test from the site,an experimental study was made in Theological and dynamic disturbance.After that,on the basis of variable rock mechanics parameters from the experimental data,numerical simulation was used to analyze the vertical stress distribution of goaf,vertical displacement and plastic area of roof in the "deterioration" caused by Theological and blasting,which shows that the mechanics properties of the rock were greatly influenced by Theological,and dynamic disturbance.The results of the experimental study and numerical simulation show that the mechanics properties of rock are greatly influenced by Theological and dynamic disturbance.As a result,the stability of goaf is greatly reduced.Finally,by comparing golf monitoring results with the analysis of theoretical calculation,it was found that the results were approximately the same,which testifies the reliability of the method.This method provides a new way of studying the stability of goaf as well as laying a basic foundation for future safety management.

Effect of vertical load difference on cracking behaviors in multistory masonry buildings and numerical simulation

To investigate the causes qf cracks in multistory masonry buildings, the effect of vertical load difference on cracking behaviors was investigated experimentally by testing and measuring the displacements at the testing points of a large sized real masonry U-shaped model. Additionally, the cracking behaviors in U-shaped model were analyzed with shear stress and numerical simulated with ANSYS software. The experimental results show that the deformation increases with the increase of the vertical load. The vertical load results in different deformation between the bearing wall and non-bearing wall, which leads to cracking on the non-beating wall. The rapid deformation happens at 160 kN and cracks occur firstly at the top section of non-bearing wall near to the bearing wall. New cracks are observed and the previous cracks are enlarged and developed with the increase of vertical load. The maximum crack opening reaches 12 mm, and the non-bearing wall is about to collapse when the vertical load arrives at 380 kN. Theoretical analysis indicates that the shear stress reaches the maximum value at the top section of the non-bearing wall, and thus cracks tend to happen at the top section of the non-bearing wall. Numerical simulation results about the cracking behaviors are in good agreement with experiments results.

Pull-out tests and slope stability analyses of nailing systems comprising single and multi rebars with grouted cement

The pull-out capacities for soil nailing systems comprising of one single 29 mm diameter(type A) and four 16 mm diameter(type B) rebars with grouted cement were examined.A field test and numerical analysis for the type A and type B systems were carried out to investigate the pull-out capacities and the slope stability reinforcement efficiency in soil and rock slopes.The results of the pull-out tests show the mobilized shear force and load transfer characteristics with respect to soil depth.The load-displacement relationship was examined for both type A and type B systems.Slope stability analyses were carried out to study the relationships between soil and nail reinforcement and bending stiffness as well as combined axial tension and shear forces.Factors of safety were calculated in relation to the number of nails and their outside diameters.Both soil and rock slopes were included in this evaluation.

Wave-induced seepage and its possible contribution to the formation of pockmarks in the Huanghe(Yellow) River delta

Wave-induced seepage and its possible contribution to the formation of pockmarks in the Huanghe(Yellow) River delta were investigated experimentally and numerically. Laboratory experiments were carried out to explore the response of a layered silty seabed with various saturation conditions under cyclic wave loads,in which the pore pressure and seepage-related phenomena were particularly monitored. Numerical models to simulate wave-induced seepage in the seabed were presented and evaluated,then applied to the Huanghe River delta. The experimental results show that the excess pore pressure decreases more rapidly at the surface layer,while the seepage-related phenomena are more pronounced when large cyclic loads are applied and the underlying layer is less saturated. The proposed numerical models were verified by comparing with the experiments. The calculated seepage depth agreed well with the depth of the pockmarks in the Huanghe River delta. The experimental and numerical results and the existing insitu investigations indicate that the wave-induced seepage may be a direct cause of the pockmarks in the Huanghe River delta. Extreme storm waves and the dual-layered structure of hard surface layer and weak underlying layer are essential external and internal factors,respectively. Wave- or current-induced scour and transport are possible contributors to the reformation of pockmarks at a later stage.

Numerical Analysis of Masonry Bell-Towers under Dynamic Loading

The seismic behavior of masonry bell-towers located in Corfu, Greece, is considered since constructions of this kind are particularly vulnerable to seismic loadings. Dynamics analyses for different total heights of Corfu bell-towers and according to the seismic data of this region are carried out. A diagram of displacements is determined. The impact of the percentages of damping in the distribution of displacements is examined.

Evaluation of coal pillar loads during longwall extraction using the numerical method and its application

It is very difficult to reasonably evaluate the loads acting on coal pillars in longwall panels during the planning of a new pillar system. The application of empirical equations is a common practice in calculating coal pillar loads while designing a new pillar. This paper proposes numerical models for evaluating coal pillar loads. The key of building a successful numerical model for calculating coal pillar loads lies in the fact that the model should represent the redistribution of stress all over the longwall panels and the surrounding areas, and it is especially important to include the characteristics of the stress rebuilding process in the gob areas, which are crucial for the building process of coal pillar loads. Based on the geo-mechanical background of the Baoshan Coal Mine, this paper details the procedures of applying numerical models to the evaluation of coal pillar loads and their local practices. The study results show it is feasible and reasonable to use numerical models to evaluate coal pillar loads.

Response of air gap and slamming loadings on semi-submersible platform

The present research deals with the numerical prediction of the air gap within the 6th generation of deepwater drilling floating semi-submersible platform and the experimental studies on the slamming loadings onto the structure. The survivability of the floating model with a mooring system was tested under extreme wave of 10-year return period. In the numerical simulation of the Gaussian method,the narrow band model was applied to obtain the first-order wave surface equation and the modified second-order wave surface equation. The hydrodynamic responses of the floating body,i.e. radiation damping,added mass,second-order wave excitation force and drifting force,were computed by using the potential flow theory based on higher order boundary element method in frequent domain. In the experimental analysis,high-frequency sensors were installed at the lower deck to measure the wave slamming loads. Equivalent truncated mooring system was applied to make sure position of the floating body in the wave tank. The comparison between the numerical and experimental results showed the numerical model underestimated the air gap of the floating body. Nevertheless,the predictions of the high risk spots underneath the floating deck that is prone to wave slamming obtained from both models were agreeable to each other. The experimental results also revealed that the wave slamming events often occurred at the connection point between the rear columns and the lower deck.

Assessment of Numerical Models for Live Load Distribution in a Road Slab Bridge

This paper presents comparison of numerical models used in an analysis of a road bridge deck. The models were adapted for computing the live load distribution coefficients in composite concrete bridge deck. The load distribution method was chosen for assessment of the usability of different numerical model in slab bridge deck analysis. The goal of the study is to determine a simplest but still accurate numerical model to estimate live load effects on composite slab bridge. In the analysis, the well-established grillage approach was adapted for representation of the bridge deck as a basic model as well as more sophisticated three-dimensional models which was supposed to better represent the real behavior of the deck under concentrated wheel loads. The bridge deck was effectively modeled using beam and shell elements. The grillage method compares well with the finite-element method. This finding is allowed to establish simplification in numerical modeling of slab bridge decks for live load effect computations.

Determination of T-stress using finite element analysis

A method that uses finite element analysis to determine the non-singular stress （T-stress） at a crack tip is proposed in this study. T-stress includes two components： the Tx-stress parallel to the tangent of the crack at its tip and the Ty-stress perpendicular to this tangent. The effects of contact and friction on both the Tx- and Ty-stresses on the crack flanks are considered in the method. Because the method uses a single standard elastic finite element analysis derived directly from the equation of the stress fields around the crack tip and does not require any assumptions or simplification, it can be used to determine the T-stress for any given geometry and loading condition. Theoretical results are used to calibrate the results, which exhibited good agreement and to discuss the T-stress computational methodology. Furthermore, the Tx- and Ty-stresses in center-cracked Brazilian disc （CCBD） specimens subjected to diametrical or partially distributed compression were numerically computed, and the effects of contact and friction on the Tx- and Ty-stresses are discussed.

Numerical simulation of non-Gaussian wind load

Compared with Gaussian wind loads, there is a higher probability of strong suction fluctuations occurrence for non-Gaussian wind pressures. These instantaneous and intermittent fluctuations are the initial cause of local damage to roof structures, par- ticularly at the edges and comers of long-span roofs. Thus, comparative errors would occur if a Gaussian model is used to de- scribe a non-Gaussian wind load, and structural security would not be guaranteed. This paper presents a simplified method based on the inverse fast Fourier transform （IFFT）, in which the amplitude spectrum is established via a target power spectrum. Also, the phase spectrum is constructed by introducing the exponential peak generation （EPG） model. Finally, a random pro- cess can be generated via IFFT that meets the specified power spectral density （PSD）, skewness and kurtosis. In contrast to a wind tunnel experiment, this method can avoid the coupled relation between the non-Gaussian and the power spectrum char- acteristics, and lead to the desired computational efficiency. Its fitting accuracy is not affected by phase spectrum. Moreover, the fitting precision of the kurtosis and PSD parameters can be guaranteed. In a few cases, the fitting precision of the skewness parameter is fairly poor, but kurtosis is more important than skewness in the description of the non-Gaussian characteristics. Above all, this algorithm is simple and stable and would be an effective method to simulate a non-Gaussian signal.

Unsteady Shock-Flow Characteristics in an Over-Expanded Rocket Nozzle

A numerical investigation of transient side-loads in an axisymmetric over-expanded thrust optimized contour nozzle is presented.These nozzles experience side-loads during start-up and shut-down operations,because of the flow separation at nozzle walls.Two types of flow separations such as FSS and RSS shock structure occur.A two-dimension numerical simulation has been carried out over an axisymmetric TOC nozzle to validate present results and investigate oscillatory flow characteristics for start-up processes.Reynolds Averaged Navier-Stokes equations are numerically solved using a fully implicit finite volume scheme.Governing equations are solved by coupled implicit scheme.Reynolds Stress turbulence model is selected.Present computed pressure at the nozzle wall closely matched with experiment data.A hysteresis phenomenon has been observed between these two shock structures.The transition from FSS to RSS pattern during start-up process has shown maximum nozzle wall pressure.Nozzle wall pressure and shear stress values have shown fluctuations during the FSS to RSS transition. The oscillatory pressure has been observed on the nozzle wall for high pressure ratio.Present results have shown that magnitude of the nozzle wall pressure variation is high for the oscillatory phenomenon.

Numerical study on explosion-induced fractures of reinforced concrete structure by beam-particle model

In the field of disaster prevention mitigation and protection engineering,it is important to identify the mechanical behaviors of reinforced concrete(RC)under explosive load by simulation.A three dimensional beam-particle model(BPM),which is suitable to simulate the fracture process of RC under explosive load,has been developed in the frame of discrete element method (DEM).In this model,only the elastic deformations of beams between concrete particles were considered.The matrix displacement method(MDM)was employed to describe the relationship between the deformation and forces of the beam.A fracture criterion expressed by stress was suggested to identify the state of the beam.A BPM for steel bar,which can simulate the deformation of steel bar under high loading rate,was also developed based on the Cowper-Symonds theory.A program has been coded using C++language.Experiments of RC slab under explosive load were carried out using the program.Good agreement was achieved between the experimental and simulated results.It is indicated that the proposed theoretical model can well simulate the fracture characteristics of RC slab under explosive load such as blasting pit formation,cracks extension, spallation formation,etc.