共查询到20条相似文献,搜索用时 15 毫秒
1.
D. J. Bernasconi P. M. Richelsen 《International journal for numerical methods in engineering》2008,73(5):706-727
A higher‐order panel method for analysing the three‐dimensional potential flow fields around bodies and wakes is presented. The geometric surfaces are represented by continuous curved patches, with no discretization into panels. These geometric patches hold singularity distributions that have C2 continuity, and which are solved by applying Dirichlet or Neumann boundary conditions at discrete collocation points. While higher‐order methods have previously been developed for thick bodies and Dirichlet boundary conditions, this potential flow method is capable of modelling continuous geometry and singularity surfaces over thin bodies and wakes. The continuous surface method has a number of advantages over conventional constant panel methods. Firstly, as curved geometries are represented exactly, changing the order of the solution does not modify the physical shape of the configuration under investigation. Furthermore, the continuous singularity distributions allow velocities to be evaluated accurately across the entire surface rather than just at collocation points. This means that pressure distributions can be calculated exactly without interpolation, and streamlines can be constructed very close to surfaces without problems of divergence. Finally, body and wake surfaces do not exhibit the strong modelling singularities that can present difficulties with wake relaxation. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
2.
Jie Shen David Yoon 《International journal for numerical methods in engineering》2004,60(15):2441-2466
In this paper we propose a new scheme for freeform shape optimization on arbitrary polygonal or polyhedral meshes. The approach consists of three main steps: (1) surface partitioning of polygonal meshes into different patches; (2) a new freeform perturbation scheme of using the Cox–de Boor basis function over arbitrary polygonal meshes, which supports multi‐resolution shape optimization and does not require CAD information; (3) freeform shape optimization of arbitrary polygonal or polyhedral meshes. Numerical experiments indicate the effectiveness of the proposed approach. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
3.
Endel V. Iarve Mark R. Gurvich David H. Mollenhauer Cheryl A. Rose Carlos G. Dávila 《International journal for numerical methods in engineering》2011,88(8):749-773
The initiation and evolution of transverse matrix cracks and delaminations are predicted within a mesh‐independent cracking (MIC) framework. MIC is a regularized extended finite element method (x‐FEM) that allows the insertion of cracks in directions that are independent of the mesh orientation. The Heaviside step function that is typically used to introduce a displacement discontinuity across a crack surface is replaced by a continuous function approximated by using the original displacement shape functions. Such regularization allows the preservation of the Gaussian integration schema regardless of the enrichment required to model cracking in an arbitrary direction. The interaction between plies is anchored on the integration point distribution, which remains constant through the entire simulation. Initiation and propagation of delaminations between plies as well as intra‐ply MIC opening is implemented by using a mixed‐mode cohesive formulation in a fully three‐dimensional model that includes residual thermal stresses. The validity of the proposed methodology was tested against a variety of problems ranging from simple evolution of delamination from existing transverse cracks to strength predictions of complex laminates withouttextita priori knowledge of damage location or initiation. Good agreement with conventional numerical solutions and/or experimental data was observed in all the problems considered. Published 2011. This article is a US Government work and is in the public domain in the USA. 相似文献
4.
A. E. J. Bogaers S. Kok A. G. Malan 《International journal for numerical methods in engineering》2011,86(8):935-952
In this paper, we implement the method of proper orthogonal decomposition (POD) to generate a reduced order model (ROM) of an optimization‐based mesh movement scheme. In this study it is shown that POD can be used effectively to generate an ROM, that accurately reproduces the full order mesh movement algorithm, with a decrease in computational time of over 99%. We further introduce a novel training procedure whereby the POD models are generated in a fully automated fashion. The technology is applicable to any mesh movement method and enables potential reductions of up to four orders of magnitude in mesh movement related costs. The proposed model can be implemented without having to pre‐train the POD model, to any fluid–structure interaction code with an existing mesh movement scheme. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
5.
Ashok V. Kumar Premdheepak S. Periyasamy 《International journal for numerical methods in engineering》2012,91(5):472-490
Mesh independent analysis is motivated by the desire to use accurate geometric models represented as equations rather than approximated by a mesh. The trial and test functions are approximated or interpolated on a background mesh that is independent of the geometry. This background mesh is easy to generate because it does not have to conform to the geometry. Essential boundary conditions can be applied using the implicit boundary method where the trial and test functions are constructed utilizing approximate step functions such that the boundary conditions are guaranteed to be satisfied. This approach has been demonstrated for two‐dimensional (2D) and three‐dimensional (3D) structural analysis and is extended in this paper to model shell‐like structures. The background mesh consists of 3D elements that use uniform B‐spline approximations, and the shell geometry is assumed to be defined as parametric surfaces to allow arbitrarily complex shell‐like structures to be modeled. Several benchmark problems are used to study the validity of these 3D B‐spline shell elements. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
6.
Jin Qian Yongjie Zhang Wenyan Wang Alexis C. Lewis M. A. Siddiq Qidwai Andrew B. Geltmacher 《International journal for numerical methods in engineering》2010,82(11):1406-1423
This paper describes a novel approach to improve the quality of non‐manifold hexahedral meshes with feature preservation for microstructure materials. In earlier works, we developed an octree‐based isocontouring method to construct unstructured hexahedral meshes for domains with multiple materials by introducing the notion of material change edge to identify the interface between two or more materials. However, quality improvement of non‐manifold hexahedral meshes is still a challenge. In the present algorithm, all the vertices are categorized into seven groups, and then a comprehensive method based on pillowing, geometric flow and optimization techniques is developed for mesh quality improvement. The shrink set in the modified pillowing technique is defined automatically as the boundary of each material region with the exception of local non‐manifolds. In the relaxation‐based smoothing process, non‐manifold points are identified and fixed. Planar boundary curves and interior spatial curves are distinguished, and then regularized using B‐spline interpolation and resampling. Grain boundary surface patches and interior vertices are improved as well. Finally, the optimization method eliminates negative Jacobians of all the vertices. We have applied our algorithms to two beta titanium data sets, and the constructed meshes are validated via a statistics study. Finite element analysis of the 92‐grain titanium is carried out based on the improved mesh, and compared with the direct voxel‐to‐element technique. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
7.
Young‐Sam Cho Seyoung Im 《International journal for numerical methods in engineering》2006,65(4):517-547
Two‐dimensional finite ‘crack’ elements for simulation of propagating cracks are developed using the moving least‐square (MLS) approximation. The mapping from the parental domain to the physical element domain is implicitly obtained from MLS approximation, with the shape functions and their derivatives calculated and saved only at the numerical integration points. The MLS‐based variable‐node elements are extended to construct the crack elements, which allow the discontinuity of crack faces and the crack‐tip singularity. The accuracy of the crack elements is checked by calculating the stress intensity factor under mode I loading. The crack elements turn out to be very efficient and accurate for simulating crack propagations, only with the minimal amount of element adjustment and node addition as the crack tip moves. Numerical results and comparison to the results from other works demonstrate the effectiveness and accuracy of the present scheme for the crack elements. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
8.
Dragan Ribarić Gordan Jelenić 《International journal for numerical methods in engineering》2014,98(7):492-517
We employ the linked interpolation concept to develop two higher‐order nine‐node quadrilateral plate finite elements with curved sides that pass the constant bending patch test for arbitrary node positions. The linked interpolation for the plate displacements is expanded with three bubble parameters to get polynomial completeness necessary to satisfy the patch test. In contrast to some other techniques, the elements developed in this way retain a symmetric stiffness matrix at a marginal computational expense at the element level. The new elements generated using this concept are tested on several examples with curved sides or some other kind of geometric distortion. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
9.
《Advanced Composite Materials》2012,21(3):209-220
The damage initiation and growth near the edges of a circular open hole of laminated composite plates under compressive load were experimentally studied to understand the physical meaning of the open hole compression strength. Damage growth sequence was investigated mainly through optical micrographs of various sections of damaged specimens. When comparatively tough interface laminates were tested, various damages occurred following initial kink band damage in the 0° layers, such as zigzag damage composed of kink bands and transverse cracks in ±45° layers, matrix cracks in 0° layers, and delaminations above and below the 0° layers. High-speed video images showed that unstable growth of the delaminations near the surface was observed prior to the final failure, accompanying buckling of the delaminated layer. The possible sequences of the damage growth were discussed through consideration of the above damage data. 相似文献
10.
11.
David González Elías Cueto Manuel Doblaré 《International journal for numerical methods in engineering》2010,83(6):741-764
We present here a generalization of local maximum entropy (max‐ent) approximation for high orders of consistency (i.e. quadratic, cubic, …). The method is based upon the application of the de Boor's algorithm to the standard, linear local max‐ent approximation. The resulting approximation possesses some interesting properties, such as non‐negativity, ??∞ smoothness, exact interpolation on the boundary and variation diminishing (no Gibbs effect). The resulting structure has many similarities with B‐spline surfaces, but without the tensor‐product structure typical of that approximation. Examples are provided of its use in the framework of a Galerkin method showing the potential of the proposed method in solving boundary value problems. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
12.
Yong Gan Zheng Sun Zhen Chen Xiong Zhang Yu Liu 《International journal for numerical methods in engineering》2018,113(3):411-431
The material point method (MPM) enhanced with B‐spline basis functions, referred to as B‐spline MPM (BSMPM), is developed and demonstrated using representative quasi‐static and dynamic example problems. Smooth B‐spline basis functions could significantly reduce the cell‐crossing error as known for the original MPM. A Gauss quadrature scheme is designed and shown to be able to diminish the quadrature error in the BSMPM analysis of large‐deformation problems for the improved accuracy and convergence, especially with the quadratic B‐splines. Moreover, the increase in the order of the B‐spline basis function is also found to be an effective way to reduce the quadrature error and to improve accuracy and convergence. For plate impact examples, it is demonstrated that the BSMPM outperforms the generalized interpolation material point (GIMP) and convected particle domain interpolation (CPDI) methods in term of the accuracy of representing stress waves. Thus, the BSMPM could become a promising alternative to the MPM, GIMP, and CPDI in solving certain types of transient problems. 相似文献
13.
A. Rassineux P. Breitkopf P. Villon 《International journal for numerical methods in engineering》2003,57(3):371-389
We present a method to adapt a tetrahedron mesh together with a surface mesh with respect to a size criterion. The originality of our work lies in the fact that both surface and tetrahedron mesh adaptation are carried out simultaneously and that no CAD is required to adapt the surface mesh. The adaptation procedure consists of splitting or removing interior and surface edges which violate a given size criterion. The enrichment process is based on a bisection technique. In order to guarantee mesh conformity during the refinement process, all possible remeshing configurations of tetrahedra have been examined. Once the tetrahedron mesh has been adapted, surface nodes are projected on a geometrical model. The building of a surface model from discrete data has already been presented in this journal. The method is based on a mesh‐free technique called Hermite Diffuse Interpolation. Surface and volume mesh optimization procedures are carried out during the adaptation and at the end of the process to enhance the mesh. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
14.
Model and mesh generation of cracked tubular Y-joints 总被引:1,自引:0,他引:1
In this study, the methods for constructing accurate and consistent geometrical and finite element (FE) models for general cracked tubular Y-joints are described. Firstly, geometrical analysis of welded tubular joint is given and it is then extended to the modelling of general cracked Y-joints. The concept of crack surface and a simple mapping approach are suggested to model either through-thickness or surface cracks which can be of any length and located at any position along the brace-chord intersection. Secondly, the geometrical model developed will be used in the generation of consistent FE meshes. The basic concepts used for the design and generation of three-dimensional FE meshes will be described. This will include the meshing procedures for discretization of tubular joints with through-thickness and surface cracks which are frequently regarded as one of the most difficult steps in the construction of tubular joint models. Finally, some mesh generation examples for uncracked and cracked Y-joints will be presented to demonstrate the use of the purposed geometrical model and mesh generation scheme developed. 相似文献
15.
A. Gargallo‐Peiró X. Roca J. Peraire J. Sarrate 《International journal for numerical methods in engineering》2016,106(13):1100-1130
A framework to validate and generate curved nodal high‐order meshes on Computer‐Aided Design (CAD) surfaces is presented. The proposed framework is of major interest to generate meshes suitable for thin‐shell and 3D finite element analysis with unstructured high‐order methods. First, we define a distortion (quality) measure for high‐order meshes on parameterized surfaces that we prove to be independent of the surface parameterization. Second, we derive a smoothing and untangling procedure based on the minimization of a regularization of the proposed distortion measure. The minimization is performed in terms of the parametric coordinates of the nodes to enforce that the nodes slide on the surfaces. Moreover, the proposed algorithm repairs invalid curved meshes (untangling), deals with arbitrary polynomial degrees (high‐order), and handles with low‐quality CAD parameterizations (independence of parameterization). Third, we use the optimization procedure to generate curved nodal high‐order surface meshes by means of an a posteriori approach. Given a linear mesh, we increase the polynomial degree of the elements, curve them to match the geometry, and optimize the location of the nodes to ensure mesh validity. Finally, we present several examples to demonstrate the features of the optimization procedure, and to illustrate the surface mesh generation process. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
16.
G. L. Ghiringhelli 《Composites Part B》1997,28(5-6):613-626
This paper describes a method of evaluating the linear stiffness matrix of a three-dimensional beam which accounts for all the possible couplings involved in the use of composite laminations appropriate to satisfy design requirements. The starting point of the procedure is represented by a complete evaluation of the beam cross section properties by use of a two-dimensional finite element analysis. A mass matrix of the beam is also evaluated as well as the external forces, to account for overall thermal deflections.
Some examples are presented and the results are compared either with their theoretical counterparts or with numerical results obtained by means of a full three-dimensional finite element analysis. 相似文献
17.
A. Gargallo‐Peiró X. Roca J. Peraire J. Sarrate 《International journal for numerical methods in engineering》2015,103(5):342-363
We present a robust method for generating high‐order nodal tetrahedral curved meshes. The approach consists of modifying an initial linear mesh by first, introducing high‐order nodes, second, displacing the boundary nodes to ensure that they are on the computer‐aided design surface, and third, smoothing and untangling the mesh obtained after the displacement of the boundary nodes to produce a valid curved high‐order mesh. The smoothing algorithm is based on the optimization of a regularized measure of the mesh distortion relative to the original linear mesh. This means that whenever possible, the resulting mesh preserves the geometrical features of the initial linear mesh such as shape, stretching, and size. We present several examples to illustrate the performance of the proposed algorithm. Furthermore, the examples show that the implementation of the optimization problem is robust and capable of handling situations in which the mesh before optimization contains a large number of invalid elements. We consider cases with polynomial approximations up to degree ten, large deformations of the curved boundaries, concave boundaries, and highly stretched boundary layer elements. The meshes obtained are suitable for high‐order finite element analyses. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
18.
杨方鸿;杜文韬;车路;张代军 《复合材料学报》2025,42(8):4774-4784
变角度复合材料具有优异的力学性能和更强的可设计性,研究变角度复合材料结构的力学行为和损伤特性对工程应用和结构设计具有重要科学意义。单元型近场动力学(Element-based peridynamics,EBPD)理论基于积分-微分方程,处理不连续问题时不存在奇异性问题,相比于经典连续介质力学方法可以方便处理断裂等破坏问题。通过虚拟单元积分点的位置与纤维角度的对应关系,提出了EBPD变角度复合材料单层板模型。通过对具有不同边缘初始角度含中心圆孔的变角度复合材料单层板变形模拟验证了模型的可靠性,并预测不同角度单层板的裂纹扩展路径。提出的模型可以准确表征变角度复合材料单层板的变形以及破坏模式,为其强度评估及结构设计提供理论指导。 相似文献
19.
Łukasz Kaczmarczyk Mohaddeseh Mousavi Nezhad Chris Pearce 《International journal for numerical methods in engineering》2014,97(7):531-550
This paper presents a computational framework for quasi‐static brittle fracture in three‐dimensional solids. The paper sets out the theoretical basis for determining the initiation and direction of propagating cracks based on the concept of configurational mechanics, consistent with Griffith's theory. Resolution of the propagating crack by the FEM is achieved by restricting cracks to element faces and adapting the mesh to align it with the predicted crack direction. A local mesh improvement procedure is developed to maximise mesh quality in order to improve both accuracy and solution robustness and to remove the influence of the initial mesh on the direction of propagating cracks. An arc‐length control technique is derived to enable the dissipative load path to be traced. A hierarchical hp‐refinement strategy is implemented in order to improve both the approximation of displacements and crack geometry. The performance of this modelling approach is demonstrated on two numerical examples that qualitatively illustrate its ability to predict complex crack paths. All problems are three‐dimensional, including a torsion problem that results in the accurate prediction of a doubly‐curved crack. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
20.
Curvature determination in the bending test of continuous fibre reinforcements 总被引:2,自引:0,他引:2
下载免费PDF全文

Bending deformation is one of the main deformation modes in the forming of continuous fibre reinforcements. It is very specific compared to the classical continuous materials. In most textile material forming simulation models, bending stiffness is neglected, but taking into account it would give a more accurate predication, particularly the shape of wrinkles. In the deformation, fibres constituting the reinforcement can slide to each other, resulting in the bending stiffness of reinforcement is not directly related to its in‐plane tensile modulus as the classical continuous materials. Consequently, it is necessary to determine the bending stiffness by the experimental method. A cantilever bending stiffness test approach was proposed to measure the bending stiffness of continuous fibre reinforcements. A CCD camera was used to take the bending deflection shape. To calculate the curvature of bending deflection curve, uniform quartic B‐spline curve was chosen. A detailed process of curvature calculating of deflection curve was given in this paper. A single bending test gives the bending moment as a function of curvature from zero to the maximum value. Bending tests were conducted to validate the capacity of the proposed approach. A smooth curvature plot along the bending deflection curve was obtained for different types of reinforcements, and the non‐linear bending deformation behaviour can be identified with the proposed bending test method. 相似文献