A novel dynamic constitutive micromechanical model to predict the strain rate dependent mechanical behavior of glass/epoxy laminated composites

In the present research, a novel dynamic constitutive micromechanical (DCM) model was developed to predict the strain rate dependent mechanical behavior of laminated glass/epoxy composites. The present model is an integration of the generalized strain rate dependent constitutive model as a constitutive model for the neat polymer, the plasticity model of Huang as a micromechanical model, and dynamic progressive failure criteria. This model is able to predict the longitudinal and transverse tensile and in-plane shear behaviors of unidirectional glass/epoxy composites with arbitrary fiber volume fractions at arbitrary strain rates. The present model can also predict the stress-strain behavior of laminated composites with different layups and fiber volume fractions at arbitrary strain rates. A comparison between the results predicted by the present model and the available experimental data showed that the model predicts the strain rate dependent mechanical behavior of glass/epoxy composites with very good accuracy.     

Prediction of environmental stress cracking in polycarbonate by molecular modeling

Environmental stress cracking (ESC) is a premature failure of a polymer exposed to a fluid, under stress which is much less than its yield stress. Many experimental works were done before in an effort to predict experimentally the ESC potential of a fluid in different polymers. None of the previous works applied molecular modeling techniques to predict this potential so this work is a pioneering work. This study's goal was to apply atomistic molecular modeling techniques to gain a better understanding of the ESC mechanism and to predict the ESC potential of different fluids in polymers. Our model experimental system was amorphous polycarbonate (PC) with water as an ESC fluid. The computational study was expanded to include a high level ESC fluid for PC–toluene and a non ESC fluid–BD, together with the moderate ESC fluid–water. A clear distinction between ESC fluids and non ESC fluids for PC was achieved by means of molecular modeling. The experimental work approved that water is an ESC fluid for PC as predicted in the computational part. Copyright © 2014 John Wiley & Sons, Ltd.     

玻璃硅烷化处理对罗丹明6G和亚甲基蓝吸附行为的影响

利用六甲基二硅烷胺对平面玻璃光波导(高折射率透明导光薄膜介质)进行硅烷化处理, 得到水接触角大于90°的疏水表面. 然后使用时间分辨光波导分光光谱技术研究水溶液中的罗丹明6G (R6G)和亚甲基蓝(MB)分子在疏水玻璃表面的吸附行为, 并与亲水玻璃条件下测得的结果进行对比. 对利用疏水玻璃光波导测得的R6G的吸附-脱附动力学曲线进行Langmuir拟合得到了R6G的吸附速率常数, 脱附速率常数以及吸附自由能. 并且发现与亲水玻璃情况相比, 吸附速率常数增大, 脱附速率常数减小, 吸附自由能更负. 在疏水玻璃表面形成的R6G和MB吸附层的吸光度与亲水玻璃情况相比显著升高, 表明这两种分子更倾向于吸附在疏水玻璃表面. 实验结果还发现玻璃硅烷化处理能够有效抑制这两种染料分子在表面的聚合反应.     

Parameterization of Models for Vapor Solubility in Semicrystalline Polyethylene

This paper extends the previous article by the authors on the solubility of hydrocarbon vapors in semicrystalline polyethylenes produced in the gas phase process. That work demonstrates a computational model for solubility based on an activity coefficient modification of the Sako–Wu–Prausnitz equation of state. In that work, by fitting a key parameter, one related to the constraint of tie chains on polymer fluid behavior, to a single isopentane solubility isotherm, accurate predictions of hydrocarbon solubility in polymer granules over a range of temperature, pressure, and composition are reported. In the present work, additional experimental solubility data are reported, an error in the authors' previous article is corrected, and a useful parameterization method that improves the predictive capabilities of the model is demonstrated. By using the model to predict much of the authors' own experimental data, as well as that published by others in the field, it is demonstrated how the proposed parameterization method allows for accurate predictions using a limited amount of experimental measurements.     

A closed form solution for the vulcanization prediction of NR cured with sulphur and different accelerators

The paper presents a novel efficient closed form approach to determine the degree of vulcanization of natural rubber (NR) vulcanized with sulphur in presence of different accelerators. The general reaction scheme proposed by Han and co-workers for vulcanized sulphur NR is re-adapted and suitably modified taking into account the single contributions of the different accelerators, focusing in particular on some experimental data, where NR was vulcanized at different temperatures (from 150 to \(180\, ^{\circ }\hbox {C}\)) and concentrations of sulphur, using TBBS and DPG in the mixture as co-agents at variable concentrations. In the model, chain reactions initiated by the formation of macro-compounds responsible for the formation of the unmatured crosslinked polymer are accounted for. It is assumed that such reactions depend on the reciprocal concentrations of all components and their chemical nature. In presence of two accelerators, reactions are assumed to proceed in parallel, making the assumption that there is no interaction between the two accelerators. Despite there is experimental evidence that a weak process by which each accelerator affects the other, the reaction chemistry is still not well understood and therefore its effect cannot be translated into any mathematical model. In any case, even disregarding such interaction, good approximations of the rheometer curves are obtained. From the simplified kinetic scheme adopted, a closed form solution is found for the crosslink density, with the only limitation that the induction period is excluded from computations. The main capability of the model stands however in the closed form determination of kinetic constants representing the velocities of single reactions in the kinetic scheme adopted, which allows avoiding a numerically demanding least-squares best fitting on rheometer experimental data. Two series of experiments available, relying into rheometer curves at different temperatures and different concentrations of sulphur and accelerators, are utilized to evaluate the fitting capabilities of the mathematical model. Very good agreement between numerical output and experimental data is experienced in all cases analyzed.     

Recursive neural networks prediction of glass transition temperature from monomer structure: an application to acrylic and methacrylic polymers

We propose a new method based on a Recursive Neural Network (RecNN) for predicting polymer properties from their structured molecular representations. RecNN allows for a completely novel approach to QSPR analysis by direct adaptive processing of molecular graphs. This model joins the representational power of structured domains with Neural Network ability to capture underlying complex relationships in the data by a process of training from examples. To this aim, a structured representation was designed for the modelling of polymer structures. The adopted representation can account also for average macromolecule characteristics, such as degree of polymerization, stereoregularity, comonomer distribution. To begin with, this model was applied to the prediction of the glass transition temperature of (meth)acrylic polymers with different degree of main chain tacticity. The results so far obtained indicate that the proposed representation of polymer structure can convey information on both the repeating unit structure and average polymer features. The ability of the proposed RecNN method of treating this structured representation makes this method more general and flexible with respect to standard literature methods. Moreover, the same model can handle at the same time the Tg of polymer samples present in only one tacticity form together with that of polymer with different stereoregularity.     

Development of a polarizable force field for proteins via ab initio quantum chemistry: first generation model and gas phase tests

We present results of developing a methodology suitable for producing molecular mechanics force fields with explicit treatment of electrostatic polarization for proteins and other molecular system of biological interest. The technique allows simulation of realistic-size systems. Employing high-level ab initio data as a target for fitting allows us to avoid the problem of the lack of detailed experimental data. Using the fast and reliable quantum mechanical methods supplies robust fitting data for the resulting parameter sets. As a result, gas-phase many-body effects for dipeptides are captured within the average RMSD of 0.22 kcal/mol from their ab initio values, and conformational energies for the di- and tetrapeptides are reproduced within the average RMSD of 0.43 kcal/mol from their quantum mechanical counterparts. The latter is achieved in part because of application of a novel torsional fitting technique recently developed in our group, which has already been used to greatly improve accuracy of the peptide conformational equilibrium prediction with the OPLS-AA force field.1 Finally, we have employed the newly developed first-generation model in computing gas-phase conformations of real proteins, as well as in molecular dynamics studies of the systems. The results show that, although the overall accuracy is no better than what can be achieved with a fixed-charges model, the methodology produces robust results, permits reasonably low computational cost, and avoids other computational problems typical for polarizable force fields. It can be considered as a solid basis for building a more accurate and complete second-generation model.     

Statistical mechanics of worm-like polymers from a new generating function

We present a mathematical approach to the worm-like chain model of semiflexible polymers. Our method is built on a novel generating function from which all the properties of the model can be derived. Moreover, this approach satisfies the local inextensibility constraint exactly. In this paper, we focus on the lowest order contribution to the generating function and derive explicit analytical expressions for the characteristic function, polymer propagator, single chain structure factor, and mean square end-to-end distance. These analytical expressions are valid for polymers with any degree of stiffness and contour length. We find that our calculations are able to capture the fully flexible and infinitely stiff limits of the aforementioned quantities exactly while providing a smooth and approximate crossover behavior for intermediate values of the stiffness of the polymer backbone. In addition, our results are in very good quantitative agreement with the exact and approximate results of five other treatments of semiflexible polymers.     

The rheology of pseudoplastic fluids in porous media using network modeling

This paper considers the rheology of pseudoplastic (shear thinning) fluids in porous media. The central problem studied is the relationship between the viscometric behavior of the polymer solution and its observed behavior in the porous matrix. In the past, a number of macroscopic approaches have been applied, usually based on capillary bundle models of the porous medium. These simplified models have been used along with constitutive equations describing the fluid behavior (usually of power law type) to establish semiempirical macroscopic equations describing the flow of non-Newtonian fluids in porous media. This procedure has been reasonably successful in correlating experimental results on the flow of polymer solutions through both consolidated and unconsolidated porous materials. However, it does not allow an interpretation of polymer flow in porous media in terms of the flows on a microscopic scale; nor does it allow us to predict changes in macroscopic behavior resulting from variations at a microscopic level in the characteristics of the porous medium such as pore size distribution. In this work, we use a network approach to the modeling of non-Newtonian rheology, in order to understand some of the more detailed features of polymjer flow in porous media. This approach provides a mathematical bridge between the behavior of the non-Newtonian fluid in a single capillary and the macroscopic behavior as deduced from the pressure drop-flow rate relation across the whole network model. It demonstrates the importance of flow redistribution within the elements of the capillary network as the overall pressure gradient varies. As an example of a pseudoplastic fluid in a porous medium, we consider the flow of xanthan biopolymer. This polymer is important as a displacing fluid viscosifier in enhanced oil recovery applications and, for that reason, a considerable amount of experimental data has been published on the flow of xanthan solutions in various porous media.     

Nanosize and microsize clay effects on the kinetics of the thermal degradation of polylactides

Polylactide (PLA)-montmorillonite (MMT) micro- and nanocomposites based on semicrystalline and amorphous polymers and unmodified or organomodified clays at 5 wt% content were produced by melt mixing. Based on the three different test methods that were used to follow thermal degradation, different conclusions were obtained. During melt processing, thermomechanical degradation was more pronounced in the presence of all fillers, which apparently acted catalytically, but to different degrees. During isothermal degradation in air from 180 °C to 200 °C, degradation rate constants were calculated from novel equations incorporating changes in intrinsic viscosity (IV). Results show that the thermal degradation rate constants of the amorphous PLA and its composites are lower than those of the semicrystalline PLA and its composites. Due to better filler dispersion in the polymer matrix, the thermal degradation rate constants of the nanocomposites are significantly lower than those of the unfilled polymers and their microcomposites under air. As per dynamic TGA data and thermal kinetic analysis from weight losses and activation energy calculations, organomodified nanofillers have a complex effect on the polymer thermal stability; the unmodified fillers, however, reduce polymer thermal stability. These TGA data and kinetic analysis results also support the findings that the thermal stability of the amorphous PLA and its composites is higher than that of the semicrystalline polymer and its composites and the thermal stability of the nanocomposites is higher than that of the microcomposites. In general, mathematical modeling based on random thermal scission equations was satisfactory for fitting the TGA experimental data.     

Combining configurational entropy and self-concentration to describe the component dynamics in miscible polymer blends

We provide a new approach to describe the component segmental dynamics of miscible polymer blends combining the concept of chain connectivity, expressed in terms of the self-concentration, and the Adam-Gibbs model. The results show an excellent agreement between the prediction of our approach and the experimental data. The self-concentrations obtained yield length scales between 1 and 3.2 nm depending on the temperature, the flexibility of the polymer, expressed in terms of the Kuhn segment, and its concentration in the blends, at temperatures above the glass transition range of the blend.     

Single-particle electrophoresis for studying the adsorption of cationic polymers onto anionic particles

Understanding the adsorption of polymers onto particles is crucial for many technological and biomedical applications. Even though polymer adsorption on particles is a dynamic process, most experimental techniques can only study the adsorption indirectly, in equilibrium and on the ensemble level. New analysis methods are required to overcome these limitations. We investigated the use of single-particle electrophoresis to study the adsorption kinetics of cationic polymers onto anionic particles and compared the resulting data to a theoretical model. In this approach, the electrophoretic mobility of single polystyrene (PS) particles, exposed to different concentrations of poly(2-guanidinoethyl methacrylate), was measured as a function of time. The polymer adsorption leads to an electrophoretic mobility change of the PS particle over time, from the initial negative value to a positive value at equilibrium. By fitting the kinetics data to the Langmuir model, the adsorption rate, desorption rate and equilibrium constant were determined. Finally, the adsorption kinetics of several other polymers was investigated. This showed that the presented technique enables direct analysis and comparison of the kinetics of polymer adsorption on the single-particle level.     

Self consistent tight binding model for dissociable water

We report results of development of a self consistent tight binding model for water. The model explicitly describes the electrons of the liquid self consistently, allows dissociation of the water and permits fast direct dynamics molecular dynamics calculations of the fluid properties. It is parameterized by fitting to first principles calculations on water monomers, dimers, and trimers. We report calculated radial distribution functions of the bulk liquid, a phase diagram and structure of solvated protons within the model as well as ac conductivity of a system of 96 water molecules of which one is dissociated. Structural properties and the phase diagram are in good agreement with experiment and first principles calculations. The estimated DC conductivity of a computational sample containing a dissociated water molecule was an order of magnitude larger than that reported from experiment though the calculated ratio of proton to hydroxyl contributions to the conductivity is very close to the experimental value. The conductivity results suggest a Grotthuss-like mechanism for the proton component of the conductivity.     

Modeling the interaction between two- and four-ring progestin models and a silicone-based polymer model: a density functional theory study

In this paper, we introduce a relatively fast and reliable method for determining the feasibility of drug delivery from transdermal and implant materials. We are using density functional theory for modeling the interaction of progestins, that is, progesterone and six of its hydroxyl derivatives, with a silicone-based polymer. The silicone-based polymer model is a linear molecule, which consists of four dimethylsiloxane units. The progestin models are (1) complete progestin structures, which are called four-ring models, and (2) their two-ring models, which are comprised of the C and D rings of the basic steroid skeletons. We are investigating the interaction between the four- and two-ring models and the polymer model in three different interaction configurations. Altogether, 42 different equilibrium geometries of progestin-polymer model complexes and the corresponding interaction energies have been calculated. Our computational results are in very good agreement with the experimental findings reported previously in the literature, which state that the release rates and permeabilities of progestin pharmaceuticals in silicone-based drug delivery systems decrease when the number of hydroxyl groups is increased in the steroid skeleton. The four-ring models take the total interaction of the steroid into account slightly better than the two-ring models. However, the two-ring models are very good for predicting the local interactions between the steroid and the polymer model.     

Nonlocal viscosity of polymer melts approaching their glassy state

The nonlocal viscosity kernels of polymer melts have been determined by means of equilibrium molecular dynamics upon cooling toward the glass transition. Previous results for the temperature dependence of the self-diffusion coefficient and the value of the glass transition temperature are confirmed. We find that it is essential to include the attractive part of the interatomic potential in order to observe a strong glass transition. The width of the reciprocal space kernel decreases dramatically near the glass transition, being described by a deltalike function near and below the glass transition, leading to a very broad kernel in physical space. Thus, spatial nonlocality turns out to play an important role in polymeric fluids at temperatures near the glass transition temperature.     

分辨伏安分析重叠峰的研究

研究了一种处理伏安(极谱)重叠峰的数学模型。将几类具峰状的极谱电流公式归纳成一般的关系式,提出了一个通用的拟合函数,经非线性最小二乘法处理,可得到重叠组份的蜂高、峰电位和半峰宽等参数。本法适用于示差脉冲极谱、交流极谱、方波极谱、一阶导数卷积伏安法及其反向溶出伏安重叠峰的分离。已用于示差脉冲极谱和交流极谱重叠峰的分辨,得到满意结果。     

Single molecule probing of the glass transition phenomenon: simulations of several types of probes

Molecular dynamics simulations of a system of short bead-spring chains containing an additional dumbbell are presented and analyzed. This system represents a coarse-grained model for a melt of short, flexible polymers containing fluorescent probe molecules at very dilute concentration. It is shown that such a system is very well suited to study aspects of the glass transition of the undercooled polymer melt via single molecule spectroscopy, which are not easily accessed by other methods. Such aspects include data which can be extracted from a study of fluctuations along a trajectory of the single molecule, probing the rugged energy landscape of the glass-forming liquid and transitions from one metabasin of this energy landscape to the next one. Such an information can be inferred from "distance maps" constructed from trajectories characterizing the translational and orientational motion of the probe. At the same time, determining autocorrelation functions along such trajectories, it is shown for several types of probes (differing in their size and/or mass within reasonable limits) that this time-averaged information of the probe is fully compatible with ensemble averaged information on the relaxation of the glass-forming matrix, accessible from bulk measurements. The analyzed quantities include the fluorescence lifetime, linear dichroism, and also various orientational correlation functions of the probe, in order to provide guidance to experimental work. Similar to earlier findings from simulations of bulk molecular fluids, deviations from the Stokes-Einstein and Stokes-Einstein-Debye relations are observed.     

Modeling the growth processes of polyelectrolyte multilayers using a quartz crystal resonator

The layer-by-layer buildup of chitosan/hyaluronan (CH/HA) and poly(l-lysine)/hyaluronan (PLL/HA) multilayers was followed on a quartz crystal resonator (QCR) in different ionic strengths and at different temperatures. These polyelectrolytes were chosen to demonstrate the method whereby useful information is retrieved from acoustically thick polymer layers during their buildup. Surface acoustic impedance recorded in these measurements gives a single or double spiral when plotted in the complex plane. The shape of this spiral depends on the viscoelasticity of the layer material and regularity of the growth process. The polymer layer is assumed to consist of one or two zones. A mathematical model was devised to represent the separation of the layer to two zones with different viscoelastic properties. Viscoelastic quantities of the layer material and the mode and parameters of the growth process were acquired by fitting a spiral to the experimental data. In all the cases the growth process was mainly exponential as a function of deposition cycles, the growth exponent being between 0.250 and 0.275.     

Novel strategy for the hyperelastic parameter fitting procedure of polymer foam materials

The most widely used approach to model the large strain elastic response of polymer foams in a finite element (FE) simulation is the use of the Ogden–Hill compressible hyperelastic material model. This model is implemented and termed as “hyperfoam” material model in the commercial FE software ABAQUS. The hyperfoam model is able to characterize the large compressibility (in volumetric sense) of the foam material. In order to find the material parameters of the model for a particular foam specimen, we need to fit the simulated responses to the available experimental data. This task is easier for incompressible hyperelastic materials because we can use the incompressibility constraint to eliminate the transverse stretch from the stress solutions. However, this simplification cannot be used for the hyperfoam model, therefore, in the stress-strain relations, the transverse stretch is included, which makes the parameter fitting procedure more complicated. In this paper, a novel strategy is proposed for the parameter fitting task. The performance of the new algorithm is demonstrated by presenting fitted material responses for a particular polymer foam material. The major advantage of the new strategy is that it can be used with any third-party optimization solver and there is no need to write our own code.     

Relations between Glass Transition Temperatures in Miscible Polymer Blends and Composition: From Volume to Mass Fractions

The use of volume fractions in the empirical mixing laws to predict the glass transition temperatures (Tg) of polymer blends provides good agreement with experimental values, even for polymer systems with different densities. No adjustment parameter is therefore required whereas Gordon-Taylor and Kwei equations based on weight fractions need the use of a fitting parameter which has to be determined from experimental data. This assumption was validated from Tg measurements through DSC experiments conducted on PMMA /PVDF blends which have significantly different densities.