Application of the p-version of the finite-element method to global-local problems

A brief survey is given of some recent developments in finite-element analysis technology which bear upon the three main research areas under consideration in this workshop: (1) analysis methods; (2) software testing and quality assurance; and (3) parallel processing. The variational principle incorporated in a finite-element computer program, together with a particular set of input data, determines the exact solution corresponding to that input data. Most finite-element analysis computer programs are based on the principle of virtual work. In the following, researchers consider only programs based on the principle of virtual work and denote the exact displacement vector field corresponding to some specific set of input data by vector u(EX). The exact solution vector u(EX) is independent of the design of the mesh or the choice of elements. Except for very simple problems, or specially constructed test problems, vector u(EX) is not known. Researchers perform a finite-element analysis (or any other numerical analysis) because they wish to make conclusions concerning the response of a physical system to certain imposed conditions, as if vector u(EX) were known

Solution of geometrically nonlinear statics problems by the p-version of the finite element method

This project is concerned with the possibility of using computers for the simulation of structural systems with the same degree of reliability as full scale physical experiments. Reliable numerical simulation will make it possible to reduce the costs of engineering and improve the quality of engineering decisions based on computed information. An error of idealization is an error between the actual physical quantities on which engineering decisions are based (e.g., maximum principal stress, first natural frequency, etc.) and the same data corresponding to the exact solution of the mathematical model. An error of discretization is an error between the quantities of interest corresponding to the exact and approximate solutions of a mathematical model. A high degree of reliability can be achieved in numerical simulation only if both the errors of idealization and errors of discretization can be shown to be small

Hierarchic models for laminated plates

The research conducted in the formulation of hierarchic models for laminated plates is described. The work is an extension of the work done for laminated strips. The use of a single parameter, beta, is investigated that represents the degree to which the equilibrium equations of three dimensional elasticity are satisfied. The powers of beta identify members in the hierarchic sequence. Numerical examples that were analyzed with the proposed sequence of models are included. The results obtained for square plates with uniform loading and homogeneous boundary conditions are very encouraging. Several cross-ply and angle-ply laminates were evaluated and the results compared with those of the fully three dimensional model, computed using MSC/PROBE, and with previously reported work on laminated strips

Hierarchic plate and shell models based on p-extension

Formulations of finite element models for beams, arches, plates and shells based on the principle of virtual work was studied. The focus is on computer implementation of hierarchic sequences of finite element models suitable for numerical solution of a large variety of practical problems which may concurrently contain thin and thick plates and shells, stiffeners, and regions where three dimensional representation is required. The approximate solutions corresponding to the hierarchic sequence of models converge to the exact solution of the fully three dimensional model. The stopping criterion is based on: (1) estimation of the relative error in energy norm; (2) equilibrium tests, and (3) observation of the convergence of quantities of interest

Elastic-plastic models for multi-site damage

This paper presents recent developments in advanced analysis methods for the computation of stress site damage. The method of solution is based on the p-version of the finite element method. Its implementation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experimental results are available for the panel. The plasticity analysis provided the same limit load as the experimentally determined load. The results of elastic-plastic analysis were compared with the results of linear elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of net-section yielding was determined also. The results show that crack growth rate is accelerated by the presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideally-plastic material model is compared against the Ramberg-Osgood model

Solution of elastic-plastic stress analysis problems by the p-version of the finite element method

The solution of small strain elastic-plastic stress analysis problems by the p-version of the finite element method is discussed. The formulation is based on the deformation theory of plasticity and the displacement method. Practical realization of controlling discretization errors for elastic-plastic problems is the main focus. Numerical examples which include comparisons between the deformation and incremental theories of plasticity under tight control of discretization errors are presented

Refined heart failure detection algorithm for improved clinical reliability of OptiVol alerts in CRT-D recipients

   Background: The reliability of intrathoracic impedance monitoring for prediction of heart failure (HF) by implantable cardiac devices is controversial. Despite using additional device-based parameters described in the PARTNERS HF study, such as new onset of arrhythmias, abnormal autonomics, low biventricular pacing rate or patient activity level, the predictive power of device diagnostic algorithm is still in doubt. The objective of this study was to compare the device diagnostic algorithm described in the PARTNERS HF study to a newly developed algorithm applying refined diagnostic criteria. Methods: Fourty two patients were prospectively enrolled who had been implanted with an intrathoracic impedance and remote monitoring capable implantable cardiac defibrillator with a cardiac resychroniza­tion therapy (CRT-D) device in this observational study. If a remote OptiVolTM alert occurred, patients were checked for presence of HF symptoms. A new algorithm was derived from the original PARTNERS HF criteria, considering more sensitive cut-offs and changes of patterns of the device-based parameters. Results: During an average follow-up of 38 months, 722 remote transmissions were received. From the total of 128 transmissions with OptiVol alerts, 32 (25%) corresponded to true HF events. Upon multivariate discriminant analysis, low patient activity, high nocturnal heart rate, and low CRT pacing (< 90%) proved to be independent predictors of true HF events (all p < 0.01). Incorporating these three refined criteria in a new algorithm, the diagnostic yield of OptiVol was improved by increasing specific­ity from 37.5% to 86.5%, positive predictive value from 34.1% to 69.8% and area under the curve from 0.787 to 0.922 (p < 0.01), without a relevant loss in sensitivity (96.9% vs. 93.8%). Conclusions: A refined device diagnostic algorithm based on the parameters of low activity level, high nocturnal heart rate, and suboptimal biventricular pacing might improve the clinical reliability of OptiVol alerts.  

Three is the magic number -- distance measurement of NGC 3147 using SN 2021hpr and its siblings

The nearby spiral galaxy NGC 3147 hosted three Type Ia supernovae (SNe Ia) in the past decades, which have been subjects of intense follow-up observations. Simultaneous analysis of their data provides a unique opportunity for testing the different light curve fitting methods and distance estimations. The detailed optical follow-up of SN 2021hpr allows us to revise the previous distance estimations to NGC 3147, and compare the widely used light curve fitting algorithms to each other. After the combination of the available and newly published data of SN 2021hpr, its physical properties can be also estimated with higher accuracy. We present and analyse new BVgriz and Swift photometry of SN 2021hpr to constrain its general physical properties. Together with its siblings, SNe 1997bq and 2008fv, we cross-compare the individual distance estimates of these three SNe given by the SALT code, and also check their consistency with the results from the MLCS2k2 method. The early spectral series of SN 2021hpr are also fit with the radiative spectral code TARDIS in order to verify the explosion properties and constrain the chemical distribution of the outer ejecta. After combining the distance estimates for the three SNe, the mean distance to their host galaxy, NGC 3127, is 42.5 $\pm$ 1.0 Mpc, which matches with the distance inferred by the most up-to-date LC fitters, SALT3 and BayeSN. We confirm that SN~2021hpr is a Branch-normal Type Ia SN that ejected $\sim 1.12 \pm 0.28$ M$_\odot$ from its progenitor white dwarf, and synthesized $\sim 0.44 \pm 0.14$ M$_\odot$ of radioactive $^{56}$Ni.Comment: 16 pages, 17 figures, 11 tables; accepted for publication in A&

Recruitment and inhibitory action of hippocampal axo-axonic cells during behavior.

The axon initial segment of hippocampal pyramidal cells is a key subcellular compartment for action potential generation, under GABAergic control by the "chandelier" or axo-axonic cells (AACs). Although AACs are the only cellular source of GABA targeting the initial segment, their in vivo activity patterns and influence over pyramidal cell dynamics are not well understood. We achieved cell-type-specific genetic access to AACs in mice and show that AACs in the hippocampal area CA1 are synchronously activated by episodes of locomotion or whisking during rest. Bidirectional intervention experiments in head-restrained mice performing a random foraging task revealed that AACs inhibit CA1 pyramidal cells, indicating that the effect of GABA on the initial segments in the hippocampus is inhibitory in vivo. Finally, optogenetic inhibition of AACs at specific track locations induced remapping of pyramidal cell place fields. These results demonstrate brain-state-specific dynamics of a critical inhibitory controller of cortical circuits