• Computers and Concrete
• /
• 제10권3호
• /
• pp.307-330
• /
• 2012

This work intends to contribute for the improvement of the procedure suggested by Brazilian Technical Code that takes into account the cracked concrete stiffness in the estimative of the displacement of reinforced concrete beams submitted to service loads. A damage constitutive model accounting for induced anisotropy, plastic deformations and bimodular elastic response is used in order to simulate the concrete behaviour, while an elastoplastic behaviour is admitted for the reinforcement. The constitutive models were implemented in a program for bars structures analysis with layered finite elements. Initially, the damage model is briefly presented as well as the parametric identification of the materials that have been used in the reinforced concrete beams. After that, beams with different geometries and reinforcement area are analyzed and a statistical method (ANOVA) is employed in order to identify the main variables in the problem. Soon after, the same procedure is used with another resistance of concrete, where the compression strength is changed. The numerical responses are compared with the ones obtained by Brazilian Technical Code and experimental tests in order to validate the use of the damage model. Finally, some remarks are discussed based on responses presented in this work.

• Computers and Concrete
• /
• 제15권5호
• /
• pp.735-758
• /
• 2015

This work deals with unilateral effect of quasi-brittle materials, such as concrete. For this propose, a two-dimensional meso-scale model is presented. The material is considered as a three-phase material consisting of interface zone, matrix and inclusions - each constituent modeled by an appropriate constitutive model. The Representative Volume Element (RVE) consists of inclusions idealized as circular shapes randomly placed into the specimen. The interface zone is modeled by means of cohesive contact finite elements developed here in order to capture the effects of phase debonding and interface crack closure/opening. As an initial approximation, the inclusion is modeled as linear elastic as well as the matrix. Our main goal here is to show a computational homogenization-based approach as an alternative to complex macroscopic constitutive models for the mechanical behavior of the quasi-brittle materials using a finite element procedure within a purely kinematical multi-scale framework. A set of numerical examples, involving the microcracking processes, is provided. It illustrates the performance of the proposed model. In summary, the proposed homogenization-based model is found to be a suitable tool for the identification of macroscopic mechanical behavior of quasi-brittle materials dealing with unilateral effect.

• 전산구조공학
• /
• 제6권1호
• /
• pp.91-98
• /
• 1993

굴삭기 붐의 변위와 응력 제한조건에 대해 각 판의 두께를 설계변수로 선정하여 자중을 최소화 하였다. 변위와 응력 제한조건식을 구성하기 위해 붐을 판 복합 구조물로 모델링한 후 3절점 삼각형 판요소의 유한요소 해석 프로그램으로 해석하여 상용 구조해석 프로그램인 ANSYS결과와 비교 검토하였다. 유한요소 해석 결과로부터 구성된 변위와 응력 제한조건을 직접 미분법으로 민감도를 해석하고 차분결과를 기준으로 검토하였다. 최종적으로 민감도 해석 프로그램과 최적화 알고리즘을 결합하여 판의 최적설계 프로그램을 구성하고 균일응력의 외팔보 문제로써 해외 정확도를 검증하였다. 굴삭기 붐에 대한 자중 최소화를 수행한 결과, 붐의 반쪽 모델의 초기무게가 453kgf이었고 최적설계 결과가 331kgf로서 약 27%의 자중 감소 효과를 가져왔다.

• Journal of Mechanical Science and Technology
• /
• 제15권4호
• /
• pp.413-421
• /
• 2001

A simple procedure to add and remove material simultaneously along the boundary is developed to optimize the shape of a two dimensional elastic problems and to minimize the maximum von Mises stress. The results for the two dimensional infinite plate with a hole, are close to the theoretical results of an elliptical boundary and the stress concentration is reduced by half for the fillet problem. The proposed shape optimization method, when compared with existing derivative based shape optimization methods has many features such as simplicity, applicability, flexibility, computational efficiency and a much better control on stresses on the design boundary.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 춘계 학술대회논문집
• /
• pp.83-90
• /
• 2004

Substructuring methods are usually used in finite element structural analyses. In this study a multi-level substructuring algorithm is developed and proposed as a possible candidate for incompressible fluid solves. Finite element formulation for incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et.al.[5]. The present algorithm consists of four stages such as a gathering stage, a condensing stage, a solving stage and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At highest level, each subdomain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each subdomain has been replaced by a sequential static condensation. The global algebraic system arising feom the assembly of each subdomains is solved using Conjugate Gradient Squared(CGS) method. In this case, pre-conditioning techniques usually accompanied by iterative solvers are not needed.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1999년도 봄 학술발표회 논문집
• /
• pp.121-128
• /
• 1999

This paper presents a time domain method for soil-structure interaction analysis for seismic loadings. It is based on the finite element formulation incorporating analytical frequency-dependent infinite elements for the far field soil. The dynamic stiffness matrices of the far field region formulated using the present method in frequency domain can be easily transformed into the corresponding matrices in time domain. At first, the equivalent earthquake forces are evaluated along the interface between the near and the far fields from the free-field response analysis carried out in frequency domain, and the results are transformed into the time domain. An efficient procedure is developed for the convolution integrals to evaluate the interaction force along the interface, which depends on the response on the interface at the past time instances as well as the concurrent instance. Then, the dynamic responses are obtained for the equivalent earthquake force and the interaction force using Newmark direct integration technique. Since the response analysis is carried out in time domain, it can be easily extended to the nonlinear analysis. Example analysis has been carried out to verify the present method in a multi-layered half-space.

• 한국전산구조공학회논문집
• /
• 제22권5호
• /
• pp.475-485
• /
• 2009

본 논문에서는 철근콘크리트 구조물의 3차원 디지털기술 기반 설계에 필요한 CAD/CAE 데이터 연동을 위하여 3D CAD 데이터로부터 철근 요소를 생성하는 2단계 알고리즘을 제안하였다. 1단계 과정은 STEP 표준형식의 3D CAD 데이터를 분석하여 여러 형태의 철근 객체 정보를 추출하는 알고리즘이며, 2단계 과정은 추출된 객체 형식의 철근 정보로부터 유한요소해석에 필요한 1차원 철근요소를 생성하는 알고리즘이다. 개발된 알고리즘을 구현하여 프로그램을 작성하였으며, 이를 사용하여 실제 3D CAD 소프트웨어로부터 생성된 철근콘크리트 구조물 데이터를 유한요소해석 데이터로 성공적으로 변환하여 제안된 알고리즘이 유효함을 보였다.

• Advances in aircraft and spacecraft science
• /
• 제1권2호
• /
• pp.125-143
• /
• 2014

This paper provides a new technique for solving the static analysis of arbitrarily shaped composite plates by using Strong Formulation Finite Element Method (SFEM). Several papers in literature by the authors have presented the proposed technique as an extension of the classic Generalized Differential Quadrature (GDQ) procedure. The present methodology joins the high accuracy of the strong formulation with the versatility of the well-known Finite Element Method (FEM). The continuity conditions among the elements is carried out by the compatibility or continuity conditions. The mapping technique is used to transform both the governing differential equations and the compatibility conditions between two adjacent sub-domains into the regular master element in the computational space. The numerical implementation of the global algebraic system obtained by the technique at issue is easy and straightforward. The main novelty of this paper is the application of the stress and strain recovery once the displacement parameters are evaluated. Computer investigations concerning a large number of composite plates have been carried out. SFEM results are compared with those presented in literature and a perfect agreement is observed.

• International Journal of Aeronautical and Space Sciences
• /
• 제14권4호
• /
• pp.310-323
• /
• 2013

In this paper, the aeroelastic static response of flexible wings with arbitrary cross-section geometry via a coupled CUF-XFLR5 approach is presented. Refined structural one-dimensional (1D) models, with a variable order of expansion for the displacement field, are developed on the basis of the Carrera Unified Formulation (CUF), taking into account cross-sectional deformability. A three-dimensional (3D) Panel Method is employed for the aerodynamic analysis, providing more accuracy with respect to the Vortex Lattice Method (VLM). A straight wing with an airfoil cross-section is modeled as a clamped beam, by means of the finite element method (FEM). Numerical results present the variation of wing aerodynamic parameters, and the equilibrium aeroelastic response is evaluated in terms of displacements and in-plane cross-section deformation. Aeroelastic coupled analyses are based on an iterative procedure, as well as a linear coupling approach for different free stream velocities. A convergent trend of displacements and aerodynamic coefficients is achieved as the structural model accuracy increases. Comparisons with 3D finite element solutions prove that an accurate description of the in-plane cross-section deformation is provided by the proposed 1D CUF model, through a significant reduction in computational cost.

• 전산구조공학
• /
• 제9권1호
• /
• pp.151-159
• /
• 1996

구조용 콘크리트의 비선형 거동을 예측하기 위하여, 압축강도 연화현상, 거시적 및 회전균열모델등의 내용을 포함하고 있는 압축장 응력장 이론(CFT)에 근거한 유한요소법이 개발/제시되었다. 또한, 이와 관련하여 CFT가 암시하는 탄젠트 및 세칸트 재료강성이 반복계산해법의 관점에서 정의/논의되었다. 최종적으로 계산상의 효율성 증대 및 최대하중 이후의 거동 포착에 주안점을 두어 초기재료 강성을 채택한 변위증분법 논리 및 빠른 수렴을 위한 Over-Relaxtion방법이 Isoparametric계의 8-Node요소에 포함/유도되었다. 이와 같이하여 제시된 비선형 해석 프로그램 NASCOM은 응력 혼돈지역에 위치하는 콘크리트 평면요소의 하중 지지능력, 탄성범위 이후의 변형 특성, 균열양상 및 보강근의 항복범위등의 예측을 가능하게 하였다. NASCOM의 제한된 검증을 위하여, Cervenka의 판넬 시험결과에 대한 하중지지능력 및 변형이력등을 예측한 결과가 전체적인 의미에서 실험결과와 상응하는 일치를 나타내었다.