In this paper, a finite element analysis based on the local approach concept to fracture in the continuum damage mechanics is performed to analyze ductile fracture in two dimensional quasi-static state. First an isotropic damage model based on the generalized concept of effective stress is proposed for structural materials in the context of large deformation. In this model, the stiffness degradation is taken as a measure of damage and so, the fracture phenomenon can be explained as the critical deterioration of stiffness at a material point. The modified Riks' continuation technique is used to solve incremental iterative equations. Crack propagation is achieved by removing critically damaged elements. The mesh size sensitivity analysis and the simulation of the well known shearing mode failure in plane strain state are carried out to verify the present formulation. As numerical examples, an edge cracked plate and the specimen with a circular hole under plane stress are taken. Load-displacement curves and successively fractured shapes are shown. From the results, it can be concluded that the proposed model based on the local approach concept in the continuum damage mechanics may be stated as a reasonable tool to explain ductile fracture initiation and crack propagation.
보강된 판 및 쉘구조의 안정성 및 후좌굴을 포함하는 기하학적 비선형 해석을 수행하기 위하여, total Lagrangian formulation에 근거한 연속체의 증분평형방정식으로부터 변형된 쉘요소인 유한요소이론을 제시하였다. 쉘구조의 곡률이 불연속적으로 변하거나 쉘부재들이 유한한 각도로 만나는 보강된 판 및 쉘구조의 비선형 해석이 가능하도록 주부재와 보강재 간의 연결점에 대한 일반적인 변환관계를 제시하였으며 좌굴해석 및 기하학적 비선형해석의 경우에 해의 정확성 및 수렴성을 개선시키기 위하여 접선강도행렬 산정시 회전각의 2차항을 포함시켰다. 또한, shear locking 현상을 극복하기 위하여 감차적분을 적용하였고 쉘구조의 좌굴해석에서는 power method를 적용하여 해석의 효율을 높였으며, 후좌굴해석에서는 변위 및 하중증분법을 적절히 결합시켜 보강된 쉘구조의 후좌굴 거동추적이 용이하였다. 또한, 입력자료를 손쉽게 준비하고 좌굴모드 및 후좌굴거동을 효율적으로 분석하기 위하여 전, 후 처리 프로그램을 개발하였고 다양한 해석예제를 통하여 다른 문헌의 해석결과를 비교함으로써 본 연구에서 개발된 유한요소 해석프로그램의 타당성 및 정확성을 입증하였다.
Valido, Anibal J.;Sousa, Luis G.;Cardoso, J. Barradas
Structural Engineering and Mechanics
/
v.4
no.1
/
pp.25-35
/
1996
This paper describes a continuum variational formulation for design optimization of nonlinear structures in the elastic-plastic domain, where unloading and reloading of the structures are allowed to occur. The Total Lagrangian procedure is used for the description of the structural deformation. The direct differentiation approach is used to derive the sensitivities of the various structural response measures with respect to the design parameters. Since the material goes into the inelastic range and unloading and reloading of the structure are allowed to occur, the structural response is path dependent and an additional step is needed to integrate the constitutive equations. It can be shown, consequently, that design sensitivity analysis is also path-dependent. The theory has been discretized by the finite element technique and implemented in a structural analysis code. Mathematical programming approach is used for the optimization process. Numerical applications on trusses are performed, where cross-sectional areas and nodal point coordinates are treated as design variables. Optimal designs have been obtained and compared by using two different strategies: a two level strategy where the levels are defined accordingly the type of design variables, cross sectional areas or node coordinates, and optimizing simultaneously with respect to both types of design variables.
Transactions of the Korean Society of Mechanical Engineers A
/
v.30
no.12
s.255
/
pp.1518-1525
/
2006
Finite element analyses of structures made of the fiber reinforced composites require an adequate method to characterize the high anisotropic behavior induced by one or several layers of fiber cords with different spatial orientation embedded in a rubber matrix. This paper newly proposes a continuum based rebar element considering change of the orientation of the fiber during deformation of the composite. The mechanical behavior of the embedded fiber is modeled using two-node bar elements in order to consider the relative deformation and spatial orientation of the embedded fiber. For improvement of the analysis accuracy, the load-displacement curve of fiber is applied to the stiffness matrix of fiber. A finite element program is constructed based on the total Lagrangian formulation considering both geometric and material nonlinearity. Finite element analyses of the tensile test are carried out in order to evaluate the validity of the proposed method. Analysis results obtained with the proposed method provides realistic representation of the fiber reinforced rubber composite compared to results of other two models by the Halpin-Tsai equation and a rebar element in ABAQUS/Standard.
Journal of the Korean Society for Precision Engineering
/
v.15
no.9
/
pp.117-126
/
1998
For investigating rubber pad sheet metal forming process, the rubber pad deformation characteristics as well as the contact problem of rubber pad-sheet metal has been analyzed. In this paper, the behavior of the rubber deformation is represented by hyper-elastic constitutive relations based on a generalized Mooney-Rivlin model. Finite element procedures for the two-dimensional responses, employing total Lagrangian formulations are implemented in an implicit form. The volumetric incompressibility condition of the rubber deformation is included in the formulation by using penalty method. The sheet metal is characterized by elasto-plastic material with strain hardening effect and analyzed by a commercial code. The contact procedure and interface program between rubber pad and sheet metal are implemented. Inflation experiment of circular rubber pad identifies the behaviour of the rubber pad deformation during the process. The various form dies and scaled down apparatus of the rubber-pad forming process are fabricated for simulating realistic forming process. The obtaining experimental data and FEM solutions were compared. The numerical solutions illustrate fair agreement with experimental results. The forming pressure distribution according to the dimensions of sheet metal and rubber pads, various rubber models and rubber material are also compared and discussed.
This paper extends the use of the hierarchic degenerated shell element to geometric non-linear analysis of composite laminated skew plates by the p-version of the finite element method. For the geometric non-linear analysis, the total Lagrangian formulation is adopted with moderately large displacement and small strain being accounted for in the sense of von Karman hypothesis. The present model is based on equivalent-single layer laminate theory with the first order shear deformation including a shear correction factor of 5/6. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. A wide variety of linear and non-linear results obtained by the p-version finite element model are presented for the laminated skew plates as well as laminated square plates. A numerical analysis is made to illustrate the influence of the geometric non-linear effect on the transverse deflections and the stresses with respect to width/depth ratio (a/h), skew angle (${\beta}$), and stacking sequence of layers. The present results are in good agreement with the results in literatures.
This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.
Rad, Mohammad Hossein Ghadiri;Shahabian, Farzad;Hosseini, Seyed Mahmoud
Steel and Composite Structures
/
v.35
no.1
/
pp.77-92
/
2020
The present paper outlined a procedure for geometrically nonlinear dynamic analysis of functionally graded graphene platelets-reinforced (GPLR-FG) nanocomposite cylinder subjected to mechanical shock loading. The governing equation of motion for large deformation problems is derived using meshless local Petrov-Galerkin (MLPG) method based on total lagrangian approach. In the MLPG method, the radial point interpolation technique is employed to construct the shape functions. A micromechanical model based on the Halpin-Tsai model and rule of mixture is used for formulation the nonlinear functionally graded distribution of GPLs in polymer matrix of composites. Energy dissipation in analyses of the structure responding to dynamic loads is considered using the Rayleigh damping. The Newmark-Newton/Raphson method which is an incremental-iterative approach is implemented to solve the nonlinear dynamic equations. The results of the proposed method for homogenous material are compared with the finite element ones. A very good agreement is achieved between the MLPG and FEM with very fine meshing. In addition, the results have demonstrated that the MLPG method is more effective method compared with the FEM for very large deformation problems due to avoiding mesh distortion issues. Finally, the effect of GPLs distribution on strength, stiffness and dynamic characteristics of the cylinder are discussed in details. The obtained results show that the distribution of GPLs changed the mechanical properties, so a classification of different types and volume fraction exponent is established. Indeed by comparing the obtained results, the best compromise of nanocomposite cylinder is determined in terms of mechanical and dynamic properties for different load patterns. All these applications have shown that the present MLPG method is very effective for geometrically nonlinear analyses of GPLR-FG nanocomposite cylinder because of vanishing mesh distortion issue in large deformation problems. In addition, since in proposed method the distributed nodes are used for discretization the problem domain (rather than the meshing), modeling the functionally graded media yields to more accurate results.
Magazine of the Korean Society of Agricultural Engineers
/
v.42
no.5
/
pp.151-159
/
2000
Dynamic loading of structures often causes excursions of stresses will into the inelastic range and the influence of geometry changes on the response is also significant in may cases. In general , the shell structures designed according to quasi-Static analysis may collapse under condition of dynamic loading. Therefore, for a more realistic prediction on the lad carrying capacity of these shell. both material and geometric nonlinear effects should be considered. In this study , the material nonlinearity effect on the dynamic response is formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a Total Lagrangian formulation. the reinforcing bars are modeled by the equivalent steel layer at the location of reinforcements, and Von Mises yield criteria is adopted for the steel layer behavior. Also, Drucker-Prager yield criteria is applied for the behavior of concrete. the shape imperfection of dome is assumed as 'dimple type' which can be expressed Wd1=Wd0(1-(r-a)m)n while the shape imperfection of wall is assumed as sinusoidal curve which is Wwi =Wwo sin(n $\pi$y/l). In numerical test, three cases of shape imperfection of 0.0 -5.0cm(opposite direction to loading ; inner shape imperfection)and 5cm (direction to loading : outward shape imperfection) and thickness of steel layer determined by steel ratio of 0,3, and 5% were analyzed. The effect of shape imperfection and steel ratio and behavior characteristics of perfect shape shell and imperfect shape shell are identified through analysis of above mentioned numerical test. Dynamic behaviors of dome and wall according toe combination of shape imperfection and steel ratio are also discussed in this paper.
Magazine of the Korean Society of Agricultural Engineers
/
v.40
no.5
/
pp.91-99
/
1998
The widespread use of thin shell structures has created a need for a systematic method of analysis which can adequately account for arbitrary geometric form and boundary conditions as well as arbitrary general type of loading. Therefore, the stress and analysis of thin shell has been one of the more challenging areas of structural mechanics. A wide variety of numerical methods have been applied to the governing differential equations for spherical and cylindrical structures with a few results applicable to practice. The analysis of axisymmetric spherical shell is almost an every day occurrence in many industrial applications. A reliable and accurate finite element analysis procedure for such structures was needed. Dynamic loading of structures often causes excursions of stresses well into the inelastic range and the influence of geometry changes on the response is also significant in many cases. Therefore both material and geometric nonlinear effects should be considered. In general, the shell structures designed according to quasi-static analysis may fail under conditions of dynamic loading. For a more realistic prediction on the load carrying capacity of these shell, in addition to the dynamic effect, consideration should also include other factors such as nonlinearities in both material and geometry since these factors, in different manner, may also affect the magnitude of this capacity. The objective of this paper is to demonstrate the dynamic characteristics of spherical shell. For these purposes, the spherical shell subjected to uniformly distributed step load was analyzed for its large displacements elasto-viscoplastic static and dynamic response. Geometrically nonlinear behaviour is taken into account using a Total Lagrangian formulation and the material behaviour is assumed to elasto-viscoplastic model highly corresponding to the real behaviour of the material. The results for the dynamic characteristics of spherical shell in the cases under various conditions of base-radius/central height(a/H) and thickness/shell radius(t/R) were summarized as follows : The dynamic characteristics with a/H. 1) AS the a/H increases, the amplitude of displacement in creased. 2) The values of displacement dynamic magnification factor (DMF) were ranges from 2.9 to 6.3 in the crown of shell and the values of factor in the mid-point of shell were ranged from 1.8 to 2.6. 3) As the a/H increases, the values of DMF in the crown of shell is decreased rapidly but the values of DMF in mid-point shell is increased gradually. 4) The values of DMF of hoop-stresses were range from 3.6 to 6.8 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.3 to 2.6, and the values of DMF of stress were larger than that of displacement. The dynamic characteristics with t/R. 5) With the thickness of shell decreases, the amplitude of the displacement and the period increased. 6) The values of DMF of the displacement were ranged from 2.8 to 3.6 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.1 to 2.2.
본 웹사이트에 게시된 이메일 주소가 전자우편 수집 프로그램이나
그 밖의 기술적 장치를 이용하여 무단으로 수집되는 것을 거부하며,
이를 위반시 정보통신망법에 의해 형사 처벌됨을 유념하시기 바랍니다.
[게시일 2004년 10월 1일]
이용약관
제 1 장 총칙
제 1 조 (목적)
이 이용약관은 KoreaScience 홈페이지(이하 “당 사이트”)에서 제공하는 인터넷 서비스(이하 '서비스')의 가입조건 및 이용에 관한 제반 사항과 기타 필요한 사항을 구체적으로 규정함을 목적으로 합니다.
제 2 조 (용어의 정의)
① "이용자"라 함은 당 사이트에 접속하여 이 약관에 따라 당 사이트가 제공하는 서비스를 받는 회원 및 비회원을
말합니다.
② "회원"이라 함은 서비스를 이용하기 위하여 당 사이트에 개인정보를 제공하여 아이디(ID)와 비밀번호를 부여
받은 자를 말합니다.
③ "회원 아이디(ID)"라 함은 회원의 식별 및 서비스 이용을 위하여 자신이 선정한 문자 및 숫자의 조합을
말합니다.
④ "비밀번호(패스워드)"라 함은 회원이 자신의 비밀보호를 위하여 선정한 문자 및 숫자의 조합을 말합니다.
제 3 조 (이용약관의 효력 및 변경)
① 이 약관은 당 사이트에 게시하거나 기타의 방법으로 회원에게 공지함으로써 효력이 발생합니다.
② 당 사이트는 이 약관을 개정할 경우에 적용일자 및 개정사유를 명시하여 현행 약관과 함께 당 사이트의
초기화면에 그 적용일자 7일 이전부터 적용일자 전일까지 공지합니다. 다만, 회원에게 불리하게 약관내용을
변경하는 경우에는 최소한 30일 이상의 사전 유예기간을 두고 공지합니다. 이 경우 당 사이트는 개정 전
내용과 개정 후 내용을 명확하게 비교하여 이용자가 알기 쉽도록 표시합니다.
제 4 조(약관 외 준칙)
① 이 약관은 당 사이트가 제공하는 서비스에 관한 이용안내와 함께 적용됩니다.
② 이 약관에 명시되지 아니한 사항은 관계법령의 규정이 적용됩니다.
제 2 장 이용계약의 체결
제 5 조 (이용계약의 성립 등)
① 이용계약은 이용고객이 당 사이트가 정한 약관에 「동의합니다」를 선택하고, 당 사이트가 정한
온라인신청양식을 작성하여 서비스 이용을 신청한 후, 당 사이트가 이를 승낙함으로써 성립합니다.
② 제1항의 승낙은 당 사이트가 제공하는 과학기술정보검색, 맞춤정보, 서지정보 등 다른 서비스의 이용승낙을
포함합니다.
제 6 조 (회원가입)
서비스를 이용하고자 하는 고객은 당 사이트에서 정한 회원가입양식에 개인정보를 기재하여 가입을 하여야 합니다.
제 7 조 (개인정보의 보호 및 사용)
당 사이트는 관계법령이 정하는 바에 따라 회원 등록정보를 포함한 회원의 개인정보를 보호하기 위해 노력합니다. 회원 개인정보의 보호 및 사용에 대해서는 관련법령 및 당 사이트의 개인정보 보호정책이 적용됩니다.
제 8 조 (이용 신청의 승낙과 제한)
① 당 사이트는 제6조의 규정에 의한 이용신청고객에 대하여 서비스 이용을 승낙합니다.
② 당 사이트는 아래사항에 해당하는 경우에 대해서 승낙하지 아니 합니다.
- 이용계약 신청서의 내용을 허위로 기재한 경우
- 기타 규정한 제반사항을 위반하며 신청하는 경우
제 9 조 (회원 ID 부여 및 변경 등)
① 당 사이트는 이용고객에 대하여 약관에 정하는 바에 따라 자신이 선정한 회원 ID를 부여합니다.
② 회원 ID는 원칙적으로 변경이 불가하며 부득이한 사유로 인하여 변경 하고자 하는 경우에는 해당 ID를
해지하고 재가입해야 합니다.
③ 기타 회원 개인정보 관리 및 변경 등에 관한 사항은 서비스별 안내에 정하는 바에 의합니다.
제 3 장 계약 당사자의 의무
제 10 조 (KISTI의 의무)
① 당 사이트는 이용고객이 희망한 서비스 제공 개시일에 특별한 사정이 없는 한 서비스를 이용할 수 있도록
하여야 합니다.
② 당 사이트는 개인정보 보호를 위해 보안시스템을 구축하며 개인정보 보호정책을 공시하고 준수합니다.
③ 당 사이트는 회원으로부터 제기되는 의견이나 불만이 정당하다고 객관적으로 인정될 경우에는 적절한 절차를
거쳐 즉시 처리하여야 합니다. 다만, 즉시 처리가 곤란한 경우는 회원에게 그 사유와 처리일정을 통보하여야
합니다.
제 11 조 (회원의 의무)
① 이용자는 회원가입 신청 또는 회원정보 변경 시 실명으로 모든 사항을 사실에 근거하여 작성하여야 하며,
허위 또는 타인의 정보를 등록할 경우 일체의 권리를 주장할 수 없습니다.
② 당 사이트가 관계법령 및 개인정보 보호정책에 의거하여 그 책임을 지는 경우를 제외하고 회원에게 부여된
ID의 비밀번호 관리소홀, 부정사용에 의하여 발생하는 모든 결과에 대한 책임은 회원에게 있습니다.
③ 회원은 당 사이트 및 제 3자의 지적 재산권을 침해해서는 안 됩니다.
제 4 장 서비스의 이용
제 12 조 (서비스 이용 시간)
① 서비스 이용은 당 사이트의 업무상 또는 기술상 특별한 지장이 없는 한 연중무휴, 1일 24시간 운영을
원칙으로 합니다. 단, 당 사이트는 시스템 정기점검, 증설 및 교체를 위해 당 사이트가 정한 날이나 시간에
서비스를 일시 중단할 수 있으며, 예정되어 있는 작업으로 인한 서비스 일시중단은 당 사이트 홈페이지를
통해 사전에 공지합니다.
② 당 사이트는 서비스를 특정범위로 분할하여 각 범위별로 이용가능시간을 별도로 지정할 수 있습니다. 다만
이 경우 그 내용을 공지합니다.
제 13 조 (홈페이지 저작권)
① NDSL에서 제공하는 모든 저작물의 저작권은 원저작자에게 있으며, KISTI는 복제/배포/전송권을 확보하고
있습니다.
② NDSL에서 제공하는 콘텐츠를 상업적 및 기타 영리목적으로 복제/배포/전송할 경우 사전에 KISTI의 허락을
받아야 합니다.
③ NDSL에서 제공하는 콘텐츠를 보도, 비평, 교육, 연구 등을 위하여 정당한 범위 안에서 공정한 관행에
합치되게 인용할 수 있습니다.
④ NDSL에서 제공하는 콘텐츠를 무단 복제, 전송, 배포 기타 저작권법에 위반되는 방법으로 이용할 경우
저작권법 제136조에 따라 5년 이하의 징역 또는 5천만 원 이하의 벌금에 처해질 수 있습니다.
제 14 조 (유료서비스)
① 당 사이트 및 협력기관이 정한 유료서비스(원문복사 등)는 별도로 정해진 바에 따르며, 변경사항은 시행 전에
당 사이트 홈페이지를 통하여 회원에게 공지합니다.
② 유료서비스를 이용하려는 회원은 정해진 요금체계에 따라 요금을 납부해야 합니다.
제 5 장 계약 해지 및 이용 제한
제 15 조 (계약 해지)
회원이 이용계약을 해지하고자 하는 때에는 [가입해지] 메뉴를 이용해 직접 해지해야 합니다.
제 16 조 (서비스 이용제한)
① 당 사이트는 회원이 서비스 이용내용에 있어서 본 약관 제 11조 내용을 위반하거나, 다음 각 호에 해당하는
경우 서비스 이용을 제한할 수 있습니다.
- 2년 이상 서비스를 이용한 적이 없는 경우
- 기타 정상적인 서비스 운영에 방해가 될 경우
② 상기 이용제한 규정에 따라 서비스를 이용하는 회원에게 서비스 이용에 대하여 별도 공지 없이 서비스 이용의
일시정지, 이용계약 해지 할 수 있습니다.
제 17 조 (전자우편주소 수집 금지)
회원은 전자우편주소 추출기 등을 이용하여 전자우편주소를 수집 또는 제3자에게 제공할 수 없습니다.
제 6 장 손해배상 및 기타사항
제 18 조 (손해배상)
당 사이트는 무료로 제공되는 서비스와 관련하여 회원에게 어떠한 손해가 발생하더라도 당 사이트가 고의 또는 과실로 인한 손해발생을 제외하고는 이에 대하여 책임을 부담하지 아니합니다.
제 19 조 (관할 법원)
서비스 이용으로 발생한 분쟁에 대해 소송이 제기되는 경우 민사 소송법상의 관할 법원에 제기합니다.
[부 칙]
1. (시행일) 이 약관은 2016년 9월 5일부터 적용되며, 종전 약관은 본 약관으로 대체되며, 개정된 약관의 적용일 이전 가입자도 개정된 약관의 적용을 받습니다.