Efficient random variable generation: ratio of uniforms and polar rejection sampling

Monte Carlo techniques, which require the generation of samples from some target density, are often the only alternative for performing Bayesian inference. Two classic sampling techniques to draw independent samples are the ratio of uniforms (RoU) and rejection sampling (RS). An efficient sampling algorithm is proposed combining the RoU and polar RS (i.e. RS inside a sector of a circle using polar coordinates). Its efficiency is shown in drawing samples from truncated Cauchy and Gaussian random variables, which have many important applications in signal processing and communications. RESUMEN. Método eficiente para generar algunas variables aleatorias de uso común en procesado de señal y comunicaciones (por ejemplo, Gaussianas o Cauchy truncadas) mediante la combinación de dos técnicas: "ratio of uniforms" y "rejection sampling".

Moment transfer and influence of transverse beams in interior waffle flat plate-column connections under lateral loading

This paper addresses two aspects of the behavior of interior reinforced concrete waffle flat plate?column connections under lateral loads: the share of the unbalanced moment between flexure and excentric shear, and the effect of the transverse beams. A non-linear finite element model (benchmark model) was developed and calibrated with the results of quasi-static cyclic tests conducted on a 3/5 scale specimen. First, from this numerical model, the portion cv of the unbalanced moment transferred by the excentricity of shear about the centroid of the critical sections defined by Eurocode 2 (EC-2) and by ACI 318-11 was calculated and compared with the share-out prescribed by these codes. It is found that while the critical section of EC-2 is consistent with the cv provided by this code, in the case of ACI 318-11, the value assigned to cv is far below (about 50% smaller) the actual one obtained with the numerical simulations. Second, from the benchmark model, seven additional models were developed by varying the depth D of the transverse beam over the thickness h of the plate. It was found that the ductility of the connection and the effective width of the plate can respectively be increased up to 50% and 10% by raising D/h to 2 and 1.5.

Methodology for assessing regional economic impacts of charging HGVs in Spain: an integrated approach through a random utility based multiregional input-output and a road transportation network model

An integrated approach composed of a random utility-based multiregional input-output model and a road transport network model was developed for evaluating the application of a fee to heavy-goods vehicles (HGVs) in Spain. For this purpose, a distance-based charge scenario (in euros per vehicle kilometer) for HGVs was evaluated for a selected motorway network in Spain. Although the aim of this charging policy was to increase the efficiency of transport, the approach strongly identified direct and indirect impacts on the regional economy. Estimates of the magnitude and extent of indirect effects on aggregated macroeconomic indicators (employment and gross domestic product) are provided. The macroeconomic effects of the charging policy were found to be positive for some regions and negative for other regions.

A random utility based multiregional input-output approach to estimate the impact of a new highway on the regional distribution of jobs

This paper addresses the economic impact assessment of the construction of a new road on the regional distribution of jobs. The paper summarizes different existing model approaches considered to assess economic impacts through a literature review. Afterwards, we present the development of a comprehensive approach for analyzing the interaction of new transport infrastructure and the economic impact through an integrated model. This model has been applied to the construction of the motorway A-40 in Spain (497 Km.) which runs across three regions without passing though Madrid City. This may in turn lead to the relocation of labor and capital due to the improvement of accessibility of markets or inputs. The result suggests the existence of direct and indirect effects in other regions derived from the improvement of the transportation infrastructure, and confirms the relevance of road freight transport in some regions. We found that the changes in regional employment are substantial for some regions (increasing or decreasing jobs), but a t the same time negligible in other regions. As a result,the approach provides broad guidance to national governments and other transport-related parties about the impacts of this transport policy.

Mechanical and thermal behaviour of isotactic polypropylene reinforced with inorganic fullerene-like WS2 nanoparticles: Effect of filler loading and temperature

The thermal and mechanical behaviour of isotactic polypropylene (iPP) nanocomposites reinforced with different loadings of inorganic fullerene-like tungsten disulfide (IF-WS2) nanoparticles was investigated. The IF-WS2 noticeably enhanced the polymer stiffness and strength, ascribed to their uniform dispersion, the formation of a large nanoparticle?matrix interface combined with a nucleating effect on iPP crystallization. Their reinforcement effect was more pronounced at high temperatures. However, a drop in ductility and toughness was found at higher IF-WS2 concentrations. The tensile behaviour of the nanocomposites was extremely sensitive to the strain rate and temperature, and their yield strength was properly described by the Eyring s equation. The activation energy increased while the activation volume decreased with increasing nanoparticle loading, indicating a reduction in polymer chain motion. The nanoparticles improved the thermomechanical properties of iPP: raised the glass transition and heat deflection temperatures while decreased the coefficient of thermal expansion. The nanocomposites also displayed superior flame retardancy with longer ignition time and reduced peak heat release rate. Further, a gradual rise in thermal conductivity was found with increasing IF-WS2 loading both in the glassy and rubbery states. The results presented herein highlight the benefits and high potential of using IF-nanoparticles for enhancing the thermomechanical properties of thermoplastic polymers compared to other nanoscale fillers.

Behaviour of Concrete Structural Members Subjected to Air Blast Loading

Numerical analysis is a suitable tool in the design of complex reinforced concrete structures under extreme impulsive loadings such as impacts or explosions at close range. Such events may be the result of terrorist attacks. Reinforced concrete is commonly used for buildings and infrastructures. For this reason, the ability to accurately run numerical simulations of concrete elements subjected to blast loading is needed. In this context, reliable constitutive models for concrete are of capital importance. In this research numerical simulations using two different constitutive models for concrete (Continuous Surface Cap Model and Brittle Damage Model) have been carried out using LS-DYNA. Two experimental benchmark tests have been taken as reference. The results of the numerical simulations with the aforementioned constitutive models show different abilities to accurately represent the structural response of the reinforced concrete elements studied.

The Cohesive Crack Model Applied To Notched PMMA Specimens Obeying a Non Linear Behaviour Under Torsion Loading

This article presents a new material model developed with the aim of analyzing failure of blunt notched components made of nonlinear brittle materials. The model, which combines the cohesive crack model with Hencky's theory of total deformations, is used to simulate an experimental benchmark carried out previously by the authors. Such combination is achieved through the embedded crack approach concept. In spite of the unavailability of precise material data, the numerical predictions obtained show good agreement with the experimental results.

A review of load models for the prediction of footfall-induced vibration in structures

Crowd induced dynamic loading in large structures, such as gymnasiums or stadium, is usually modelled as a series of harmonic loads which are defined in terms of their Fourier coefficients. Different values of these coefficients that were obtained from full scale measurements can be found in codes. Recently, an alternative has been proposed, based on random generation of load time histories that take into account phase lag among individuals inside the crowd. This paper presents the testing done on a structure designed to be a gymnasium. Two series of dynamic test were performed on the gym slab. For the first test an electrodynamic shaker was placed at several locations and during the second one people located inside a marked area bounced and jumped guided by different metronome rates. A finite element model (FEM) is presented and a comparison of numerically predicted and experimentally observed vibration modes and frequencies has been used to assess its validity. The second group of measurements will be compared with predictions made using the FEM model and three alternatives for crowd induced load modelling.

Model of human rhythmic activities on structures

Crowd induced dynamic loading in large structures, such as gymnasiums or stadiums, is usually modelled as a series of harmonic loads which are defined in terms of their Fourier coefficients. Different values of these Fourier coefficients that were obtained from full scale measurements can be found in codes. Recently, an alternative has been proposed, based on random generation of load time histories that take into account phase lags among individuals inside the crowd. Generally the testing is performed on platforms or structures that can be considered rigid because their natural frequencies are higher than the excitation frequencies associated with crowd loading. In this paper we shall present the testing done on a structure designed to be a gymnasium, which has natural frequencies within that range. In this test the gym slab was instrumented with acceleration sensors and different people jumped on a force plate installed on the floor. Test results have been compared with predictions based on the two abovementioned load modelling alternatives and a new methodology for modelling jumping loads has been proposed in order to reduce the difference between experimental and numerical results at high frequency range.

Evaluation of damping properties of structural glass panes under impact loading

The latest technology and architectural trends have significantly improved the use of a large variety of glass products in construction which, in function of their own characteristocs, allow to design and calculate structural glass elements under safety conditions. This paper presents the evaluation and analysis of the damping properties of rectangular laminated glass plates of 1.938 m x 0.876 m with different thickness depending on the number of PVB interlayers arranged. By means of numerical simulation and experimental verification, using modal analysis, natural frequencies and damping of the glass plates were calculated, both under free boundary conditions and operational conditions for the impact test equipment used in the experimental program, as the European standard UNE-EN 12600:2003 specifies.

On the loading-rate dependence of the Al 7017-T73 fracture-initiation toughness

While static fracture toughness is a widely studied and standardised parameter, its dynamic counterpart has not been exhaustively examined. Therefore, in this research a series of quasi-static and different loading-rate dynamic tests were carried out to determine the evolution of fracture toughness with the velocity of the application of the load on aluminium 7017-T73 alloy. Three-point bending tests of pre-fatigued standard specimens (ASTM E399) at four loading-rates were carried out. The experiments were conducted by employing the subsequent apparatus ordered from lowest to highest load application velocity: a servo-hydraulic universal testing machine, a free-drop tower, a modified Split Hopkinson Pressure Bar and an explosive load testing device. In order to perform the dynamic fracture toughness tests, it was necessary to design and develop some experimental devices. The fracture-initiation toughness of the aluminium 7017-T73 alloy did not exhibit a significant variation for the studied cases. As a conclusion, the research showed that fracture-initiation toughness remained constant regardless of the velocity at which the load was applied.

Random positions of dendritic spines in human cerebral cortex

Dendritic spines establish most excitatory synapses in the brain and are located in Purkinje cell’s dendrites along helical paths, perhaps maximizing the probability to contact different axons. To test whether spine helixes also occur in neocortex, we reconstructed >500 dendritic segments from adult human cortex obtained from autopsies. With Fourier analysis and spatial statistics, we analyzed spine position along apical and basal dendrites of layer 3 pyramidal neurons from frontal, temporal, and cingulate cortex. Although we occasionally detected helical positioning, for the great majority of dendrites we could not reject the null hypothesis of spatial randomness in spine locations, either in apical or basal dendrites, in neurons of different cortical areas or among spines of different volumes and lengths. We conclude that in adult human neocortex spine positions are mostly random. We discuss the relevance of these results for spine formation and plasticity and their functional impact for cortical circuits.

Neurite, a finite difference large scale parallel program for the simulation of the electrical signal propagation in neurites under mechanical loading

With the growing body of research on traumatic brain injury and spinal cord injury, computational neuroscience has recently focused its modeling efforts on neuronal functional deficits following mechanical loading. However, in most of these efforts, cell damage is generally only characterized by purely mechanistic criteria, function of quantities such as stress, strain or their corresponding rates. The modeling of functional deficits in neurites as a consequence of macroscopic mechanical insults has been rarely explored. In particular, a quantitative mechanically based model of electrophysiological impairment in neuronal cells has only very recently been proposed (Jerusalem et al., 2013). In this paper, we present the implementation details of Neurite: the finite difference parallel program used in this reference. Following the application of a macroscopic strain at a given strain rate produced by a mechanical insult, Neurite is able to simulate the resulting neuronal electrical signal propagation, and thus the corresponding functional deficits. The simulation of the coupled mechanical and electrophysiological behaviors requires computational expensive calculations that increase in complexity as the network of the simulated cells grows. The solvers implemented in Neurite-explicit and implicit-were therefore parallelized using graphics processing units in order to reduce the burden of the simulation costs of large scale scenarios. Cable Theory and Hodgkin-Huxley models were implemented to account for the electrophysiological passive and active regions of a neurite, respectively, whereas a coupled mechanical model accounting for the neurite mechanical behavior within its surrounding medium was adopted as a link between lectrophysiology and mechanics (Jerusalem et al., 2013). This paper provides the details of the parallel implementation of Neurite, along with three different application examples: a long myelinated axon, a segmented dendritic tree, and a damaged axon. The capabilities of the program to deal with large scale scenarios, segmented neuronal structures, and functional deficits under mechanical loading are specifically highlighted.

Random positioning of dendritic spines in the human cerebral cortex

Dendritic spines establish most excitatory synapses in the brain and are located in Purkinje cell?s dendrites along helical paths, perhaps maximizing the probability to contact different axons. To test whether spine helixes also occur in neocortex, we reconstructed ?500 dendritic segments from adult human cortex obtained from autopsies. With Fourier analysis and spatial statistics, we analyzed spine position along apical and basal dendrites of layer 3 pyramidal neurons from frontal, temporal, and cingulate cortex. Although we occasionally detected helical positioning, for the great majority of dendrites we could not reject the null hypothesis of spatial randomness in spine locations, either in apical or basal dendrites, in neurons of different cortical areas or among spines of different volumes and lengths. We conclude that in adult human neocortex spine positions are mostly random. We discuss the relevance of these results for spine formation and plasticity and their functional impact for cortical circuits.

Computer simulation of random parckings for self-similar particle size distributions in soil and granular materials: porosity and pore size distribution

A 2D computer simulation method of random packings is applied to sets of particles generated by a self-similar uniparametric model for particle size distributions (PSDs) in granular media. The parameter p which controls the model is the proportion of mass of particles corresponding to the left half of the normalized size interval [0,1]. First the influence on the total porosity of the parameter p is analyzed and interpreted. It is shown that such parameter, and the fractal exponent of the associated power scaling, are efficient packing parameters, but this last one is not in the way predicted in a former published work addressing an analogous research in artificial granular materials. The total porosity reaches the minimum value for p = 0.6. Limited information on the pore size distribution is obtained from the packing simulations and by means of morphological analysis methods. Results show that the range of pore sizes increases for decreasing values of p showing also different shape in the volume pore size distribution. Further research including simulations with a greater number of particles and image resolution are required to obtain finer results on the hierarchical structure of pore space.