scispace - formally typeset
Search or ask a question

Showing papers in "Acta Mechanica Sinica in 2013"


Journal ArticleDOI
TL;DR: A novel implementation of the asymptotic homogenization method is developed, which has rigorous mathematical foundation of the AH method, and also simplicity as the RVE method.
Abstract: Representative volume element (RVE) method and asymptotic homogenization (AH) method are two widely used methods in predicting effective properties of periodic materials. This paper develops a novel implementation of the AH method, which has rigorous mathematical foundation of the AH method, and also simplicity as the RVE method. This implementation can be easily realized using commercial software as a black box, and can use all kinds of elements available in commercial software to model unit cells with rather complicated microstructures, so the model may remain a fairly small scale. Several examples were carried out to demonstrate the simplicity and effectiveness of the new implementation.

134 citations


Journal ArticleDOI
TL;DR: In this article, the authors studied the transition of a flared cone with zero angle of attack in a Mach 6 quiet wind tunnel (M6QT) via wall pressure measurement and flow visualization.
Abstract: Transition on a flared cone with zero angle of attack was studied in our newly established Mach 6 quiet wind tunnel (M6QT) via wall pressure measurement and flow visualization. High-frequency pressure transducers were used to measure the second-mode waves’ amplitudes and frequencies. Using pulsed schlieren diagnostic and Rayleigh scattering technique, we got a clear evolution of the second-mode disturbances. The second-mode waves exist for a long distance, which means that the second-mode waves grow linearly in a large region. Strong Mach waves are radiated from the edge of the boundary layer. With further development, the second-mode waves reach their maximum magnitude and harmonics of the second-mode instability appear. Then the disturbances grow nonlinearly. The second modes become weak and merge with each other. Finally, the nonlinear interaction of disturbance leads to a relatively quiet zone, which further breaks down, resulting in the transition of the boundary layer. Our results show that transition is determined by the second mode. The quiet zone before the final breakdown is observed in flow visualization for the first time. Eventual transition requires the presence of a quiet zone generated by nonlinear interactions.

79 citations


Journal ArticleDOI
TL;DR: In this article, a uniform magnetic field is applied vertically to the flow direction and the governing equations are reduced to non-linear coupled partial differential equations and solved by means of homotopy analysis method.
Abstract: In this paper, we study the unsteady coupled heat and mass transfer of two-dimensional MHD fluid over a moving oscillatory stretching surface with Soret and Dufour effects. Viscous dissipation effects are adopted in the energy equation. A uniform magnetic field is applied vertically to the flow direction. The governing equations are reduced to non-linear coupled partial differential equations and solved by means of homotopy analysis method (HAM). The effects of some physical parameters such as magnetic parameter, Dufour number, Soret number, the Prandtl number and the ratio of the oscillation frequency of the sheet to its stretching rate on the flow and heat transfer characteristics are illustrated and analyzed.

64 citations


Journal ArticleDOI
TL;DR: In this article, the dislocation density tensor computed as the curl of plastic distortion is regarded as a new constitutive variable in crystal plasticity, and the dependence of the free energy function on the dislocalization tensor is explored starting from a quadratic ansatz.
Abstract: The dislocation density tensor computed as the curl of plastic distortion is regarded as a new constitutive variable in crystal plasticity. The dependence of the free energy function on the dislocation density tensor is explored starting from a quadratic ansatz. Rank one and logarithmic dependencies are then envisaged based on considerations from the statistical theory of dislocations. The relevance of the presented free energy potentials is evaluated from the corresponding analytical solutions of the periodic two-phase laminate problem under shear where one layer is a single crystal material undergoing single slip and the second one remains purely elastic.

57 citations


Journal ArticleDOI
TL;DR: In this paper, a comprehensive review on the research and development in multi-scale numerical modeling and simulation of PEM fuel cells is provided, including fluid flow and species transport, electron and proton transport, heat transfer and thermal management, liquid water transport and water management, transient response behaviors, and cold start processes.
Abstract: This paper provides a comprehensive review on the research and development in multi-scale numerical modeling and simulation of PEM fuel cells. An overview of recent progress in PEM fuel cell modeling has been provided. Fundamental transport phenomena in PEM fuel cells and the corresponding mathematical formulation of macroscale models are analyzed. Various important issues in PEM fuel cell modeling and simulation are examined in detail, including fluid flow and species transport, electron and proton transport, heat transfer and thermal management, liquid water transport and water management, transient response behaviors, and cold-start processes. Key areas for further improvements have also been discussed.

46 citations


Journal ArticleDOI
TL;DR: In this paper, a real CRH2 model with a leading car, a middle car and a trailing car included was studied, and the approaches of unsteady Reynold-averaged Navier-Stokes (URANS) and detached eddy simulation (DES) were utilized, respectively.
Abstract: Intensive turbulence exists in the wakes of high speed trains, and the aerodynamic performance of the trailing car could deteriorate rapidly due to complicated features of the vortices in the wake zone. As a result, the safety and amenity of high speed trains would face a great challenge. This paper considers mainly the mechanism of vortex formation and evolution in the train flow field. A real CRH2 model is studied, with a leading car, a middle car and a trailing car included. Different running speeds and cross wind conditions are considered, and the approaches of unsteady Reynold-averaged Navier-Stokes (URANS) and detached eddy simulation (DES) are utilized, respectively. Results reveal that DES has better capability of capturing small eddies compared to URANS. However, for large eddies, the effects of two approaches are almost the same. In conditions without cross winds, two large vortex streets stretch from the train nose and interact strongly with each other in the wake zone. With the reinforcement of the ground, a complicated wake vortex system generates and becomes strengthened as the running speed increases. However, the locations of flow separations on the train surface and the separation mechanism keep unchanged. In conditions with cross winds, three large vortices develop along the leeward side of the train, among which the weakest one has no obvious influence on the wake flow while the other two stretch to the tail of the train and combine with the helical vortices in the train wake. Thus, optimization of the aerodynamic performance of the trailing car should be aiming at reducing the intensity of the wake vortex system.

45 citations


Journal ArticleDOI
TL;DR: In this paper, a combination of a tracking algorithm and a resistor network approach is developed to predict the connectivity and effective conductivity for the various densified structures of solid oxide fuel cells (SOFC) and lithium-ion batteries (LIB).
Abstract: Optimization of composition and microstructure is important to enhance performance of solid oxide fuel cells (SOFC) and lithium-ion batteries (LIB). For this, the porous electrode structures of both SOFC and LIB are modeled as a binary mixture of electronic and ionic conducting particles to estimate effective transport properties. Particle packings of 10 000 spherical, binary sized and randomly positioned particles are created numerically and densified considering the different manufacturing processes in SOFC and LIB: the sintering of SOFC electrodes is approximated geometrically, whereas the calendering process and volume change due to intercalation in LIB are modeled physically by a discrete element approach. A combination of a tracking algorithm and a resistor network approach is developed to predict the connectivity and effective conductivity for the various densified structures. For SOFC, a systematic study of the influence of morphology on connectivity and conductivity is performed on a large number of assemblies with different compositions and particle size ratios between 1 and 10. In comparison to percolation theory, an enlarged percolation area is found, especially for large size ratios. It is shown that in contrast to former studies the percolation threshold correlates to varying coordination numbers. The effective conductivity shows not only an increase with volume fraction as expected but also with size ratio. For LIB, a general increase of conductivity during the intercalation process was observed in correlation with increasing contact forces. The positive influence of calendering on the percolation threshold and the effective conductivity of carbon black is shown. The anisotropy caused by the calendering process does not influence the carbon black phase.

44 citations


Journal ArticleDOI
TL;DR: In this paper, the Galerkin method and its convergence for the response of a Timoshenko beam supported by a nonlinear foundation is studied for the first time, where the effects of shear deformable beams and the shear deformation of foundations are considered at the same time.
Abstract: The present paper investigates the dynamic response of finite Timoshenko beams resting on a sixparameter foundation subjected to a moving load. It is for the first time that the Galerkin method and its convergence are studied for the response of a Timoshenko beam supported by a nonlinear foundation. The nonlinear Pasternak foundation is assumed to be cubic. Therefore, the effects of the shear deformable beams and the shear deformation of foundations are considered at the same time. The Galerkin method is utilized for discretizing the nonlinear partial differential governing equations of the forced vibration. The dynamic responses of Timoshenko beams are determined via the fourth-order Runge-Kutta method. Moreover, the effects of different truncation terms on the dynamic responses of a Timoshenko beam resting on a complex foundation are discussed. The numerical investigations shows that the dynamic response of Timoshenko beams supported by elastic foundations needs super high-order modes. Furthermore, the system parameters are compared to determine the dependence of the convergences of the Galerkin method.

42 citations


Journal ArticleDOI
TL;DR: In this article, a simplified interplanetary kite-craft accelerated by radiation of the Sun (IKAROS) model is established by using the absolute-coordinate-based (ACB) method that combines the natural coordinate formulation (NCF) describing the central rigid hub and the absolute nodal coordinate formulation describing flexible parts.
Abstract: The spinning solar sail of large scale has been well developed in recent years. Such a solar sail can be considered as a rigid-flexible multibody system mainly composed of a spinning central rigid hub, a number of flexible thin tethers, sail membranes, and tip masses. A simplified interplanetary kite-craft accelerated by radiation of the Sun (IKAROS) model is established in this study by using the absolute-coordinate-based (ACB) method that combines the natural coordinate formulation (NCF) describing the central rigid hub and the absolute nodal coordinate formulation (ANCF) describing flexible parts. The initial configuration of the system in the second-stage deployment is determined through both dynamic and static analyses. The huge set of stiff equations of system dynamics is solved by using the generalized-alpha method, and thus the deployment dynamics of the system can be well understood.

41 citations


Journal ArticleDOI
TL;DR: In this article, a hybrid method of nonlinear acoustic solver (NLAS) and Fowcs Williams-Hawkings (FW-H) acoustic analogy is used to predict the aerodynamic noise of pantograph system in this speed range.
Abstract: Pantograph system of high-speed trains become significant source of aerodynamic noise when travelling speed exceeds 300 km/h. In this paper, a hybrid method of non-linear acoustic solver (NLAS) and Ffowcs Williams-Hawkings (FW-H) acoustic analogy is used to predict the aerodynamic noise of pantograph system in this speed range. When the simulation method is validated by a benchmark problem of flows around a cylinder of finite span, we calculate the near flow field and far acoustic field surrounding the pantograph system. And then, the frequency spectra and acoustic attenuation with distance are analyzed, showing that the pantograph system noise is a typical broadband one with most acoustic power restricted in the medium-high frequency range from 200 Hz to 5 kHz. The aerodynamic noise of pantograph systems radiates outwards in the form of spherical waves in the far field. Analysis of the overall sound pressure level (OASPL) at different speeds exhibits that the acoustic power grows approximately as the 4th power of train speed. The comparison of noise reduction effects for four types of pantograph covers demonstrates that only case 1 can lessen the total noise by about 3 dB as baffles on both sides can shield sound wave in the spanwise direction. The covers produce additional aerodynamic noise themselves in the other three cases and lead to the rise of OASPLs.

40 citations


Journal ArticleDOI
TL;DR: In this article, the authors developed a general method to analytically tackle a kind of movable boundary problem from the viewpoint of energy variation, and grouped the adhesion of a micro-beam, droplet and carbon nanotube (CNT) ring on a substrate into one framework.
Abstract: In this study, we developed a general method to analytically tackle a kind of movable boundary problem from the viewpoint of energy variation. Having grouped the adhesion of a micro-beam, droplet and carbon nanotube (CNT) ring on a substrate into one framework, we used the developed line of reasoning to investigate the adhesion behaviors of these systems. Based upon the derived governing equations and transversality conditions, explicit solutions involving the critical parameters and morphologies for the three systems are successfully obtained, and then the parameter analogies and common characteristics of them are thoroughly investigated. The presented method has been verified via the concept of energy release rate in fracture mechanics. Our analyses provide a new approach for exploring the mechanism of different systems with similarities as well as for understanding the unity of nature. The analysis results may be beneficial for the design of nano-structured materials, and hold potential for enhancing their mechanical, chemical, optical and electronic properties.

Journal ArticleDOI
Sean R. Bishop1
TL;DR: The importance of oxygen non-stoichiometry induced expansion, known as chemical expansion, for the mechanical properties of solid oxide fuel cells (SOFCs) is discussed in this paper.
Abstract: The importance of oxygen non-stoichiometry induced expansion, known as chemical expansion, for the mechanical properties of solid oxide fuel cells (SOFCs) is discussed. The methods used to measure chemical expansion and the defects responsible for its existence are introduced. Recent work demonstrating the origin of chemical expansion in fluorite structured oxides for SOFCs is presented. Models used to predict stress induced by chemical expansion in SOFCs, highlighting the necessity of considering electro-chemo-mechanical coupling relationships, are discussed.

Journal ArticleDOI
TL;DR: In this article, a sessile water droplet evaporation on polydimethylsiloxane (PDMS) surfaces with three different curing ratios (5: 1, 10: 1 and 20: 1) was experimentally investigated.
Abstract: Evaporation of sessile water droplet on polydimethylsiloxane (PDMS) surfaces with three different curing ratios (5: 1, 10: 1, and 20: 1) was experimentally investigated in this paper. We show that the constant contact radius (CCR) evaporation on surface with high curing ratio lasts longer than that with low curing ratio. We also measured Young’s moduli of PDMS films by using atomic force microscopy (AFM) and simulated surface deformation of PDMS films induced by sessile water droplet. With increasing curing ratio of PDMS film, Young’s modulus of PDMS film is getting lower, and then there will be larger surface deformation and more elastic stored energy. Since such energy acts as a barrier to keep the three-phase contact line pinned, thus it will result in longer CCR evaporation on PDMS surface with higher curing ratio.

Journal ArticleDOI
TL;DR: A framework dedicated to the stochastic modeling of random properties is first introduced and a probabilistic model for matrix-valued second-order random fields with symmetry propertries, recently proposed in the literature, is further reviewed.
Abstract: This paper is concerned with the modeling of randomness in multiscale analysis of heterogeneous materials. More specifically, a framework dedicated to the stochastic modeling of random properties is first introduced. A probabilistic model for matrix-valued second-order random fields with symmetry propertries, recently proposed in the literature, is further reviewed. Algorithms adapted to the Monte Carlo simulation of the proposed representation are also provided. The derivations and calibration procedure are finally exemplified through the modeling of the apparent properties associated with an elastic porous microstructure containing stochastic interphases.

Journal ArticleDOI
TL;DR: In this article, the linearized Poisson-Boltzmann equation and the Cauchy momentum equation have been solved in semi-analytical solutions for transient electroosmotic and pressure-driven flows (EOF/PDF) between microparallel plates.
Abstract: By method of the Laplace transform, this article presents semi-analytical solutions for transient electroosmotic and pressure-driven flows (EOF/PDF) of two-layer fluids between microparallel plates. The linearized Poisson-Boltzmann equation and the Cauchy momentum equation have been solved in this article. At the interface, the Maxwell stress is included as the boundary condition. By numerical computations of the inverse Laplace transform, the effects of dielectric constant ratio ɛ, density ratio ρ, pressure ratio p, viscosity ratio µ of layer II to layer I, interface zeta potential difference \(\Delta \bar \psi\), interface charge density jump Q, the ratios of maximum electro-osmotic velocity to pressure velocity α, and the normalized pressure gradient B on transient velocity amplitude are presented. We find the velocity amplitude becomes large with the interface zeta potential difference and becomes small with the increase of the viscosity. The velocity will be large with the increases of dielectric constant ratio; the density ratio almost does not influence the EOF velocity. Larger interface charge density jump leads to a strong jump of velocity at the interface. Additionally, the effects of the thickness of fluid layers (h1 and h2) and pressure gradient on the velocity are also investigated.

Journal ArticleDOI
TL;DR: In this article, the elastic behavior of saturated porous materials under undrained freezing is investigated by using a poromechanical approach, where thermodynamic equilibria are used to describe the crystallization process of the partially frozen solution in bulk state and confined state in pores.
Abstract: The elastic behavior of saturated porous materials under undrained freezing is investigated by using a poromechanical approach. Thermodynamic equilibria are used to describe the crystallization process of the partially frozen solution in bulk state and confined state in pores. By phase transition at freezing, fusion energy, thermal contraction of solid, solution and ice crystals, volume changes of crystallization build up remarkable pore pressure that induces expansion or shrinkage of solid matrix. Owing to the lower chemical potential when pore water mixes with salts, fewer ice forms in pores. Penetration of ice into the porous materials increases the capillary pressure, but limits effect on the pore liquid pressure and the strain of solid matrix. On the contrary, the pore pressure induced by solution density rises as salt concentration increases and causes significant shrinkage of solid matrix.

Journal ArticleDOI
TL;DR: In this paper, a non-probabilistic robust reliability method for LQR-based static output feedback robust control of uncertain structures is presented by treating bounded uncertain parameters as interval variables.
Abstract: Uncertainty is inherent and unavoidable in almost all engineering systems. It is of essential significance to deal with uncertainties by means of reliability approach and to achieve a reasonable balance between reliability against uncertainties and system performance in the control design of uncertain systems. Nevertheless, reliability methods which can be used directly for analysis and synthesis of active control of structures in the presence of uncertainties remain to be developed, especially in non-probabilistic uncertainty situations. In the present paper, the issue of vibration control of uncertain structures using linear quadratic regulator (LQR) approach is studied from the viewpoint of reliability. An efficient non-probabilistic robust reliability method for LQR-based static output feedback robust control of uncertain structures is presented by treating bounded uncertain parameters as interval variables. The optimal vibration controller design for uncertain structures is carried out by solving a robust reliability-based optimization problem with the objective to minimize the quadratic performance index. The controller obtained may possess optimum performance under the condition that the controlled structure is robustly reliable with respect to admissible uncertainties. The proposed method provides an essential basis for achieving a balance between robustness and performance in controller design of uncertain structures. The presented formulations are in the framework of linear matrix inequality and can be carried out conveniently. Two numerical examples are provided to illustrate the effectiveness and feasibility of the present method.

Journal ArticleDOI
TL;DR: In this article, the behavior of a bubble near a rigid cylinder is studied experimentally as the positions of bubble induction change, and several cylinders with different diameters are used in the experiment.
Abstract: In this paper, the behavior of a bubble near a rigid cylinder is studied experimentally as the positions of bubble induction change, and several cylinders with different diameters are used in the experiment. The main results are as follows. The behavior of a bubble near a rigid cylinder is distinct from that near a rigid plate. When the cylinders are laid in deep water, there will occur three kinds of typical bubble shapes as the distance between bubble and cylinder increases. And the bubble shapes are different as the diameter of cylinder varies. When the cylinders are laid near a free surface, the behaviors of bubble near cylinders with different diameters are similar. For a certain distance between bubble and free surface, as the distance between bubble and cylinder increases, “double jet”, “inclined jet” and “downward jet” will take place successively.

Journal ArticleDOI
TL;DR: In this paper, the penetration of a circular indenter with a flat base into a soft functionally graded elastic layer is considered, where the elastic properties of a functionally graded layer arbitrarily vary with depth and the foundation is assumed to be elastic, yet much harder than a layer.
Abstract: The paper addresses a contact problem of the theory of elasticity, i.e., the penetration of a circular indenter with a flat base into a soft functionally graded elastic layer. The elastic properties of a functionally graded layer arbitrarily vary with depth, and the foundation is assumed to be elastic, yet much harder than a layer. Approximated analytical solution is constructed, and it is shown that the solutions are asymptotically exact both for large and small values of characteristic dimensionless geometrical parameter of the problem. Numerical examples are analyzed for the cases of monotonic and nonmonotonic variations of elastic properties. Numerical results for the case of homogeneous layer are compared with the results for nondeformable foundation.

Journal ArticleDOI
TL;DR: In this paper, an exact solution of the forced vibration of circular arches subjected to subsurface denotation forces is obtained, with the introduction of an interfacial damping between the structure element and the surrounding soil into the equation of motion.
Abstract: Due to the wide applications of arches in underground protective structures, dynamic analysis of circular arches including soil-structure interactions is important. In this paper, an exact solution of the forced vibration of circular arches subjected to subsurface denotation forces is obtained. The dynamic soil-structure interaction is considered with the introduction of an interfacial damping between the structure element and the surrounding soil into the equation of motion. By neglecting the influences of shear, rotary inertia and tangential forces and assuming the arch incompressible, the equations of motion of the buried arches were set up. Analytical solutions of the dynamic responses of the protective arches were deduced by means of modal superposition. Arches with different opening angles, acoustic impedances and rise-span ratios were analyzed to discuss their influences on an arch. The theoretical analysis suggests blast loads for elastic designs and predicts the potential failure modes for buried protective arches.

Journal ArticleDOI
TL;DR: In this article, the authors demonstrate that strain engineering can serve as an effective way to control the surface wettability by using a graphene sheet as a sample material and molecular dynamic simulations.
Abstract: Reversible control of surface wettability has wide applications in lab-on-chip systems, tunable optical lenses, and microfluidic tools. Using a graphene sheet as a sample material and molecular dynamic simulations, we demonstrate that strain engineering can serve as an effective way to control the surface wettability. The contact angles θ of water droplets on a graphene vary from 72.5° to 106° under biaxial strains ranging from −10% to 10% that are applied on the graphene layer. For an intrinsic hydrophilic surface (at zero strain), the variation of θ upon the applied strains is more sensitive, i.e., from 0° to 74.8°. Overall the cosines of the contact angles exhibit a linear relation with respect to the strains. In light of the inherent dependence of the contact angle on liquid-solid interfacial energy, we develop an analytic model to show the cos θ as a linear function of the adsorption energy E ads of a single water molecule over the substrate surface. This model agrees with our molecular dynamic results very well. Together with the linear dependence of E ads on biaxial strains, we can thus understand the effect of strains on the surface wettability. Thanks to the ease of reversibly applying mechanical strains in micro/nano-electromechanical systems, we believe that strain engineering can be a promising means to achieve the reversibly control of surface wettability.

Journal ArticleDOI
TL;DR: Wang et al. as discussed by the authors used the total Lagrangian smoothed particle hydrodynamics (SPH) method equation and moving least square (MLS) function to improve the SPH method regarding the stability and consistency.
Abstract: With the incorporation of total Lagrangian smoothed particle hydrodynamics (SPH) method equation and moving least square (MLS) function, the traditional SPH method is improved regarding the stability and consistency. Based on Mindlin-Ressiner plate theory, the SPH method simulating dynamic behavior via one layer of particles is applied to plate’s mid-plane, i.e., a SPH shell model is constructed. Finally, through comparative analyses on the dynamic response of square, stiffened shells and cylindrical shells under various strong impact loads with common finite element software, the feasibility, validity and numerical accuracy of the SPH shell method are verified. Consequently, further researches on SPH shell may well pave the way towards solving problems involving dynamic plastic damage, tearing or even crushing.

Journal ArticleDOI
TL;DR: The numerical optimization framework developed herein can be applied directly to cell and pack design and enable the quantification of the tradeoffs between energy and power density, and the effect of optimizing the electrode thickness and porosity.
Abstract: In this study, the effects of discharge rate and LiMn2O4 cathode properties (thickness, porosity, particle size, and solid-state diffusivity and conductivity) on the gravimetric energy and power density of a lithium-ion battery cell are analyzed simultaneously using a cell-level model. Surrogate-based analysis tools are applied to simulation data to construct educed-order models, which are in turn used to perform global sensitivity analysis to compare the relative importance of cathode properties. Based on these results, the cell is then optimized for several distinct physical scenarios using gradient-based methods. The complementary nature of the gradient- and surrogate-based tools is demonstrated by establishing proper bounds and constraints with the surrogate model, and then obtaining accurate optimized solutions with the gradient-based optimizer. These optimal solutions enable the quantification of the tradeoffs between energy and power density, and the effect of optimizing the electrode thickness and porosity. In conjunction with known guidelines, the numerical optimization framework developed herein can be applied directly to cell and pack design.

Journal ArticleDOI
TL;DR: In this article, the influence of the acoustic impedance and shock strength on the jet formation in shock-heavy gas bubble interaction is numerically studied. And the authors conclude that the pressure perturbation and baroclinic vorticity deposition are the two dominant factors for jet formation.
Abstract: The influences of the acoustic impedance and shock strength on the jet formation in shock-heavy gas bubble interaction are numerically studied in this work. The process of a shock interacting with a krypton or a SF6 bubble is studied by the numerical method VAS2D. As a validation, the experiments of a SF6 bubble accelerated by a planar shock were performed. The results indicate that, due to the mismatch of acoustic impedance, the way of jet formation in heavy gas bubble with different species is diversified under the same initial condition. With respect to the same bubble, the manner of jet formation is also distinctly different under different shock strengths. The disparities of the acoustic impedance result in different effects of shock focusing in the bubble, and different behaviors of shock wave inside and outside the bubble. The analyses of the wave pattern and the pressure variation indicate that the jet formation is closely associated with the pressure perturbation. Moreover, the analysis of the vorticity deposition, and comparisons of circulation and baroclinic torque show that the baroclinic vorticity also contributes to the jet formation. It is concluded that the pressure perturbation and baroclinic vorticity deposition are the two dominant factors for the jet formation in shock-heavy gas bubble interaction.

Journal ArticleDOI
TL;DR: In this article, the authors derived the nonlinear governing equations of motion for the system with the effect of electrostatic force, intermolecular tractions and base rotation using extended Hamilton principle.
Abstract: The nonlinear dynamic and static deflection of a micro/nano gyroscope under DC voltages and base rotation are investigated. The gyroscope undertakes two coupled bending motions along the drive and sense directions and subjected to electrostatic actuations and intermolecular forces. The nonlinear governing equations of motion for the system with the effect of electrostatic force, intermolecular tractions and base rotation are derived using extended Hamilton principle. Under constant voltage, the gyroscope finds the preformed shape. First, the deflection of the micro/nano gyroscope under electrostatic forces is obtained by static and dynamic analyses. Furthermore, the static and dynamic instability of the system are investigated. Afterward the oscillatory behavior of the pre-deformed micro/nano gyroscope around equilibrium is studied. The effects of intermolecular and nonlinear parameters on the static and dynamic deflection, natural frequencies and instability of the micro/nano gyroscope are studied. The presented model can be used to exactly determine static and the dynamic behavior of vibratory micro/nano gyroscopes.

Journal ArticleDOI
TL;DR: In this paper, Tan et al. studied the propagation of a cellular bar in the framework of continuum solids by adopting two idealized material models, viz. the dynamic rigid, perfectly plastic, locking (D-R-PP-L) model and the dynamic linear hardening plastic locking model, both considering the effects of strain-rate on the material properties.
Abstract: The propagation of shock waves in a cellular bar is systematically studied in the framework of continuum solids by adopting two idealized material models, viz. the dynamic rigid, perfectly plastic, locking (D-R-PP-L) model and the dynamic rigid, linear hardening plastic, locking (D-R-LHP-L) model, both considering the effects of strain-rate on the material properties. The shock wave speed relevant to these two models is derived. Consider the case of a bar made of one of such material with initial length L0 and initial velocity vi impinging onto a rigid target. The variations of the stress, strain, particle velocity, specific internal energy across the shock wave and the cease distance of shock wave are all determined analytically. In particular the “energy conservation condition” and the “kinematic existence condition” as proposed by Tan et al. (2005) is re-examined, showing that the “energy conservation condition” and the consequent “critical velocity”, i.e. the shock can only be generated and sustained in R-PP-L bars when the impact velocity is above this critical velocity, is incorrect. Instead, with elastic deformation, strain-hardening and strain-rate sensitivity of the cellular materials being considered, it is appropriate to redefine a first and a second critical impact velocity for the existence and propagation of shock waves in cellular solids. Starting from the basic relations for shock wave propagating in D-R-LHP-L cellular materials, a new method for inversely determining the dynamic stress-strain curve for cellular materials is proposed. By using e.g. a combination of Taylor bar and Hopkinson pressure bar impact experimental technique, the dynamic stress-strain curve of aluminum foam could be determined. Finally, it is demonstrated that this new formulation of shock theory in this one-dimensional stress state can be generalized to shocks in a one-dimensional strain state, i.e. for the case of plate impact on cellular materials, by simply making proper replacements of the elastic and plastic constants.

Journal ArticleDOI
TL;DR: In this paper, the authors proposed the concept of the "heliocentric oscillating-Kepler orbit" for near-Earth asteroid flyby trajectories, where the classical orbital elements of a flyby trajectory oscillate significantly with respect to time.
Abstract: Driven by curiosity about possible flight options for the Chang’e-2 spacecraft after it remains at the Sun-Earth L2 point, effective approaches were developed for designing preliminary fuel-optimal near-Earth asteroid flyby trajectories. The approaches include the use of modified unstable manifolds, grid search of the manifolds’ parameters, and a two-impulse maneuver for orbital phase matching and z-axis bias change, and are demonstrated to be effective in asteroid target screening and trajectory optimization. Asteroid flybys are expected to be within a distance of 2 × 107 km from the Earth owing to the constrained Earth-spacecraft communication range. In this case, the spacecraft’s orbital motion is significantly affected by the gravities of both the Sun and the Earth, and therefore, the concept of the “heliocentric oscillating-Kepler orbit” is proposed, because the classical orbital elements of the flyby trajectories referenced in the heliocentric inertial frame oscillate significantly with respect to time. The analysis and results presented in this study show that, among the asteroids whose orbits are the most accurately predicted, “Toutatis”, “2005 NZ6”, or “2010 CL19” might be encountered by Chang’e-2 in late 2012 or 2013 with total impulses less than 100m/s.

Journal ArticleDOI
TL;DR: In this article, the fracture behavior of a thermoelastic cylinder subjected to a sudden temperature change on its outer surface within the framework of non-classical heat conduction was studied.
Abstract: This paper studies the fracture behavior of a thermoelastic cylinder subjected to a sudden temperature change on its outer surface within the framework of non-classical heat conduction. The heat conduction equation is solved by separation of variable technique. Closed form solution for the temperature field and the associated thermal stress are established. The critical parameter governing the level of the transient thermal stress is identified. Exact expression for the transient stress intensity factor is obtained for a crack in the cylinder. The difference between the non-classical solutions and the classical solution are discussed. It is found that the traditional classical heat conduction considerably underestimates the transient thermal stress and thermal stress intensity factor.

Journal ArticleDOI
TL;DR: In this article, the effect of edge diffusion on the delamination process is evaluated and the edge effect is significant for active disks with a small aspect ratio, but negligible for the case of large aspect ratio.
Abstract: Progressive delamination driven by Li-ion diffusion in elastic disk-like thin film electrodes of Li-ion batteries is modeled based on the cohesive model. Axisymmetric diffusion model is considered under both galvanostatic and potentiostatic operations. The effect of edge diffusion on the delamination process is evaluated. It is found that the diffusion from edge leads to an earlier delamination initiation. The edge effect is significant for active disks with a small aspect ratio, but negligible for the case of large aspect ratio. The edge diffusion is weaker in the potentiostatic operation than in the galvanostatic operation.

Journal ArticleDOI
TL;DR: A nodal design variable-based adaptive method for topology optimization of continuum structures that allows for sufficient further topological evolution of the structural layout in higher adaptivity levels and thus essentially enables high-resolution solutions.
Abstract: For the purpose of achieving high-resolution optimal solutions this paper proposes a nodal design variable-based adaptive method for topology optimization of continuum structures. The analysis mesh-independent density field, interpolated by the nodal design variables at a given set of density points, is adaptively refined/coarsened according to a criterion regarding the gray-scale measure of local regions. New density points are added into the gray regions and redundant ones are removed from the regions occupied by purely solid/void phases for decreasing the number of design variables. A penalization factor adaptivity technique is employed to prevent premature convergence of the optimization iterations. Such an adaptive scheme not only improves the structural boundary description quality, but also allows for sufficient further topological evolution of the structural layout in higher adaptivity levels and thus essentially enables high-resolution solutions. Moreover, compared with the case with uniformly and finely distributed density points, the proposed adaptive method can achieve a higher numerical efficiency of the optimization process.