• Interaction and multiscale mechanics
• /
• 제5권3호
• /
• pp.229-265
• /
• 2012

Hygroexpansion of wood is a known and undesired characteristic in civil engineering. When wood is exposed to changing environmental humidity, it adsorbs or desorbs moisture and warps. The resulting distortions or - at restrained conditions - cracks are a major concern in timber engineering. We herein present a multiscale model for prediction of the macroscopic hygroexpansion behavior of individual pieces of softwood from their microstructure, demonstrated for spruce. By applying poromicromechanics, we establish a link between the swelling pressure, driving the hygroexpansion of wood at the nanoscale, and the resulting macroscopic dimensional changes. The model comprises six homogenization steps, which are performed by means of continuum micromechanics, the unit cell method and laminate theory, all formulated in a poromechanical framework. Model predictions for elastic properties of wood as functions of the moisture content closely approach corresponding experimental data. As for the hygroexpansion behavior, the swelling pressure has to be back-calculated from macroscopic hygroexpansion data. The good reproduction of the anisotropy of wood hygroexpansion, based on only a single scalar calibration parameter, underlines the suitability of the model. The multiscale model constitutes a valuable tool for studying the effect of microstructural features on the macroscopic behavior and for assessing the hygroexpansion behavior at smaller length scales, which are inaccessible to experiments. The model predictions deliver input parameters for the analysis of timber at the structural scale, therewith enabling to optimize the use of timber and to prevent moisture-induced damage or failure.

• Smart Structures and Systems
• /
• 제13권4호
• /
• pp.501-515
• /
• 2014

This paper presents a linear computational homogenization framework to evaluate the effective (or generalized) electromechanical coupling coefficient (EMCC) of adaptive structures with piezoelectric structural fiber (PSF) composite elements. The PSF consists of a silicon carbide (SiC) or carbon core fiber as reinforcement to a fragile piezo-ceramic shell. For the micro-scale analysis, a micromechanics model based on the variational asymptotic method for unit cell homogenization (VAMUCH) is used to evaluate the overall electromechanical properties of the PSF composites. At the macro-scale, a finite element (FE) analysis with the commercial FE code ABAQUS is performed to evaluate the effective EMCC for structures with the PSF composite patches. The EMCC is postprocessed from free-vibrations analysis under short-circuit (SC) and open-circuit (OC) electrodes of the patches. This linear two-scale computational framework may be useful for the optimal design of active structure multi-functional composites which can be used for multi-functional applications such as structural health monitoring, power harvest, vibration sensing and control, damping, and shape control through anisotropic actuation.

• Interaction and multiscale mechanics
• /
• 제1권3호
• /
• pp.369-379
• /
• 2008

A combined stochastic diffusion and mean-field model is developed for a systematic study of the grain growth in a pure single-phase polycrystalline material. A corresponding Fokker-Planck continuity equation is formulated, and the interplay/competition of stochastic and curvature-driven mechanisms is investigated. Finite difference results show that the stochastic diffusion coefficient has a strong effect on the growth of small grains in the early stage in both two-dimensional columnar and three-dimensional grain systems, and the corresponding growth exponents are ~0.33 and ~0.25, respectively. With the increase in grain size, the deterministic curvature-driven mechanism becomes dominant and the growth exponent is close to 0.5. The transition ranges between these two mechanisms are about 2-26 and 2-15 nm with boundary energy of 0.01-1 J $m^{-2}$ in two- and three-dimensional systems, respectively. The grain size distribution of a three-dimensional system changes dramatically with increasing time, while it changes a little in a two-dimensional system. The grain size distribution from the combined model is consistent with experimental data available.

• Structural Engineering and Mechanics
• /
• 제76권4호
• /
• pp.529-540
• /
• 2020

Concurrent topology optimization of macrostructure and microstructure has attracted significant interest due to its high structural performance. However, most of the existing works are carried out under deterministic conditions, the obtained design may be vulnerable or even cause catastrophic failure when the load position exists uncertainty. Therefore, it is necessary to take load position uncertainty into consideration in structural design. This paper presents a computational method for robust concurrent topology optimization with consideration of load position uncertainty. The weighted sum of the mean and standard deviation of the structural compliance is defined as the objective function with constraints are imposed to both macro- and micro-scale structure volume fractions. The Bivariate Dimension Reduction method and Gauss-type quadrature (BDRGQ) are used to quantify and propagate load uncertainty to calculate the objective function. The effective properties of microstructure are evaluated by the numerical homogenization method. To release the computation burden, the decoupled sensitivity analysis method is proposed for microscale design variables. The bi-directional evolutionary structural optimization (BESO) method is used to obtain the black-and-white designs. Several 2D and 3D examples are presented to validate the effectiveness of the proposed robust concurrent topology optimization method.

• 한국전산구조공학회논문집
• /
• 제20권3호
• /
• pp.293-301
• /
• 2007

본 연구에서는 이동최소제곱 다절점 유한요소를 이용한 새로운 전역-국부해석기법을 제시하였다. 다절점 유한요소는 요소의 변에 임의의 수 절점을 가질 수 있으므로, 여러 개의 유한요소를 요소망의 재구성 없이 동시적으로 결합시킬 수 있다. 이는 응력구배가 집중되는 곳에 유한요소망을 구성하는 데에 있어 큰 편의를 제공한다. 또한 기존의 전역-국부해석기법처럼 중첩된 요소망을 사용하거나, 지배방정식을 두 번 해석할 필요가 없기 때문에 매우 간편하고 정확하다. 제시된 방법론의 성능을 검증하기 위해 응력 집중과 관련된 다양한 다중스케일 문제를 해석하였다.

• Multiscale and Multiphysics Mechanics
• /
• 제1권2호
• /
• pp.143-156
• /
• 2016

Reinforced concrete (RC) structures require advanced analysis techniques for better estimation of their seismic responses, especially in the case of exhibiting complex three-dimensional coupling of torsional and flexural behaviors. This study focuses on validating a numerical approach for evaluating the seismic response of a three-dimensional unsymmetrical RC structure through the participation in the SMART 2013 international benchmark program. The benchmark program provides material properties, detailed drawings of the RC structure, and input ground motions for the seismic response evaluation. In this study, nonlinear constitutive models of concrete and rebar were formed and local tests were conducted to verify the constitutive models in finite element analysis. Elastic calibration of the finite element model of the SMART 2013 RC structure was performed by comparing numerical and experimental results in modal and linear time history analyses. Using the calibrated model, nonlinear earthquake analysis and seismic fragility analysis were performed to estimate the behavior and vulnerability of the RC structure with various ground motions.

• Structural Engineering and Mechanics
• /
• 제81권3호
• /
• pp.267-280
• /
• 2022

Multiscale structure has attracted significant interest due to its high stiffness/strength to weight ratios and multifunctional performance. However, most of the existing concurrent topology optimization works are carried out under deterministic load conditions. Hence, this paper proposes a robust concurrent topology optimization method based on the bidirectional evolutionary structural optimization (BESO) method for the design of structures composed of periodic microstructures subjected to uncertain dynamic loads. The robust objective function is defined as the weighted sum of the mean and standard deviation of the module of dynamic structural compliance with constraints are imposed to both macro- and microscale structure volume fractions. The polynomial chaos expansion (PCE) method is used to quantify and propagate load uncertainty to evaluate the objective function. The effective properties of microstructure is evaluated by the numerical homogenization method. To release the computation burden, the decoupled sensitivity analysis method is proposed for microscale design variables. The proposed method is a non-intrusive method, and it can be conveniently extended to many topology optimization problems with other distributions. Several numerical examples are used to validate the effectiveness of the proposed robust concurrent topology optimization method.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2006년도 정기 학술대회 논문집
• /
• pp.395-402
• /
• 2006

In analyzing the nano-scale phenomena or behaviors of nano devices or materials, it is often desirable to deal with more atoms than can be treated only with a full atomistic simulation. However, even now, it is advisable to apply the atomistic simulation to the narrow region where the deformation field changes rapidly but to apply the conventional continuum model to the region far from that region. This equivalent continuum model can be formulated by applying the Cauchy-Born rule to the exact atomistic potential as in the quasicontinuum method. To couple the atomistic model with the equivalent continuum model, continuum displacements are conformed to the molecular displacements at the discrete positions of the atoms within the bridging domain. To satisfy the coupling constraints, we apply the Lagrange multiplier method. The continuum model in the bridging model should be applied on the region where the deformation field changes gradually. Then we can make the nodal spacing in the continuum model be much larger than the atomic spacing. In the first step, we generate the atomic-resolution mesh with the nodal spacing equal to the atomic spacing, and then we eliminate the nodal degrees of freedom adaptively using the node deactivation techniques. We eliminate more DOFs as the regions are more far from the atomistic region. Computing time and computational resources can be greatly reduced by the present node deactivation technique in multi scale analysis.

• Interaction and multiscale mechanics
• /
• 제4권4호
• /
• pp.273-289
• /
• 2011

The vibration response of the bridges under the moving vehicular load is of importance for engineers to estimate the serviceability of existing bridges and to design new bridges. This paper deals with the three dimensional vibration analysis of prestressed concrete bridges under moving vehicles. The prestressed bridges are modeled by four-node isoparametric flat shell elements with the transverse shearing deformation taken into account. The usual five degrees-of-freedom (DOFs) per node of the element are appended with a drilling DOF to accommodate the transformation of the local stiffness and mass matrices to the global coordinates. The vehicle is modeled as a single or two-DOF system. A single-span prestressed Tee beam and two-span prestressed box-girder bridge are studied as the two numerical examples. The effects of prestress forces on the natural frequencies and dynamic responses of the bridges are investigated.

• Interaction and multiscale mechanics
• /
• 제4권4호
• /
• pp.257-271
• /
• 2011

The analysis of interfacial stresses in structural component has been the subject of several investigations but it still requires more effort and studies. In this study a general three-dimensional interface element has been formulated for stress and displacement analyses in the interfacial area between two adjacent plate bending element and brick element. Interface element has 16 nodes with 5 degrees of freedom (DOF) in each node adjacent to plate bending element and 3 DOF in each node adjacent to brick element. The interface element has ability to transfer three translations from each side of interface element and two rotations in the side adjacent to the plate element. Stiffness matrix of this element was formulated and implemented in three-dimensional finite element code. Application of this element to the reinforced concrete (RC) beam strengthened with fiber reinforced polymer (FRP) including variation of deflection, slip between plate and concrete, normal and shear stresses distributions in FRP plates have been verified using experimental and numerical work of strengthened RC beams carried out by some researchers. The results show that this interface element is effective and can be used for structural component with these types of interface elements.