Topology optimization for crashworthiness of thin-walled structures under axial crash considering nonlinear plastic buckling and locations of plastic hinges

Although topology optimization is established for linear static problems, more effort is required for solving nonlinear plastic problems. A new topology optimization approach with equivalent static loads (ESLs) is suggested to find the optimum topologies and locations of plastic hinges of thin-walled crash boxes by considering crash-induced deformation, the main crash energy-absorbing mechanism. Together with finite element method crashworthiness analyses, considering all nonlinearities with rate-dependent plasticity, the method was developed using an appropriate time-incremental scheme of ESLs without removing any high values of loads. Analyses show that the crash boxes with optimum topologies have energy-absorbing capabilities equivalent to the original structure. The proposed method is evaluated for two crashes: a crash box at low speed and a double cell subjected to high-speed collision. The results indicate that this method captures nonlinear crushing behaviours and accurate locations of plastic hinges where, if proper reinforcements are made, energy absorption can be enhanced.     

冲击载荷下TC18钛合金力学性能模拟研究

TC18钛合金在退火状态下的屈服强度非常高,进行动态力学性能测试存在一定困难,可优先考虑采用数值模拟方法对其动态力学性能进行研究,而数值模拟结果的精确性主要取决于合金本构方程的精确性。为了减少数值模拟与实验结果之间的差异,在Johnson-Cook模型的基础上做了调整,利用ABAQUS软件的子程序接口自定义了材料本构子程序UMAT,以代码的形式来扩展ABAQUS程序的功能。通过与ABAQUS/Standard软件分析结果以及实验结果的对比,验证了自定义材料本构子程序UMAT计算结果的精确性与可靠性。该研究可为TC18钛合金的结构设计提供参考。     

Impact of industrial virtual power plant on renewable energy integration

An industrial park is one of the typical energy consumption schemes in power systems owing to the heavy industrial loads and their abilities to respond to electricity price changes. Therefore, energy integration in the industrial sector is significant. Accordingly, the concept of industrial virtual power plant (IVPP) has been proposed to deal with such problems. This study demonstrates an IVPP model to manage resources in an eco-industrial park, including energy storage systems, demand response (DR) resources, and distributed energies. In addition, fuzzy theory is used to change the deterministic system constraints to fuzzy parameters, considering the uncertainty of renewable energy, and fuzzy chance constraints are then set based on the credibility theory. By maximizing the daily benefits of the IVPP owners in day-ahead markets, DR and energy storage systems can be scheduled economically. Therefore, the energy between the grid and IVPP can flow in both directions: the surplus renewable electricity of IVPP can be sold in the market; when the electricity generated inside IVPP is not enough for its use, IVPP can also purchase power through the market. Case studies based on three wind-level scenarios demonstrate the efficient synergies between IVPP resources. The validation results indicate that IVPP can optimize the supply and demand resources in industrial parks, thereby decarbonizing the power systems.     

The potential of energy flexibility of space heating and cooling in Portugal

This study assesses the potential of energy flexibility of space heating and cooling for a typical household under different geographical conditions in Portugal. The proposed approach modifies the demand through the optimization of the thermostat settings using a genetic algorithm to reduce either operational costs or interaction with the grid. The results show that the used energy flexibility indicator expresses the available potential and that flexibility depends on several factors, namely: i) thermal inertia of the archetypical household; ii) the time of use electricity tariffs; iii) users’ comfort boundaries; and iv) the geographical location of the houses.     

Ballistic behavior of plain and reinforced concrete slabs under high velocity impact

This work presents a numerical simulation of ballistic penetration and high velocity impact behavior of plain and reinforced concrete slabs. In this paper, we focus on the comparison of the performance of the plain and reinforced concrete slabs of unconfined compressive strength 41 MPa under ballistic impact. The concrete slab has dimensions of 675 mm × 675 mm × 200 mm, and is meshed with 8-node hexahedron solid elements in the impact and outer zones. The ogive-nosed projectile is considered as rigid element that has a mass of 0.386 kg and a length of 152 mm. The applied velocities vary between 540 and 731 m/s. 6 mm of steel reinforcement bars were used in the reinforced concrete slabs. The constitutive material modeling of the concrete and steel reinforcement bars was performed using the Johnson-Holmquist-2 damage and the Johnson-Cook plasticity material models, respectively. The analysis was conducted using the commercial finite element package Abaqus/Explicit. Damage diameters and residual velocities obtained by the numerical model were compared with the experimental results and effect of steel reinforcement and projectile diameter were studies. The validation showed good agreement between the numerical and experimental results. The added steel reinforcements to the concrete samples were found efficient in terms of ballistic resistance comparing to the plain concrete sample.     

Environmental contours for circular‐linear variables based on the direct sampling method

Environmental contours are often used in the design of engineering structures to identify extreme environmental conditions that may give rise to extreme loads and responses. The perhaps most common application of environmental contours is for wave climate variables such as significant wave height and wave period. However, for the design of wind energy installations, the joint distribution of wind speed and wind direction may be equally important. In this case, joint modelling of linear (wind speed) and circular (wind direction) variables are needed, and methods for establishing environmental contours for circular‐linear variables will be required. In this paper, different ways of establishing environmental contours for circular‐linear variables will be presented and applied to a joint distribution model for wind speed and wind direction. In particular, the direct sampling approach to environmental contours will be modified to the case where one of the variables is cyclic. In addition, contours based on exceedance planes in polar coordinates will be established, and circular‐linear contours will also be calculated based on the inverse FORM (I‐FORM) approach.     

Experimentelle und analytische Untersuchung des out‐of‐plane‐Verhaltens von unbewehrten Mauerwerkswänden unter Erdbebenbelastung

Experimental and analytical investigation of the seismic out‐of‐plane behavior of unreinforced masonry walls In addition to the vertical and horizontal load‐bearing in‐plane, masonry must also withstand out‐of‐plane loads that occur in earthquake scenarios. The out‐of‐plane behavior of unreinforced masonry walls depends on a variety of parameters and is very complex due to the strong non‐linearity. Current design methods in German codes and various international codes have not been explicitly developed for out‐of‐plane behavior and contain considerable conservatism. In the present work, shaking‐table experiments with heat‐insulating masonry walls have been conducted to investigate the out‐of‐plane behavior of vertical spanning unreinforced masonry walls. As shown in previous numerical investigations, important parameters are neglected in existing design and analysis models and the out‐of‐plane capacity is underestimated significantly. In the conducted experiments the results of these numerical investigations are verified. Furthermore, the development of an analytical design model to determine the force‐displacement relationship and the out‐of‐plane load‐bearing capacity considering all significant parameters is presented.     

Dynamic response of a beam to the train of moving forces driven by an Erlang renewal process

Dynamic response of a beam to a random train of moving forces moving with the same velocity is considered. Unlike a widely used Poisson process model, a more adequate Erlang renewal process is used as a process driving the train of forces. Normal-mode approach is used to convert the problem into that of a renewal driven train of general pulses. Consequently the modal responses are the filtered renewal processes and are expressed as integrals with respect to the response to a single pulse (passage of a force) and to the increments of the counting renewal process. The expressions for the mean values and cross-correlations of modal responses are obtained as single and double integrals, respectively. The results are presented in terms of the renewal density of the underlying Erlang counting process. Mean value and variance of the mid-span deflection of the beam are determined by numerical evaluation of the pertinent integrals. Numerical analysis is carried out for the Erlang processes with integer parameter $k=2$, $k=3$ and 4 and, for comparison, also for a Poisson process. Different traffic conditions such as the velocity and the mean arrival rate of vehicles are taken into account.     

A new direct version of the Cook-Mayne method for wind pressure probabilities in temperate storms

The original Cook-Mayne (CM) method for obtaining a wind pressure or load of a prescribed probability is reviewed. A new direct calculation method is introduced. This does not require either Monte Carlo simulations or the assumption that extreme pressures conform to the ultimate Fisher-Tippett Type I asymptote. The required value is calculated directly with no intermediate stages. The new method still relies on the original assumptions that both the pressure coefficient and the wind dynamic head have extremes which conform to the ultimate FT1 asymptote.Violation of these two assumptions is studied in some cases where calculations based on these assumptions are compared with exact results. The direct CM method is found to retain accuracy in spite of departures from the assumptions and therefore can provide a robust design tool.