A Thermal-Solid–Fluid Method for Topology Optimization of Structures with Design-Dependent Pressure Load

For the topology optimization of structures with design-dependent pressure,an intuitive way is to directly describe the loading boundary of the structure,and then update the load on it.However,boundary recognition is usually cumbersome and inaccurate.Furthermore,the pressure is always loaded either outside or inside the structure,instead of both.Hence,the inner enclosed and outer open spaces should be distinguished to recognize the loading surfaces.To handle the above issues,a thermal-solid–fluid method for topology optimization with design-dependent pressure load is proposed in this paper.In this method,the specific void phase is defined to be an incompressible hydrostatic fluid,through which the pressure load can be transferred without any needs for special loading surface recognition.The nonlinear-virtual thermal method(N-VTM)is used to distinguish the enclosed and open voids by the temperature difference between the enclosed(with higher temperature)and open(with lower temperature)voids,where the solid areas are treated as the thermal insulation material,and other areas are filled with the self-heating highly thermally conductive material.The mixed displacement–pressure formulation is used to model this solid–fluid problem.The method is easily implemented in the standard density approach and its effectiveness is verified and illustrated by several typical examples at the end of the paper.

Topology optimization of 3D structures with design-dependent loads

Topology optimization of continuum structures with design-dependent loads has long been a challenge. In this paper, the topology optimization of 3D structures subjected to design-dependent loads is investigated. A boundary search scheme is proposed for 3D problems, by means of which the load surface can be identified effectively and efficiently, and the difficulties arising in other approaches can be overcome. The load surfaces are made up of the boundaries of finite elements and the loads can be directly applied to corresponding element nodes, which leads to great convenience in the application of this method. Finally, the effectiveness and efficiency of the proposed method is validated by several numerical examples.

LOADING CAPACITY ANALYSIS OF THE MISALIGNED CONICAL-CYLINDRICAL BEARING WITH NON NEWTONIAN LUBRICANTS

A novel numerical method to lubricate a conventional finite diameterconical-cylindrical bearing with a non-Newtonian lubricant, while adhering to the power-law model,is presented. The elastic deformation of bearing and varied viscosity of lubrication due to thepressure distribution of film thickness are also considered. Simulation results indicate that thenormal load carrying capacity is more pronounced for higher values of flow behavior index n, highereccentricity ratios and larger misalignment factors. It is found that the viscosity-pressure to theeffect of lubricant viscosity is significant.

Numerical Analyses of Caisson Breakwaters on Soft Foundations Under Wave Cyclic Loading

A caisson breakwater is built on soft foundations after replacing the upper soft layer with sand. This paper presents a dynamic finite element method to investigate the strength degradation and associated pore pressure development of the intercalated soft layer under wave cyclic loading. By combining the undrained shear strength with the empirical formula of overconsolidation clay produced by unloading and the development model of pore pressure, the dynamic degradation law that describes the undrained shear strength as a function of cycle number and stress level is derived. Based on the proposed dynamic degradation law and M-C yield criterion, a dynamic finite element method is numerically implemented to predict changes in undrained shear strength of the intercalated soft layer by using the general-purpose FEM software ABAQUS, and the accuracy of the method is verified. The effects of cycle number and amplitude of the wave force on the degradation of the undrained shear strength of the intercalated soft layer and the associated excess pore pressure response are investigated by analyzing an overall distribution and three typical sections underneath the breakwater. By comparing the undrained shear strength distributions obtained by the static method and the quasi-static method with the undrained shear strength distributions obtained by the dynamic finite element method in the three typical sections, the superiority of the dynamic finite element method in predicting changes in undrained shear strength is demonstrated.

Fatigue Properties of Plain Concrete Under Triaxial Compressive Cyclic Loading

Experiments were made on plain concrete subjected to triaxial static loading and constant amplitude compressive fatigue loading with a constant lateral pressure in two directions. The initial confining pressure was 0, 0.1 f c , 0.25 f c and 0.4 f c , respectively, for the static test, and 0.1 f c and 0.25 f c for the fatigue test. Based on the triaxial compressive constitutive behavior of concrete, the inflexion of confining pressure evolution was chosen to be the fatigue damage criterion during the test. The rule of evolution of longitudinal maximum and minimum strains, longitudinal cyclic modulus and damage were recorded and analyzed. According to the Fardis Chen criterion model and the concept of equivalent fatigue life and equivalent stress level, a unified S N curve for multi axial compressive fatigue loading was proposed. Thus, the fatigue strength factors for different fatigue loading cases can be obtained. The present investigation provides information for the fatigue design of concrete structures.

实现柔性机翼后缘形状变化的综合优化(英文)

Adaptive wings have long used smooth morphing technique of compliant leading and trailing edge to improve their aerodynamic characteristics. This paper introduces a systematic approach to design compliant structures to carry out required shape changes under distributed pressure loads. In order to minimize the deviation of the deformed shape from the target shape, this method uses MATLAB and ANSYS to optimize the distributed compliant mechanisms by way of the ground approach and genetic algorithm （GA） to remove the elements possessive of very low stresses. In the optimization process, many factors should be considered such as airloads, input displacements, and geometric nonlinearities. Direct search method is used to locally optimize the dimension and input displacement after the GA optimization. The resultant structure could make its shape change from 0 to 9.3 degrees. The experimental data of the model confirms the feasibility of this approach.

海上水下爆炸的爆炸压力测量

Blast pressure measurements of a controlled underwater explosion in the sea were carried out.An explosive of 25-kg trinitro-toluene(TNT)equivalent was detonated,and the blast pressures were recorded by eight diferent high-performance pressure sensors that work at the nonresonant high-voltage output in adverse underwater conditions.Recorded peak pressure values are used to establish a relationship in the well-known form of empirical underwater explosion(UNDEX)loading formula.Constants of the formula are redetermined by employing the least-squares method in two diferent forms for best ftting to the measured data.The newly determined constants are found to be only slightly diferent from the generally accepted ones.

Experimental Study on Seepage Characteristics of a Soil-Rock Mixture in a Fault Zone

A mixture of fault gouge and rubble taken out from a fault zone is used to prepare a S-RM(Soil-Rock Mixture)sample with rock block proportions of 20%,30%,40%,50%,60%and 70%,respectively.A GDS triaxial test system is used accordingly to measure the seepage characteristics of such samples under different loading and unloading confining pressures in order to determine the variation law of the permeability coefficient.The test results show that:(1)The permeability coefficient of the S-RM samples decreases as the pressure increases,and the decrease rate of this coefficient in the initial stage of confining pressure loading is obviously higher than in the semi-late period;(2)The permeability coefficient at different confining pressure levels presents a common trend as the rock block proportion is increased,i.e.,it decreases first then it increases(the permeability coefficient of the sample with rock block proportion 40%being the smallest,70%the largest);(3)In the stage of confining pressure unloading,the recovery degree of the permeability coefficient grows with the increase of rock block proportion(the recovery rate of S-RM sample with rock block proportion 70%reaches 50.2%);(4)In the stage of confining pressure loading and unloading,the sensitivity of the permeability coefficient to the rock block proportion displays the inverse“Z”variation rule(when rock block proportion reaches 60%,the sensitivity is highest);(5)In the stage of confining pressure loading,the relationship between the permeability coefficient and confining pressure can be described by an exponential relationship.

Dynamic rock tests using split Hopkinson （Kolsky） bar system - A review

Dynamic properties of rocks are important in a variety of rock mechanics and rock engineering problems. Due to the transient nature of the loading, dynamic tests of rock materials are very different from and much more challenging than their static counterparts. Dynamic tests are usually conducted using the split Hopkinson bar or Kolsl^j bar systems, which include both split Hopkinson pressure bar （SHPB） and split Hopkinson tension bar （SHTB） systems. Significant progress has been made on the quantification of various rock dynamic properties, owing to the advances in the experimental techniques of SHPB system. This review aims to fully describe and critically assess the detailed procedures and principles of tech- niques for dynamic rock tests using split Hopkinson bars. The history and principles of SHPB are outlined, followed by the key loading techniques that are useful for dynamic rock tests with SHPB （i.e. pulse shaping, momentum-trap and multi-axial loading techniques）. Various measurement techniques for rock tests in SHPB （i.e. X-ray micro computed tomography （CT）, laser gap gauge （LGG）, digital image corre- lation （DIC）, Moir~ method, caustics method, photoelastic coating method, dynamic infrared thermog- raphy） are then discussed. As the main objective of the review, various dynamic measurement techniques for rocks using SHPB are described, including dynamic rock strength measurements （i.e. dynamic compression, tension, bending and shear tests）, dynamic fracture measurements （i.e. dynamic imitation and propagation fracture toughness, dynamic fracture energy and fracture velocity）, and dy- namic techniques for studying the influences of temperature and pore water.

Analysis of flow response to fluctuation of rotational speed in a radial impeller

By discretizing the convection terms with AUSM+-up scheme in the rotating coordinate system,a finite volume analysis code based on multi-block structured grids was developed independently so as to realize the numerical solving of internal flow fields of turbomachineries.Taking an unshrouded radial impeller with the working fluid of water vapour as the research object,the flow response to the fluctuation of rotational speed was calculated.By comparing the surface pressure profiles and velocity contours calculated by the code and commercial software respectively,the accuracy of flow solver was verified.The analysis of flow response data indicates that,as the working condition shifts closer towards the surge boundary,the response of flow parameters such as mass flow and aerodynamic torque will be more nonsynchronous with the fluctuation of rotational speed,and also the influence of density variation on mass flow variation will be smaller.Moreover,the transient variation region of working condition performance will deviate farther away from the steady performance curve as the working condition approaches the surge boundary.Compared to the working conditions with small mass flows,the distribution characteristics of pressure difference load on the blade surface vary little under large mass flow conditions.The reduction of fluctuation amplitude of rotational speed exerts no influence on abating the hysteresis of flow response.

Degree One Loading by Pressure Variations at the CMB

Hemispherical asymmetry in core dynamics induces degree-1 pressure variations at the core mantle boundary （CMB）, which in turn deforms the overlaying elastic mantle, at the same time keeps center of mass of the whole Earth stationary in space. We develop a systematic procedure to deal with the degree-1 CMB pressure loading. We find by direct calculation a surprisingly negative load Love number h1=-1.425 for vertical displacement. Further analysis indicates that the negative hi corre- sponds to thickening above the positive load that defies intuition that pressure inflation pushes over- laying material up and thins the enveloping shell. We also redefine the pressure load Love numbers in general to enable comparison between the surface mass load and the CMB pressure load for the whole spectrum of harmonic degrees. We find that the gravitational perturbations from the two kinds of loads at degrees n〉l are very similar in amplitude but opposite in sign. In particular, if the CMB pressure variation at degree 2 is at the level of -1 hpa/yr （1 cm water height per year）, it would perturb the variation of Earth＇s oblateness, known as the J2, at the observed level.

Micro-scale behavior of granular materials during cyclic loading

This study presents the micro-scale behavior of granular materials under biaxial cyclic loading for differ- ent confining pressures using the two-dimensional （2D） discrete element method （DEM）. Initially, 8450 ovals were generated in a rectangular frame without any overlap. Four dense samples having confining pressures of 15, 25, 50, and 100 kPa were prepared from the initially generated sparse sample. Numeri- cal simulations were performed under biaxial cyclic loading using these isotropically compressed dense samples. The numerical results depict stress-strain-dilatancy behavior that was similar to that observed in experimental studies. The relationship between the stress ratio and dilatancy rate is almost indepen- dent of confining pressures during loading but significantly dependent on the confining pressures during unloading. The evolution of the coordination number, effective coordination number and slip coordina- tion number depends on both the confining pressures and cyclic loading. The cyclic loading significantly affects the microtopology of the granular assembly. The contact fabric and the fabric-related anisotropy are reported, as well. A strong correlation between the stress ratio and the fabric related to contact normals is observed during cyclic loading, irrespective of confining pressures.

Experimental investigation into transient pressure pulses during pneumatic conveying of fine powders using Shannon entropy

This paper presents the results of an ongoing investigation into transient pressure pulses using Shan- non entropy. Pressure fluctuations （produced by gas-solid two-phase flow during fluidized dense-phase conveying） are recorded by pressure transducers installed at strategic locations along a pipeline. This work validates previous work on identifying the flow mode from pressure signals （Mittal, Mallick, ＆ Wypych, 2014）. Two different powders, namely fly ash （median particle diameter 45 μm, particle den- sity 1950 kg/m3. loosely poured bulk density 950 kg/m3） and cement （median particle diameter 15 p,m, particle density 3060 kg/m3, loosely poured bulk density 1070 kg/m3）, are conveyed through different pipelines （51 mm I.D. × 70 m length and 63 mm I.D. × 24 m length）. The transient nature of pressure fluc- tuations （instead of steady-state behavior） is considered in investigating flow characteristics. Shannon entropy is found to increase along straight pipe sections for both solids and both pipelines. However, Shannon entropy decreases after a bend. A comparison of Shannon entropy among different ranges of superficial air velocity reveals that high Shannon entropy corresponds to very low velocities （i.e. 3-5 m/s） and very high velocities （i.e. 11-14 m/s） while low Shannon entropy corresponds to mid-range velocities （i.e. 6-8 m/s）.

The development and application of an inviscid inverse method

A two-dimensional inviscid inverse method is developed,verified and applied inthis paper.The method solves the Euler equation in absolute reference frame by a cell-centeredfinite volume method,and the hybrid Runge-Kutta method ls used for time integration.Different from the direct method,the inverse method imposes a unique"transpiration"boundary condition on the blade surfaces.The inputs of inverse method are pressure loadingand blade tangential thickness distribution along the blade chord.During the time marchingprocess,the blade shape is periodically updated.When the solution is converged,the bladeshape will be stabled.In the paper,the principle of the inverse method is described in detail.Then the developed inverse method is verified against a consistence test:recover an axialcompressor cascade from a different start.Finally,to demonstrate the powerful capability ofthe method,it is used to redesign the cascade,and final results give an improved aerodynamicperformance.

Mass transport in a thin layer of power-law fluid in an Eulerian coordinate system

The mass transport velocity in a thin layer of muddy fluid is studied theoretically. The mud motion is driven by a periodic pressure load on the free surface, and the mud is described by a power-law model. Based on the key assumptions of the shallowness and the small deformation, a perturbation analysis is conducted up to the second order to find the mean Eulerian velocity in an Eulerian coordinate system. The numerical iteration method is adopted to solve these non-linear equations of the leading order. From the numerical results, both the first-order flow fields and the second-order mass transport velocities are examined. The verifications are made by comparing the numerical results with experimental results in the literature, and a good agreement is confirmed.