Aziz Laghdir; Yves Fortin; Carmen Mariella De la Cruz; Hernández, Roger E.
2008-01-01
The objective of this study was to develop an accurate and simple method for measuring the engineering coefficients of the 3D elasticity tensor of wood. A method using a semi-ring extensometer (SRE) and a compression specimen (6-specimen technique) is proposed. The SRE is made of a semi-ring stainless steel blade pin-jointed to two aluminum fixing plates, and two resistance strain gauges bonded to the top and bottom faces of the blade at mid-span position. Groups of five matched compression s...
3D reconstruction of tensors and vectors
Defrise, Michel; Gullberg, Grant T.
2005-02-17
Here we have developed formulations for the reconstruction of 3D tensor fields from planar (Radon) and line-integral (X-ray) projections of 3D vector and tensor fields. Much of the motivation for this work is the potential application of MRI to perform diffusion tensor tomography. The goal is to develop a theory for the reconstruction of both Radon planar and X-ray or line-integral projections because of the flexibility of MRI to obtain both of these type of projections in 3D. The development presented here for the linear tensor tomography problem provides insight into the structure of the nonlinear MRI diffusion tensor inverse problem. A particular application of tensor imaging in MRI is the potential application of cardiac diffusion tensor tomography for determining in vivo cardiac fiber structure. One difficulty in the cardiac application is the motion of the heart. This presents a need for developing future theory for tensor tomography in a motion field. This means developing a better understanding of the MRI signal for diffusion processes in a deforming media. The techniques developed may allow the application of MRI tensor tomography for the study of structure of fiber tracts in the brain, atherosclerotic plaque, and spine in addition to fiber structure in the heart. However, the relations presented are also applicable to other fields in medical imaging such as diffraction tomography using ultrasound. The mathematics presented can also be extended to exponential Radon transform of tensor fields and to other geometric acquisitions such as cone beam tomography of tensor fields.
Overall Dynamic Properties of 3-D periodic elastic composites
Srivastava, Ankit
2011-01-01
A method for the homogenization of 3-D periodic elastic composites is presented. It allows for the evaluation of the averaged overall frequency dependent dynamic material constitutive tensors relating the averaged dynamic ?eld variable tensors of velocity, strain, stress, and linear momentum. The formulation is based on micromechanical modeling of a representative unit cell of a composite proposed by Nemat-Nasser & Hori (1993), Nemat-Nasser et. al. (1982) and Mura (1987) and is the 3-D generalization of the 1-D elastodynamic homogenization scheme presented by Nemat-Nasser & Srivastava (2011). We show that for 3-D periodic composites the overall compliance (stiffness) tensor is hermitian, irrespective of whether the corresponding unit cell is geometrically or materially symmetric.Overall mass density is shown to be a tensor and, like the overall compliance tensor, always hermitian. The average strain and linear momentum tensors are, however, coupled and the coupling tensors are shown to be each others'...
3D Shape Similarity Using Vectors of Locally Aggregated Tensors
Tabia, Hedi; Picard, David; Laga, Hamid; Gosselin, Philippe-Henri
2013-01-01
IEEE International Conference on Image Processing 2013 International audience In this paper, we present an efficient 3D object retrieval method invariant to scale, orientation and pose. Our approach is based on the dense extraction of discriminative local descriptors extracted from 2D views. We aggregate the descriptors into a single vector signature using tensor products. The similarity between 3D models can then be efficiently computed with a simple dot product. Experiments on the SHR...
Causal Dynamical Triangulation of 3D Tensor Model
Kawabe, Hiroshi
2016-01-01
We extend the string field theory of the two dimensional (2D) generalized causal dynamical triangulation (GCDT) with the Ishibashi-Kawai (IK-) type interaction formulated by the matrix model, to the three dimensional (3D) model of the surface field theory. Based on the loop gas model, we construct a tensor model for the discretized surface field and then apply it the stochastic quantization method. In the double scaling limit, the model is characterized by two scaling dimensions $D$ and $D_N$, the power indices of the minimal length as the scaling parameter. The continuum GCDT model with the IK-type interaction is realized with the similar restriction in the $D_N$-$D$ space, to the 2D model. The distinct property in the 3D model is that the quantum effect contains the IK-type interaction only, while the ordinary splitting interaction is excluded.
Local Tensor Radiation Conditions For Elastic Waves
Krenk, S.; Kirkegaard, Poul Henning
2001-01-01
A local boundary condition is formulated, representing radiation of elastic waves from an arbitrary point source. The boundary condition takes the form of a tensor relation between the stress at a point on an arbitrarily oriented section and the velocity and displacement vectors at the point....... The tensor relation generalizes the traditional normal incidence impedance condition by accounting for the angle between wave propagation and the surface normal and by including a generalized stiffness term due to spreading of the waves. The effectiveness of the local tensor radiation condition...
General Expression of Elastic Tensor for Anisotropic Materials
HUANG Bo
2005-01-01
In order to formulate a general expression of elastic tensor for anisotropic materials, a method of tensor derivative is used for determining relationship between fourth-order elastic tensor and second-order structure tensor that has satisfied material symmetrical conditions. From this general expression of elastic tensor, specific expressions of elastic tensor for different anisotropic materials, such as isotropic materials, transverse isotropic materials and orthogonal-anisotropic materials, can be deduced. This expression underlies the scalar description of anisotropic factors, which are used for classifying and analyzing anisotropic materials. Cubic crystals are analyzed macroscopically by means of the general expression and anisotropic factor.
3-D joint inversion of the magnetotelluric phase tensor and vertical magnetic transfer functions
Tietze, Kristina; Ritter, Oliver; Egbert, Gary D.
2015-11-01
With advancing computational resources, 3-D inversion techniques have become feasible in recent years and are now a more widely used tool for magnetotelluric (MT) data interpretation. Galvanic distortion caused by small-scale near-surface inhomogeneities remains an obstacle for 3-D MT inversion which so far has experienced little attention. If not considered properly, the effect on 3-D inversion can be immense and result in erroneous subsurface models and interpretations. To tackle the problem we implemented inversion of the distortion-free phase tensor into the ModEM inversion package. The dimensionless phase tensor components describe only variations of the conductivity structure. When inverting these data, particular care has to be taken of the conductivity structure in the a priori model, which provides the reference frame when transferring the information from phase tensors into absolute conductivity values. Our results obtained with synthetic data show that phase tensor inversion can recover the regional conductivity structure in presence of galvanic distortion if the a priori model provides a reasonable assumption for the regional resistivity average. Joint inversion of phase tensor data and vertical magnetic transfer functions improves recovery of the absolute resistivity structure and is less dependent on the prior model. We also used phase tensor inversion for a data set of more than 250 MT sites from the central San Andreas fault, California, where a number of sites showed significant galvanic distortion. We find the regional structure of the phase tensor inversion results compatible with previously obtained models from impedance inversion. In the vicinity of distorted sites, phase tensor inversion models exhibit more homogeneous/smoother conductivity structures.
Elasticity tensor and ultrasonic velocities for anisotropic cubic polycrystal
2008-01-01
The orientation distribution of crystallites in a polycrystal can be described by the orientation distribution function(ODF) . The ODF can be expanded under the Wigner D-bases. The expanded coefficients in the ODF are called the texture coefficients. In this paper,we use the Clebsch-Gordan expression to derive an explicit expression of the elasticity tensor for an anisotropic cubic polycrystal. The elasticity tensor contains three material constants and nine texture coefficients. In order to measure the nine texture coefficients by ultrasonic wave,we give relations between the nine texture coefficients and ultrasonic propagation velocities. We also give a numerical example to check the relations.
Second order Method for Solving 3D Elasticity Equations with Complex and Sharp Interfaces
Wang , Bao; Xia, Kelin; Wei, Guowei
2014-01-01
Elastic materials are ubiquitous in nature and indispensable components in man-made devices and equipments. When a device or equipment involves composite or multiple elastic materials, elasticity interface problems come into play. The solution of three dimensional (3D) elasticity interface problems is significantly more difficult than that of elliptic counterparts due to the coupled vector components and cross derivatives in the governing elasticity equation. This work introduces the matched ...
Centroid Moment Tensor Inversion in a 3D heterogeneous Earth: Application to the Australasian region
Hejrani, B.; Tkalcic, H.; Fichtner, A.
2015-12-01
Australia is surrounded by active complex tectonic belts causing significant seismicity. The recent expansion of permanent seismic networks in the Australasian region provides great opportunity to study Earth structure and a great variety of physical mechanisms responsible for earthquakes.On one hand, a better understanding of the Australasian lithosphere, which is now available through tomographic images from full waveform modelling (Fichtner et al. 2010), provides a powerful tool to scrutinize the determination of earthquake source parameters. Even at relatively long periods (40-200s), the 3D effects of regional structure were shown to significantly alter the global centroid moment tensor solutions (Hingee et al. 2012). Thus, we can now explore other types of uncertainties and test the accuracy of global centroid moment tensor (GCMT) solution for the earthquakes in the Australasian region while checking for the systematic inconsistencies in the solutions. This has a significant bearing on tectonic interpretations. For example, azimuth and plunge of fault planes can be investigated in search for systematic biases.On the other hand, the time has come to take a full advantage of the 3D Earth structural model and embrace ongoing advances in computational power and storage. We develop a semi-automated procedure to calculate the Centroid Moment Tensors in a 3D heterogeneous Earth. We utilize the reciprocity theorem to create Green's functions for point sources covering seismogenic zones of Australasia. We focus on improving the capacity of the method to fully complement the existing monitoring tools at Geosciences Australia. Furthermore, we investigate the effects of detailed velocity structure on Centroid location and double-couple percentages. Moreover Azimuth and Plunge of focal mechanisms in GCMT (Global CMT), were investigated in search for any systematic bias.References: Fichtner, A., Kennett, B.L.N., Igel, H., Bunge, H.-P., 2010. Full waveform tomography for
Motion-induced phase error estimation and correction in 3D diffusion tensor imaging.
Van, Anh T; Hernando, Diego; Sutton, Bradley P
2011-11-01
A multishot data acquisition strategy is one way to mitigate B0 distortion and T2∗ blurring for high-resolution diffusion-weighted magnetic resonance imaging experiments. However, different object motions that take place during different shots cause phase inconsistencies in the data, leading to significant image artifacts. This work proposes a maximum likelihood estimation and k-space correction of motion-induced phase errors in 3D multishot diffusion tensor imaging. The proposed error estimation is robust, unbiased, and approaches the Cramer-Rao lower bound. For rigid body motion, the proposed correction effectively removes motion-induced phase errors regardless of the k-space trajectory used and gives comparable performance to the more computationally expensive 3D iterative nonlinear phase error correction method. The method has been extended to handle multichannel data collected using phased-array coils. Simulation and in vivo data are shown to demonstrate the performance of the method. PMID:21652284
In this paper, we derive the Green tensor of anisotropic gradient elasticity with separable weak non-locality, a special version of Mindlin's form II anisotropic gradient elasticity theory with up to six independent length scale parameters. The framework models materials where anisotropy is twofold, namely the bulk material anisotropy and a weak non-local anisotropy relevant at the nano-scale. In contrast with classical anisotropic elasticity, it is found that both the Green tensor and its gradient are non-singular at the origin, and that they rapidly converge to their classical counterparts away from the origin. Therefore, the Green tensor of Mindlin's anisotropic gradient elasticity with separable weak non-locality can be used as a physically-based regularization of the classical Green tensor for materials with strong anisotropy. - Highlights: • Theory of Mindlin's anisotropic gradient elasticity with separable weak non-locality is presented. • The non-singular (3D) Green tensor is given. • The gradient of the non-singular Green tensor is calculated
Eshelby tensor for a crack in an orthotropic elastic medium
Gruescu, Cosmin; Monchiet, Vincent; Kondo, Djimedo
2005-06-01
In the present Note, we provide new analytical expressions of the components of Hill tensor P (or equivalently the Eshelby tensor S) associated to an arbitrarily oriented crack in orthotropic elastic medium. The crack is modelled as an infinite cylinder along a symmetry axis of the matrix, with low aspect ratio. The three dimensional results obtained show explicitly the interaction between the primary (structural) anisotropy and the crack-induced anisotropy. They are validated by comparison with existing results in the case where the crack is in a symmetry plane. To cite this article: C. Gruescu et al., C. R. Mecanique 333 (2005).
A BOUNDARY INTEGRAL METHOD FOR COMPUTING ELASTIC MOMENT TENSORS FOR ELLIPSES AND ELLIPSOIDS
Habib Ammari; Hyeonbae Kang; Hyundae Lee
2007-01-01
The concept of elastic moment tensor occurs in several interesting contexts, in particular in imaging small elastic inclusions and in asymptotic models of dilute elastic composites.In this paper, we compute the elastic moment tensors for ellipses and ellipsoids by using a systematic method based on layer potentials. Our computations reveal an underlying elegant relation between the elastic moment tensors and the single layer potential.
Proton chemical shift tensors determined by 3D ultrafast MAS double-quantum NMR spectroscopy
Zhang, Rongchun; Mroue, Kamal H.; Ramamoorthy, Ayyalusamy, E-mail: ramamoor@umich.edu [Biophysics and Department of Chemistry, The University of Michigan, Ann Arbor, Michigan 48109-1055 (United States)
2015-10-14
Proton NMR spectroscopy in the solid state has recently attracted much attention owing to the significant enhancement in spectral resolution afforded by the remarkable advances in ultrafast magic angle spinning (MAS) capabilities. In particular, proton chemical shift anisotropy (CSA) has become an important tool for obtaining specific insights into inter/intra-molecular hydrogen bonding. However, even at the highest currently feasible spinning frequencies (110–120 kHz), {sup 1}H MAS NMR spectra of rigid solids still suffer from poor resolution and severe peak overlap caused by the strong {sup 1}H–{sup 1}H homonuclear dipolar couplings and narrow {sup 1}H chemical shift (CS) ranges, which render it difficult to determine the CSA of specific proton sites in the standard CSA/single-quantum (SQ) chemical shift correlation experiment. Herein, we propose a three-dimensional (3D) {sup 1}H double-quantum (DQ) chemical shift/CSA/SQ chemical shift correlation experiment to extract the CS tensors of proton sites whose signals are not well resolved along the single-quantum chemical shift dimension. As extracted from the 3D spectrum, the F1/F3 (DQ/SQ) projection provides valuable information about {sup 1}H–{sup 1}H proximities, which might also reveal the hydrogen-bonding connectivities. In addition, the F2/F3 (CSA/SQ) correlation spectrum, which is similar to the regular 2D CSA/SQ correlation experiment, yields chemical shift anisotropic line shapes at different isotropic chemical shifts. More importantly, since the F2/F1 (CSA/DQ) spectrum correlates the CSA with the DQ signal induced by two neighboring proton sites, the CSA spectrum sliced at a specific DQ chemical shift position contains the CSA information of two neighboring spins indicated by the DQ chemical shift. If these two spins have different CS tensors, both tensors can be extracted by numerical fitting. We believe that this robust and elegant single-channel proton-based 3D experiment provides useful atomistic
Towards an Anisotropic Whole Mantle 3D Elastic Velocity Model
Panning, M. P.; Romanowicz, B.; Gung, Y.
2001-12-01
Many studies have documented the existence of anisotropy in the earth's upper mantle, concentrated in the top 200 km. This evidence comes from the study of surface waves as well as shear wave splitting. There is also evidence for shear wave splitting in D", at least in well sampled regions. There are some hints of anisotropy at the base of the transition zone. Tomographic models of the upper mantle have been developed with simplifying assumptions about the nature of the anisotropy, in order to minimize the number of free parameters in the inversions. Some assume transverse isotropy (e.g Ekström and Dziewonski, 1997), others include additional degrees of freedom with some realistic constraints on mineralogy (e.g. Montagner and Tanimoto, 1991). Our goal is to investigate anisotropy in the whole mantle, using the framework of waveform inversion, and the nonlinear asymptotic mode coupling theory (NACT), previously developed and applied to the construction of whole-mantle SH velocity models (Li and Romanowicz, 1996; Mégnin and Romanowicz, 2000). For this we require a 3 component dataset, and we have extended our automatic transverse (T) component wavepicking procedures to the vertical (Z) and longitudinal (L) component - a non-trivial task given the large number of phases present in the coupled P-SV system. A useful initial assumption, for which the theory has been readily adapted, is that of transverse isotropy. As a first step towards this, we have been investigating inversions using T component and Z,L component data separately. In particular, this allows us to explore the sampling that can be achieved with Z,L component data alone in the deepest part of the mantle. Indeed, D" is in general much better sampled in SH than in SV, owing to the availability of SHdiff at large distances, while SVdiff decays more rapidly due to mantle-core coupling. We present the results of our resolution experiments and discuss the differences between the 3D SV model obtained in well
Liang, Yun; Zhu, Hui; Friedman, Morton H.
2010-11-01
Intravascular ultrasound (IVUS) elastography is a promising tool for studying atherosclerotic plaque composition and assessing plaque vulnerability. Current IVUS elastography techniques can measure the 1D or 2D strain of the vessel wall using various motion tracking algorithms. Since biological soft tissue tends to deform non-uniformly in 3D, measurement of the complete 3D strain tensor is desirable for more rigorous analysis of arterial wall mechanics. In this paper, we extend our previously developed method of 2D arterial wall strain measurement based on non-rigid image registration into 3D strain measurement. The new technique registers two image volumes acquired from the same vessel segment under different levels of luminal pressure and longitudinal stress. The 3D displacement field obtained from the image registration is used to calculate the local 3D strain tensor. From the 3D strain tensor, radial, circumferential and longitudinal strain distributions can be obtained and displayed. This strain tensor measurement method is validated and evaluated using IVUS images of healthy porcine carotid arteries subjected to a luminal pressure increase and longitudinal stretch. The ability of the algorithm to overcome systematic noise was tested, as well as the consistency of the results under different longitudinal frame resolutions.
3D Elastic Registration of Ultrasound Images Based on Skeleton Feature
LI Dan-dan; LIU Zhi-Yan; SHEN Yi
2005-01-01
In order to eliminate displacement and elastic deformation between images of adjacent frames in course of 3D ultrasonic image reconstruction, elastic registration based on skeleton feature was adopt in this paper. A new automatically skeleton tracking extract algorithm is presented, which can extract connected skeleton to express figure feature. Feature points of connected skeleton are extracted automatically by accounting topical curvature extreme points several times. Initial registration is processed according to barycenter of skeleton. Whereafter, elastic registration based on radial basis function are processed according to feature points of skeleton. Result of example demonstrate that according to traditional rigid registration, elastic registration based on skeleton feature retain natural difference in shape for organ's different part, and eliminate slight elastic deformation between frames caused by image obtained process simultaneously. This algorithm has a high practical value for image registration in course of 3D ultrasound image reconstruction.
ELASTIC BEHAVIOR ANALYSIS OF 3D ANGLE-INTERLOCK WOVEN CERAMIC COMPOSITES
Chang Yanjun; Jiao Guiqiong; Wang Bo; Liu Wei
2006-01-01
A micromechanical model for elastic behavior analysis of angle-interlock woven ceramic composites is proposed in this paper. This model takes into account the actual fabric structure by considering the fiber undulation and continuity in space, the cavities between adjacent yarns and the actual cross-section geometry of the yarn. Based on the laminate theory, the elastic properties of 3D angle-interlock woven ceramic composites are predicted. Different numbers of interlaced wefts have almost the same elastic moduli. The thickness of ceramic matrix has little effect on elastic moduli. When the undulation ratio increases longitudinal modulus decreases and the other Young's moduli increase. Good agreement between theoretical predictions and experimental results demonstrates the feasibility of the proposed model in analyzing the elastic properties of3D angle-interlock woven ceramic composites. The results of this paper verify the fact that the method of analyzing polyester matrix composites is suitable for woven ceramic composites.
A comprehensive statistical framework for elastic shape analysis of 3D faces
Kurtek, Sebastian; Drira, Hassen
2015-01-01
We develop a comprehensive statistical framework for analyzing shapes of 3D faces. In particular, we adapt a recent elastic shape analysis framework to the case of hemispherical surfaces, and explore its use in a number of processing applications. This framework provides a parameterization-invariant, elastic Riemannian metric, which allows the development of mathematically rigorous tools for statistical analysis. Specifically, this paper describes methods for registration, comparison and defo...
Laplace-domain waveform modeling and inversion for the 3D acoustic-elastic coupled media
Shin, Jungkyun; Shin, Changsoo; Calandra, Henri
2016-06-01
Laplace-domain waveform inversion reconstructs long-wavelength subsurface models by using the zero-frequency component of damped seismic signals. Despite the computational advantages of Laplace-domain waveform inversion over conventional frequency-domain waveform inversion, an acoustic assumption and an iterative matrix solver have been used to invert 3D marine datasets to mitigate the intensive computing cost. In this study, we develop a Laplace-domain waveform modeling and inversion algorithm for 3D acoustic-elastic coupled media by using a parallel sparse direct solver library (MUltifrontal Massively Parallel Solver, MUMPS). We precisely simulate a real marine environment by coupling the 3D acoustic and elastic wave equations with the proper boundary condition at the fluid-solid interface. In addition, we can extract the elastic properties of the Earth below the sea bottom from the recorded acoustic pressure datasets. As a matrix solver, the parallel sparse direct solver is used to factorize the non-symmetric impedance matrix in a distributed memory architecture and rapidly solve the wave field for a number of shots by using the lower and upper matrix factors. Using both synthetic datasets and real datasets obtained by a 3D wide azimuth survey, the long-wavelength component of the P-wave and S-wave velocity models is reconstructed and the proposed modeling and inversion algorithm are verified. A cluster of 80 CPU cores is used for this study.
3D inversion of full gravity gradient tensor data using SL0 sparse recovery
Meng, Zhaohai
2016-04-01
We present a new method dedicated to the interpretation of full gravity gradient tensor data, based on SL0 sparse recovery inversion. The SL0 sparse recovery method aims to find out the minimum value of the objective function to fit the data function and to solve the non-zero solution to the objective function. Based on continuous iteration, we can easily obtain the final global minimum (namely the property and space attribute of the inversion target). We consider which type of tensor data combination produces the best inversion results based on the inversion results of different full gravity gradient tensor data combinations (separate tensor data and combined tensor data). We compare the recovered models obtained by inverting the different combinations of different gravity gradient tensor components to understand how different component combinations contribute to the resolution of the recovered model. Based on the comparison between the SL0 sparse recovery inversion results and the smoothest and focusing inversion results of the full gravity gradient tensor data, we show that SL0 sparse recovery inversion can obtain more stable and efficient inversion results with relatively sharp edge information, and that this method can also produce a stable solution of the inverse problem for complex geological structures. This new method to resolve very large full gravity gradient tensor datasets has the considerable advantage of being highly efficient; the full gravity gradient tensor inversion requires very little time. This new method is very effective in explaining the full gravity tensor which is very sensitive to small changes in local anomaly. The numerical simulation and inversion results of the compositional model indicates that including multiple components for inversion increases the resolution of the recovered density model and improves the structure delineation. We apply our inversion method to invert the gravity gradient tensor survey data from the Vinton salt
On Prediction of 3d Stress State in Elastic Shell by Higher-order Shell Formulations
Brank, Boštjan; Ibrahimbegović, Adnan; Bohinc, Uroš
2008-01-01
In this work we study the accuracy of modem higher-order shell finite element formulations in computation of 3d stress state in elastic shells. In that sense we compare three higher-order shell models: (i) with seven dislacement-like kinematic parameters, and (ii, iii) with six displacement-like kinematic parameters plus one strain-like kinematic parameter introduced by two different versions of enhanced assumed strain (EAS) concept. The finite element approximations of all shell models are b...
Freedman, Daniel A.; Roundy, D.; Arias, T. A.
2008-01-01
We present a study of the local strain effects associated with vacancy defects in strontium titanate and report the first calculations of elastic dipole tensors and chemical strains for point defects in perovskites. The combination of local and long-range results will enable determination of x-ray scattering signatures that can be compared with experiments. We find that the oxygen vacancy possesses a special property -- a highly anisotropic elastic dipole tensor which almost vanishes upon ave...
X-ray scattering in the elastic regime as source for 3D imaging reconstruction technique
Kocifaj, Miroslav; Mego, Michal
2015-11-01
X-ray beams propagate across a target object before they are projected onto a regularly spaced array of detectors to produce a routine X-ray image. A 3D attenuation coefficient distribution is obtained by tomographic reconstruction where scattering is usually regarded as a source of parasitic signals which increase the level of electromagnetic noise that is difficult to eliminate. However, the elastically scattered radiation could be a valuable source of information, because it can provide a 3D topology of electron densities and thus contribute significantly to the optical characterization of the scanned object. The scattering and attenuation data form a complementary base for concurrent retrieval of both electron density and attenuation coefficient distributions. In this paper we developed the 3D reconstruction method that combines both data inputs and produces better image resolution compared to traditional technology.
Moussawi, Ali
2015-02-24
Summary: The post-treatment of (3D) displacement fields for the identification of spatially varying elastic material parameters is a large inverse problem that remains out of reach for massive 3D structures. We explore here the potential of the constitutive compatibility method for tackling such an inverse problem, provided an appropriate domain decomposition technique is introduced. In the method described here, the statically admissible stress field that can be related through the known constitutive symmetry to the kinematic observations is sought through minimization of an objective function, which measures the violation of constitutive compatibility. After this stress reconstruction, the local material parameters are identified with the given kinematic observations using the constitutive equation. Here, we first adapt this method to solve 3D identification problems and then implement it within a domain decomposition framework which allows for reduced computational load when handling larger problems.
崔洁; 江权; 冯夏庭; 李邵军; 高红; 李帅军
2016-01-01
Discontinuities constitute an integral part of rock mass and inherently affect its anisotropic deformation behavior. This work focuses on the equivalent elastic deformation of rock mass with multiple persistent joint sets. A new method based on the space geometric and mechanical properties of the modified crack tensor is proposed, providing an analytical solution for the equivalent elastic compliance tensor of rock mass. A series of experiments validate the capability of the compliance tensor to accurately represent the deformation of rock mass with multiple persistent joint sets, based on conditions set by the basic hypothesis. The spatially varying rules of the equivalent elastic parameters of rock mass with a single joint set are analyzed to reveal the universal law of the stratified rock mass.
Numerical study of elastic turbulence in a 3D curvilinear micro-channel
Zhang, Hongna; Kunugi, Tomoaki; Li, Fengchen
2012-11-01
Elastic turbulence is an intriguing phenomenon of viscoelastic fluid flow, and dominated by the strong nonlinear elasticity due to the existence of flexible microstructures. It implies the possibility to generate a turbulent state (so-called an elastic turbulence) in the micro-scale devices by introducing the viscoelastic fluids, which could significantly enhance the mixing efficiency therein. Several experiments have been carried out to study its characteristics and underlying physics. However, the difficulty in measuring the flow information and behaviors of the microstructures, especially in the cross section normal to the mean flow direction, limits our current understanding and controlling. In the present study, the nondimensionalization method in which the characteristic velocity is defined as the ratio of the solution viscosity to the width of the channel was adopted to simulate the elastic turbulence in the micro-scale devices. And the elastic turbulent flow was obtained numerically in the 3D curvilinear micro-channel. Therein, the characteristics of the velocity field and polymer's behavior are discussed. Moreover, the energy transfer between the kinetic energy and the polymer's elastic energy is also investigated to understand its physical mechanism. Supported by the Japan Society for the Promotion of Science research fellowship and the Ministry of Education, Culture, Sports, Science and Technology via `Energy Science in the Age of Global Warming' of Global Center of Excellence (G-COE) program (J-051).
Ultrasound elastic tensor imaging: comparison with MR diffusion tensor imaging in the myocardium
Lee, Wei-Ning; Larrat, Benoît; Pernot, Mathieu; Tanter, Mickaël
2012-08-01
We have previously proven the feasibility of ultrasound-based shear wave imaging (SWI) to non-invasively characterize myocardial fiber orientation in both in vitro porcine and in vivo ovine hearts. The SWI-estimated results were in good correlation with histology. In this study, we proposed a new and robust fiber angle estimation method through a tensor-based approach for SWI, coined together as elastic tensor imaging (ETI), and compared it with magnetic resonance diffusion tensor imaging (DTI), a current gold standard and extensively reported non-invasive imaging technique for mapping fiber architecture. Fresh porcine (n = 5) and ovine (n = 5) myocardial samples (20 × 20 × 30 mm3) were studied. ETI was firstly performed to generate shear waves and to acquire the wave events at ultrafast frame rate (8000 fps). A 2.8 MHz phased array probe (pitch = 0.28 mm), connected to a prototype ultrasound scanner, was mounted on a customized MRI-compatible rotation device, which allowed both the rotation of the probe from -90° to 90° at 5° increments and co-registration between two imaging modalities. Transmural shear wave speed at all propagation directions realized was firstly estimated. The fiber angles were determined from the shear wave speed map using the least-squares method and eigen decomposition. The test myocardial sample together with the rotation device was then placed inside a 7T MRI scanner. Diffusion was encoded in six directions. A total of 270 diffusion-weighted images (b = 1000 s mm-2, FOV = 30 mm, matrix size = 60 × 64, TR = 6 s, TE = 19 ms, 24 averages) and 45 B0 images were acquired in 14 h 30 min. The fiber structure was analyzed by the fiber-tracking module in software, MedINRIA. The fiber orientation in the overlapped myocardial region which both ETI and DTI accessed was therefore compared, thanks to the co-registered imaging system. Results from all ten samples showed good correlation (r2 = 0.81, p 0.05, unpaired, one-tailed t-test, N = 10). In
Ultrasound elastic tensor imaging: comparison with MR diffusion tensor imaging in the myocardium
We have previously proven the feasibility of ultrasound-based shear wave imaging (SWI) to non-invasively characterize myocardial fiber orientation in both in vitro porcine and in vivo ovine hearts. The SWI-estimated results were in good correlation with histology. In this study, we proposed a new and robust fiber angle estimation method through a tensor-based approach for SWI, coined together as elastic tensor imaging (ETI), and compared it with magnetic resonance diffusion tensor imaging (DTI), a current gold standard and extensively reported non-invasive imaging technique for mapping fiber architecture. Fresh porcine (n = 5) and ovine (n = 5) myocardial samples (20 × 20 × 30 mm3) were studied. ETI was firstly performed to generate shear waves and to acquire the wave events at ultrafast frame rate (8000 fps). A 2.8 MHz phased array probe (pitch = 0.28 mm), connected to a prototype ultrasound scanner, was mounted on a customized MRI-compatible rotation device, which allowed both the rotation of the probe from −90° to 90° at 5° increments and co-registration between two imaging modalities. Transmural shear wave speed at all propagation directions realized was firstly estimated. The fiber angles were determined from the shear wave speed map using the least-squares method and eigen decomposition. The test myocardial sample together with the rotation device was then placed inside a 7T MRI scanner. Diffusion was encoded in six directions. A total of 270 diffusion-weighted images (b = 1000 s mm−2, FOV = 30 mm, matrix size = 60 × 64, TR = 6 s, TE = 19 ms, 24 averages) and 45 B0 images were acquired in 14 h 30 min. The fiber structure was analyzed by the fiber-tracking module in software, MedINRIA. The fiber orientation in the overlapped myocardial region which both ETI and DTI accessed was therefore compared, thanks to the co-registered imaging system. Results from all ten samples showed good correlation (r2 = 0.81, p 0.05, unpaired, one-tailed t-test, N = 10
2D Riemann-Christoffel curvature tensor via a 3D space using a specialized permutation scheme
Omerbashich, Mensur
2009-01-01
When a space in which Christoffel symbols of the second kind are symmetrical in lower indices exists, it makes for a supplement to the standard procedure when a 2D surface is normally induced from the geometry of the surrounding 3D space in which the surface is embedded. There it appears appropriate to use a scheme for straightforward permutation of indices of Gkij, when such a space would make this transformation possible, so as to obtain the components of the 2D Riemann-Christoffel tensor (here expressed in geodetic coordinates for an ellipsoid of revolution, of use in geophysics). By applying my scheme I find the corresponding indices in 2D and 3D supplement-spaces, and I compute components of the Riemann-Christoffel tensor. By operating over the elements of the projections alone, the all-known value of 1/MN for the Gaussian curvature on an ellipsoid of revolution is obtained. To further validate my scheme, I show that in such a 3D space the tangent vector to a PHI-curve for LAM=const1 would be parallel to...
Mathematical Description and Finite Element Equation of 3D Coupled Thermo-elastic Contact Problem
Shi Yu; Xiao Yougang; Chen Guoxin
2006-01-01
Through defining slide yield function and floating potential function of thermo-contact surface, the complementary equation of thermo-contact boundary has been reached, the fundamental equations to solve 3D thermo-contact coupled problem have been listed. On this foundation, the finite element equation and definite solution condition of contact heat transfer have been given out. Based on virtual work principle and contact element technology, the finite element equation of 3D elastic contact system has been deduced under the effect of thermal stress. The pseudo load brought by contact gap have been introduced into this equation in order to reflect the contact state change. During iteration, once contact rigidity matrix is formed, it won't change,which will make calculation reduce greatly.
2D Riemann-Christoffel curvature tensor via a 3D space using a specialized permutation scheme
Omerbashich, Mensur
2009-01-01
When a space in which Christoffel symbols of the second kind are symmetrical in lower indices exists, it makes for a supplement to the standard procedure when a 2D surface is normally induced from the geometry of the surrounding 3D space in which the surface is embedded. There it appears appropriate to use a scheme for straightforward permutation of indices of Gkij, when such a space would make this transformation possible, so as to obtain the components of the 2D Riemann-Christoffel tensor (...
3D PATTERN OF BRAIN ABNORMALITIES IN WILLIAMS SYNDROME VISUALIZED USING TENSOR-BASED MORPHOMETRY
Chiang, Ming-Chang; Reiss, Allan L.; Lee, Agatha D.; Bellugi, Ursula; Galaburda, Albert M.; Korenberg, Julie R.; Mills, Debra L.; Toga, Arthur W.; Thompson, Paul M.
2007-01-01
Williams syndrome (WS) is a neurodevelopmental disorder associated with deletion of ~20 contiguous genes in chromosome band 7q11.23. Individuals with WS exhibit mild to moderate mental retardation, but are relatively more proficient in specific language and musical abilities. We used tensor-based morphometry (TBM) to visualize the complex pattern of gray/white matter reductions in WS, based on fluid registration of structural brain images.
Jammed elastic shells - a 3D experimental soft frictionless granular system
Jose, Jissy; Blab, Gerhard A.; van Blaaderen, Alfons; Imhof, Arnout
2015-03-01
We present a new experimental system of monodisperse, soft, frictionless, fluorescent labelled elastic shells for the characterization of structure, universal scaling laws and force networks in 3D jammed matter. The interesting fact about these elastic shells is that they can reversibly deform and therefore serve as sensors of local stress in jammed matter. Similar to other soft particles, like emulsion droplets and bubbles in foam, the shells can be packed to volume fractions close to unity, which allows us to characterize the contact force distribution and universal scaling laws as a function of volume fraction, and to compare them with theoretical predictions and numerical simulations. However, our shells, unlike other soft particles, deform rather differently at large stresses. They deform without conserving their inner volume, by forming dimples at contact regions. At each contact one of the shells buckled with a dimple and the other remained spherical, closely resembling overlapping spheres. We conducted 3D quantitative analysis using confocal microscopy and image analysis routines specially developed for these particles. In addition, we analysed the randomness of the process of dimpling, which was found to be volume fraction dependent.
Panebianco, Valeria, E-mail: valeria.panebianco@gmail.com [Department of Radiological Sciences, Oncology and Anatomical Pathology, Sapienza University of Rome, Viale Regina Elena, 324, 00161 Rome (Italy); Barchetti, Flavio [Department of Radiological Sciences, Oncology and Anatomical Pathology, Sapienza University of Rome, Viale Regina Elena, 324, 00161 Rome (Italy); Sciarra, Alessandro [Department of Urology, Sapienza University of Rome (Italy); Marcantonio, Andrea; Zini, Chiara [Department of Radiological Sciences, Oncology and Anatomical Pathology, Sapienza University of Rome, Viale Regina Elena, 324, 00161 Rome (Italy); Salciccia, Stefano [Department of Urology, Sapienza University of Rome (Italy); Collettini, Federico [Department of Radiology, Charité, Berlin (Germany); Gentile, Vincenzo [Department of Urology, Sapienza University of Rome (Italy); Hamm, Bernard [Department of Radiology, Charité, Berlin (Germany); Catalano, Carlo [Department of Radiological Sciences, Oncology and Anatomical Pathology, Sapienza University of Rome, Viale Regina Elena, 324, 00161 Rome (Italy)
2013-10-01
Objectives: To evaluate if Diffusion Tensor Imaging technique (DTI) can improve the visualization of periprostatic nerve fibers describing the location and distribution of entire neurovascular plexus around the prostate in patients who are candidates for prostatectomy. Materials and methods: Magnetic Resonance Imaging (MRI), including a 2D T2-weighted FSE sequence in 3 planes, 3D T2-weighted and DTI using 16 gradient directions and b = 0 and 1000, was performed on 36 patients. Three out of 36 patients were excluded from the analysis due to poor image quality (blurring N = 2, artifact N = 1). The study was approved by local ethics committee and all patients gave an informed consent. Images were evaluated by two radiologists with different experience in MRI. DTI images were analyzed qualitatively using dedicated software. Also 2D and 3D T2 images were independently considered. Results: 3D-DTI allowed description of the entire plexus of the periprostatic nerve fibers in all directions, while 2D and 3D T2 morphological sequences depicted part of the fibers, in a plane by plane analysis of fiber courses. DTI demonstrated in all patients the dispersion of nerve fibers around the prostate on both sides including the significant percentage present in the anterior and anterolateral sectors. Conclusions: DTI offers optimal representation of the widely distributed periprostatic plexus. If validated, it may help guide nerve-sparing radical prostatectomy.
Objectives: To evaluate if Diffusion Tensor Imaging technique (DTI) can improve the visualization of periprostatic nerve fibers describing the location and distribution of entire neurovascular plexus around the prostate in patients who are candidates for prostatectomy. Materials and methods: Magnetic Resonance Imaging (MRI), including a 2D T2-weighted FSE sequence in 3 planes, 3D T2-weighted and DTI using 16 gradient directions and b = 0 and 1000, was performed on 36 patients. Three out of 36 patients were excluded from the analysis due to poor image quality (blurring N = 2, artifact N = 1). The study was approved by local ethics committee and all patients gave an informed consent. Images were evaluated by two radiologists with different experience in MRI. DTI images were analyzed qualitatively using dedicated software. Also 2D and 3D T2 images were independently considered. Results: 3D-DTI allowed description of the entire plexus of the periprostatic nerve fibers in all directions, while 2D and 3D T2 morphological sequences depicted part of the fibers, in a plane by plane analysis of fiber courses. DTI demonstrated in all patients the dispersion of nerve fibers around the prostate on both sides including the significant percentage present in the anterior and anterolateral sectors. Conclusions: DTI offers optimal representation of the widely distributed periprostatic plexus. If validated, it may help guide nerve-sparing radical prostatectomy
Spatial Parallelism of a 3D Finite Difference, Velocity-Stress Elastic Wave Propagation Code
MINKOFF,SUSAN E.
1999-12-09
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately. finite difference simulations for 3D elastic wave propagation are expensive. We model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MP1 library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speed up. Because i/o is handled largely outside of the time-step loop (the most expensive part of the simulation) we have opted for straight-forward broadcast and reduce operations to handle i/o. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ''ghost cells''. When this communication is balanced against computation by allocating subdomains of reasonable size, we observe excellent scaled speed up. Allocating subdomains of size 25 x 25 x 25 on each node, we achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
Simulation domain size requirements for elastic response of 3D polycrystalline materials
Ozturk, Tugce; Stein, Clayton; Pokharel, Reeju; Hefferan, Christopher; Tucker, Harris; Jha, Sushant; John, Reji; Lebensohn, Ricardo A.; Kenesei, Peter; Suter, Robert M.; Rollett, Anthony D.
2016-01-01
A fast Fourier transform (FFT) based spectral algorithm is used to compute the full field mechanical response of polycrystalline microstructures. The field distributions in a specific region are used to determine the sensitivity of the method to the number of surrounding grains through quantification of the divergence of the field values from the largest simulation domain, as successively smaller surrounding volumes are included in the simulation. The analysis considers a mapped 3D structure where the location of interest is taken to be a particular pair of surface grains that enclose a small fatigue crack, and synthetically created statistically representative microstructures to further investigate the effect of anisotropy, loading condition, loading direction, and texture. The synthetic structures are generated via DREAM3D and the measured material is a cyclically loaded, Ni-based, low solvus high refractory (LSHR) superalloy that was characterized via 3D high energy x-ray diffraction microscopy (HEDM). Point-wise comparison of distributions in the grain pairs shows that, in order to obtain a Pearson correlation coefficient larger than 99%, the domain must extend to at least the third nearest neighbor. For an elastic FFT calculation, the stress-strain distributions are not sensitive to the shape of the domain. The main result is that convergence can be specified in terms of the number of grains surrounding a region of interest.
Irreducible decomposition of strain gradient tensor in isotropic strain gradient elasticity
Lazar, Markus
2016-01-01
In isotropic strain gradient elasticity, we decompose the strain gradient tensor into its irreducible pieces under the n-dimensional orthogonal group O(n). Using the Young tableau method for traceless tensors, four irreducible pieces (n>2), which are canonical, are obtained. In three dimensions, the strain gradient tensor can be decomposed into four irreducible pieces with 7+5+3+3 independent components whereas in two dimensions, the strain gradient tensor can be decomposed into three irreducible pieces with 2+2+2 independent components. The knowledge of these irreducible pieces is extremely useful when setting up constitutive relations and strain energy.
Green-Naghdi rate of the Kirchhoff stress and deformation rate: the elasticity tensor
Bellini, Chiara; Federico, Salvatore
2015-06-01
The elasticity tensor providing the power-conjugation of the Green-Naghdi rate of the Kirchhoff stress and the deformation rate is required, e.g. by the commercially available Finite Element package ABAQUS/Standard for the material user subroutine UMAT, used to input material behaviours other than those included in the libraries of the package. This elasticity tensor had been studied in the literature, but its symmetries have only been briefly discussed, and only its component form in Cartesian coordinates was known. In this work, we derived a covariant, component-free expression of this elasticity tensor and thoroughly studied its symmetries. We found that, although symmetry on both pair of feet (indices) has been deemed to be desirable in the literature, the expression of the tensor available to-date in fact possesses only symmetry on the first pair of feet (indices), whereas the second pair lacks symmetry, and therefore carries a skew-symmetric contribution. This contribution is unnecessary, as it is automatically filtered in the contraction of the elasticity tensor with the symmetric deformation rate tensor. In order to avoid carrying this unnecessary skew-symmetric contribution in the computations, we employ a tensor identity that naturally symmetrises the second pair of feet of the elasticity tensor. We demonstrated the validity and robustness of the implementation of the user-defined material based on this tensor representation by simulating a benchmark problem consisting in biaxial tests of porcine and human atrial tissue, with material properties taken from previously performed experiments. We compared the results obtained by means of our user-defined material and those obtained through an equivalent built-in material, and obtained identical results.
Compact Vectors of Locally Aggregated Tensors for 3D shape retrieval
Tabia, Hedi; Picard, David; Laga, Hamid; Gosselin, Philippe-Henri
2013-01-01
International audience During the last decade, a significant attention has been paid, by the computer vision and the computer graphics communities, to three dimensional (3D) object retrieval. Shape retrieval methods can be divided into three main steps: the shape descriptors extraction, the shape signatures and their associated similarity measures, and the machine learning relevance functions. While the first and the last points have vastly been addressed in recent years, in this paper, we...
GREEN＇S FUNCTION AND EFFECTIVE ELASTIC STIFFNESS TENSOR FOR ARBITRARY AGGREGATES OF CUBIC CRYSTALS
HuangMojia; ZhengChaomei
2004-01-01
A closed but approximate formula of Green's function for an arbitrary aggregate of cubic crystallites is given to derive the effective elastic stiffness tensor of the polycrystal. This formula, which includes three elastic constants of single cubic crystal and five texture coefficients,accounts for the effects of the orientation distribution function (ODF) up to terms linear in the texture coefficients. Thus it is expected that our formula would be applicable to arbitrary aggregates with weak texture or to materials such as aluminum whose single crystal has weak anisotropy.Three examples are presented to compare predictions from our formula with those from Nishioka and Lothe's formula and Synge's contour integral through numerical integration. As an application of Green's function, we briefly describe the procedure of deriving the effective elastic stiffness tensor for an orthorhombic aggregate of cubic crystallites. The comparison of the computational results given by the finite element method and our effective elastic stiffness tensor is made by an example.
Spatio-Temporal Video Object Segmentation via Scale-Adaptive 3D Structure Tensor
Hai-Yun Wang
2004-06-01
Full Text Available To address multiple motions and deformable objects' motions encountered in existing region-based approaches, an automatic video object (VO segmentation methodology is proposed in this paper by exploiting the duality of image segmentation and motion estimation such that spatial and temporal information could assist each other to jointly yield much improved segmentation results. The key novelties of our method are (1 scale-adaptive tensor computation, (2 spatial-constrained motion mask generation without invoking dense motion-field computation, (3 rigidity analysis, (4 motion mask generation and selection, and (5 motion-constrained spatial region merging. Experimental results demonstrate that these novelties jointly contribute much more accurate VO segmentation both in spatial and temporal domains.
Yu eCai
2011-10-01
Full Text Available High resolution diffusion tensor imaging (DTI can provide important information on brain development, yet it is challenging in live neonatal rats due to the small size of neonatal brain and motion-sensitive nature of DTI. Imaging in live neonatal rats has clear advantages over fixed brain scans, as longitudinal and functional studies would be feasible to understand neuro-developmental abnormalities. In this study, we developed imaging strategies that can be used to obtain high resolution 3D DTI images in live neonatal rats at postnatal day 5 (PND5 and postnatal day 14 (PND14, using only 3 hours of imaging acquisition time. An optimized 3D DTI pulse sequence and appropriate animal setup to minimize physiological motion artifacts are the keys to successful high resolution 3D DTI imaging. Thus, a 3D RARE DTI sequence with twin navigator echoes was implemented to accelerate imaging acquisition time and minimize motion artifacts. It has been suggested that neonatal mammals possess a unique ability to tolerate mild to moderate hypothermia and hypoxia without long term impact. Thus, we additionally utilized this ability to minimize motion artifacts in MR images by carefully suppressing the respiratory rate to around 15/min for PND5 and 30/min for PND14 using mild to moderate hypothermia. These imaging strategies have been successfully implemented to study how the effect of cocaine exposure in dams might affect brain development in their rat pups. Image quality resulting from this in vivo DTI study was comparable to ex vivo scans. FA values were also similar between the live and fixed brain scans. The capability of acquiring high quality in vivo DTI imaging offers a valuable opportunity to study many neurological disorders in brain development in an authentic living environment.
Mollemans, W.; Schutyser, F.; Nadjmi, N.; Maes, F.; Suetens, P.
2006-03-01
In this paper we present an extensive quantitative validation on 3D facial soft tissue simulation for maxillofacial surgery planning. The study group contained 10 patients. In previous work we presented a new Mass Tensor Model to simulate the new facial appearance after maxillofacial surgery in a fast way. 10 patients were preoperatively CT-scanned and the surgical intervention was planned. 4 months after surgery, a post-operative control CT was acquired. In this study, the simulated facial outlook is compared with post-operative image data. After defining corresponding points between the predicted and actual post-operative facial skin surface, using a variant of the non-rigid TPS-RPM algorithm, distances between these correspondences are quantified and visualized in 3D. As shown, the average median distance measures only 0.60 mm and the average 90% percentile stays below 1.5 mm. We can conclude that our model clearly provides an accurate prediction of the real post-operative outcome and is therefore suitable for use in clinical practice.
Ice-shelf forced vibrations modelled with a full 3-D elastic model
Y. V. Konovalov
2014-12-01
Full Text Available Ice-shelf forced vibrations modelling was performed using a full 3-D finite-difference elastic model, which takes into account sub-ice seawater flow. The sub-ice seawater flow was described by the wave equation, so the ice-shelf flexures result from the hydrostatic pressure perturbations in sub-ice seawater layer. The numerical experiments were performed for idealized ice-shelf geometry, which was considered in the numerical experiments in Holdsworth and Glynn (1978. The ice-plate vibrations were modelled for harmonic ingoing pressure perturbations and for a wide spectrum of the ocean swell periodicities, ranging from infragravity wave periods down to periods of a few seconds (0.004–0.2 Hz. The spectrums for the vibration amplitudes were obtained in this range and are published in this manuscript. The spectrums contain distinct resonant peaks, which corroborate the ability of resonant-like motion in suitable conditions of the forcing. The impact of local irregularities in the ice-shelf geometry to the amplitude spectrums was investigated for idealized sinusoidal perturbations of the ice surface and the sea bottom. The results of the numerical experiments presented in this manuscript, are approximately in agreement with the results obtained by the thin-plate model in the research carried out by Holdsworth and Glynn (1978. In addition, the full model allows to observe 3-D effects, for instance, vertical distribution of the stress components in the plate. In particular, the model reveals the increasing in shear stress, which is neglected in the thin-plate approximation, from the terminus towards the grounding zone with the maximum at the grounding line in the case of considered high-frequency forcing. Thus, the high-frequency forcing can reinforce the tidal impact to the ice-shelf grounding zone additionally exciting the ice fracture there.
Clochard, V.
1998-12-02
3D VSP imaging is nowadays a strategic requirement by petroleum companies. It is used to precise in details the geology close to the well. Because of the lack of redundancy and limited coverage in the data. this kind of technology is more restrictive than surface seismic which allows an investigation at a higher scale. Our contribution was to develop an elastic quantitative imagine (GRT migration) which can be applied to 3 components borehole dataset. The method is similar to the Kirchhoff migration using sophistical weighting of the seismic amplitudes. In reality. GRT migration uses pre-calculated Green functions (travel time. amplitude. polarization). The maps are obtained by 3D ray tracing (wavefront construction) in the velocity model. The migration algorithm works with elementary and independent tasks. which is useful to process different kind of dataset (fixed or moving geophone antenna). The study has been followed with validations using asymptotic analytical solution. The ability of reconstruction in 3D borehole survey has been tested in the Overthrust synthetic model. The application to a real circular 3D VSP shows various problems like velocity model building, anisotropy factor and the preprocessing (deconvolution. wave mode separation) which can destroy seismic amplitudes. An isotropic 3 components preprocessing of the whole dataset allows a better lateral reconstruction. The choice of a big migration aperture can help the reconstruction of strong geological dip in spite of migration smiles. Finally, the methodology can be applied to PS converted waves. (author)
Petit, J.L.
1997-07-21
This thesis is devoted to the inversion of VSP (vertical seismic profile) seismic data in order to determine the elastic properties of horizontally stratified media. The VSP records are computed using the full wave elastic modelling in isotropic and transversely isotropic media using Hankel transform, a finite difference scheme and an inverse Hankel transform algorithm, and the propagation equations are determined and numerically solved; the importance of considering a 3D wave propagation model instead of a 1 D one is emphasized. The theoretical VSP inverse problem is then considered, with the seismic waveform inversion set as a least-squares problem, consisting in recovering the distribution of physical parameters which minimize the misfit between calculated and observed VSP. The corresponding problem requires the knowledge of the source function
Park, Kwan Soo; Lee, Sam Sun; Huh, Kyung Hoe; Yi, Wan Jin; Heo, Min Suk; Choi, Soon Chul [Department of Oral and Maxillofacial Radiology, School of Dentistry, Seoul National University, Seoul (Korea, Republic of)
2008-06-15
To investigate the relationship between 3D bone architectural parameters and direction-related elastic moduli of cancellous bone of mandibular condyle. Two micro-pigs (Micro-pigR, PWG Genetics Korea) were used. Each pig was about 12 months old and weighing around 44 kg. 31 cylindrical bone specimen were obtained from cancellous bone of condyles for 3D analysis and measured by micro-computed tomography. Six parameters were trabecular thickness (Tb.Th), bone specific surface (BS/BV), percent bone volume (BV/TV), structure model index (SMI), degree of anisotropy (DA) and 3-dimensional fractal dimension (3DFD). Elastic moduli of three orthogonal directions (superiorinferior (SI), medial-lateral (ML), andterior-posterior (AP) direction) were calculated through finite element analysis. Elastic modulus of superior-inferior direction was higher than those of other directions. Elastic moduli of 3 orthogonal directions showed different correlation with 3D architectural parameters. Elastic moduli of SI and ML directions showed significant strong to moderate correlation with BV/TV, SMI and 3DFD. Elastic modulus of cancellous bone of pig mandibular condyle was highest in the SI direction and it was supposed that the change into plate-like structure of trabeculae was mainly affected by increase of trabeculae of SI and ML directions.
Embedding SAS approach into conjugate gradient algorithms for asymmetric 3D elasticity problems
Chen, Hsin-Chu; Warsi, N.A. [Clark Atlanta Univ., GA (United States); Sameh, A. [Univ. of Minnesota, Minneapolis, MN (United States)
1996-12-31
In this paper, we present two strategies to embed the SAS (symmetric-and-antisymmetric) scheme into conjugate gradient (CG) algorithms to make solving 3D elasticity problems, with or without global reflexive symmetry, more efficient. The SAS approach is physically a domain decomposition scheme that takes advantage of reflexive symmetry of discretized physical problems, and algebraically a matrix transformation method that exploits special reflexivity properties of the matrix resulting from discretization. In addition to offering large-grain parallelism, which is valuable in a multiprocessing environment, the SAS scheme also has the potential for reducing arithmetic operations in the numerical solution of a reasonably wide class of scientific and engineering problems. This approach can be applied directly to problems that have global reflexive symmetry, yielding smaller and independent subproblems to solve, or indirectly to problems with partial symmetry, resulting in loosely coupled subproblems. The decomposition is achieved by separating the reflexive subspace from the antireflexive one, possessed by a special class of matrices A, A {element_of} C{sup n x n} that satisfy the relation A = PAP where P is a reflection matrix (symmetric signed permutation matrix).
Barth, Andreas
2016-04-01
On January 6, 2016 the Democratic People's Republic of Korea (DPRK) carried out an announced nuclear test, which was the fourth after tests conducted in 2006, 2009, and 2013. An important task in discriminating a man-made explosion and a natural tectonic earthquake is the analysis of seismic waveforms. To determine the isotropic and non-isotropic characteristics of the detonation source, I invert long-period seismic data for the full seismic moment tensor to match the observed seismic signals by synthetic waveforms based on a 3D earth model. Here, I show that the inversion of long-period seismic data of the 2016 test reveals a clear explosive (isotropic) component combined with a significant release of shear energy by the double-couple part of the moment tensor. The short- and long-period waveforms of the recent test are very similar to the previous ones. First data show that the energy release of the recent event on long periods greater than 10 s is enlarged by 20-30% compared to the nuclear test in 2013. As shown previously, the double-couple part of the 2009 event was lower by a factor of 0.55 compared to the explosion in 2013, while the isotropic parts of the nuclear tests in 2009 and 2013 were similar (Barth, 2014). However, the recent test again shows a rather small double-couple part, indicating a lower amount of shear-energy radiation than in 2013. This highlights the importance of considering the release of shear energy in understanding near source damaging effects and the containment of nuclear explosions.
Chu, Chunlei
2009-01-01
We present two Lax‐Wendroff type high‐order time stepping schemes and apply them to solving the 3D elastic wave equation. The proposed schemes have the same format as the Taylor series expansion based schemes, only with modified temporal extrapolation coefficients. We demonstrate by both theoretical analysis and numerical examples that the modified schemes significantly improve the stability conditions.
The influence of 3d-metal alloy additions on the elastic and thermodynamic properties of CuPd3
Huang Shuo; Zhang Chuan-Hui; Sun Jing; Shen Jiang
2013-01-01
Embedded-atom method (EAM) potentials are used to investigate the effects of alloying (e.g.3d-metals) on the trends of elastic and thermodynamic properties for CuPd3 alloy.Our calculated lattice parameter,cohesive energy,and elastic constants of CuPd3 are consistent with the available experimental and theoretical data.The results of elastic constants indicate that all these alloys are mechanically stable.Further mechanical behavior analysis shows that the additions of Cr,Fe,Co,and Ni could improve the hardness of CuPd3 while V could well increase its ductility.Moreover,in order to evaluate the thermodynamic contribution of 3d-metals,the Debye temperature,phonon density of states,and vibrational entropy for CuMPd6 alloy are also investigated.
Revealing the z~2.5 Cosmic Web With 3D Lyman-Alpha Forest Tomography: A Deformation Tensor Approach
Lee, Khee-Gan
2016-01-01
Studies of cosmological objects should take into account their positions within the cosmic web of large-scale structure. Unfortunately, the cosmic web has only been extensively mapped at low-redshifts ($z1$, the required galaxy densities are inaccessible for the foreseeable future, but 3D reconstructions of Lyman-$\\alpha$ forest absorption in closely-separated background QSOs and star-forming galaxies already offer a detailed window into $z\\sim2-3$ large-scale structure. We quantify the utility of such maps for studying the cosmic web by using realistic $z=2.5$ Ly$\\alpha$ forest simulations matched to observational properties of upcoming surveys. A deformation tensor-based analysis is used to classify voids, sheets, filaments and nodes in the flux, which is compared to those determined from the underlying dark matter field. We find an extremely good correspondence, with $70\\%$ of the volume in the flux maps correctly classified relative to the dark matter web, and $99\\%$ classified to within 1 eigenvalue. Thi...
Simultaneous elastic parameter inversion in 2-D/3-D TTI medium combined later arrival times
Bai, Chao-ying; Wang, Tao; Yang, Shang-bei; Li, Xing-wang; Huang, Guo-jiao
2016-04-01
Traditional traveltime inversion for anisotropic medium is, in general, based on a "weak" assumption in the anisotropic property, which simplifies both the forward part (ray tracing is performed once only) and the inversion part (a linear inversion solver is possible). But for some real applications, a general (both "weak" and "strong") anisotropic medium should be considered. In such cases, one has to develop a ray tracing algorithm to handle with the general (including "strong") anisotropic medium and also to design a non-linear inversion solver for later tomography. Meanwhile, it is constructive to investigate how much the tomographic resolution can be improved by introducing the later arrivals. For this motivation, we incorporated our newly developed ray tracing algorithm (multistage irregular shortest-path method) for general anisotropic media with a non-linear inversion solver (a damped minimum norm, constrained least squares problem with a conjugate gradient approach) to formulate a non-linear inversion solver for anisotropic medium. This anisotropic traveltime inversion procedure is able to combine the later (reflected) arrival times. Both 2-D/3-D synthetic inversion experiments and comparison tests show that (1) the proposed anisotropic traveltime inversion scheme is able to recover the high contrast anomalies and (2) it is possible to improve the tomographic resolution by introducing the later (reflected) arrivals, but not as expected in the isotropic medium, because the different velocity (qP, qSV and qSH) sensitivities (or derivatives) respective to the different elastic parameters are not the same but are also dependent on the inclination angle.
Geodynamic background of the 2008 Wenchuan earthquake based on 3D visco-elastic numerical modelling
Liu, Chang; Zhu, Bojing; Yang, Xiaolin; Shi, Yaolin
2016-03-01
The 2008 Wenchuan earthquake (Mw7.9) occurred in the Longmen Shan fault zone. The stress change and crustal deformation during the accumulation period is computed using 3D finite element modelling assuming visco-elastic rheology. Our results support that the eastward movement of the Tibetan Plateau resulting from the India-Eurasia collision is obstructed at the Longmen Shan fault zone by the strong Yangtze craton. In response, the Tibetan ductile crust thickens and accumulates at the contact between the Tibetan Plateau and the Sichuan Basin. This process implies a strong uplift with the rate of about 1.8 mm/a of the upper crust and induces a stress concentration nearly at the bottom of the Longmen Shan fault zone. We believe that the stress concentration in the Longmen Shan fault zone provides a very important geodynamic background of the 2008 Wenchuan earthquake. Using numerical experiments we find that the key factor controlling this stress concentration process is the large viscosity contrast in the middle and lower crusts between the Tibetan Plateau and the Sichuan Basin. The results show that large viscosity contrast in the middle and lower crusts accelerates the stress concentration in the Longmen Shan fault zone. Fast moving lower crustal flow accelerates this stress accumulation process. During the inter-seismic period, spatially the maximum stress accumulation rate of the eastern margin of the Tibetan Plateau is located nearly at the bottom of the brittle upper crust of the Longmen Shan fault zone. The spatial distribution of the stress accumulation along the strike of the Longmen Shan fault zone is as follows: the normal stress decreases while the shear stress increases from southwest to northeast along the Longmen Shan fault zone. This stress distribution explains the thrust motion in the SW and strike-slip motion in the NE during the 2008 Wenchuan earthquake.
Levashov, Valentin A
2015-01-01
It is possible to associate with every atom or molecule in a liquid its own atomic stress tensor. These atomic stress tensors can be used to describe liquids' structures and to investigate the connection between structural and dynamic properties. In particular, atomic stresses allow to address atomic scale correlations relevant to the Green-Kubo expression for viscosity. Previously correlations between the atomic stresses of different atoms were studied using the Cartesian representation of the stress tensors or the representation based on spherical harmonics. In this paper we address structural correlations in a model 3D binary liquid using the eigenvalues and eigenvectors of the atomic stress tensors. Thus correlations relevant to the Green-Kubo expression for viscosity are interpreted in a simple geometric way. On decrease of temperature the changes in the relevant stress correlation function between different atoms are significantly more pronounced than the changes in the pair density function. We demonst...
Saxena, Nishank; Mavko, Gary
2016-03-01
Estimation of elastic rock moduli using 2D plane strain computations from thin sections has several numerical and analytical advantages over using 3D rock images, including faster computation, smaller memory requirements, and the availability of cheap thin sections. These advantages, however, must be weighed against the estimation accuracy of 3D rock properties from thin sections. We present a new method for predicting elastic properties of natural rocks using thin sections. Our method is based on a simple power-law transform that correlates computed 2D thin section moduli and the corresponding 3D rock moduli. The validity of this transform is established using a dataset comprised of FEM-computed elastic moduli of rock samples from various geologic formations, including Fontainebleau sandstone, Berea sandstone, Bituminous sand, and Grossmont carbonate. We note that using the power-law transform with a power-law coefficient between 0.4-0.6 contains 2D moduli to 3D moduli transformations for all rocks that are considered in this study. We also find that reliable estimates of P-wave (Vp) and S-wave velocity (Vs) trends can be obtained using 2D thin sections.
Studying conformally flat spacetimes with an elastic stress energy tensor using 1+3 formalism
Brito, I
2012-01-01
Conformally flat spacetimes with an elastic stress energy tensor given by a diagonal trace-free anisotropic pressure tensor are investigated using 1+3 formalism. We show how the null tetrad Ricci components are related to the pressure components and energy density. The 1+3 Bianchi and Jacobi identities and Einstein field equations are written for this particular case. In general the commutators must be considered since they supply potentially new information on higher order derivatives of the 1+3 quantities. We solve the system for the non rotating case which consist of ODEs of a spatial coordinate.
Tensor Analyzing Powers for Quasi-Elastic Electron Scattering from Deuterium
Z.-L. Zhou; M. Bouwhuis; M. Ferro-Luzzi; E. Passchier; R. Alarcon; M. Anghinolfi; H. Arenhoevel; R. van Bommel; T. Botto; J.F.J. van den Brand; H.J. Bulten; S. Choi; J. Comfort; S.M. Dolfini; R. Ent; C. Gaulard; D.W. Higinbotham; C.W. de Ja ger; E. Konstantinov; J. Lang; W. Leidemann; D.J. de Lange; M.A. Miller; D. Niko lenko; N. Papadakis; I. Passchier; H.R. Poolman; S.G. Popov; I. Rachek; M. Ripan i; E. Six; J.J.M. Steijger; M. Taiuti; O. Unal; N. Vodinas; H. de Vries
1999-01-01
We report on a first measurement of tensor analyzing powers in quasi-elastic electron-deuteron scattering at an average three-momentum transfer of 1.7 fm{sup -1}. Data sensitive to the spin-dependent nucleon density in the deuteron were obtained for missing momenta up to 150 MeV/c with a tensor polarized {sup 2}H target internal to an electron storage ring. The data are well described by a calculation that includes the effects of final-state interaction, meson-exchange and isobar currents, and leading-order relativistic contributions.
Integral identities for a semi-infinite interfacial crack in 2D and 3D elasticity
Piccolroaz, Andrea
2011-01-01
The paper is concerned with the problem of a semi-infinite crack at the interface between two dissimilar elastic half-spaces, loaded by a general asymmetrical system of forces distributed along the crack faces. On the basis of the weight function approach and the fundamental reciprocal identity (Betti formula), we formulate the elasticity problem in terms of singular integral equations relating the applied loading and the resulting crack opening. Such formulation is fundamental in the theory of elasticity and extensively used to solve several problems in linear elastic fracture mechanics (for instance various classic crack problems in homogeneous and heterogeneous media). This formulation is also crucial in important recent multiphysics applications, where the elastic problem is coupled with other concurrent physical phenomena. A paradigmatic example is hydraulic fracturing, where the elasticity equations are coupled with fluid dynamics.
ELATE: an open-source online application for analysis and visualization of elastic tensors
Gaillac, Romain; Pullumbi, Pluton; Coudert, François-Xavier
2016-07-01
We report on the implementation of a tool for the analysis of second-order elastic stiffness tensors, provided with both an open-source Python module and a standalone online application allowing the visualization of anisotropic mechanical properties. After describing the software features, how we compute the conventional elastic constants and how we represent them graphically, we explain our technical choices for the implementation. In particular, we focus on why a Python module is used to generate the HTML web page with embedded Javascript for dynamical plots.
ELATE: an open-source online application for analysis and visualization of elastic tensors.
Gaillac, Romain; Pullumbi, Pluton; Coudert, François-Xavier
2016-07-13
We report on the implementation of a tool for the analysis of second-order elastic stiffness tensors, provided with both an open-source Python module and a standalone online application allowing the visualization of anisotropic mechanical properties. After describing the software features, how we compute the conventional elastic constants and how we represent them graphically, we explain our technical choices for the implementation. In particular, we focus on why a Python module is used to generate the HTML web page with embedded Javascript for dynamical plots. PMID:27199239
Levine, Lyle; Okoro, Chukwudi; Lee, I-Fang; Phan, Thien; Kassner, Michael; Xu, Ruqing; Tischler, Jon; Liu, Wenjun
2014-01-01
The use of depth resolved, submicrometer X-ray beams for studying deformation microstructures in plastically deformed metals has come a long way over the past 5 years. We can identify phases, measure crystallographic orientations, and measure lattice constants from buried, submicrometer sample volumes throughout extended sample regions within single crystal and polycrystalline samples. In special cases, we can also measure both deviatoric and complete elastic strain tensors with reliable unce...
无
2010-01-01
The squirrel-cage elastic support is one of the most important components of an aero-engine rotor system.A proper structural design will favor the static and dynamic performances of the system.In view of the deficiency of the current shape optimization techniques,a new mapping approach is proposed to define shape design variables based on the parametric equations of 3D curves and surfaces.It is then applied for the slot shape optimization of a squirrel-cage elastic support.To this end,an automatic design procedure that integrates the Genetic Algorithm (GA) is developed to solve the problem.Two typical examples with different shape constraints are considered.Numerical results provide reasonable optimum designs for the improvement of stiffness and strength of the squirrel-cage elastic support.
Elasticity of 3D networks with rigid filaments and compliant crosslinks
Heidemann, Knut M; Rehfeldt, Florian; Schmidt, Christoph F; Wardetzky, Max
2014-01-01
Disordered filamentous networks with compliant crosslinks exhibit a low linear elastic shear modulus at small strains, but stiffen dramatically at high strains. Experiments have shown that the elastic modulus can increase by up to three orders of magnitude while the networks withstand relatively large stresses without rupturing. Here, we perform an analytical and numerical study on model networks in three dimensions. Our model consists of a collection of randomly oriented rigid filaments connected by flexible crosslinks that are modeled as wormlike chains. Due to zero probability of filament intersection in three dimensions, our model networks are by construction prestressed in terms of initial tension in the crosslinks. We demonstrate how the linear elastic modulus can be related to the prestress in these network. Under the assumption of affine deformations in the limit of infinite crosslink density, we show analytically that the nonlinear elastic regime in 1- and 2-dimensional networks is characterized by p...
Scheunert, M.; Ullmann, A.; Afanasjew, M.; Börner, R.-U.; Siemon, B.; Spitzer, K.
2016-06-01
We present an inversion concept for helicopter-borne frequency-domain electromagnetic (HEM) data capable of reconstructing 3-D conductivity structures in the subsurface. Standard interpretation procedures often involve laterally constrained stitched 1-D inversion techniques to create pseudo-3-D models that are largely representative for smoothly varying conductivity distributions in the subsurface. Pronounced lateral conductivity changes may, however, produce significant artifacts that can lead to serious misinterpretation. Still, 3-D inversions of entire survey data sets are numerically very expensive. Our approach is therefore based on a cut-&-paste strategy whereupon the full 3-D inversion needs to be applied only to those parts of the survey where the 1-D inversion actually fails. The introduced 3-D Gauss-Newton inversion scheme exploits information given by a state-of-the-art (laterally constrained) 1-D inversion. For a typical HEM measurement, an explicit representation of the Jacobian matrix is inevitable which is caused by the unique transmitter-receiver relation. We introduce tensor quantities which facilitate the matrix assembly of the forward operator as well as the efficient calculation of the Jacobian. The finite difference forward operator incorporates the displacement currents because they may seriously affect the electromagnetic response at frequencies above 100. Finally, we deliver the proof of concept for the inversion using a synthetic data set with a noise level of up to 5%.
Shiann-Jong Lee
2010-01-01
Full Text Available Moment tensor inversion is a routine procedure to obtain information on an earthquake source for moment magnitude and focal mechanism. However, the inversion quality is usually controlled by factors such as knowledge of an earthquake location and the suitability of a 1-D velocity model used. Here we present an improved method to invert the moment tensor solution for local earthquakes. The proposed method differs from routine centroid-moment-tensor inversion of the Broadband Array in Taiwan for Seismology in three aspects. First, the inversion is repeated in the neighborhood of an _ hypocenter on a grid basis. Second, it utilizes _ functions based on a true three-dimensional velocity model. And third, it incorporates most of the input waveforms from strong-motion records. The proposed grid-based moment tensor inversion is applied to a local earthquake that occurred near the Taipei basin on 23 October 2004 to demonstrate its effectiveness and superiority over methods used in previous studies. By using the grid-based moment tensor inversion technique and 3-D _ functions, the earthquake source parameters, including earthquake location, moment magnitude and focal mechanism, are accurately found that are sufficiently consistent with regional ground motion observations up to a frequency of 1.0 Hz. This approach can obtain more precise source parameters for other earthquakes in or near a well-modeled basin and crustal structure.
Levashov, V. A.
2016-03-01
It is possible to associate with every atom or molecule in a liquid its own atomic stress tensor. These atomic stress tensors can be used to describe liquids' structures and to investigate the connection between structural and dynamic properties. In particular, atomic stresses allow to address atomic scale correlations relevant to the Green-Kubo expression for viscosity. Previously correlations between the atomic stresses of different atoms were studied using the Cartesian representation of the stress tensors or the representation based on spherical harmonics. In this paper we address structural correlations in a 3D model binary liquid using the eigenvalues and eigenvectors of the atomic stress tensors. This approach allows to interpret correlations relevant to the Green-Kubo expression for viscosity in a simple geometric way. On decrease of temperature the changes in the relevant stress correlation function between different atoms are significantly more pronounced than the changes in the pair density function. We demonstrate that this behaviour originates from the orientational correlations between the eigenvectors of the atomic stress tensors. We also found correlations between the eigenvalues of the same atomic stress tensor. For the studied system, with purely repulsive interactions between the particles, the eigenvalues of every atomic stress tensor are positive and they can be ordered: λ1 ≥ λ2 ≥ λ3 ≥ 0. We found that, for the particles of a given type, the probability distributions of the ratios (λ2/λ1) and (λ3/λ2) are essentially identical to each other in the liquids state. We also found that λ2 tends to be equal to the geometric average of λ1 and λ3. In our view, correlations between the eigenvalues may represent "the Poisson ratio effect" at the atomic scale.
Oh, Ju-Won
2016-07-04
Multiparameter full waveform inversion (FWI) applied to an elastic orthorhombic model description of the subsurface requires in theory a nine-parameter representation of each pixel of the model. Even with optimal acquisition on the Earth surface that includes large offsets, full azimuth, and multicomponent sensors, the potential for trade-off between the elastic orthorhombic parameters are large. The first step to understanding such trade-off is analysing the scattering potential of each parameter, and specifically, its scattering radiation patterns. We investigate such radiation patterns for diffraction and for scattering from a horizontal reflector considering a background isotropic model. The radiation patterns show considerable potential for trade-off between the parameters and the potentially limited resolution in their recovery. The radiation patterns of C11, C22, and C33 are well separated so that we expect to recover these parameters with limited trade-offs. However, the resolution of their recovery represented by recovered range of model wavenumbers varies between these parameters. We can only invert for the short wavelength components (reflection) of C33 while we can mainly invert for the long wavelength components (transmission) of the elastic coefficients C11 and C22 if we have large enough offsets. The elastic coefficients C13, C23, and C12 suffer from strong trade-offs with C55, C44, and C66, respectively. The trade-offs between C13 and C55, as well as C23 and C44, can be partially mitigated if we acquire P–SV and SV–SV waves. However, to reduce the trade-offs between C12 and C66, we require credible SH–SH waves. The analytical radiation patterns of the elastic constants are supported by numerical gradients of these parameters.
Chu, Chunlei
2009-01-01
We analyze the dispersion properties and stability conditions of the high‐order convolutional finite difference operators and compare them with the conventional finite difference schemes. We observe that the convolutional finite difference method has better dispersion properties and becomes more efficient than the conventional finite difference method with the increasing order of accuracy. This makes the high‐order convolutional operator a good choice for anisotropic elastic wave simulations on rotated staggered grids since its enhanced dispersion properties can help to suppress the numerical dispersion error that is inherent in the rotated staggered grid structure and its efficiency can help us tackle 3D problems cost‐effectively.
Force sensing using 3D displacement measurements in linear elastic bodies
Feng, Xinzeng; Hui, Chung-Yuen
2016-07-01
In cell traction microscopy, the mechanical forces exerted by a cell on its environment is usually determined from experimentally measured displacement by solving an inverse problem in elasticity. In this paper, an innovative numerical method is proposed which finds the "optimal" traction to the inverse problem. When sufficient regularization is applied, we demonstrate that the proposed method significantly improves the widely used approach using Green's functions. Motivated by real cell experiments, the equilibrium condition of a slowly migrating cell is imposed as a set of equality constraints on the unknown traction. Our validation benchmarks demonstrate that the numeric solution to the constrained inverse problem well recovers the actual traction when the optimal regularization parameter is used. The proposed method can thus be applied to study general force sensing problems, which utilize displacement measurements to sense inaccessible forces in linear elastic bodies with a priori constraints.
A 3D Orthotropic Strain-Rate Dependent Elastic Damage Material Model.
English, Shawn Allen
2014-09-01
A three dimensional orthotropic elastic constitutive model with continuum damage and cohesive based fracture is implemented for a general polymer matrix composite lamina. The formulation assumes the possibility of distributed (continuum) damage followed b y localized damage. The current damage activation functions are simply partially interactive quadratic strain criteria . However, the code structure allows for changes in the functions without extraordinary effort. The material model formulation, implementation, characterization and use cases are presented.
3D printed elastic honeycombs with graded density for tailorable energy absorption
Bates, Simon R. G.; Farrow, Ian R.; Trask, Richard S.
2016-04-01
This work describes the development and experimental analysis of hyperelastic honeycombs with graded densities, for the purpose of energy absorption. Hexagonal arrays are manufactured from thermoplastic polyurethane (TPU) via fused filament fabrication (FFF) 3D printing and the density graded by varying cell wall thickness though the structures. Manufactured samples are subject to static compression tests and their energy absorbing potential analysed via the formation of energy absorption diagrams. It is shown that by grading the density through the structure, the energy absorption profile of these structures can be manipulated such that a wide range of compression energies can be efficiently absorbed.
Abbott, D
2000-01-01
Tensor polarization observables (t20, t21 and t22) have been measured in elastic electron-deuteron scattering for six values of momentum transfer between 0.66 and 1.7 (GeV/c)^2. The experiment was performed at the Jefferson Laboratory in Hall C using the electron HMS Spectrometer, a specially designed deuteron magnetic channel and the recoil deuteron polarimeter POLDER. The new data determine to much larger Q^2 the deuteron charge form factors G_C and G_Q. They are in good agreement with relativistic calculations and disagree with pQCD predictions.
Garçon, M; Ahmidouch, A; Anklin, H; Arvieux, J; Ball, J; Beedoe, S; Beise, E J; Bimbot, L; Böglin, W; Breuer, H; Carlini, R; Chant, N S; Danagulyan, S; Dow, K; Ducret, J E; Dunne, J; Ewell, L A; Eyraud, L; Furget, C; Gilman, R; Glashausser, C; Gueye, P; Gustafsson, K K; Hafidi, K; Honegger, A; Jourdan, J; Kox, S; Kumbartzki, G; Lü, L; Lung, A; Mack, D; Markowitz, P; McIntyre, J; Meekins, D; Merchez, F; Mitchell, J; Möhring, R H; Mtingwa, S; Mrktchyan, H; Pitz, D; Qin, L; Ransome, R; Real, J S; Roos, P G; Rutt, P; Schmidt, W; Sawafta, R; Stepanyan, S; Stephenson, E J; Tieulent, R; Tomasi-Gustafsson, E; Turchinetz, W E; Vansyoc, K; Volmer, J; Voutier, E; Vulcan, W; Williamson, C; Wood, S A; Yan, C; Zhao, J; Zhao, W
1999-01-01
In elastic electron-deuteron scattering, the tensor polarization moments t sub 2 sub 0 , t sub 2 sub 1 and t sub 2 sub 2 , together with the unpolarized cross-sections, have been measured up to a momentum transfer of 1.8 (GeV/c) sup 2 , or 6.8 fm sup - sup 1. The experiment was performed at Jefferson Laboratory using the recoil deuteron polarimeter POLDER. Preliminary results are presented and discussed, especially in view of their significance concerning the applicability of perturbative QCD to this exclusive process.
Energy dependence of the tensor polarisation t20 in π+-d elastic scattering
The tensor polarisation t20 of the recoil deuterons from π+-d elastic scattering has been measured between Tsub(π)=110 and 150 MeV. An excitation curve of t20 shows a clear peak near the mass 2.14 GeV with a structure 20 MeV wide. The angular distributions at 117, 134 and 151 MeV confirm the rapidly changing behaviour with energy. A phase-shift analysis of all π+-d scattering data in this energy range has been performed. The strong energy and angular dependence is believed to manifest evidence for non-conventional dynamics. (author)
SHELL-5, Elastic Stress Analysis of 3-D Thin Shells Using Finite Elements Method
1 - Description of problem or function: SHELL5 performs an elastic stress analysis of smoothly curved, arbitrarily shaped, three- dimensional thin shells with any desired distributions of material properties, boundary constraints, and mechanical, thermal, and displacement loading conditions. 2 - Method of solution: SHELL5 uses the finite-element method to formulate a triangular plate element whose membrane displacement fields are linear polynomial functions and bending displacement field is a cubic polynomial function. The shell surface is approximated by a network of these plate elements of arbitrary orientation. Five degrees of freedom (3 displacements and 2 bending rotations) are obtained at each nodal point. The direct solution of the resulting system of equilibrium equations is obtained by the segment block tridiagonal Gaussian elimination method. 3 - Restrictions on the complexity of the problem - Maxima of: 1000 nodal points, 2000 elements, The maximum difference allowed for coupled nodal point indexes is 23
Pathak, Himanshu; Singh, Akhilendra; Singh, Indra Vir
2016-06-01
In this work, finite element method (FEM) and element free Galerkin method (EFGM) are coupled for solving 3D crack domains subjected to cyclic thermal load of constant amplitude. Crack growth contours and fatigue life have been obtained for each of the considered numerical examples. Thermo-elastic problems are decoupled into thermal and elastic problems . Firstly, the unknown temperature field is obtained by solving heat conduction equation, then, it is used as the input load in the elastic problem to calculate the displacement and stress fields. The geometrical discontinuity across crack surface is modelled by extrinsically enriched EFGM and the remaining part of the domain is approximated by standard finite element method. At the crack interface, a ramp function based interpolation scheme has been implemented. This coupled approach combines the advantages of both EFGM and FEM. A linear successive crack increment approach is used to model crack growth. The growing crack surface is traced by level set function. Standard Paris law is used for life estimation of the three-dimensional crack models. Different cases of planar and non-planar crack problems have been solved and their results are compared with the results obtained using extended finite element method to check accuracy, efficiency and robustness of the coupled FE-EFG approach implemented in this study.
Lyle E. Levine
2015-11-01
Full Text Available Nondestructive measurements of the full elastic strain and stress tensors from individual dislocation cells distributed along the full extent of a 50 µm-long polycrystalline copper via in Si is reported. Determining all of the components of these tensors from sub-micrometre regions within deformed metals presents considerable challenges. The primary issues are ensuring that different diffraction peaks originate from the same sample volume and that accurate determination is made of the peak positions from plastically deformed samples. For these measurements, three widely separated reflections were examined from selected, individual grains along the via. The lattice spacings and peak positions were measured for multiple dislocation cell interiors within each grain and the cell-interior peaks were sorted out using the measured included angles. A comprehensive uncertainty analysis using a Monte Carlo uncertainty algorithm provided uncertainties for the elastic strain tensor and stress tensor components.
Elastic and inelastic neutron scattering studies on 3d and 4f magnetic compounds
First, some theoretical aspects of neutron scattering techniques are given, and the cyrogenic equipment and the neutron spectrometers employed are described. Experiments on a 3-d Ising system are described, performed at very low temperatures and in a magnetic field. Experimental proof has been obtained for the theoretical prediction that the critical behaviour of a d-dimensional Ising system in a transverse magnetic field near T=0 is identical to that of a d+1 dimensional Ising system as a function of temperature in zero field. Experiments are described on a Ni2+ compound which represents a good example of a 1-d antiferromagnetic Heisenberg (HAF), spin s=1, system. The results give evidence for the so called 'Haldane conjecture', a theory which predicts that the ground state of HAF systems with integer spin is a nonmagnetic many-body singlet. The excited states are separated from the ground state by an energy gap. Contrastingly, half-integer spin systems are predicted to have no such gap. A short introduction is given to phenomena in rare earth, 4f compunds, like the Kondo effect and heavy fermion behaviour. Experimental results on the RE hexaborides are reported, among which CeB6, a typical Kondo system with complex magnetic orderings. Furthermore, inelastic neutron scattering experiments on NdB6 and CeB6, performed in order to get insight in the various reaction mechanisms, are presented. Finally a report is given on magnetic correlations and excitations in two nonmagnetically ordered heavy fermion compounds, CeCu6 and CeRu2Si2 and their interpretation in the light of existing theories. 201 refs.; 61 figs.; 4 tabs
Assi, Hisham
2016-01-01
Numerical simulation of wave propagation in an infinite medium is made possible by surrounding a finite region by a perfectly matched layer (PML). Using this approach a generalized three-dimensional (3D) formulation is proposed for time-domain modeling of elastic wave propagation in an unbounded lossless anisotropic medium. The formulation is based on a second-order approach that has the advantages of, physical relationship to the underlying equations, and amenability to be implemented in common numerical schemes. Specifically, our formulation uses three second-order equations of the displacement field and nine auxiliary equations, along with the three time histories of the displacement field. The properties of the PML, which are controlled by a complex two-parameter stretch function, are such that it acts as near perfect absorber. Using finite element method (FEM) 3D numerical results are presented for a highly anisotropic medium. An extension of the formulation to the particular case of a Kelvin-Vogit visco...
Wang, S.; De Hoop, M. V.; Xia, J.; Li, X.
2011-12-01
We consider the modeling of elastic seismic wave propagation on a rectangular domain via the discretization and solution of the inhomogeneous coupled Helmholtz equation in 3D, by exploiting a parallel multifrontal sparse direct solver equipped with Hierarchically Semi-Separable (HSS) structure to reduce the computational complexity and storage. In particular, we are concerned with solving this equation on a large domain, for a large number of different forcing terms in the context of seismic problems in general, and modeling in particular. We resort to a parsimonious mixed grid finite differences scheme for discretizing the Helmholtz operator and Perfect Matched Layer boundaries, resulting in a non-Hermitian matrix. We make use of a nested dissection based domain decomposition, and introduce an approximate direct solver by developing a parallel HSS matrix compression, factorization, and solution approach. We cast our massive parallelization in the framework of the multifrontal method. The assembly tree is partitioned into local trees and a global tree. The local trees are eliminated independently in each processor, while the global tree is eliminated through massive communication. The solver for the inhomogeneous equation is a parallel hybrid between multifrontal and HSS structure. The computational complexity associated with the factorization is almost linear with the size of the Helmholtz matrix. Our numerical approach can be compared with the spectral element method in 3D seismic applications.
Revealing the z~2.5 Cosmic Web With 3D Lyman-Alpha Forest Tomography: A Deformation Tensor Approach
Lee, Khee-Gan; White, Martin
2016-01-01
Studies of cosmological objects should take into account their positions within the cosmic web of large-scale structure. Unfortunately, the cosmic web has only been extensively mapped at low-redshifts ($z1$, the required galaxy densities are inaccessible for the foreseeable future, but 3D reconstructions of Lyman-$\\alpha$ forest absorption in closely-separated background QSOs and star-forming galaxies already offer a detailed window into $z\\sim2-3$ large-scale structure. We quantify the utili...
Borisov, Dmitry; Singh, Satish C.; Fuji, Nobuaki
2015-09-01
Seismic full waveform inversion is an objective method to estimate elastic properties of the subsurface and is an important area of research, particularly in seismic exploration community. It is a data-fitting approach, where the difference between observed and synthetic data is minimized iteratively. Due to a very high computational cost, the practical implementation of waveform inversion has so far been restricted to a 2-D geometry with different levels of physics incorporated in it (e.g. elasticity/viscoelasticity) or to a 3-D geometry but using an acoustic approximation. However, the earth is three-dimensional, elastic and heterogeneous and therefore a full 3-D elastic inversion is required in order to obtain more accurate and valuable models of the subsurface. Despite the recent increase in computing power, the application of 3-D elastic full waveform inversion to real-scale problems remains quite challenging on the current computer architecture. Here, we present an efficient method to perform 3-D elastic full waveform inversion for time-lapse seismic data using a finite-difference injection method. In this method, the wavefield is computed in the whole model and is stored on a surface above a finite volume where the model is perturbed and localized inversion is performed. Comparison of the final results using the 3-D finite-difference injection method and conventional 3-D inversion performed within the whole volume shows that our new method provides significant reductions in computational time and memory requirements without any notable loss in accuracy. Our approach shows a big potential for efficient reservoir monitoring in real time-lapse experiments.
Martin, O. J. F.; Gay-Balmaz, P.
2002-01-01
We describe a library to compute various types of Green's tensor for three-dimensional electromagnetic scattering calculations. This library includes the retarded and non-retarded (quasi-static) Green's tensors for infinite homogeneous space and the non-retarded Green's tensor associated with a surface. Both standard and filtered Green's tensor can be computed. Filtered Green's tensor can be used to accurately investigate high permittivity scatterers with the coupled-dipole approximation. (C)...
Liu, Xiaolong; Eckert, Andreas; Connolly, Peter
2016-06-01
Buckle folds of sedimentary strata commonly feature a variety of different fracture sets. Some fracture sets including outer arc tensile fractures and inner arc shear fractures at the fold hinge zones are well understood by the extensional and compressional strain/stress pattern. However, other commonly observed fracture sets, including tensile fractures parallel to the fold axis, tensile fractures cutting through the limb, extensional faults at the fold hinge, and other shear fractures of various orientations in the fold limb, fail to be intuitively explained by the strain/stress regimes during the buckling process. To obtain a better understanding of the conditions for the initiation of the various fractures sets associated with single-layer cylindrical buckle folds, a 3D finite element modeling approach using a Maxwell visco-elastic rheology is utilized. The influences of three model parameters with significant influence on fracture initiation are considered: burial depth, viscosity, and permeability. It is concluded that these parameters are critical for the initiation of major fracture sets at the hinge zone with varying degrees. The numerical simulation results further show that the buckling process fails to explain most of the fracture sets occurring in the limb unless the process of erosional unloading as a post-fold phenomenon is considered. For fracture sets that only develop under unrealistic boundary conditions, the results demonstrate that their development is realistic for a perclinal fold geometry. In summary, a more thorough understanding of fractures sets associated with buckle folds is obtained based on the simulation of in-situ stress conditions during the structural development of buckle folds.
Full text: Due to its low yield, the mb 3.9 North Korean nuclear test on 9 October 2006 could not be verified undoubtedly using solely seismic methods. The low signal-to-noise ratio (SNR) limits conventional discrimination methods such as P/S amplitude-ratios, and thus we perform a moment tensor inversion (MTI) to estimate the explosive (isotropic) component of the event. Long-period recordings are mainly controlled by the seismic source mechanism, while at shorter periods variations of the crustal structure have a stronger influence on the waveforms. On the other hand, long periods have a lower SNR than shorter periods. To obtain a good SNR for source-sensitive periods, we use an intermediate pass-band (T = 10-50 s) and include a regional 3D velocity structure deduced from the Crust 2.0 model to calculate synthetic waveforms. We study the influence of the crustal model on the inversion result using regional recordings with source-receiver distances up to 3000 km and compare the results with a tectonic mb 4.0 event that occurred in a distance of 34 km on 16 April 2002. We observe a higher isotropic part (Miso) for the nuclear test compared to the tectonic event. However, the absolute value of Miso cannot be determined because a trade-off with the inversion damping exists. To study the influence of seismic noise on Miso, we add artificial noise to the recordings of the tectonic event. Compared with the MTI of the original data we also find an increased Miso, but with other dependencies on the damping compared to the nuclear explosion data. We show trade-off curves of damping versus model length (as a function of MW) that can be a helpful tool for a timely verification of nuclear tests. (author)
Gao, Yurui; Parvathaneni, Prasanna; Schilling, Kurt G.; Wang, Feng; Stepniewska, Iwona; Xu, Zhoubing; Choe, Ann S.; Ding, Zhaohua; Gore, John C.; Chen, Li min; Landman, Bennett A.; Anderson, Adam W.
2016-03-01
Modern magnetic resonance imaging (MRI) brain atlases are high quality 3-D volumes with specific structures labeled in the volume. Atlases are essential in providing a common space for interpretation of results across studies, for anatomical education, and providing quantitative image-based navigation. Extensive work has been devoted to atlas construction for humans, macaque, and several non-primate species (e.g., rat). One notable gap in the literature is the common squirrel monkey - for which the primary published atlases date from the 1960's. The common squirrel monkey has been used extensively as surrogate for humans in biomedical studies, given its anatomical neuro-system similarities and practical considerations. This work describes the continued development of a multi-modal MRI atlas for the common squirrel monkey, for which a structural imaging space and gray matter parcels have been previously constructed. This study adds white matter tracts to the atlas. The new atlas includes 49 white matter (WM) tracts, defined using diffusion tensor imaging (DTI) in three animals and combines these data to define the anatomical locations of these tracks in a standardized coordinate system compatible with previous development. An anatomist reviewed the resulting tracts and the inter-animal reproducibility (i.e., the Dice index of each WM parcel across animals in common space) was assessed. The Dice indices range from 0.05 to 0.80 due to differences of local registration quality and the variation of WM tract position across individuals. However, the combined WM labels from the 3 animals represent the general locations of WM parcels, adding basic connectivity information to the atlas.
Petrov, P.; Newman, G. A.
2010-12-01
-Fourier domain we had developed 3D code for full-wave field simulation in the elastic media which take into account nonlinearity introduced by free-surface effects. Our approach is based on the velocity-stress formulation. In the contrast to conventional formulation we defined the material properties such as density and Lame constants not at nodal points but within cells. This second order finite differences method formulated in the cell-based grid, generate numerical solutions compatible with analytical ones within the range errors determinate by dispersion analysis. Our simulator will be embedded in an inversion scheme for joint seismic- electromagnetic imaging. It also offers possibilities for preconditioning the seismic wave propagation problems in the frequency domain. References. Shin, C. & Cha, Y. (2009), Waveform inversion in the Laplace-Fourier domain, Geophys. J. Int. 177(3), 1067- 1079. Shin, C. & Cha, Y. H. (2008), Waveform inversion in the Laplace domain, Geophys. J. Int. 173(3), 922-931. Commer, M. & Newman, G. (2008), New advances in three-dimensional controlled-source electromagnetic inversion, Geophys. J. Int. 172(2), 513-535. Newman, G. A., Commer, M. & Carazzone, J. J. (2010), Imaging CSEM data in the presence of electrical anisotropy, Geophysics, in press.
De Silva, R. T.; Pasbakhsh, Pooria; Goh, K. L.;
2014-01-01
A real-structure based 3-D micromechanical computational model of poly (lactic acid) nanocomposites reinforced by randomly oriented halloysite nanotubes (HNTs) was developed and compared with an idealized model (conventional model) and experimental results. The developed idealized model consists of...
3D tensor tractography (DTT) has been applied to central nervous system (CNS) diseases to depict neuronal fibers. In this study with 3 tesla MRI, we have evaluated DTT to predict outcome of motor function in patients with lacunar infarcts. Fifteen patients with New lacunar infarcts, underwent DTTs with at least one in the acute (mean 1.4 days) and another in the subacute phase (mean 18.7 days). Patients were separated to 2 groups, recovery and non-recovery. Patients in former group had almost complete recovery in motor function 3 month later, while those to latter had a residual hemiparesis. Motor function was assessed with MMT score, which was uniquely stratefied into 12 levels by a modified MMT (manual muscle testing) protocol. DTT was implemented with 3 tesla MRI (Signa Excite; GE) and analyzed with dTV. IISR which was produced by the Department of Radiology, Tokyo University. The pyramidal tract was delineated by setting each region of interest (ROI), with the cerebral peduncle as the seed point and the motor cortex as the target point. The number of pyramidal fibers was identified as drawn lines obtained from a result display. The ratio of the number of fibers (RF) was calculated based on the number of fibers in the injured side relative to the number of fibers in the intact side x 100. In acute phase mean RFs the recovery (70.8±21.6%) and non-recovery (63.5±23.4%) groups were not significantly different. RF of recovery group in subacute phase was 100.5±28.3%, which was significantly higher with that in acute phase, meanwhile there was no significance difference between RFs of non-recovery group in two phases. In addition there was a significant correlation (R2=0.89) between MMT score 3 month later and RF in subacute phase in all patients group. There seems to be a correlation between long-term recovery of motor function and increased numbers in pyramidal fibers defected by DTT. Therefore, DTT may have a potential use in predicting the outcome of patients
Morita, Yasuyuki; Kawase, Naoki; Ju, Yang; Yamauchi, Takashi
2016-07-01
Cells maintain homeostasis and perform various functions by interacting mechanically with a cell-adhesive matrix. Regarding cellular differentiation, it has been found that matrix elasticity can determine the differentiation lineage of mesenchymal stem cells (MSCs). Direct quantitative measurements of the mechanical interaction between MSCs and matrix for differentiation, however, have yet to be reported. Herein, the displacement field of the cell-adhesive matrix was observed quantitatively using a digital volume correlation (DVC) method. Maximum displacement and cellular traction stress were analyzed when the MSC differentiated into a neuron-like cell or an osteoblast-like cell on a soft or hard elastic matrix, respectively. The function of non-muscle myosin II (NMM II), which plays an important role in intracellular cytoskeletal dynamics, was investigated during cellular differentiation. The mechanical interaction (maximum displacement and subjected area of the matrix) between the cell and matrix was dependent on matrix elasticity. It has also been shown that the mechanical interaction between the intracellular cytoskeleton and cell-adhesion matrix is indispensable for cellular differentiation. This work provides the first quantitative visualization of the mechanical interaction between MSCs and the cell-adhesion matrix for differentiation. PMID:26945874
Levashov, Valentin A.
2015-01-01
It is possible to associate with every atom or molecule in a liquid its own atomic stress tensor. These atomic stress tensors can be used to describe liquids' structures and to investigate the connection between structural and dynamic properties. In particular, atomic stresses allow to address atomic scale correlations relevant to the Green-Kubo expression for viscosity. Previously correlations between the atomic stresses of different atoms were studied using the Cartesian representation of t...
Obermann, Anne; Planès, Thomas; Hadziioannou, Céline; Campillo, Michel
2016-07-01
In the context of seismic monitoring, recent studies made successful use of seismic coda waves to locate medium changes on the horizontal plane. Locating the depth of the changes, however, remains a challenge. In this paper, we use 3-D wavefield simulations to address two problems: firstly, we evaluate the contribution of surface and body wave sensitivity to a change at depth. We introduce a thin layer with a perturbed velocity at different depths and measure the apparent relative velocity changes due to this layer at different times in the coda and for different degrees of heterogeneity of the model. We show that the depth sensitivity can be modelled as a linear combination of body- and surface-wave sensitivity. The lapse-time dependent sensitivity ratio of body waves and surface waves can be used to build 3-D sensitivity kernels for imaging purposes. Secondly, we compare the lapse-time behavior in the presence of a perturbation in horizontal and vertical slabs to address, for instance, the origin of the velocity changes detected after large earthquakes.
This technical note demonstrates the relevance of the isotropic 3D T2 turbo-spin-echo (TSE) sequence with short-term inversion recovery (STIR) and variable flip angle RF excitations (SPACE: Sampling Perfection with Application optimized Contrasts using different flip angle Evolutions) for high-resolution brachial plexus imaging. The sequence was used in 11 patients in the diagnosis of brachial plexus pathologies involving primary and secondary tumors, and in six volunteers. We show that 3D STIR imaging is not only a reliable alternative to 2D STIR imaging, but it also better evaluates the anatomy, nerve site compression and pathology of the plexus, especially to depict space-occupying tumors along its course. Finally, due to its appropriate contrast we describe how 3D-STIR can be used as a high-resolution mask to be fused with fraction of anisotropy (FA) maps calculated from diffusion tensor imaging (DTI) data of the plexus. (orig.)
Viallon, M.; Vargas, M.I.; Jlassi, H.; Loevblad, K.O.; Delavelle, J. [University Hospital of Geneva, Department of Radiology, Geneva (Switzerland)
2008-05-15
This technical note demonstrates the relevance of the isotropic 3D T2 turbo-spin-echo (TSE) sequence with short-term inversion recovery (STIR) and variable flip angle RF excitations (SPACE: Sampling Perfection with Application optimized Contrasts using different flip angle Evolutions) for high-resolution brachial plexus imaging. The sequence was used in 11 patients in the diagnosis of brachial plexus pathologies involving primary and secondary tumors, and in six volunteers. We show that 3D STIR imaging is not only a reliable alternative to 2D STIR imaging, but it also better evaluates the anatomy, nerve site compression and pathology of the plexus, especially to depict space-occupying tumors along its course. Finally, due to its appropriate contrast we describe how 3D-STIR can be used as a high-resolution mask to be fused with fraction of anisotropy (FA) maps calculated from diffusion tensor imaging (DTI) data of the plexus. (orig.)
Bending analysis of a general cross-ply laminate using 3D elasticity solution and layerwise theory
Yazdani Sarvestani, H.; Naghashpour, A.; Heidari-Rarani, M.
2015-12-01
In this study, the analytical solution of interlaminar stresses near the free edges of a general (symmetric and unsymmetric layups) cross-ply composite laminate subjected to pure bending loading is presented based on Reddy's layerwise theory (LWT) for the first time. First, the reduced form of displacement field is obtained for a general cross-ply composite laminate subjected to a bending moment by elasticity theory. Then, first-order shear deformation theory of plates and LWT is utilized to determine the global and local deformation parameters appearing in the displacement fields, respectively. One of the main advantages of the developed solution based on the LWT is exact prediction of interlaminar stresses at the boundary layer regions. To show the accuracy of this solution, three-dimensional elasticity bending problem of a laminated composite is solved for special set of boundary conditions as well. Finally, LWT results are presented for edge-effect problems of several symmetric and unsymmetric cross-ply laminates under the bending moment. The obtained results indicate high stress gradients of interlaminar stresses near the edges of laminates.
The source mechanisms of microseismics in hydraulic fracturing present guiding significance to the research on source types, crustal stress analysis and crack prediction. Numerical simulations based on various source mechanisms can be used to investigate the stress characteristics and response characteristics of different source types. In this paper, a method based on the seismic moment tensor (SMT) and elastic wave equation (EWE) was presented for forward modeling. Additionally, we have given the expressions of nine couples of force which can be combined into different kinds of source types. The calculations of wave fields and records with three basic types of sources showed the features in homogeneous isotropic and anisotropic media by the finite-difference (FD) method. Lastly, analysis of the relationship between the polarizing angle and incident angle provided us with some evidence to distinguish the type of media in single media. The work offers methods of instruction for identification and interpretation in microseismic monitoring. (paper)
Gras, Renaud
2015-03-01
Performing a single but complex mechanical test on small structures rather than on coupons to probe multiple strain states/histories for identification purposes is nowadays possible thanks to full-field measurements. The aim is to identify many parameters thanks to the heterogeneity of mechanical fields. Such an approach is followed herein, focusing on a blade root made of 3D woven composite. The performed test, which is analyzed using global Digital Image Correlation (DIC), provides heterogeneous kinematic fields due to the particular shape of the sample. This displacement field is further processed to identify the four in-plane material parameters of the macroscopic equivalent orthotropic behavior. The key point, which may limit the ability to draw reliable conclusions, is the presence of acquisition noise in the original images that has to be tracked along the DIC/identification processing to provide uncertainties on the identified parameters. A further regularization based on a priori knowledge is finally introduced to compensate for possible lack of experimental information needed for completing the identification.
Ha, Jiho; Shin, Sungryul; Shin, Changsoo; Chung, Wookeen
2015-05-01
Because complex mixed waves are typically generated in elastic media, wavefield decomposition is required for such media to obtain migration images accurately. In isotropic media, this is achieved according to the Helmholtz decomposition theorem; in particular, the divergence operator is commonly applied to P-wavefield decomposition. In this study, two types of elastic reverse-time migration algorithms are proposed for decomposition of the P-wavefield without requiring the divergence operator. The first algorithm involves formulation of the stress tensor by spatially differentiated displacement according to the stress-strain relationship and is utilized to construct an imaging condition for the decomposed P-wavefield. We demonstrate this approach through numerical testing. The second algorithm allows us to obtain emphasized interfaces through the application of the absolute value function to decomposed wavefield in imaging condition. Because reverse-time migration can be defined by a zero-lag cross-correlation relationship between the partial-derivative wavefield and the observed wavefield data, we derive the virtual source to construct the partial-derivative wavefield based on a 2D staggered-grid finite-difference modeling method in the time domain. The explicitly computed partial-derivative wavefield from virtual sources with the stress tensor is in agreement with the partial-derivative wavefield directly computed from residual by between with and without a perturbation point in the subsurface. Moreover, the back-propagation technique is used to enhance the computational efficiency. To validate our two types of imaging conditions, numerical tests are conducted. The migration images created according to our imaging conditions can represent the subsurface structure accurately. Thus, we can confirm the feasibility of obtaining migration images of the decomposed P-wavefield without requiring the application of the divergence operator.
Ju Wei; Sun Weifeng; Ma Xiaojing; Jiang Hui
2016-07-01
Future earthquake potential in the Bohai–Zhangjiakou Seismotectonic Zone (BZSZ) in North Chinadeserves close attention. Tectonic stress accumulation state is an important indicator for earthquakes;therefore, this study aims to analyse the stress accumulation state in the BZSZ via three-dimensionalvisco-elastic numerical modelling. The results reveal that the maximum shear stress in the BZSZ increasesgradually as the depth increases, and the stress range is wider in the lower layer. In the upper layer, themaximum shear stress is high in the Zhangjiakou area, whereas in the lower layer, relatively high valuesoccur in the Penglai–Yantai area, which may be affected by the depth of the Moho surface. Besides,weak fault zones will be easily fractured when the maximum shear stress is not sufficiently high due totheir low strengths, resulting in earthquakes. Therefore, based on the modelling results, the upper layerof the Zhangjiakou area and the lower layer of the Penglai–Yantai area in the BZSZ in North China aremore likely to experience earthquakes.
Bignardi, S.; Mantovani, A.; Abu Zeid, N.
2016-08-01
OpenHVSR is a computer program developed in the Matlab environment, designed for the simultaneous modeling and inversion of large Horizontal-to-Vertical Spectral Ratio (HVSR or H/V) datasets in order to construct 2D/3D subsurface models (topography included). The program is designed to provide a high level of interactive experience to the user and still to be of intuitive use. It implements several effective and established tools already present in the code ModelHVSR by Herak (2008), and many novel features such as: -confidence evaluation on lateral heterogeneity -evaluation of frequency dependent single parameter impact on the misfit function -relaxation of Vp/Vs bounds to allow for water table inclusion -a new cost function formulation which include a slope dependent term for fast matching of peaks, which greatly enhances convergence in case of low quality HVSR curves inversion -capability for the user of editing the subsurface model at any time during the inversion and capability to test the changes before acceptance. In what follows, we shall present many features of the program and we shall show its capabilities on both simulated and real data. We aim to supply a powerful tool to the scientific and professional community capable of handling large sets of HSVR curves, to retrieve the most from their microtremor data within a reduced amount of time and allowing the experienced scientist the necessary flexibility to integrate into the model their own geological knowledge of the sites under investigation. This is especially desirable now that microtremor testing has become routinely used. After testing the code over different datasets, both simulated and real, we finally decided to make it available in an open source format. The program is available by contacting the authors.
Ichimura, Tsuyoshi; Agata, Ryoichiro; Hori, Takane; Hirahara, Kazuro; Hashimoto, Chihiro; Hori, Muneo; Fukahata, Yukitoshi
2016-07-01
As a result of the accumulation of high-resolution observation data, 3-D high-fidelity crustal structure data for large domains are becoming available. However, it has been difficult to use such data to perform elastic/viscoelastic crustal deformation analyses in large domains with quality assurance of the numerical simulation that guarantees convergence of the numerical solution with respect to the discretization size because the costs of analysis are significantly high. This paper proposes a method of constructing a high-fidelity crustal structure finite element (FE) model using high-fidelity crustal structure data and fast FE analysis to reduce the costs of analysis (based on automatic FE model generation for parallel computation, OpenMP/MPI hybrid parallel computation on distributed memory computers, a geometric multigrid, variable preconditioning and multiple precision arithmetic). Using the proposed methods, we construct 10 billion degree-of-freedom high-fidelity crustal structure FE models for the entire Japan, and conduct elastic/viscoelastic crustal deformation analysis using this model with enough high accuracy of the numerical simulation.
Jechumtálová, Zuzana; Bulant, P.
Prague : Charles University, 2012, s. 17-27 R&D Projects: GA ČR GAP210/10/1728 Institutional support: RVO:67985530 Keywords : ray tracing * 3-D velocity model * earthquake mechanism * amplitude inversion Subject RIV: DC - Siesmology, Volcanology, Earth Structure
Jechumtálová, Zuzana; Bulant, P.
Prague : Charles University, 2013, s. 253-269 R&D Projects: GA ČR GAP210/10/1728 Institutional support: RVO:67985530 Keywords : ray tracing * 3-D velocity model * earthquake mechanism * amplitude inversion Subject RIV: DC - Siesmology, Volcanology, Earth Structure
Jechumtálová, Zuzana; Bulant, P.
2014-01-01
Roč. 18, č. 3 (2014), s. 511-531. ISSN 1383-4649 R&D Projects: GA ČR GAP210/10/1728 EU Projects: European Commission(XE) 230669 - AIM Institutional support: RVO:67985530 Keywords : ray tracing * 1-D and 3-D velocity models * earthquake mechanism * amplitude inversion * Dobrá Voda earthquakes Subject RIV: DC - Siesmology, Volcanology, Earth Structure Impact factor: 1.386, year: 2014
We present a newly developed combined scanning particle tracking velocimetry (SPTV) and scanning laser-induced fluorescence (SLIF) technique. The new method allows for the first time to measure the full velocity gradient tensor and the three-dimensional density field simultaneously in a refractive index matched environment. The data thus obtained will be valuable in investigating the interaction between turbulence, i.e. the dynamics of vorticity and strain, and the density field. In this study, we describe the implementation of the measurement system in detail. As a showcase, the new technique is applied to a gravity current flow and the measured data are validated by imposing various checks. Further results are presented that illustrate the capability of the SPTV/SLIF technique in the investigation of interface dynamics. (paper)
3D Face Recognition Using Facial Curves Elastic Matching%基于面部曲线弹性匹配的三维人脸识别方法
潘仁林; 达飞鹏; 邹红艳; 王朝阳
2014-01-01
提出一种基于面部径向曲线弹性匹配的三维人脸识别方法。使用人脸曲面上的多条曲线表征人脸曲面，提取三维人脸上从鼻尖点发射的多条面部径向曲线，对其进行分层弹性匹配和点距对应匹配，根据人脸不同部位受表情影响程度不同，对不同曲线识别相似度赋予不同权重进行加权融合作为总相似度用于识别。测试结果表明该方法具有很好的识别性能，并且对表情、遮挡和噪声具有较好的鲁棒性。%A 3D face recognition algorithm based on facial curves elastic matching is presented in this paper. Facial surfaces are approximated with indexed sets of level curves. Firstly, the radial curves emanating from the nose tip are extracted after 3D faces are preprocessed. Secondly, the curves of the probe face match the curves of gallery face using the improved hierarchical matching method and the two curves are matched exploiting the distance between the point and nose tip. Finally, the similarities of different curves are weighted fused according to the different degree influenced by expression as sum similarity which is used for matching. Experiment results verify that the proposed method can obtain better recognition rate, and is robust to occlusion, spikes and expression.
Zhou, Li; Zhang, Wei; Shen, Yang; Chen, Xiaofei; Zhang, Jie
2016-06-01
With dense seismic arrays and advanced imaging methods, regional three-dimensional (3D) Earth models have become more accurate. It is now increasingly feasible and advantageous to use a 3D Earth model to better locate earthquakes and invert their source mechanisms by fitting synthetics to observed waveforms. In this study, we develop an approach to determine both the earthquake location and source mechanism from waveform information. The observed waveforms are filtered in different frequency bands and separated into windows for the individual phases. Instead of picking the arrival times, the traveltime differences are measured by cross-correlation between synthetic waveforms based on the 3D Earth model and observed waveforms. The earthquake location is determined by minimizing the cross-correlation traveltime differences. We then fix the horizontal location of the earthquake and perform a grid search in depth to determine the source mechanism at each point by fitting the synthetic and observed waveforms. This new method is verified by a synthetic test with noise added to the synthetic waveforms and a realistic station distribution. We apply this method to a series of M W3.4-5.6 earthquakes in the Longmenshan fault (LMSF) zone, a region with rugged topography between the eastern margin of the Tibetan plateau and the western part of the Sichuan basin. The results show that our solutions result in improved waveform fits compared to the source parameters from the catalogs we used and the location can be better constrained than the amplitude-only approach. Furthermore, the source solutions with realistic topography provide a better fit to the observed waveforms than those without the topography, indicating the need to take the topography into account in regions with rugged topography.
Applications of tensor analysis
McConnell, A J
2011-01-01
Standard work applies tensorial methods to subjects within realm of advanced college mathematics. Text explains fundamental ideas and notation of tensor theory; covers geometrical treatment of tensor algebra; introduces theory of differentiation of tensors; and applies mathematics to dynamics, electricity, elasticity and hydrodynamics. 685 exercises, most with answers.
Erba, Alessandro
2016-05-18
Symmetry features of the internal-strain tensor of crystals (whose components are mixed second-energy derivatives with respect to atomic displacements and lattice strains) are formally presented, which originate from translational-invariance, atomic equivalences, and atomic invariances. A general computational scheme is devised, and implemented into the public Crystal program, for the quantum-mechanical evaluation of the internal-strain tensor of crystals belonging to any space-group, which takes full-advantage of the exploitation of these symmetry-features. The gain in computing time due to the full symmetry exploitation is documented to be rather significant not just for high-symmetry crystalline systems such as cubic, hexagonal or trigonal, but also for low-symmetry ones such as monoclinic and orthorhombic. The internal-strain tensor is used for the evaluation of the nuclear relaxation term of the fourth-rank elastic and third-rank piezoelectric tensors of crystals, where, apart from a reduction of the computing time, the exploitation of symmetry is documented to remarkably increase the numerical precision of computed coefficients. PMID:27150599
Combinatorial 3D Mechanical Metamaterials
Coulais, Corentin; Teomy, Eial; de Reus, Koen; Shokef, Yair; van Hecke, Martin
2015-03-01
We present a class of elastic structures which exhibit 3D-folding motion. Our structures consist of cubic lattices of anisotropic unit cells that can be tiled in a complex combinatorial fashion. We design and 3d-print this complex ordered mechanism, in which we combine elastic hinges and defects to tailor the mechanics of the material. Finally, we use this large design space to encode smart functionalities such as surface patterning and multistability.
Deuteron tensor and vector analyzing powers have been measured for the first time in p-d vector elastic scattering at 800 MeV proton laboratory energy for a-t range between 0.032 and 1.038 (GeV/c)2. The data provide a very sensitive test of the multiple scattering theory at intermediate energies. The inclusion of the noneikonal correction in the Glauber diffraction theory is essential in the description of the data, especially for-t approx. >= 0.2(GeV/c)2, where discrepancies of up to 150% are observed with the theory based on the eikonal approximation. (orig.)
Soutas-Little, Robert William
2010-01-01
According to the author, elasticity may be viewed in many ways. For some, it is a dusty, classical subject . . . to others it is the paradise of mathematics."" But, he concludes, the subject of elasticity is really ""an entity itself,"" a unified subject deserving comprehensive treatment. He gives elasticity that full treatment in this valuable and instructive text. In his preface, Soutas-Little offers a brief survey of the development of the theory of elasticity, the major mathematical formulation of which was developed in the 19th century after the first concept was proposed by Robert Hooke
Ivankina, T. I.; Kern, H.; Lokajíček, Tomáš; Nikitin, A. N.
s. l. : European Seismological Comission, 2010. ---. [European Seismological Comission General Assembly /32./. 06.09.2010-10.09.2010, Montpellier] R&D Projects: GA AV ČR IAA300130906 Institutional research plan: CEZ:AV0Z30130516 Keywords : Fine-grained biotite gneiss * neutron diffraction * elastic anisotropy lattice preferred orientation * shape preferred orientation Subject RIV: DC - Siesmology, Volcanology, Earth Structure
Tensorial analysis of Eshelby stresses in 3D supercooled liquids
Lemaître, Anaël
2015-10-01
It was recently proposed that the local rearrangements governing relaxation in supercooled liquids impress on the liquid medium long-ranged (Eshelby) stress fluctuations that accumulate over time. From this viewpoint, events must be characterized by elastic dipoles, which are second order tensors, and Eshelby fields are expected to show up in stress and stress increment correlations, which are fourth order tensor fields. We construct here an analytical framework that permits analyzing such tensorial correlations in isotropic media in view of accessing Eshelby fields. Two spherical bases are introduced, which correspond to Cartesian and spherical coordinates for tensors. We show how they can be used to decompose stress correlations and thus test such properties as isotropy and power-law scalings. Eshelby fields and the predicted stress correlations in an infinite medium are shown to belong to an algebra that can conveniently be described using the spherical tensor bases. Using this formalism, we demonstrate that the inherent stress field of 3D supercooled liquids is power law correlated and carries the signature of Eshelby fields, thus supporting the idea that relaxation events give rise to Eshelby stresses that accumulate over time.
Tensorial analysis of Eshelby stresses in 3D supercooled liquids.
Lemaître, Anaël
2015-10-28
It was recently proposed that the local rearrangements governing relaxation in supercooled liquids impress on the liquid medium long-ranged (Eshelby) stress fluctuations that accumulate over time. From this viewpoint, events must be characterized by elastic dipoles, which are second order tensors, and Eshelby fields are expected to show up in stress and stress increment correlations, which are fourth order tensor fields. We construct here an analytical framework that permits analyzing such tensorial correlations in isotropic media in view of accessing Eshelby fields. Two spherical bases are introduced, which correspond to Cartesian and spherical coordinates for tensors. We show how they can be used to decompose stress correlations and thus test such properties as isotropy and power-law scalings. Eshelby fields and the predicted stress correlations in an infinite medium are shown to belong to an algebra that can conveniently be described using the spherical tensor bases. Using this formalism, we demonstrate that the inherent stress field of 3D supercooled liquids is power law correlated and carries the signature of Eshelby fields, thus supporting the idea that relaxation events give rise to Eshelby stresses that accumulate over time. PMID:26520535
Zhang, Mian; Huang, Cheng-li
2012-08-01
Generalized spherical harmonics (GSH) are usually applied on the problems where the Earth model is elliptical and elastic stress tensor is involved in, as stress tensor can’t be represented in vector spherical harmonics. However, the divergence of the te ns or and a vector dot - product with the tensor are only needed on computation rotation modes of the Earth which can be written in the vector spherical harmonics. We extend the equations on the spherical Earth to asymmetric 3D model by means of linear operator method. This method doesn’t use the complicated generalized spherical harmonics nor Wigner 3 - j symbol. As a validation of this method, the practical calculation of rotational modes of 3D Earth will be made and discussed.
Ichimura, Tsuyoshi; Agata, Ryoichiro; Hori, Takane; Hirahara, Kazuro; Hashimoto, Chihiro; Hori, Muneo; Fukahata, Yukitoshi
2016-04-01
As a result of the accumulation of high-resolution observation data, three-dimensional high-fidelity crustal structure data for large domains are becoming available. However, it has been difficult to use such data to perform elastic/viscoelastic crustal deformation analyses in large domains with quality assurance of the numerical simulation that guarantees convergence of the numerical solution with respect to the discretisation size, because the costs of analysis are significantly high. This paper proposes a method of constructing a high-fidelity crustal structure finite element (FE) model using high-fidelity crustal structure data and fast FE analysis to reduce the costs of analysis (based on automatic FE model generation for parallel computation, OpenMP/MPI hybrid parallel computation on distributed memory computers, a geometric multigrid, variable preconditioning, and multiple precision arithmetic). Using the proposed methods, we construct 10 billion degree-of-freedom high-fidelity crustal structure FE models for the entire Japan, and conduct elastic/viscoelastic crustal deformation analysis using this model with enough high accuracy of the numerical simulation.
Thizy, C.; Eliot, F.; Ballhause, D.; Olympio, K. R.; Kluge, R.; Shannon, A.; Laduree, G.; Logut, D.; Georges, M. P.
2013-04-01
Thermo-elastic distortions of composite structures have been measured by a holographic camera using a BSO photorefractive crystal as the recording medium. The first test campaign (Phase 1) was performed on CFRP struts with titanium end-fittings glued to the tips of the strut. The samples were placed in a vacuum chamber. The holographic camera was located outside the chamber and configured with two illuminations to measure the relative out-of-plane and in-plane (in one direction) displacements. The second test campaign (Phase 2) was performed on a structure composed of a large Silicon Carbide base plate supported by 3 GFRP struts with glued Titanium end-fittings. Thermo-elastic distortions have been measured with the same holographic camera used in phase 1, but four illuminations, instead of two, have been used to provide the three components of displacement. This technique was specially developed and validated during the phase 2 in CSL laboratory. The system has been designed to measure an object size of typically 250x250 mm2; the measurement range is such that the sum of the largest relative displacements in the three measurement directions is maximum 20 μm. The validation of the four-illuminations technique led to measurement uncertainties of 120 nm for the relative in-plane and out-of-plane displacements, 230 nm for the absolute in-plane displacement and 400 nm for the absolute out-of-plane displacement. For both campaigns, the test results have been compared to the predictions obtained by finite element analyses and the correlation of these results was good.
Ando, R.
2014-12-01
The boundary integral equation method formulated in the real space and time domain (BIEM-ST) has been used as a powerful tool to analyze the earthquake rupture dynamics on non-planar faults. Generally, BIEM is more accurate than volumetric methods such as the finite difference method and the finite difference method. With the recent development of the high performance computing environment, the earthquake rupture simulation studies have been conducted considering three dimensional realistic fault geometry models. However, the utility of BIEM-ST has been limited due to its heavy computational demanding increased depending on square of time steps (N2), which was needed to evaluate the historic integration. While BIEM can be efficient with the spectral domain formulation, the applications of such a method are limited to planar fault cases. In this study, we propose a new method to reduce the calculation time of BIEM-ST to linear of time step (N) without degrading the accuracy in the 3 dimensional modeling space. We extends the method proposed earlier for the case of the 2 dimensional framework, applying the asymptotic expressions of the elasto-dynamic Green's functions. This method uses the physical nature of the stress Green's function as dividing the causality cone according to the distances from the wave-fronts. The scalability of this method is shown on the parallel computing environment of the distributed memory. We demonstrate the applicability to analyses of subduction earthquake cases, suffering long time from the numerical limitations of previously available BIEMs. We analyze the dynamic rupture processes on dipping reverse faults embed in a three dimensional elastic half space.
DYNA3D, 3-D Finite Elements for Dynamic Response of Inelastic Solids
1 - Description of program or function: DYNA3D is an explicit, three- dimensional, finite element program for analyzing the large deformation dynamic response of inelastic solids and structures. DYNA3D contain 30 material models and 10 equations of state (EOS) to cover a wide range of material behavior. The material models implemented are: elastic, orthotropic elastic, kinematic/isotropic plasticity, thermo-elastoplastic, soil and crushable foam, linear viscoelastic, Blatz-Ko rubber, high explosive burn, hydrodynamic without deviatoric stresses, elastoplastic hydrodynamic, temperature dependent elastoplastic, isotropic elastoplastic, isotropic elastoplastic with failure, soil and crushable foam with failure, Johnson/Cook plasticity model, pseudo TENSOR geological model, elastoplastic with fracture, power law isotropic plasticity, strain rate dependent plasticity, rigid, thermal orthotropic, composite damage model, thermal orthotropic with 12 curves, piecewise linear isotropic plasticity, and inviscid two invariant geologic cap, orthotropic crushable model, Moonsy-Rivlin rubber, resultant plasticity, closed form update shell plasticity, and Frazer-Nash rubber model. The IBM 3090 version does not contain the last two models mentioned. The hydrodynamic material models determine only the deviatoric stresses. Pressure is determined by one of ten equations of state including linear polynomial, JWL high explosive, Sack 'Tuesday' high explosive, Gruneisen, ratio of polynomials, linear polynomial with energy deposition, ignition and growth of reaction in HE, tabulated compaction, tabulated, and TENSOR pore collapse. DYNA3D generates three binary output databases. One contains information for complete states at infrequent intervals; 50 to 100 states is typical. The second contains information for a subset of nodes and elements at frequent intervals; 1,000 to 10,000 states is typical. The last contains interfaces data for contact surfaces. 2 - Method of solution: A contact
Lucas, Laurent; Loscos, Céline
2013-01-01
While 3D vision has existed for many years, the use of 3D cameras and video-based modeling by the film industry has induced an explosion of interest for 3D acquisition technology, 3D content and 3D displays. As such, 3D video has become one of the new technology trends of this century.The chapters in this book cover a large spectrum of areas connected to 3D video, which are presented both theoretically and technologically, while taking into account both physiological and perceptual aspects. Stepping away from traditional 3D vision, the authors, all currently involved in these areas, provide th
Beane, Andy
2012-01-01
The essential fundamentals of 3D animation for aspiring 3D artists 3D is everywhere--video games, movie and television special effects, mobile devices, etc. Many aspiring artists and animators have grown up with 3D and computers, and naturally gravitate to this field as their area of interest. Bringing a blend of studio and classroom experience to offer you thorough coverage of the 3D animation industry, this must-have book shows you what it takes to create compelling and realistic 3D imagery. Serves as the first step to understanding the language of 3D and computer graphics (CG)Covers 3D anim
Modelling elasticity in solids using active cubes - application to simulated operations
Bro-Nielsen, Morten
The paper describes an approach to elastic modelling of human tissue based on the use of 3D solid active models-active cubes (M. Bro-Nielsen, 1994)-and a shape description based on the metric tensor in a solid. Active cubes are used because they provide a natural parameterization of the surface and....... The application of this approach to modelling the elastic deformation of human tissue in response to movement of bones is demonstrated...
D. Pletinckx
2012-09-01
Full Text Available The current 3D hype creates a lot of interest in 3D. People go to 3D movies, but are we ready to use 3D in our homes, in our offices, in our communication? Are we ready to deliver real 3D to a general public and use interactive 3D in a meaningful way to enjoy, learn, communicate? The CARARE project is realising this for the moment in the domain of monuments and archaeology, so that real 3D of archaeological sites and European monuments will be available to the general public by 2012. There are several aspects to this endeavour. First of all is the technical aspect of flawlessly delivering 3D content over all platforms and operating systems, without installing software. We have currently a working solution in PDF, but HTML5 will probably be the future. Secondly, there is still little knowledge on how to create 3D learning objects, 3D tourist information or 3D scholarly communication. We are still in a prototype phase when it comes to integrate 3D objects in physical or virtual museums. Nevertheless, Europeana has a tremendous potential as a multi-facetted virtual museum. Finally, 3D has a large potential to act as a hub of information, linking to related 2D imagery, texts, video, sound. We describe how to create such rich, explorable 3D objects that can be used intuitively by the generic Europeana user and what metadata is needed to support the semantic linking.
This book explains modeling of solid works 3D and application of 3D CAD/CAM. The contents of this book are outline of modeling such as CAD and 2D and 3D, solid works composition, method of sketch, writing measurement fixing, selecting projection, choosing condition of restriction, practice of sketch, making parts, reforming parts, modeling 3D, revising 3D modeling, using pattern function, modeling necessaries, assembling, floor plan, 3D modeling method, practice floor plans for industrial engineer data aided manufacturing, processing of CAD/CAM interface.
A strain tensor that couples to the Madelung stress tensor
Delphenich, D H
2013-01-01
Ordinarily, the stress tensor that one derives for a Madelung fluid is not regarded as being coupled to a strain tensor, which is consistent with the fluid hypothesis. However, based upon earlier work regarding the geometric nature of the quantum potential, one can, in fact, define a strain tensor, which is not, however, due to a deformation of a spatial region, but to a deformation of a frame field on that region. When one expresses the Madelung stress tensor as a function of the strain tensor and its derivatives, one then defines a constitutive law for the Madelung medium that might lead to a more detailed picture of its elementary structure. It is pointed out that the resulting constitutive law is strongly analogous to laws that were presented by Kelvin and Tait for the bending and torsion of elastic wires and plates, as well as the Einstein equations for gravitation if one takes the viewpoint of metric elasticity.
Felician ALECU
2010-01-01
Full Text Available Many professionals and 3D artists consider Blender as being the best open source solution for 3D computer graphics. The main features are related to modeling, rendering, shading, imaging, compositing, animation, physics and particles and realtime 3D/game creation.
3d-3d correspondence revisited
Chung, Hee-Joong; Dimofte, Tudor; Gukov, Sergei; Sułkowski, Piotr
2016-04-01
In fivebrane compactifications on 3-manifolds, we point out the importance of all flat connections in the proper definition of the effective 3d {N}=2 theory. The Lagrangians of some theories with the desired properties can be constructed with the help of homological knot invariants that categorify colored Jones polynomials. Higgsing the full 3d theories constructed this way recovers theories found previously by Dimofte-Gaiotto-Gukov. We also consider the cutting and gluing of 3-manifolds along smooth boundaries and the role played by all flat connections in this operation.
Brdnik, Lovro
2015-01-01
Diplomsko delo analizira trenutno stanje 3D tiskalnikov na trgu. Prikazan je razvoj in principi delovanja 3D tiskalnikov. Predstavljeni so tipi 3D tiskalnikov, njihove prednosti in slabosti. Podrobneje je predstavljena zgradba in delovanje koračnih motorjev. Opravljene so meritve koračnih motorjev. Opisana je programska oprema za rokovanje s 3D tiskalniki in komponente, ki jih potrebujemo za izdelavo. Diploma se oklepa vprašanja, ali je izdelava 3D tiskalnika bolj ekonomična kot pa naložba v ...
Meulien Ohlmann, Odile
2013-02-01
Today the industry offers a chain of 3D products. Learning to "read" and to "create in 3D" becomes an issue of education of primary importance. 25 years professional experience in France, the United States and Germany, Odile Meulien set up a personal method of initiation to 3D creation that entails the spatial/temporal experience of the holographic visual. She will present some different tools and techniques used for this learning, their advantages and disadvantages, programs and issues of educational policies, constraints and expectations related to the development of new techniques for 3D imaging. Although the creation of display holograms is very much reduced compared to the creation of the 90ies, the holographic concept is spreading in all scientific, social, and artistic activities of our present time. She will also raise many questions: What means 3D? Is it communication? Is it perception? How the seeing and none seeing is interferes? What else has to be taken in consideration to communicate in 3D? How to handle the non visible relations of moving objects with subjects? Does this transform our model of exchange with others? What kind of interaction this has with our everyday life? Then come more practical questions: How to learn creating 3D visualization, to learn 3D grammar, 3D language, 3D thinking? What for? At what level? In which matter? for whom?
DYNA3D2000*, Explicit 3-D Hydrodynamic FEM Program
1 - Description of program or function: DYNA3D2000 is a nonlinear explicit finite element code for analyzing 3-D structures and solid continuum. The code is vectorized and available on several computer platforms. The element library includes continuum, shell, beam, truss and spring/damper elements to allow maximum flexibility in modeling physical problems. Many materials are available to represent a wide range of material behavior, including elasticity, plasticity, composites, thermal effects and rate dependence. In addition, DYNA3D has a sophisticated contact interface capability, including frictional sliding, single surface contact and automatic contact generation. 2 - Method of solution: Discretization of a continuous model transforms partial differential equations into algebraic equations. A numerical solution is then obtained by solving these algebraic equations through a direct time marching scheme. 3 - Restrictions on the complexity of the problem: Recent software improvements have eliminated most of the user identified limitations with dynamic memory allocation and a very large format description that has pushed potential problem sizes beyond the reach of most users. The dominant restrictions remain in code execution speed and robustness, which the developers constantly strive to improve
Tournay, Bruno; Rüdiger, Bjarne
2006-01-01
3d digital model af Arkitektskolens gård med virtuel udstilling af afgangsprojekter fra afgangen sommer 2006. 10 s.......3d digital model af Arkitektskolens gård med virtuel udstilling af afgangsprojekter fra afgangen sommer 2006. 10 s....
Roberto Rinaldi
2014-12-01
Full Text Available After an experimental phase of many years, 3D filming is now effective and successful. Improvements are still possible, but the film industry achieved memorable success on 3D movie’s box offices due to the overall quality of its products. Special environments such as space (“Gravity” and the underwater realm look perfect to be reproduced in 3D. “Filming in space” was possible in “Gravity” using special effects and computer graphic. The underwater realm is still difficult to be handled. Underwater filming in 3D was not that easy and effective as filming in 2D, since not long ago. After almost 3 years of research, a French, Austrian and Italian team realized a perfect tool to film underwater, in 3D, without any constrains. This allows filmmakers to bring the audience deep inside an environment where they most probably will never have the chance to be.
Valenza, Enrico
2015-01-01
This book is aimed at the professionals that already have good 3D CGI experience with commercial packages and have now decided to try the open source Blender and want to experiment with something more complex than the average tutorials on the web. However, it's also aimed at the intermediate Blender users who simply want to go some steps further.It's taken for granted that you already know how to move inside the Blender interface, that you already have 3D modeling knowledge, and also that of basic 3D modeling and rendering concepts, for example, edge-loops, n-gons, or samples. In any case, it'
Hundebøl, Jesper
wave of new building information modelling tools demands further investigation, not least because of industry representatives' somewhat coarse parlance: Now the word is spreading -3D digital modelling is nothing less than a revolution, a shift of paradigm, a new alphabet... Research qeustions. Based...... on empirical probes (interviews, observations, written inscriptions) within the Danish construction industry this paper explores the organizational and managerial dynamics of 3D Digital Modelling. The paper intends to - Illustrate how the network of (non-)human actors engaged in the promotion (and arrest) of 3......D Modelling (in Denmark) stabilizes - Examine how 3D Modelling manifests itself in the early design phases of a construction project with a view to discuss the effects hereof for i.a. the management of the building process. Structure. The paper introduces a few, basic methodological concepts...
3D, 9-C anisotropic seismic modeling and inversion
Rusmanugroho, Herurisa
The most complete representation of an elastic medium consists of an elastic tensor with 21 independent moduli. All 21 can be estimated from compressional and shear wave polarization and slowness vectors corresponding to wide apertures of polar and azimuth angles. In isotropic media, when seismic source and receiver components have the same orientation (such as XX and YY), the reflection amplitude contours align approximately perpendicular to the particle motions. The mixed components (such as XY and YX) have amplitude patterns that are in symmetrical pairs of either the same, or of opposite, polarity on either side of the diagonal of the 9-C response matrix. In anisotropic media, amplitude variations with azimuth show the same basic patterns and symmetries as for isotropic, but with a superimposed tendency for alignment parallel to the strike of the vertical cracks. Solutions for elastic tensor elements from synthetic slowness and polarization data calculated directly from the Christoffel equation are more sensitive to the polar angle aperture than to the azimuth aperture. Nine-component synthetic elastic vertical seismic profile data for a model with triclinic symmetry calculated by finite-differencing allows estimation of the elastic 21 tensor elements in the vicinity of a three-component borehole receiver. Wide polar angle and azimuth apertures are needed for accurately estimating the elastic tensor elements. The tensor elements become less independent as the data apertures decrease. Results obtained by extracting slowness and polarization data from the corresponding synthetic seismograms show similar results. The inversion algorithm has produced good results from field vertical seismic profile data set from the Weyburn Field in Southern Saskatchewan in Canada. Synthetic nine-component seismograms calculated from the extracted tensor are able to explain most of the significant features in the field data. The inverted stiffness elastic tensor shows orthorhombic
Lively, Michael
2010-01-01
Professional Papervision3D describes how Papervision3D works and how real world applications are built, with a clear look at essential topics such as building websites and games, creating virtual tours, and Adobe's Flash 10. Readers learn important techniques through hands-on applications, and build on those skills as the book progresses. The companion website contains all code examples, video step-by-step explanations, and a collada repository.
Alaoui-Ismaili, N.; Guy, P.; Chassignole, B.
2014-02-01
The aim of this work is to measure the complex elastic tensor and Euler angles in very complex anisotropic media like austenitic steel welds, by inverse problem resolution from experimental data. The obtained experimental characteristics of the anisotropic material will be injected in a FE code developed by EDF enabling the simulation of an actual ultrasonic NDE of welds. The present work aims to provide reliable input data to the 3D future development of the code. In particular, this complex elastic tensor will allow to predict by modeling beam skewing ant attenuation in an austenitic weld. The investigation of such anisotropic media is very complex because of the directional dependency of the elastic stiffness tensor. Then we will discuss the use of a hybrid genetic algorithm to overcome this difficulty. The identification method is based on waveforms spectra reconstruction associated to a physical model describing wave propagation in plates, during underwater measurements. The entire procedure is qualified and validated using simulated data. Moreover, a comparison of the estimated elastic coefficients with literature values and ultrasonic measurements obtained in transmission is also given, at the end of the paper.
Alaoui-Ismaili, N. [INSA-Lyon, MATEIS, UMR5510 Villeurbanne, F-69621 (France); Guy, P. [INSA-Lyon, LVA, EA677 Villeurbanne, F-69621 (France); Chassignole, B. [EDF R and D, Moret sur Loing, F77818 (France)
2014-02-18
The aim of this work is to measure the complex elastic tensor and Euler angles in very complex anisotropic media like austenitic steel welds, by inverse problem resolution from experimental data. The obtained experimental characteristics of the anisotropic material will be injected in a FE code developed by EDF enabling the simulation of an actual ultrasonic NDE of welds. The present work aims to provide reliable input data to the 3D future development of the code. In particular, this complex elastic tensor will allow to predict by modeling beam skewing ant attenuation in an austenitic weld. The investigation of such anisotropic media is very complex because of the directional dependency of the elastic stiffness tensor. Then we will discuss the use of a hybrid genetic algorithm to overcome this difficulty. The identification method is based on waveforms spectra reconstruction associated to a physical model describing wave propagation in plates, during underwater measurements. The entire procedure is qualified and validated using simulated data. Moreover, a comparison of the estimated elastic coefficients with literature values and ultrasonic measurements obtained in transmission is also given, at the end of the paper.
3D Spectroscopic Instrumentation
Bershady, Matthew A
2009-01-01
In this Chapter we review the challenges of, and opportunities for, 3D spectroscopy, and how these have lead to new and different approaches to sampling astronomical information. We describe and categorize existing instruments on 4m and 10m telescopes. Our primary focus is on grating-dispersed spectrographs. We discuss how to optimize dispersive elements, such as VPH gratings, to achieve adequate spectral resolution, high throughput, and efficient data packing to maximize spatial sampling for 3D spectroscopy. We review and compare the various coupling methods that make these spectrographs ``3D,'' including fibers, lenslets, slicers, and filtered multi-slits. We also describe Fabry-Perot and spatial-heterodyne interferometers, pointing out their advantages as field-widened systems relative to conventional, grating-dispersed spectrographs. We explore the parameter space all these instruments sample, highlighting regimes open for exploitation. Present instruments provide a foil for future development. We give an...
Halskov, Kim; Johansen, Stine Liv; Bach Mikkelsen, Michelle
2014-01-01
Three-dimensional projection installations are particular kinds of augmented spaces in which a digital 3-D model is projected onto a physical three-dimensional object, thereby fusing the digital content and the physical object. Based on interaction design research and media studies, this article...... contributes to the understanding of the distinctive characteristics of such a new medium, and identifies three strategies for designing 3-D projection installations: establishing space; interplay between the digital and the physical; and transformation of materiality. The principal empirical case, From...... Fingerplan to Loop City, is a 3-D projection installation presenting the history and future of city planning for the Copenhagen area in Denmark. The installation was presented as part of the 12th Architecture Biennale in Venice in 2010....
Francisco R. Feito Higueruela
2010-04-01
Full Text Available Applications of Geographical Information Systems on several Archeology fields have been increasing during the last years. Recent avances in these technologies make possible to work with more realistic 3D models. In this paper we introduce a new paradigm for this system, the GIS Thetrahedron, in which we define the fundamental elements of GIS, in order to provide a better understanding of their capabilities. At the same time the basic 3D characteristics of some comercial and open source software are described, as well as the application to some samples on archeological researchs
Iliesiu, Luca; Kos, Filip; Poland, David; Pufu, Silviu S.; Simmons-Duffin, David; Yacoby, Ran
2016-03-01
We study the conformal bootstrap for a 4-point function of fermions in 3D. We first introduce an embedding formalism for 3D spinors and compute the conformal blocks appearing in fermion 4-point functions. Using these results, we find general bounds on the dimensions of operators appearing in the ψ × ψ OPE, and also on the central charge C T . We observe features in our bounds that coincide with scaling dimensions in the GrossNeveu models at large N . We also speculate that other features could coincide with a fermionic CFT containing no relevant scalar operators.
Ms. Swapnali R. Ghadge
2013-01-01
In today’s ever-shifting media landscape, it can be a complex task to find effective ways to reach your desired audience. As traditional media such as television continue to lose audience share, one venue in particular stands out for its ability to attract highly motivated audiences and for its tremendous growth potential the 3D Internet. The concept of '3D Internet' has recently come into the spotlight in the R&D arena, catching the attention of many people, and leading to a lot o...
Kotek, L.
2015-01-01
This paper is about 3D scan of plaster dental casts. The main aim of the work is a hardware and software proposition of 3D scan system for scanning of dental casts. There were used camera, projector and rotate table for this scanning system. Surface triangulation was used, taking benefits of projections of structured light on object, which is being scanned. The rotate table is controlled by PC. The camera, projector and rotate table are synchronized by PC. Controlling of stepper motor is prov...
Villaume, René Domine; Ørstrup, Finn Rude
2002-01-01
Projektet undersøger potentialet for interaktiv 3D design via Internettet. Arkitekt Jørn Utzons projekt til Espansiva blev udviklet som et byggesystem med det mål, at kunne skabe mangfoldige planmuligheder og mangfoldige facade- og rumudformninger. Systemets bygningskomponenter er digitaliseret som...... 3D elementer og gjort tilgængelige. Via Internettet er det nu muligt at sammenstille og afprøve en uendelig række bygningstyper som systemet blev tænkt og udviklet til....
M.M. Voormolen
2007-01-01
textabstractThree dimensional (3D) echocardiography has recently developed from an experimental technique in the â€™90 towards an imaging modality for the daily clinical practice. This dissertation describes the considerations, implementation, validation and clinical application of a unique
Hejlesen, Aske K.; Ovesen, Nis
2012-01-01
This paper presents an experimental approach to teaching 3D modelling techniques in an Industrial Design programme. The approach includes the use of tangible free form models as tools for improving the overall learning. The paper is based on lecturer and student experiences obtained through...
Stenholt, Rasmus; Madsen, Claus B.
2011-01-01
Enabling users to shape 3-D boxes in immersive virtual environments is a non-trivial problem. In this paper, a new family of techniques for creating rectangular boxes of arbitrary position, orientation, and size is presented and evaluated. These new techniques are based solely on position data...
Klusoň, Jindřich
2010-01-01
Computer animation has a growing importance and application in the world. With expansion of technologies increases quality of the final animation as well as number of 3D animation software. This thesis is currently mapped animation software for creating animation in film, television industry and video games which are advisable users requirements. Of them were selected according to criteria the best - Autodesk Maya 2011. This animation software is unique with tools for creating special effects...
Andringa, Roel; de Roo, Mees; Hohm, Olaf; Sezgin, Ergin; Townsend, Paul K
2009-01-01
We construct the N=1 three-dimensional supergravity theory with cosmological, Einstein-Hilbert, Lorentz Chern-Simons, and general curvature squared terms. We determine the general supersymmetric configuration, and find a family of supersymmetric adS vacua with the supersymmetric Minkowski vacuum as a limiting case. Linearizing about the Minkowski vacuum, we find three classes of unitary theories; one is the supersymmetric extension of the recently discovered `massive 3D gravity'. Another is a `new topologically massive supergravity' (with no Einstein-Hilbert term) that propagates a single (2,3/2) helicity supermultiplet.
Andringa, Roel; Bergshoeff, Eric A; De Roo, Mees; Hohm, Olaf [Centre for Theoretical Physics, University of Groningen, Nijenborgh 4, 9747 AG Groningen (Netherlands); Sezgin, Ergin [George and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, Texas A and M University, College Station, TX 77843 (United States); Townsend, Paul K, E-mail: E.A.Bergshoeff@rug.n, E-mail: O.Hohm@rug.n, E-mail: sezgin@tamu.ed, E-mail: P.K.Townsend@damtp.cam.ac.u [Department of Applied Mathematics and Theoretical Physics, Centre for Mathematical Sciences, University of Cambridge, Wilberforce Road, Cambridge, CB3 0WA (United Kingdom)
2010-01-21
We construct the N=1 three-dimensional supergravity theory with cosmological, Einstein-Hilbert, Lorentz Chern-Simons, and general curvature squared terms. We determine the general supersymmetric configuration, and find a family of supersymmetric adS vacua with the supersymmetric Minkowski vacuum as a limiting case. Linearizing about the Minkowski vacuum, we find three classes of unitary theories; one is the supersymmetric extension of the recently discovered 'massive 3D gravity'. Another is a 'new topologically massive supergravity' (with no Einstein-Hilbert term) that propagates a single (2,3/2) helicity supermultiplet.
Ms. Swapnali R. Ghadge
2013-08-01
Full Text Available In today’s ever-shifting media landscape, it can be a complex task to find effective ways to reach your desired audience. As traditional media such as television continue to lose audience share, one venue in particular stands out for its ability to attract highly motivated audiences and for its tremendous growth potential the 3D Internet. The concept of '3D Internet' has recently come into the spotlight in the R&D arena, catching the attention of many people, and leading to a lot of discussions. Basically, one can look into this matter from a few different perspectives: visualization and representation of information, and creation and transportation of information, among others. All of them still constitute research challenges, as no products or services are yet available or foreseen for the near future. Nevertheless, one can try to envisage the directions that can be taken towards achieving this goal. People who take part in virtual worlds stay online longer with a heightened level of interest. To take advantage of that interest, diverse businesses and organizations have claimed an early stake in this fast-growing market. They include technology leaders such as IBM, Microsoft, and Cisco, companies such as BMW, Toyota, Circuit City, Coca Cola, and Calvin Klein, and scores of universities, including Harvard, Stanford and Penn State.
The structures of linear chain Fe(II) spin-crossover compounds of α,β- and α,ω-bis (tetrazol-1-yl)alkane type ligands are described in relation to their magnetic properties. The first threefold interlocked 3-D catenane Fe(II) spin-transition system, [p-tris(1,4-bis(tetrazol-1-yl)butane-N1,N1') iron(II)] bis(perchlorate), will be discussed. An analysis is made among the structures and the cooperativity of the spin-crossover behavior of polynuclear Fe(II) spin-transition materials. Refs. 42 (author)
Holographic renormalization of 3D minimal massive gravity
Alishahiha, Mohsen; Qaemmaqami, Mohammad; Naseh, Ali; Shirzad, Ahmad
2016-01-01
We study holographic renormalization of 3D minimal massive gravity using the Chern-Simons-like formulation of the model. We explicitly present Gibbons- Hawking term as well as all counterterms needed to make the action finite in terms of dreibein and spin-connection. This can be used to find correlation functions of stress tensor of holographic dual field theory.
A spin-4 analog of 3D massive gravity
Bergshoeff, Eric A.; Kovacevic, Marija; Rosseel, Jan; Townsend, Paul K.; Yin, Yihao
2011-01-01
A sixth-order, but ghost-free, gauge-invariant action is found for a fourth-rank symmetric tensor potential in a three-dimensional (3D) Minkowski spacetime. It propagates two massive modes of spin 4 that are interchanged by parity and is thus a spin-4 analog of linearized 'new massive gravity'. Also
NIKE3D, Static and Dynamic Response of 3-D Solids
1 - Description of program or function: NIKE3D is a vectorized, fully implicit, three-dimensional, finite-element program for analyzing the finite-strain, static and dynamic response of inelastic solids, shells, and beams. Capabilities currently available include sliding interfaces, body force loads due to base acceleration, body force loads due to spinning (geometry dependent), concentrated nodal loads, pressure boundary conditions (geometry dependent), displacement boundary conditions, thermal stresses, plot-file generation, and problem restart. 2 - Method of solution: Spatial discretization is accomplished using eight-node solid elements that are integrated using a two-point Gauss quadrature rule, four-node shell elements that use a 2x2 Gauss integration in the surface and from 1 to 5 integration points through the thickness, and two-node thick beam elements using from 1 to 5 integration points in each local direction of the cross section at the center of the beam. Elastic, orthotropic elastic, elastic-plastic, thermo-elastic-plastic, soil and crushable foam, linear viscoelastic, thermo-orthotropic elastic, and thermo-elastic- creep material models are implemented. An incremental-interactive numerical algorithm, based on the Green-Naghdi and Jaumann stress rate formulations, is implemented in NIKE3D. Five iterative schemes are available including the BFGS quasi-Newton method, Broyden's quasi-Newton method, modified Newton-Raphson, full Newton, and full Newton with line search
Hausman, Kalani Kirk
2014-01-01
Get started printing out 3D objects quickly and inexpensively! 3D printing is no longer just a figment of your imagination. This remarkable technology is coming to the masses with the growing availability of 3D printers. 3D printers create 3-dimensional layered models and they allow users to create prototypes that use multiple materials and colors. This friendly-but-straightforward guide examines each type of 3D printing technology available today and gives artists, entrepreneurs, engineers, and hobbyists insight into the amazing things 3D printing has to offer. You'll discover methods for
Szkandera, Jan
2009-01-01
Tato bakalářská práce se zabývá návrhem a realizací systému, který umožní obraz scény zobrazovaný na ploše vnímat prostorově. Prostorové vnímání 2D obrazové informace je umožněno jednak stereopromítáním a jednak tím, že se obraz mění v závislosti na poloze pozorovatele. Tato práce se zabývá hlavně druhým z těchto problémů. This Bachelor's thesis goal is to design and realize system, which allows user to perceive 2D visual information as three-dimensional. 3D visual preception of 2D image i...
Mobile tomographs often have the problem that high spatial resolution is impossible owing to the position or setup of the tomograph. While the tree tomograph developed by Messrs. Isotopenforschung Dr. Sauerwein GmbH worked well in practice, it is no longer used as the spatial resolution and measuring time are insufficient for many modern applications. The paper shows that the mechanical base of the method is sufficient for 3D CT measurements with modern detectors and X-ray tubes. CT measurements with very good statistics take less than 10 min. This means that mobile systems can be used, e.g. in examinations of non-transportable cultural objects or monuments. Enhancement of the spatial resolution of mobile tomographs capable of measuring in any position is made difficult by the fact that the tomograph has moving parts and will therefore have weight shifts. With the aid of tomographies whose spatial resolution is far higher than the mechanical accuracy, a correction method is presented for direct integration of the Feldkamp algorithm
X3D: Extensible 3D Graphics Standard
Daly, Leonard; Brutzman, Don
2007-01-01
The article of record as published may be located at http://dx.doi.org/10.1109/MSP.2007.905889 Extensible 3D (X3D) is the open standard for Web-delivered three-dimensional (3D) graphics. It specifies a declarative geometry definition language, a run-time engine, and an application program interface (API) that provide an interactive, animated, real-time environment for 3D graphics. The X3D specification documents are freely available, the standard can be used without paying any royalties,...
3D game environments create professional 3D game worlds
Ahearn, Luke
2008-01-01
The ultimate resource to help you create triple-A quality art for a variety of game worlds; 3D Game Environments offers detailed tutorials on creating 3D models, applying 2D art to 3D models, and clear concise advice on issues of efficiency and optimization for a 3D game engine. Using Photoshop and 3ds Max as his primary tools, Luke Ahearn explains how to create realistic textures from photo source and uses a variety of techniques to portray dynamic and believable game worlds.From a modern city to a steamy jungle, learn about the planning and technological considerations for 3D modelin
Aboufadel, Edward F.
2014-01-01
The purpose of this short paper is to describe a project to manufacture a regular octohedron on a 3D printer. We assume that the reader is familiar with the basics of 3D printing. In the project, we use fundamental ideas to calculate the vertices and faces of an octohedron. Then, we utilize the OPENSCAD program to create a virtual 3D model and an STereoLithography (.stl) file that can be used by a 3D printer.
Higher derivative extensions of 3 d Chern–Simons models: conservation laws and stability
Kaparulin, D. S.; Karataeva, I. Yu.; Lyakhovich, S. L.
2015-01-01
We consider the class of higher derivative $3d$ vector field models with the field equation operator being a polynomial of the Chern-Simons operator. For $n$-th order theory of this type, we provide a general receipt for constructing $n$-parameter family of conserved second rank tensors. The family includes the canonical energy-momentum tensor, which is unbounded, while there are bounded conserved tensors that provide classical stability of the system for certain combinations of the parameter...
Rodrigues, Marcos; Robinson, Alan; Alboul, Lyuba; Brink, Willie
2006-01-01
3D face recognition is an open field. In this paper we present a method for 3D facial recognition based on Principal Components Analysis. The method uses a relatively large number of facial measurements and ratios and yields reliable recognition. We also highlight our approach to sensor development for fast 3D model acquisition and automatic facial feature extraction.
3-D contextual Bayesian classifiers
Larsen, Rasmus
distribution for the pixel values as well as a prior distribution for the configuration of class variables within the cross that is made of a pixel and its four nearest neighbours. We will extend these algorithms to 3-D, i.e. we will specify a simultaneous Gaussian distribution for a pixel and its 6 nearest 3......-D neighbours, and generalise the class variable configuration distributions within the 3-D cross given in 2-D algorithms. The new 3-D algorithms are tested on a synthetic 3-D multivariate dataset....
Taming Supersymmetric Defects in 3d-3d Correspondence
Gang, Dongmin; Romo, Mauricio; Yamazaki, Masahito
2015-01-01
We study knots in 3d Chern-Simons theory with complex gauge group $SL(N,\\mathbb{C})$, in the context of its relation with 3d $\\mathcal{N}=2$ theory (the so-called 3d-3d correspondence). The defect has either co-dimension 2 or co-dimension 4 inside the 6d $(2,0)$ theory, which is compactified on a 3-manifold $\\hat{M}$. We identify such defects in various corners of the 3d-3d correspondence, namely in 3d $SL(N,\\mathbb{C})$ Chern-Simons theory, in 3d $\\mathcal{N}=2$ theory, in 5d $\\mathcal{N}=2$ super Yang-Mills theory, and in the M-theory holographic dual. We can make quantitative checks of the 3d-3d correspondence by computing partition functions at each of these theories. This Letter is a companion to a longer paper, which contains more details and more results.
Burgos, Gaël.; Capdeville, Yann; Guillot, Laurent
2016-06-01
We investigate the effect of small-scale heterogeneities close to a seismic explosive source, at intermediate periods (20-50 s), with an emphasis on the resulting nonisotropic far-field radiation. First, using a direct numerical approach, we show that small-scale elastic heterogeneities located in the near-field of an explosive source, generate unexpected phases (i.e., long period S waves). We then demonstrate that the nonperiodic homogenization theory applied to 2-D and 3-D elastic models, with various pattern of small-scale heterogeneities near the source, leads to accurate waveforms at a reduced computational cost compared to direct modeling. Further, it gives an interpretation of how nearby small-scale features interact with the source at low frequencies, through an explicit correction to the seismic moment tensor. In 2-D simulations, we find a deviatoric contribution to the moment tensor, as high as 21% for near-source heterogeneities showing a 25% contrast of elastic values (relative to a homogeneous background medium). In 3-D this nonisotropic contribution reaches 27%. Second, we analyze intermediate-periods regional seismic waveforms associated with some underground nuclear explosions conducted at the Nevada National Security Site and invert for the full moment tensor, in order to quantify the relative contribution of the isotropic and deviatoric components of the tensor. The average value of the deviatoric part is about 35%. We conclude that the interactions between an explosive source and small-scale local heterogeneities of moderate amplitude may lead to a deviatoric contribution to the seismic moment, close to what is observed using regional data from nuclear test explosions.
Rubber Impact on 3D Textile Composites
Heimbs, Sebastian; Van Den Broucke, Björn; Duplessis Kergomard, Yann; Dau, Frederic; Malherbe, Benoit
2012-06-01
A low velocity impact study of aircraft tire rubber on 3D textile-reinforced composite plates was performed experimentally and numerically. In contrast to regular unidirectional composite laminates, no delaminations occur in such a 3D textile composite. Yarn decohesions, matrix cracks and yarn ruptures have been identified as the major damage mechanisms under impact load. An increase in the number of 3D warp yarns is proposed to improve the impact damage resistance. The characteristic of a rubber impact is the high amount of elastic energy stored in the impactor during impact, which was more than 90% of the initial kinetic energy. This large geometrical deformation of the rubber during impact leads to a less localised loading of the target structure and poses great challenges for the numerical modelling. A hyperelastic Mooney-Rivlin constitutive law was used in Abaqus/Explicit based on a step-by-step validation with static rubber compression tests and low velocity impact tests on aluminium plates. Simulation models of the textile weave were developed on the meso- and macro-scale. The final correlation between impact simulation results on 3D textile-reinforced composite plates and impact test data was promising, highlighting the potential of such numerical simulation tools.
3D Cell Culture in Alginate Hydrogels
Therese Andersen
2015-03-01
Full Text Available This review compiles information regarding the use of alginate, and in particular alginate hydrogels, in culturing cells in 3D. Knowledge of alginate chemical structure and functionality are shown to be important parameters in design of alginate-based matrices for cell culture. Gel elasticity as well as hydrogel stability can be impacted by the type of alginate used, its concentration, the choice of gelation technique (ionic or covalent, and divalent cation chosen as the gel inducing ion. The use of peptide-coupled alginate can control cell–matrix interactions. Gelation of alginate with concomitant immobilization of cells can take various forms. Droplets or beads have been utilized since the 1980s for immobilizing cells. Newer matrices such as macroporous scaffolds are now entering the 3D cell culture product market. Finally, delayed gelling, injectable, alginate systems show utility in the translation of in vitro cell culture to in vivo tissue engineering applications. Alginate has a history and a future in 3D cell culture. Historically, cells were encapsulated in alginate droplets cross-linked with calcium for the development of artificial organs. Now, several commercial products based on alginate are being used as 3D cell culture systems that also demonstrate the possibility of replacing or regenerating tissue.
Harvesting vibrations via 3D phononic isolators
Psarobas, Ioannis E.; Yannopapas, Vassilios; Matikas, Theodore E.
2016-05-01
We report on the existence of unidirectional phononic band gaps that may span over extended regions of the Brillouin zone and can find application in trapping elastic (acoustic) waves in properly designed multilayered 3D structures. Phononic isolators operate as a result of asymmetrical wave transmission through a slab of a crystallographic phononic structure with broken mirror symmetry. Due to the use of lossless materials in the crystal, the absorption rate is dramatically enhanced when the proposed isolator is placed next to a vibrational harvesting cell. xml:lang="fr"
Bjoerck, A.; Thor, S.E. [Aeronautical Research Inst. of Sweden, Bromma (Sweden)
1996-12-01
The JOULE II project `Dynamic stall and 3D effects` started in January 1994 and was completed in September 1995. The objective of the project has been to increase the understanding of the three-dimensional and unsteady aerodynamics of stall controlled HAWT`s. The objectives have also been to develop `engineering models` suitable for inclusion into aero-elastic codes. The project included the participation of 13 parties within Europe. This paper describes an overview of the work carried out within the project and key results. 3 refs, 4 figs
Asymptotic expansion of the multi-orientable random tensor model
Fusy, Eric
2014-01-01
Three-dimensional random tensor models are a natural generalization of the celebrated matrix models. The associated tensor graphs, or 3D maps, can be classified with respect to a particular integer or half-integer, the degree of the respective graph. In this paper we analyze the general term of the asymptotic expansion in N, the size of the tensor, of a particular random tensor model, the multi-orientable tensor model. We perform their enumeration and we establish which are the dominant configurations of a given degree.
Three-dimensional conjugate gradients inversion of magnetotelluric impedance tensor
Complete text of publication follows. We have developed a three-dimensional conjugate gradients inversion algorithm of magnetotelluric (MT) impedance tensor. This inversion algorithm can be used to invert 3D MT area data and 2D MT profile data. The importance of including the on-diagonal impedance tensor terms, Zxx and Zyy, in 3D inversion is considered in the synthetic example. The results from the 3D inversion of synthetic 3D area data indicate that the results using all impedance tensor elements, Zxx, Zxy, Zyx and Zyy, in the inversion are better than those without using on-diagonal elements. Therefore, 3D inversion should include all tensor elements if possible. The results from the 3D inversion of synthetic 2D profile data show that reasonable images for structures beneath and near a single MT profile can be obtained if all tensor elements are used in the inversion of the profile data. Currently, most of MT data are still collected on 2D profiles. This inversion algorithm of impedance tensor can promote the practical application of 3D MT inversion. We also use this inversion algorithm to invert the real data in the kayabe area in Japan. From the trial inversion with the synthetic and real data, the validity and practicability of this inversion algorithm is verified.
On microcontinuum field theories: the Eshelby stress tensor and incompatibility conditions
Lazar, Markus; Maugin, Gerard A.
2015-01-01
We investigate linear theories of incompatible micromorphic elasticity, incompatible microstretch elasticity, incompatible micropolar elasticity and the incompatible dilatation theory of elasticity (elasticity with voids). The incompatibility conditions and Bianchi identities are derived and discussed. The Eshelby stress tensor (static energy momentum) is calculated for such inhomogeneous media with microstructure. Its divergence gives the driving forces for dislocations, disclinations, point...
Tensor Target Polarization at TRIUMF
Smith, G
2014-10-27
The first measurements of tensor observables in $\\pi \\vec{d}$ scattering experiments were performed in the mid-80's at TRIUMF, and later at SIN/PSI. The full suite of tensor observables accessible in $\\pi \\vec{d}$ elastic scattering were measured: $T_{20}$, $T_{21}$, and $T_{22}$. The vector analyzing power $iT_{11}$ was also measured. These results led to a better understanding of the three-body theory used to describe this reaction. %Some measurements were also made in the absorption and breakup channels. A direct measurement of the target tensor polarization was also made independent of the usual NMR techniques by exploiting the (nearly) model-independent result for the tensor analyzing power at 90$^\\circ _{cm}$ in the $\\pi \\vec{d} \\rightarrow 2p$ reaction. This method was also used to check efforts to enhance the tensor polarization by RF burning of the NMR spectrum. A brief description of the methods developed to measure and analyze these experiments is provided.
3D Printing Functional Nanocomposites
Leong, Yew Juan
2016-01-01
3D printing presents the ability of rapid prototyping and rapid manufacturing. Techniques such as stereolithography (SLA) and fused deposition molding (FDM) have been developed and utilized since the inception of 3D printing. In such techniques, polymers represent the most commonly used material for 3D printing due to material properties such as thermo plasticity as well as its ability to be polymerized from monomers. Polymer nanocomposites are polymers with nanomaterials composited into the ...
3D IBFV : Hardware-Accelerated 3D Flow Visualization
Telea, Alexandru; Wijk, Jarke J. van
2003-01-01
We present a hardware-accelerated method for visualizing 3D flow fields. The method is based on insertion, advection, and decay of dye. To this aim, we extend the texture-based IBFV technique for 2D flow visualization in two main directions. First, we decompose the 3D flow visualization problem in a
3D Elevation Program—Virtual USA in 3D
Lukas, Vicki; Stoker, J.M.
2016-01-01
The U.S. Geological Survey (USGS) 3D Elevation Program (3DEP) uses a laser system called ‘lidar’ (light detection and ranging) to create a virtual reality map of the Nation that is very accurate. 3D maps have many uses with new uses being discovered all the time.
The 3D model: explaining densification and deformation mechanisms by using 3D parameter plots.
Picker, Katharina M
2004-04-01
The aim of the study was to analyze very differently deforming materials using 3D parameter plots and consequently to gain deeper insights into the densification and deformation process described with the 3D model in order to define an ideal tableting excipient. The excipients used were dicalcium phosphate dihydrate (DCPD), sodium chloride (NaCl), microcrystalline cellulose (MCC), xylitol, mannitol, alpha-lactose monohydrate, maltose, hydroxypropyl methylcellulose (HPMC), sodium carboxymethylcellulose (NaCMC), cellulose acetate (CAC), maize starch, potato starch, pregelatinized starch, and maltodextrine. All of the materials were tableted to graded maximum relative densities (rhorel, max) using an eccentric tableting machine. The data which resulted, namely force, displacement, and time, were analyzed by the application of 3D modeling. Different particle size fractions of DCPD, CAC, and MCC were analyzed in addition. Brittle deforming materials such as DCPD exhibited a completely different 3D parameter plot, with low time plasticity, d, and low pressure plasticity, e, and a strong decrease in omega values when densification increased, in contrast to the plastically deforming MCC, which had much higher d, e, and omega values. e and omega values changed only slightly when densification increased for MCC. NaCl showed less of a decrease in omega values than DCPD did, and the d and e values were between those of MCC and DCPD. The sugar alcohols, xylitol and mannitol, behaved in a similar fashion to sodium chloride. This is also valid for the crystalline sugars, alpha-lactose monohydrate, and maltose. However, the sugars are more brittle than the sugar alcohols. The cellulose derivatives, HPMC, NaCMC, and CAC, are as plastic as MCC, however, their elasticity depends on substitution indicated by lower (more elastic) or higher (less elastic) omega values. The native starches, maize starch and potato starch, are very elastic, and pregelatinized starch and maltodextrine are
Connell, Ellery
2011-01-01
Helping graphic designers expand their 2D skills into the 3D space The trend in graphic design is towards 3D, with the demand for motion graphics, animation, photorealism, and interactivity rapidly increasing. And with the meteoric rise of iPads, smartphones, and other interactive devices, the design landscape is changing faster than ever.2D digital artists who need a quick and efficient way to join this brave new world will want 3D for Graphic Designers. Readers get hands-on basic training in working in the 3D space, including product design, industrial design and visualization, modeling, ani
Griffey, Jason
2014-01-01
As the maker movement continues to grow and 3-D printers become more affordable, an expanding group of hobbyists is keen to explore this new technology. In the time-honored tradition of introducing new technologies, many libraries are considering purchasing a 3-D printer. Jason Griffey, an early enthusiast of 3-D printing, has researched the marketplace and seen several systems first hand at the Consumer Electronics Show. In this report he introduces readers to the 3-D printing marketplace, covering such topics asHow fused deposition modeling (FDM) printing workBasic terminology such as build
Larsen, Rasmus
1997-01-01
. This includes the specification of a Gaussian distribution for the pixel values as well as a prior distribution for the configuration of class variables within the cross that is m ade of a pixel and its four nearest neighbours. We will extend this algorithm to 3-D, i.e. we will specify a simultaneous Gaussian...... distr ibution for a pixel and its 6 nearest 3-D neighbours, and generalise the class variable configuration distribution within the 3-D cross. The algorithm is tested on a synthetic 3-D multivariate dataset....
3D Bayesian contextual classifiers
Larsen, Rasmus
2000-01-01
We extend a series of multivariate Bayesian 2-D contextual classifiers to 3-D by specifying a simultaneous Gaussian distribution for the feature vectors as well as a prior distribution of the class variables of a pixel and its 6 nearest 3-D neighbours.......We extend a series of multivariate Bayesian 2-D contextual classifiers to 3-D by specifying a simultaneous Gaussian distribution for the feature vectors as well as a prior distribution of the class variables of a pixel and its 6 nearest 3-D neighbours....
Interactive 3D multimedia content
Cellary, Wojciech
2012-01-01
The book describes recent research results in the areas of modelling, creation, management and presentation of interactive 3D multimedia content. The book describes the current state of the art in the field and identifies the most important research and design issues. Consecutive chapters address these issues. These are: database modelling of 3D content, security in 3D environments, describing interactivity of content, searching content, visualization of search results, modelling mixed reality content, and efficient creation of interactive 3D content. Each chapter is illustrated with example a
The physical property tensors of one-dimensional quasicrystals
Li Cui-Lian; Liu You-Yan
2004-01-01
According to the group representation theory, we derive the character formulae of representation-matrices of the physical property tensors for the one-dimensional (1D) quasicrystals. Based on this, we have calculated the numbers of independent components of representation-matrices for thermal expansion coefficient tensors, piezoelectric coefficient tensors and elastic constant tensors under 31 point-groups for the 1D quasicrystals. Moreover, we have deduced the particular matrix forms of these tensors under the 31 point-groups. This is an important complement of quasicrystal physical property.
3-D Experimental Fracture Analysis at High Temperature
John H. Jackson; Albert S. Kobayashi
2001-09-14
T*e, which is an elastic-plastic fracture parameter based on incremental theory of plasticity, was determined numerically and experimentally. The T*e integral of a tunneling crack in 2024-T3 aluminum, three point bend specimen was obtained through a hybrid analysis of moire interferometry and 3-D elastic-plastic finite element analysis. The results were verified by the good agreement between the experimentally and numerically determined T*e on the specimen surface.
Goda, Ibrahim; Ganghoffer, Jean-François
2015-11-01
The purpose of this paper is to develop a homogeneous, orthotropic couple-stress continuum model as a substitute of the 3D periodic heterogeneous cellular solid model of vertebral trabecular bone. Vertebral trabecular bone is modeled as a porous material with an idealized periodic structure made of 3D open cubic cells, which is effectively orthotropic. The chosen architecture is based on studies of samples taken from the central part of vertebral bodies. The effective properties are obtained based on the response of the representative volume element under prescribed boundary conditions. Mixed boundary conditions comprising both traction and displacement boundary conditions are applied on the structure boundaries. In this contribution, the effective mechanical constants of the effective couple-stress continuum are deduced by an equivalent strain energy method. The characteristic lengths for bending and torsion are identified from the resulting homogenized orthotropic moduli. We conduct this study computationally using a finite element approach. Vertebral trabecular bone is modeled either as a cellular solid or as a two-phase material consisting of bone tissue (stiff phase) forming a trabecular network, and a surrounding soft tissue referring to the bone marrow present in the pores. Both the bone tissue forming the network and the pores are assumed to be homogeneous linear elastic, and isotropic media. The scale effects on the predicted couple stress moduli of these networks are investigated by varying the size of the bone specimens over which the boundary conditions are applied. The analysis using mixed boundary conditions gives results that are independent of unit cell size when computing the first couple stress tensor, while it is dependent on the cell size as to the second couple stress tensor moduli. This study provides overall guidance on how the size of the trabecular specimen influence couple stresses elastic moduli of cellular materials, with focus on bones
无
2003-01-01
The disadvantage remaining in 3D scanning system and its reasons are discussed. A new host-and-slave structure with high speed image acquisition and processing system is proposed to quicken the image processing and improve the performance of 3D scanning system.
ECT Team, Purdue
2015-01-01
Building Bytes, by Brian Peters, is a project that uses desktop 3D printers to print bricks for architecture. Instead of using an expensive custom-made printer, it uses a normal standard 3D printer which is available for everyone and makes it more accessible and also easier for fabrication.
MT3D was first developed by Chunmiao Zheng in 1990 at S.S. Papadopulos & Associates, Inc. with partial support from the U.S. Environmental Protection Agency (USEPA). Starting in 1990, MT3D was released as a pubic domain code from the USEPA. Commercial versions with enhanced capab...
Strictly nonnegative tensors and nonnegative tensor partition
HU ShengLong; HUANG ZhengHai; QI LiQun
2014-01-01
We introduce a new class of nonnegative tensors—strictly nonnegative tensors.A weakly irreducible nonnegative tensor is a strictly nonnegative tensor but not vice versa.We show that the spectral radius of a strictly nonnegative tensor is always positive.We give some necessary and su？cient conditions for the six wellconditional classes of nonnegative tensors,introduced in the literature,and a full relationship picture about strictly nonnegative tensors with these six classes of nonnegative tensors.We then establish global R-linear convergence of a power method for finding the spectral radius of a nonnegative tensor under the condition of weak irreducibility.We show that for a nonnegative tensor T,there always exists a partition of the index set such that every tensor induced by the partition is weakly irreducible;and the spectral radius of T can be obtained from those spectral radii of the induced tensors.In this way,we develop a convergent algorithm for finding the spectral radius of a general nonnegative tensor without any additional assumption.Some preliminary numerical results show the feasibility and effectiveness of the algorithm.
Wood, Jo; Kirschenbauer, Sabine; Döllner, Jürgen;
2005-01-01
The notion of three-dimensionality is applied to five stages of the visualization pipeline. While 3D visulization is most often associated with the visual mapping and representation of data, this chapter also identifies its role in the management and assembly of data, and in the media used...... to display 3D imagery. The extra cartographic degree of freedom offered by using 3D is explored and offered as a motivation for employing 3D in visualization. The use of VR and the construction of virtual environments exploit navigational and behavioral realism, but become most usefil when combined...... with abstracted representations embedded in a 3D space. The interactions between development of geovisualization, the technology used to implement it and the theory surrounding cartographic representation are explored. The dominance of computing technologies, driven particularly by the gaming industry...
Finding Apparent Horizons in Dynamic 3D Numerical Spacetimes
Anninos, P.; Camarda, K.; Libson, J.; Masso, J.; Seidel, E; Suen, W.
1996-01-01
We have developed a general method for finding apparent horizons in 3D numerical relativity. Instead of solving for the partial differential equation describing the location of the apparent horizons, we expand the closed 2D surfaces in terms of symmetric trace--free tensors and solve for the expansion coefficients using a minimization procedure. Our method is applied to a number of different spacetimes, including numerically constructed spacetimes containing highly distorted axisymmetric blac...
Reeb Graph for Automatic 3D Cephalometry
Mestiri Makram
2014-03-01
Full Text Available The purpose of this study is to present a method of three-dimensional computed tomographic (3D-CT cephalometrics and its use to study cranio/maxilla-facial malformations. We propose a system for automatic localization of cephalometric landmarks using reeb graphs. Volumetric images of a patient were reconstructed into 3D mesh. The proposed method is carried out in three steps: we begin by applying 3d mesh skull simplification, this mesh was reconstructed from a head volumetric medical image, and then we extract a reeb graph. Reeb graph mesh extraction represents a skeleton composed in a number of nodes and arcs. We are interested in the node position; we noted that some reeb nodes could be considered as cephalometric landmarks under specific conditions. The third step is to identify these nodes automatically by using elastic mesh registration using “thin plate” transformation and clustering. Preliminary results show a landmarks recognition rate of more than 90%, very close to the manually provided landmarks positions made by a medical stuff.
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITHOUT TURB3D)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITH TURB3D)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
ADT-3D Tumor Detection Assistant in 3D
Jaime Lazcano Bello
2008-12-01
Full Text Available The present document describes ADT-3D (Three-Dimensional Tumor Detector Assistant, a prototype application developed to assist doctors diagnose, detect and locate tumors in the brain by using CT scan. The reader may find on this document an introduction to tumor detection; ADT-3D main goals; development details; description of the product; motivation for its development; result’s study; and areas of applicability.
Unassisted 3D camera calibration
Atanassov, Kalin; Ramachandra, Vikas; Nash, James; Goma, Sergio R.
2012-03-01
With the rapid growth of 3D technology, 3D image capture has become a critical part of the 3D feature set on mobile phones. 3D image quality is affected by the scene geometry as well as on-the-device processing. An automatic 3D system usually assumes known camera poses accomplished by factory calibration using a special chart. In real life settings, pose parameters estimated by factory calibration can be negatively impacted by movements of the lens barrel due to shaking, focusing, or camera drop. If any of these factors displaces the optical axes of either or both cameras, vertical disparity might exceed the maximum tolerable margin and the 3D user may experience eye strain or headaches. To make 3D capture more practical, one needs to consider unassisted (on arbitrary scenes) calibration. In this paper, we propose an algorithm that relies on detection and matching of keypoints between left and right images. Frames containing erroneous matches, along with frames with insufficiently rich keypoint constellations, are detected and discarded. Roll, pitch yaw , and scale differences between left and right frames are then estimated. The algorithm performance is evaluated in terms of the remaining vertical disparity as compared to the maximum tolerable vertical disparity.
无
2002-01-01
This paper presents a principle and a method of col or 3D laser scanning measurement. Based on the fundamental monochrome 3D measureme nt study, color information capture, color texture mapping, coordinate computati on and other techniques are performed to achieve color 3D measurement. The syste m is designed and composed of a line laser light emitter, one color CCD camera, a motor-driven rotary filter, a circuit card and a computer. Two steps in captu ring object's images in the measurement process: Firs...
3-D neutron transport benchmarks
A set of 3-D neutron transport benchmark problems proposed by the Osaka University to NEACRP in 1988 has been calculated by many participants and the corresponding results are summarized in this report. The results of Keff, control rod worth and region-averaged fluxes for the four proposed core models, calculated by using various 3-D transport codes are compared and discussed. The calculational methods used were: Monte Carlo, Discrete Ordinates (Sn), Spherical Harmonics (Pn), Nodal Transport and others. The solutions of the four core models are quite useful as benchmarks for checking the validity of 3-D neutron transport codes
Puntar, Matej
2012-01-01
The purpose of this thesis is the presentation of already established and new technologies of displaying 3D content in a web browser. The thesis begins with a short presentation of the history of 3D content available on the internet and its development together with advantages and disadvantages of individual technologies. The latter two are described in detail as well is their use and the differences among them. Special emphasis has been given to WebGL, the newest technology of 3D conte...
Grutle, Øyvind Kallevik
2015-01-01
3D printers have in recent years become extremely popular. Even though 3D printing technology have existed since the late 1980's, it is now considered one of the most significant technological breakthroughs of the twenty-first century. Several different 3D printing processes have been invented during the years. But it is the fused deposition modeling (FDM) which was one of the first invented that is considered the most popular today. Even though the FDM process is the most popular, it still s...
Garrou , Philip; Ramm , Peter
2014-01-01
Edited by key figures in 3D integration and written by top authors from high-tech companies and renowned research institutions, this book covers the intricate details of 3D process technology.As such, the main focus is on silicon via formation, bonding and debonding, thinning, via reveal and backside processing, both from a technological and a materials science perspective. The last part of the book is concerned with assessing and enhancing the reliability of the 3D integrated devices, which is a prerequisite for the large-scale implementation of this emerging technology. Invaluable reading fo
Lin, Zeyu
2014-01-01
3D printing technology is introduced and defined in this Thesis. Some methods of 3D printing are illustrated and their principles are explained with pictures. Most of the essential parts are presented with pictures and their effects are explained within the whole system. Problems on Up! Plus 3D printer are solved and a DIY product is made with this machine. The processes of making product are recorded and the items which need to be noticed during the process are the highlight in this th...
Tuotekehitysprojekti: 3D-tulostin
Pihlajamäki, Janne
2011-01-01
Opinnäytetyössä tutustuttiin 3D-tulostamisen teknologiaan. Työssä käytiin läpi 3D-tulostimesta tehty tuotekehitysprojekti. Sen lisäksi esiteltiin yleisellä tasolla tuotekehitysprosessi ja syntyneiden tulosten mahdollisia suojausmenetelmiä. Tavoitteena tässä työssä oli kehittää markkinoilta jo löytyvää kotitulostin-tasoista 3D-laiteteknologiaa lähemmäksi ammattilaistason ratkaisua. Tavoitteeseen pyrittiin keskittymällä parantamaan laitteella saavutettavaa tulostustarkkuutta ja -nopeutt...
Gravitational Superenergy Tensor
Mashhoon, Bahram; McClune, James C.; Quevedo, Hernando
1996-01-01
We provide a physical basis for the local gravitational superenergy tensor. Furthermore, our gravitoelectromagnetic deduction of the Bel-Debever-Robinson superenergy tensor permits the identification of the gravitational stress-energy tensor. This {\\it local} gravitational analog of the Maxwell stress-energy tensor is illustrated for a plane gravitational wave.
The 3S1-3D1 mixing parameter ε1 and the 3PJ phase shifts in the N-N system at low energy are investigated in the three nucleon system. The sensitivity of polarization transfer coefficients in p-d elastic scattering is determined by rigorous Faddeev calculations using Paris and Bonn potentials. The experimental results of Kysup(y') favours the tensor force component of the Bonn A potential. (orig.)
Main: TATCCAYMOTIFOSRAMY3D [PLACE
Full Text Available TATCCAYMOTIFOSRAMY3D S000256 01-August-2006 (last modified) kehi TATCCAY motif found in rice (O. ... otif and G motif (see S000130) are responsible for sugar ... repression (Toyofuku et al. 1998); GATA; amylase; ...
Yang, Xu; Zhang, Yong; Yang, Chenghua; Xu, Lu; Wang, Qiang; Zhao, Yuan
2016-06-01
Conventional three dimensional (3D) ghost imaging measures range of target based on pulse fight time measurement method. Due to the limit of data acquisition system sampling rate, range resolution of the conventional 3D ghost imaging is usually low. In order to take off the effect of sampling rate to range resolution of 3D ghost imaging, a heterodyne 3D ghost imaging (HGI) system is presented in this study. The source of HGI is a continuous wave laser instead of pulse laser. Temporal correlation and spatial correlation of light are both utilized to obtain the range image of target. Through theory analysis and numerical simulations, it is demonstrated that HGI can obtain high range resolution image with low sampling rate.
Conducting polymer 3D microelectrodes
Sasso, Luigi; Vazquez, Patricia; Vedarethinam, Indumathi;
2010-01-01
Conducting polymer 3D microelectrodes have been fabricated for possible future neurological applications. A combination of micro-fabrication techniques and chemical polymerization methods has been used to create pillar electrodes in polyaniline and polypyrrole. The thin polymer films obtained...
Hess, Siegfried
2015-01-01
This book presents the science of tensors in a didactic way. The various types and ranks of tensors and the physical basis is presented. Cartesian Tensors are needed for the description of directional phenomena in many branches of physics and for the characterization the anisotropy of material properties. The first sections of the book provide an introduction to the vector and tensor algebra and analysis, with applications to physics, at undergraduate level. Second rank tensors, in particular their symmetries, are discussed in detail. Differentiation and integration of fields, including generalizations of the Stokes law and the Gauss theorem, are treated. The physics relevant for the applications in mechanics, quantum mechanics, electrodynamics and hydrodynamics is presented. The second part of the book is devoted to tensors of any rank, at graduate level. Special topics are irreducible, i.e. symmetric traceless tensors, isotropic tensors, multipole potential tensors, spin tensors, integration and spin-...
Ahmed, Zeeshan
2010-01-01
In this short paper I briefly discuss 3D war Game based on artificial intelligence concepts called AI WAR. Going in to the details, I present the importance of CAICL language and how this language is used in AI WAR. Moreover I also present a designed and implemented 3D War Cybug for AI WAR using CAICL and discus the implemented strategy to defeat its enemies during the game life.
Lading, Brian; Larsen, Rasmus; Åström, Kalle
2006-01-01
We build a 3d face shape model, including inter- and intra-shape variations, derive the analytical jacobian of its resulting 2d rendered image, and show example of its fitting performance with light, pose, id, expression and texture variations.}......We build a 3d face shape model, including inter- and intra-shape variations, derive the analytical jacobian of its resulting 2d rendered image, and show example of its fitting performance with light, pose, id, expression and texture variations.}...
Lading, Brian; Larsen, Rasmus; Astrom, K
2006-01-01
We build a 3D face shape model, including inter- and intra-shape variations, derive the analytical Jacobian of its resulting 2D rendered image, and show example of its fitting performance with light, pose, id, expression and texture variations......We build a 3D face shape model, including inter- and intra-shape variations, derive the analytical Jacobian of its resulting 2D rendered image, and show example of its fitting performance with light, pose, id, expression and texture variations...
Brox, Thomas; Weickert, Joachim; Burgeth, Bernhard; Mrázek, Pavel
2004-01-01
In this article we introduce nonlinear versions of the popular structure tensor, also known as second moment matrix. These nonlinear structure tensors replace the Gaussian smoothing of the classical structure tensor by discontinuity-preserving nonlinear diffusions. While nonlinear diffusion is a well-established tool for scalar and vector-valued data, it has not often been used for tensor images so far. Two types of nonlinear diffusion processes for tensor data are studied: an isotropic one w...
Design of manufacturable 3D extremal elastic microstructure
Andreassen, Erik; Lazarov, Boyan Stefanov; Sigmund, Ole
2014-01-01
, and the deviations are due to the imposed manufacturing constraints. The designs are generated using topology optimization. Due to high resolution and the imposed robustness requirement they are manufacturable without any need for post-processing. This has been validated by the manufacturing of an...
3D viscous-spring artificial boundary in time domain
Liu Jingbo; Du Yixin; Du Xiuli; Wang Zhenyu; Wu Jun
2006-01-01
After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary (VSAB) in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.
Pomeron as a Reggeized Tensor Glueball
MA Wei-Xing; A.W.Thomas; SHEN Peng-Nian; ZHOU Li-Juan
2001-01-01
We study gluonic content of the pomeron and propose that the pomeron could be a reggeized tensor glueball ζ(2230) with quantum numbers IG JPc = 0+2++.This conjecture is examined in high energy proton-proton elastic scattering,and the calculations lend a favorable support to our physical idea.``
Dima, M.; Farisato, G.; Bergomi, M.; Viotto, V.; Magrin, D.; Greggio, D.; Farinato, J.; Marafatto, L.; Ragazzoni, R.; Piazza, D.
2014-08-01
In the last few years 3D printing is getting more and more popular and used in many fields going from manufacturing to industrial design, architecture, medical support and aerospace. 3D printing is an evolution of bi-dimensional printing, which allows to obtain a solid object from a 3D model, realized with a 3D modelling software. The final product is obtained using an additive process, in which successive layers of material are laid down one over the other. A 3D printer allows to realize, in a simple way, very complex shapes, which would be quite difficult to be produced with dedicated conventional facilities. Thanks to the fact that the 3D printing is obtained superposing one layer to the others, it doesn't need any particular work flow and it is sufficient to simply draw the model and send it to print. Many different kinds of 3D printers exist based on the technology and material used for layer deposition. A common material used by the toner is ABS plastics, which is a light and rigid thermoplastic polymer, whose peculiar mechanical properties make it diffusely used in several fields, like pipes production and cars interiors manufacturing. I used this technology to create a 1:1 scale model of the telescope which is the hardware core of the space small mission CHEOPS (CHaracterising ExOPlanets Satellite) by ESA, which aims to characterize EXOplanets via transits observations. The telescope has a Ritchey-Chrétien configuration with a 30cm aperture and the launch is foreseen in 2017. In this paper, I present the different phases for the realization of such a model, focusing onto pros and cons of this kind of technology. For example, because of the finite printable volume (10×10×12 inches in the x, y and z directions respectively), it has been necessary to split the largest parts of the instrument in smaller components to be then reassembled and post-processed. A further issue is the resolution of the printed material, which is expressed in terms of layers
YouDash3D: exploring stereoscopic 3D gaming for 3D movie theaters
Schild, Jonas; Seele, Sven; Masuch, Maic
2012-03-01
Along with the success of the digitally revived stereoscopic cinema, events beyond 3D movies become attractive for movie theater operators, i.e. interactive 3D games. In this paper, we present a case that explores possible challenges and solutions for interactive 3D games to be played by a movie theater audience. We analyze the setting and showcase current issues related to lighting and interaction. Our second focus is to provide gameplay mechanics that make special use of stereoscopy, especially depth-based game design. Based on these results, we present YouDash3D, a game prototype that explores public stereoscopic gameplay in a reduced kiosk setup. It features live 3D HD video stream of a professional stereo camera rig rendered in a real-time game scene. We use the effect to place the stereoscopic effigies of players into the digital game. The game showcases how stereoscopic vision can provide for a novel depth-based game mechanic. Projected trigger zones and distributed clusters of the audience video allow for easy adaptation to larger audiences and 3D movie theater gaming.
Remote 3D Medical Consultation
Welch, Greg; Sonnenwald, Diane H.; Fuchs, Henry; Cairns, Bruce; Mayer-Patel, Ketan; Yang, Ruigang; State, Andrei; Towles, Herman; Ilie, Adrian; Krishnan, Srinivas; Söderholm, Hanna M.
Two-dimensional (2D) video-based telemedical consultation has been explored widely in the past 15-20 years. Two issues that seem to arise in most relevant case studies are the difficulty associated with obtaining the desired 2D camera views, and poor depth perception. To address these problems we are exploring the use of a small array of cameras to synthesize a spatially continuous range of dynamic three-dimensional (3D) views of a remote environment and events. The 3D views can be sent across wired or wireless networks to remote viewers with fixed displays or mobile devices such as a personal digital assistant (PDA). The viewpoints could be specified manually or automatically via user head or PDA tracking, giving the remote viewer virtual head- or hand-slaved (PDA-based) remote cameras for mono or stereo viewing. We call this idea remote 3D medical consultation (3DMC). In this article we motivate and explain the vision for 3D medical consultation; we describe the relevant computer vision/graphics, display, and networking research; we present a proof-of-concept prototype system; and we present some early experimental results supporting the general hypothesis that 3D remote medical consultation could offer benefits over conventional 2D televideo.
Dagiuklas, Tasos
2015-01-01
This book describes recent innovations in 3D media and technologies, with coverage of 3D media capturing, processing, encoding, and adaptation, networking aspects for 3D Media, and quality of user experience (QoE). The contributions are based on the results of the FP7 European Project ROMEO, which focuses on new methods for the compression and delivery of 3D multi-view video and spatial audio, as well as the optimization of networking and compression jointly across the future Internet. The delivery of 3D media to individual users remains a highly challenging problem due to the large amount of data involved, diverse network characteristics and user terminal requirements, as well as the user’s context such as their preferences and location. As the number of visual views increases, current systems will struggle to meet the demanding requirements in terms of delivery of consistent video quality to fixed and mobile users. ROMEO will present hybrid networking solutions that combine the DVB-T2 and DVB-NGH broadcas...
Dagiuklas, Tasos
2014-01-01
This book describes recent innovations in 3D media and technologies, with coverage of 3D media capturing, processing, encoding, and adaptation, networking aspects for 3D Media, and quality of user experience (QoE). The main contributions are based on the results of the FP7 European Projects ROMEO, which focus on new methods for the compression and delivery of 3D multi-view video and spatial audio, as well as the optimization of networking and compression jointly across the Future Internet (www.ict-romeo.eu). The delivery of 3D media to individual users remains a highly challenging problem due to the large amount of data involved, diverse network characteristics and user terminal requirements, as well as the user’s context such as their preferences and location. As the number of visual views increases, current systems will struggle to meet the demanding requirements in terms of delivery of constant video quality to both fixed and mobile users. ROMEO will design and develop hybrid-networking solutions that co...
Materialedreven 3d digital formgivning
Hansen, Flemming Tvede
2010-01-01
Formålet med forskningsprojektet er for det første at understøtte keramikeren i at arbejde eksperimenterende med digital formgivning, og for det andet at bidrage til en tværfaglig diskurs om brugen af digital formgivning. Forskningsprojektet fokuserer på 3d formgivning og derved på 3d digital...... formgivning og Rapid Prototyping (RP). RP er en fællesbetegnelse for en række af de teknikker, der muliggør at overføre den digitale form til 3d fysisk form. Forskningsprojektet koncentrerer sig om to overordnede forskningsspørgsmål. Det første handler om, hvordan viden og erfaring indenfor det keramiske...... fagområde kan blive udnyttet i forhold til 3d digital formgivning. Det andet handler om, hvad en sådan tilgang kan bidrage med, og hvordan den kan blive udnyttet i et dynamisk samspil med det keramiske materiale i formgivningen af 3d keramiske artefakter. Materialedreven formgivning er karakteriseret af en...
DYADIC METHOD FOR TENSOR FUNCTIONS
黄永念; 鲁昊
2002-01-01
In this paper, we discuss tensor functions by dyadic representation of tensor. Two different cases of scalar invariants and two different cases of tensor invariants are calculated. It is concluded that there are six independent scale invariants for a symmetrical tensor and an antisymmetrical tensor, and there are twelve invariants for two symmetrical tensors and an antisymmetrical tensor. And we present a new list of tensor invariants for the tensor-valued isotropic function.
Modification of 3D milling machine to 3D printer
Halamíček, Lukáš
2015-01-01
Tato práce se zabývá přestavbou gravírovací frézky na 3D tiskárnu. V první části se práce zabývá možnými technologiemi 3D tisku a možností jejich využití u přestavby. Dále jsou popsány a vybrány vhodné součásti pro přestavbu. V další části je realizováno řízení ohřevu podložky, trysky a řízení posuvu drátu pomocí softwaru TwinCat od společnosti Beckhoff na průmyslovém počítači. Výsledkem práce by měla být oživená 3D tiskárna. This thesis deals with rebuilding of engraving machine to 3D pri...
3D Imager and Method for 3D imaging
Kumar, P.; Staszewski, R.; Charbon, E.
2013-01-01
3D imager comprising at least one pixel, each pixel comprising a photodetectorfor detecting photon incidence and a time-to-digital converter system configured for referencing said photon incidence to a reference clock, and further comprising a reference clock generator provided for generating the re
Becker, T; Roth, M M; Becker, Thomas; Fabrika, Sergei; Roth, Martin M.
2003-01-01
The quantitative spectroscopy of stellar objects in complex environments is mainly limited by the ability of separating the object from the background. Standard slit spectroscopy, restricting the field of view to one dimension, is obviously not the proper technique in general. The emerging Integral Field (3D) technique with spatially resolved spectra of a two-dimensional field of view provides a great potential for applying advanced subtraction methods. In this paper an image reconstruction algorithm to separate point sources and a smooth background is applied to 3D data. Several performance tests demonstrate the photometric quality of the method. The algorithm is applied to real 3D observations of a sample Planetary Nebula in M31, whose spectrum is contaminated by the bright and complex galaxy background. The ability of separating sources is also studied in a crowded stellar field in M33.
TRAB-3D is a reactor dynamics code with three-dimensional neutronics coupled to core and circuit thermal-hydraulics. The code, entirely developed at VTT, can be used in transient and accident analyses of boiling (BWR) and pressurized water (PWR) reactors with rectangular fuel bundle geometry. The validation history of TRAB-3D includes calculation of international benchmark exercises, as well as comparisons with measured data from real plant transients. The most recent validation case is a load rejection test performed at the Olkiluoto 1 nuclear power plant in 1998 in connection with the power uprating project. The fact that there is local power measurement data available from this test makes it a suitable case for three-dimensional core model validation. The agreement between the results of the TRAB-3D calculation and the measurements is very good. (orig.)
3D-grafiikkamoottori mobiililaitteille
Vahlman, Lauri
2014-01-01
Tässä insinöörityössä käydään läpi mobiililaitteille suunnatun yksinkertaisen 3D-grafiikkamoottorin suunnittelu ja toteutus käyttäen OpenGL ES -rajapintaa. Työssä esitellään grafiikkamoottorin toteutuksessa käytettyjä tekniikoita sekä tutustutaan moottorin rakenteeseen ja toteutuksellisiin yksityiskohtiin. Työn alkupuolella tutustutaan myös modernin 3D-grafiikan yleisiin periaatteisiin ja toimintaan sekä käydään läpi 3D-grafiikkaan liittyviä suorituskykyongelmia. Työn loppupuolella esitel...
Walder, Christian; Breidt, Martin; Bulthoff, Heinrich;
2009-01-01
We present a novel algorithm for the markerless tracking of deforming surfaces such as faces. We acquire a sequence of 3D scans along with color images at 40Hz. The data is then represented by implicit surface and color functions, using a novel partition-of-unity type method of efficiently...... combining local regressors using nearest neighbor searches. Both these functions act on the 4D space of 3D plus time, and use temporal information to handle the noise in individual scans. After interactive registration of a template mesh to the first frame, it is then automatically deformed to track...... the scanned surface, using the variation of both shape and color as features in a dynamic energy minimization problem. Our prototype system yields high-quality animated 3D models in correspondence, at a rate of approximately twenty seconds per timestep. Tracking results for faces and other objects...
Pihl, Michael Johannes
The main purpose of this PhD project is to develop an ultrasonic method for 3D vector flow imaging. The motivation is to advance the field of velocity estimation in ultrasound, which plays an important role in the clinic. The velocity of blood has components in all three spatial dimensions, yet...... conventional methods can estimate only the axial component. Several approaches for 3D vector velocity estimation have been suggested, but none of these methods have so far produced convincing in vivo results nor have they been adopted by commercial manufacturers. The basis for this project is the Transverse...... on the TO fields are suggested. They can be used to optimize the TO method. In the third part, a TO method for 3D vector velocity estimation is proposed. It employs a 2D phased array transducer and decouples the velocity estimation into three velocity components, which are estimated simultaneously based on 5...
Focused high energy ion beam micromachining is the newest of the micromachining techniques. There are about 50 scanning proton microprobe facilities worldwide, but so far only few of them showed activity in this promising field. High energy ion beam micromachining using a direct-write scanning MeV ion beam is capable of producing 3D microstructures and components with well defined lateral and depth geometry. The technique has high potential in the manufacture of 3D molds, stamps, and masks for X-ray lithography (LIGA), and also in the rapid prototyping of microcomponents either for research purposes or for components testing prior to batch production. (R.P.)
3D Computations and Experiments
Couch, R; Faux, D; Goto, D; Nikkel, D
2004-04-05
This project consists of two activities. Task A, Simulations and Measurements, combines all the material model development and associated numerical work with the materials-oriented experimental activities. The goal of this effort is to provide an improved understanding of dynamic material properties and to provide accurate numerical representations of those properties for use in analysis codes. Task B, ALE3D Development, involves general development activities in the ALE3D code with the focus of improving simulation capabilities for problems of mutual interest to DoD and DOE. Emphasis is on problems involving multi-phase flow, blast loading of structures and system safety/vulnerability studies.
3D Printing: Exploring Capabilities
Samuels, Kyle; Flowers, Jim
2015-01-01
As 3D printers become more affordable, schools are using them in increasing numbers. They fit well with the emphasis on product design in technology and engineering education, allowing students to create high-fidelity physical models to see and test different iterations in their product designs. They may also help students to "think in three…
Manos, Harry
2016-01-01
Visual aids are important to student learning, and they help make the teacher's job easier. Keeping with the "TPT" theme of "The Art, Craft, and Science of Physics Teaching," the purpose of this article is to show how teachers, lacking equipment and funds, can construct a durable 3-D model reference frame and a model gravity…
Lading, Brian; Larsen, Rasmus; Åström, Kalle
2006-01-01
We build a 3d face shape model, including inter- and intra-shape variations, derive the analytical jacobian of its resulting 2d rendered image, and show example of its fitting performance with light, pose, id, expression and texture variations.}
Lading, Brian; Larsen, Rasmus; Astrom, K
2006-01-01
We build a 3D face shape model, including inter- and intra-shape variations, derive the analytical Jacobian of its resulting 2D rendered image, and show example of its fitting performance with light, pose, id, expression and texture variations
Krajnović, Davor
2016-01-01
Thanks to a technique that reveals galaxies in 3D, astronomers can now show that many galaxies have been wrongly classified. Davor Krajnovi\\'c argues that the classification scheme proposed 85 years ago by Edwin Hubble now needs to be revised.
Pedersen, Pernille Klarskov; Strikwerda, Andrew; Wang, Tianwu;
2013-01-01
We present a characterization of THz beams generated in both a two-color air plasma and in a LiNbO3 crystal. Using a commercial THz camera, we record intensity images as a function of distance through the beam waist, from which we extract 2D beam profiles and visualize our measurements into 3D beam...
3D electrohydrodynamic simulation of electrowetting displays
The fluid dynamic behavior within a pixel of an electrowetting display (EWD) is thoroughly investigated through a 3D simulation. By coupling the electrohydrodynamic (EHD) force deduced from the Maxwell stress tensor with the laminar phase field of the oil–water dual phase, the complete switch processes of an EWD, including the break-up and the electrowetting stages in the switch-on process (with voltage) and the oil spreading in the switch-off process (without voltage), are successfully simulated. By considering the factor of the change in the apparent contact angle at the contact line, the electro–optic performance obtained from the simulation is found to agree well with its corresponding experiment. The proposed model is used to parametrically predict the effect of interfacial (e.g. contact angle of grid) and geometric (e.g. oil thickness and pixel size) properties on the defects of an EWD, such as oil dewetting patterns, oil overflow, and oil non-recovery. With the help of the defect analysis, a highly stable EWD is both experimentally realized and numerically analyzed. (paper)
Introduction to Tensor Calculus
Sochi, Taha
2016-01-01
These are general notes on tensor calculus which can be used as a reference for an introductory course on tensor algebra and calculus. A basic knowledge of calculus and linear algebra with some commonly used mathematical terminology is presumed.
Priprava 3D modelov za 3D tisk
Pikovnik, Tomaž
2015-01-01
Po mnenju nekaterih strokovnjakov bo aditivna proizvodnja (ali 3D tiskanje) spremenila proizvodnjo industrijo, saj si bo vsak posameznik lahko natisnil svoj objekt po želji. V diplomski nalogi so predstavljene nekatere tehnologije aditivne proizvodnje. V nadaljevanju diplomske naloge je predstavljena izdelava makete hiše v merilu 1:100, vse od modeliranja do tiskanja. Poseben poudarek je posvečen predelavi modela, da je primeren za tiskanje, kjer je razvit pristop za hitrejše i...
Post processing of 3D models for 3D printing
Pikovnik, Tomaž
2015-01-01
According to the opinion of some experts the additive manufacturing or 3D printing will change manufacturing industry, because any individual could print their own model according to his or her wishes. In this graduation thesis some of the additive manufacturing technologies are presented. Furthermore in the production of house scale model in 1:100 is presented, starting from modeling to printing. Special attention is given to postprocessing of the building model elements us...