1. Field of the Invention
The present invention generally relates to methods, systems and software product used in the area of computer-aided engineering analysis (e.g., finite element method (FEM), meshfree method, etc.), more particularly to methods and systems for numerically simulation structural behaviors of embedded bi-materials (e.g., particle-reinforced composites, fiber-reinforced composites, etc.).
2. Description of the Related Art
A composite material is a microscopic or macroscopic combination of two or more distinct materials with a recognizable interface between them. Composite materials are developed because no single, homogeneous structural material could be found that had all of the desired characteristics for a given application. For example, fiber-reinforced composites are developed to replace aluminum alloys to provide high strength and fairly high stiffness at low weight without corrosion and fatigue.
For predicting the structural behaviors of such composite materials, computer aided engineering analysis has been used. However, numerically simulating embedded bi-materials using prior art approaches have a number of shortcomings.
A simplified embedded bi-material model in two-dimension and a corresponding prior art FEM model 120 are shown in
Other prior art approaches have problems also. In one example, motar FEM uses a domain decomposition technique to treat mismatching meshes. However, motar FEM requires using Lagrange multipliers that sometimes violate the so-called “inf-sup” condition, which can lead to numerical instability (i.e., no solution). In another examples, “generalized FEM (GFEM)”, “the extended FEM (XFEM)”, “immersed FEM (IFEM)” have been used but the results are too expensive to achieve (e.g., very high computational requirements/costs or not clear to apply in three-dimension).
Another prior art approach is based on meshfree method. Meshfree method has become one of the focused research topics during the 1990's. Many applications of using mesh-free analysis have been achieved in the past decade.
An exemplary mesh-free model 200 is shown in
Due to the flexibility of the meshfree nodal representation 204 of the physical domain 202, a practical way to create a computer model for the meshfree method is to use the FEM model's nodal data that is readily generated from a pre-processing software package. The pre-processing software may be a stand-alone software package or a built-in portion of an engineering design or analysis computer program product package.
However, for simulating structural behaviors of an embedded bi-material, prior art meshfree methods require adding interface constraints and a set of interface nodes to ensure the visibility of the interface in the numerical simulation results. Adding the interface nodes must be done manually. Further, each added node's integration cell must also be adjusted such that the domain integration can be properly conducted. Extending the manual nodal adjustment in the three-dimensional cases, it is not only nontrivial but sometimes impossible.
It would, therefore, be desirable to have a new improved computer aided engineering analysis method and system for numerically simulation structural behaviors of embedded bi-materials.
This section is for the purpose of summarizing some aspects of the present invention and to briefly introduce some preferred embodiments. Simplifications or omissions in this section as well as in the abstract and the title herein may be made to avoid obscuring the purpose of the section. Such simplifications or omissions are not intended to limit the scope of the present invention.
The present invention discloses methods and systems for numerically simulating structural behaviors of embedded bi-materials (e.g., composites). According to one aspect of the present invention, at least first and second grid models are created independently for an embedded bi-material that contains an immersed material (i.e., second material) embedded entirely within a base material (i.e., first material).
The first and second grid models are used for defining the geometry of the domain (i.e., embedded bi-material) and used as integration cell for meshfree method. The first grid model represents the base and immersed materials. The second grid model represents the immersed material.
A set of interface nodes is determined to define the material interface between the base and immersed materials. By examining the second grid model, those nodes of the second grid model located on the outer border or edge are designated as the interfaced nodes.
First and second groups of meshfree nodes are created basing on the grid models. The first group of meshfree nodes represents the entire domain (i.e., base plus immersed materials). The second group of meshfree nodes represents the immersed or embedded material, which includes all interface nodes and nodes located within a space bordered by the material interface.
Numerical structural behaviors of the embedded bi-material are simulated using the first and second set of meshfree nodes with a meshfree method that combines two meshfree approximations. The first meshfree approximation covers the first set of meshfree nodes and is based on properties of the base material, while the second meshfree approximation covers the second set of meshfree nodes and is based on a differential between the immersed and base materials.
Other objects, features, and advantages of the present invention will become apparent upon examining the following detailed description of an embodiment thereof, taken in conjunction with the attached drawings.
These and other features, aspects, and advantages of the present invention will be better understood with regard to the following description, appended claims, and accompanying drawings as follows:
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will become obvious to those skilled in the art that the present invention may be practiced without these specific details. The descriptions and representations herein are the common means used by those experienced or skilled in the art to most effectively convey the substance of their work to others skilled in the art. In other instances, well-known methods, procedures, components, and circuitry have not been described in detail to avoid unnecessarily obscuring aspects of the present invention.
Reference herein to “one embodiment” or “an embodiment” means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the invention. The appearances of the phrase “in one embodiment” in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. Further, the order of blocks in process flowcharts or diagrams representing one or more embodiments of the invention do not inherently indicate any particular order nor imply any limitations in the invention.
Embodiments of the present invention are discussed herein with reference to
Referring first to
Process 300 starts by receiving configuration of an embedded bi-material at step 302. For example, the embedded bi-material contains base (outer) and immersed (inner) materials. The immersed material is entirely embedded inside the base material (i.e., see
Creations of the grid models are independent of each other. In other words, there is no correlation between these two grid models, for example, nodes need not to be matched or aligned. Each grid model contains a plurality of nodes that defines the mesh.
Next, at step 306, a set of interface nodes is determined and designated. The set of interface nodes defines the material interface between the base and immersed materials. On exemplary method to determine the interface nodes is to examine the grid model of the immersed material (i.e., the second computerized grid model). Nodes located on the outer border/edge are designated as the interface nodes (shown as solid black circular dots 441 in
At step 308, a first group of meshfree nodes representing the entire embedded bi-material is created. One exemplary method is to combine the nodes from the first and second grid models without any duplication. In other words, any duplicated node is removed during the combination procedure. An exemplary first group of meshfree nodes 430 (i.e., all circles—hollow, solid, filled) with a corresponding grid model 435 are shown in
At step 310, a second group of meshfree nodes 440, shown in
At step 312, at least two meshfree approximations (i.e., first and second meshfree approximations) are performed. The first meshfree approximation covers the first group of meshfree nodes and is carried out using the properties of the base material. The second meshfree approximation covers the second group of meshfree nodes and is carried out using a property differential between the immersed material and the base material. Two integrals of Eqs. (17)-(18) listed below represent these two meshfree approximations, respectively.
Finally, at step 314, numerically simulated structural behaviors (e.g., stresses, strains, displacements, etc.) are obtained by combing the first and the second meshfree approximations' result.
Problem Description and Variational Equation
We consider an elastic solid occupying a bounded and open domain Ω⊂R2 with Lipschitz boundary. Let ∂ΩD and ∂ΩN be two open subsets of boundary ∂Ω such that ∂Ω=∂ΩD∪∂ΩN and ∂ΩD∩∂ΩN0. g(x) is the prescribed displacement applied on the Dirichlet boundary ∂ΩD, and t(x)εL2(∂ΩN) is the prescribed traction applied on the Neumann boundary ∂ΩN with n0 denoting the outward unit normal to the boundary ∂ΩN. The elastic body is composed of two perfectly bounded material with zero-thickness interface Γ. The equilibrium configuration of the elastic body is characterized by the continuity of displacement and continuity of normal stress across the material interface Γ. The elasticity interface problem can be described by the following second-order elliptic boundary value problem with the associated homogenous Dirichlet and Neumann jump conditions on the interface Γ.
where function u: Ω→R2 is the displacement and ƒ: Ω→R2 is the prescribed body force over the domain Ω. Normally, ƒεL2(Ω). The notation ∇su denotes the symmetric gradient of the displacement, ∇su=½(∇u+(∇u)T). Without loss of generality, we assume the interface Γ is a smoothed and closed curve that divides the global domain Ω into two regions: Ω+ representing the matrix and Ω− denoting the immersed media or inclusion such that their union gives the global domain Ω,
[[q]](x)=q+(x)−q−(x) (2)
in which + and − denote the two sides of the interface Γ with the jump of quantity q across the interface. The body force and material constants can exhibit discontinuities across the interface Γ, but have smooth restrictions ƒ+, C+ to Ω+ and ƒ−, C− to Ω−. They are given by
The infinitesimal strain tensor ε(u) is defined by
ε(u)=∇su=½(∇u+(∇u)T) (4)
C+ and C−εL∞(Ω) are elasticity tensors with major and minor symmetries and are corresponding to different materials in Ω+ and Ω− respectively. In the case of linear isotropic elasticity, we take C+ and C− to be constants. The Cauchy stress tensor σ and strain tensor ε have the following relationship
where the positive constants μ+, μ− and λ+, λ− are Lamé constants. The Lamé constants are related to the Young's modulus E and Poisson ratio v by
The variational form of this problem is to find the displacement uεVg={vεH1(Ω): v=g on ∂ΩD} such that for all vεV
a(u,v)=l(v) (7)
where the functional space V=H01(Ω) consists of functions in Sobolev space H1(Ω) which vanish on the boundary ∂Ω and is defined by
V(Ω)={v:vεH1, v=0 on ∂ΩD} (8)
The bilinear form a(•,•) and linear functional l(•) are obtained by multiplying the test function vεV to both sides of Eq. (7) and integrating over the regions Ω+ and Ω− separately using Green's formula.
Using the fact that n+=−n−, we can rewrite the above equation to
Applying the homogenous Neumann jump condition to Eq. (10) yields
It is noted that the elasticity tensors C+ and C− are symmetric, and homogenous Neumann jump condition is enforced in the variational level. Obviously, the bilinear form a(•,•) in Eq. (12) is symmetric, bounded and coercive by Friedrichs inequality. The existence and uniqueness of the problem is ensured by the Lax-Milgram theorem
Embedded Meshfree Method for Interface Elasticity Problem
This section first devotes to the development of a new meshfree discretization that spans across the material interface on the overlapping meshes. The definition of computation domain for the proposed method and the construction of the meshfree approximation that satisfies a point-wise continuity across the interface are described subsequently. Finally, a primal problem which is equivalent to a degenerated Lagrangain-type mortar method is devised and followed by a proof of the optimality.
Embedded Meshfree Discretization and Integration Cells
For simplicity, we assume that Ω is a convex polygonal domain and only consider the case of pure displacement problem under homogenous boundary condition (∂ΩD=0). The standard meshfree Galerkin method is formulated on a finite dimensional space Vh⊂V employing the variational formulation of Eq. (11) to find uhεVh such that
a(uh,vh)=l(vh)∀vhεVh (14)
where Vh=span{ΨI: IεZI} and ZI is an index set. {ΨI(x)}Iε
We let hI=diam(rI) denote the nodal support radius and assume the following overlapping condition
∃MεR ∀xεΩ card ZI≦M (16)
Often, a shape function ΨI(x) is associated with a particle xIεR2 and particles are distinct. Using the superposition principle for the above linear system, the discrete bilinear form in Eq. (14) can be re-expressed by
Similarly, we have the discrete linear functional to be rewritten as
The expression of discrete forms in Eqs. (17) and (18) allows us to reset the meshfree computation domain by two overlapping sub-regions; namely Ω+∪Ω− and Ω−. The advantage of overlapping sub-regions is their flexibility to accommodate complex immersed structures in which different discretizations can be made easily and independently in each sub-region using finite element model. Since the computation domain Ω− is “embedded” in computation domain Ω+∪Ω− under the meshfree Galerkin framework, this method is referred to as an embedded meshfree method.
Consequently, we can define the computational sub-regions as follows: Given a bounded domain Ω⊂R2, we consider sub-regions Ω1 and Ω2 to be the overlapping sub-regions satisfying
where Ω1=Ω+∪Ω− is the computation domain containing the base matrix, and Ω2=Ω− is the computation domain containing the inclusion. We assume that each Ωi is partitioned by a finite element triangulation Th
In the sub-region Ω1, the total number of meshfree nodes consists of a set of component nodes that overlap and cover the domain Ω.
Let Z1={xl, l=l, . . . NP} be the set of distinct nodes in Ω1. NP indicates the total number of overlapping nodes in sub-region Ω1. For each xlεZ1, Ψl(x) denotes the corresponding meshfree shape function. We define the meshfree interpolant of u(x) by the formula
where uI=u(xI) is called the ‘generalized’ displacement of node I. In other words, uI(x) is considered as an interpolant of u(x) in a generalized sense. In general, conventional meshfree approximations are not interpolants, i.e, uI≠uI(xI). For this reason, special treatment is required to impose the essential boundary conditions on the global boundary ∂Ω of the model problem. We employ the generalized meshfree approximation method (GMF) to obtain the meshfree convex approximation. The first-order convex GMF approximation is constructed using the inverse tangent basis function and the cubic spline window function is chosen to be the weight function in GMF method. Giving a convex hull co(Zl) of a node set Zl={xI, I=1, . . . NP}⊂R2 defined by
the GMF method is to construct a convex approximation of a given (smooth) function u in the form of Eq. (20) such that the shape function ΨI: co(Zl)→R satisfies the following linear polynomial reproduction property
The convex approximation space constructed by GMF method is a subspace of H01(Ω) and conforms to the boundary conditions if the approximating domain is convex. Ideally this conforming meshfree approximation secures H1—-compatibility and the homogenous Dirichlet jump condition across the interface is verified automatically. On the other hand, the meshfree approximation also introduces the non-locality across the interface. This gives rise to a solution that exhibits a smearing near the interface. To remove the smearing, we invoke a second meshfree approximation in sub-region Ω2, e.g. by zero extension in the sub-region Ω1.
To be more precise, we let {tilde over (Z)}2={{tilde over (x)}l, l=1, . . . MP}⊂Z1 be the subset nodes that contain the overlapping nodes in the sub-region Ω2. MP is the total number of overlapping nodes in sub-region Ω2.
Analogously, every function ũh(x)ε{tilde over (V)}h2(Ω2) has a unique representation of the form
Since the sub-region Ω2 is assumed to be convex, the subspace {tilde over (V)}h2(Ω2) is defined by
{tilde over (V)}h2(Ω2)={v:v|Ω
Apparently, the subspace {tilde over (V)}h2(Ω2) is not a subset of subspace Vh(Ω1), i.e., {tilde over (V)}h2(Ω2)Vh(Ω1). Since the approximations in sub-region Ω1 and sub-region Ω2 are constructed independently, the meshfree shape functions {tilde over (Ψ)}l(x) and Ψl(x) of the same node xlε{tilde over (Z)}2⊂Z1 are not necessarily the same. This is true in particular when support of node xl covers the interface Γ, i.e.,
{tilde over (Ψ)}l(x)≠Ψl(x) for xlε{tilde over (Z)}2⊂Z1 and supp(xl)∩Γ≠0 (25)
The support of shape function in sub-region Ω2 is defined by
supp(xl)=supp({tilde over (Ψ)}l(x))={x|εΩ2 and {tilde over (Ψ)}l(x)≠0} (26)
As a result, it leads to a non-conforming meshfree approximation and the continuity of displacement across the interface is not ensured. A one-dimensional example is shown in
Construction of Meshfree Approximation and Modified Variational Formulation
In mortar finite element method, the non-conforming problem is caused by the non-matching meshes across the interface where no sharing nodes are established to provide the displacement continuity. Therefore an additional constraint equation corresponding to the Dirichlet jump condition is required to connect the disjointed nodes between two triangulations. In standard mortar finite element method, the constraint equation is imposed weakly across the interface and is sufficient to guarantee an approximation with a consistency error of order h if the weak solution u is smooth enough. However in embedded meshfree method, the sharing nodes are well-defined in the subset of interface-fitted nodes {tilde over ({tilde over (Z)}2={{tilde over ({tilde over (x)}l, l=1, . . . IP}. Indeed, the non-conformity of approximation in embedded meshfree method is due to the non-matching meshfree shape functions in the overlapping domain as illustrated in
Namely we define a new constrained discrete meshfree approximation space by
where the definition of subspace {tilde over (V)}h2(Ω2)⊂{circumflex over (V)}h(Ω) is given in Eq. (24) and additional subspace {tilde over (V)}h1(Ω1) is defined by
{tilde over (V)}h1(Ω1)={v:v|Ω
We note that {tilde over (V)}hi(Ωi), i=1, 2 stand for the spaces of linear conforming meshfree approximations that satisfy homogenous Dirichlet boundary conditions on ∂Ω∩∂Ωi, i=1, 2 in the sub-regions Ω1\Ω2 and Ω2 respectively. Let {tilde over (Z)}1={{tilde over (x)}l, l=1, . . . LP}⊂Z1 be the subset nodes such that
{tilde over (Z)}1=(Z1\{tilde over (Z)}2)∪{tilde over ({tilde over (Z)}2 and LP=NP−MP+IP (29)
In a one dimensional case, the new approximation is conforming since {circumflex over (V)}h(Ω)⊂V(Ω). An improved meshfree approximation is shown in
However, this is not the case in a two dimensional problem. Since sub-region Ω1\Ω2 is not convex (by assuming Ω2 is convex previously), the approximation of the interface-fitted nodes generated by the GMF method does not possess a Kronecker-delta property. In other words, the displacement continuity at the interface-fitted node does not hold. A graphic sketch of this nonconforming approximation near the interface is given in
We impose the displacement continuity weakly and point-wisely across the interface by introducing the Kronecker-delta property to the approximation of the interface-fitted nodes. This can be done by either employing a singular kernel function to the interface-fitted nodes or applying a transformation method to those nodes whose supports cover the interface nodes. We have adopted the transformation method in this study. With the introduction of Kronecker-delta property to the shape functions at the interface-fitted nodes set {tilde over ({tilde over (Z)}2={{tilde over ({tilde over (x)}l, l=1, . . . IP}, the reconstructed meshfree shape functions in sub-region Ω1 have the following form:
In general {circumflex over (Ψ)}I1(x)≠{circumflex over (Ψ)}I2(x) on interface Γ except at the interface nodes x=xIε{tilde over ({tilde over (Z)}2. A nonconforming meshfree approximation with point-wise continuity across the interface is illustrated in
Since the displacement continuity across the interface is only imposed point-wisely at the interface-fitted nodes, the weak form of Eq. (14) is reformulated based on the nonconforming meshfree approximation. We now define the embedded meshfree solution of the elasticity interface problem as a function uhε{circumflex over (V)}h satisfying
and the associated discrete (broken) energy norm is defined by
∥vh∥1,h=â(vh,vh)1/2, vhε{circumflex over (V)}h (35)
Using the Cauchy-Schwarz inequality and triangle inequality, it can be shown that the modified bilinear form â(•,•) is bounded on {circumflex over (V)}h×{circumflex over (V)}h with respect to the broken energy norm on Ω.
Lemma 1
There exists a positive constant cb such that for any uh, vhε{circumflex over (V)}h, we have
|â(uh,vh)|≦cb∥uh∥1,h∥vh∥1,h (36)
Proof.
In inequality (Eq. (37)), we have used the discrete semi-norm on space {circumflex over (V)}h defined by standard Sobolev notation as
|vh|1,h2=|vh|1,Ω
The constants γmax+ and γmax− in inequality (Eq. (37)) are the largest eigenvalues of C+ and C− respectively.
Lemma 2
There exists a positive constant cc such that for any vhε{circumflex over (V)}h, we have
â(vh,vh)≧cc∥vh∥1,h2 (39)
Proof.
Observing that â(vh,vh)=0 implies vh is constant. Since vh vanishes on global boundary ∂Ω and satisfies continuity at the interface-fitted nodes, we have vh=0 in Ω and thus ensure the coercivity of the modified bilinear form.
Theorem 1
Let uεV be the solution of the variational problem (Eq. (7)), the discretization problem (Eq. (32)) in embedded meshfree method admits a uniqueness solution uhε{circumflex over (V)}h.
A Degenerated Mortar Method and Energy-Norm Error Estimate
Since â(•,•) is coercive, we can apply the well-known second Strang's lemma for the energy-norm error estimate.
where the first term on the RHS of Eq. (40) is the best approximation error which can be obtained by the meshfree approximation error estimate and the second term is the consistency error which comes from the nonconforming of {circumflex over (V)}h. Assume the regularity on the exact solution uεH2(Ω), we have the following error estimate by Céa's inequality.
where h is the largest nodal support radius, i.e, h=supIεZ
Because the interface Γ is associated with a one dimensional triangulation, we call this one dimensional triangulation on Γ, Ξh. Each integration segment eiεΞh is a boundary edge of integration cell (finite element triangulation) Th
where the term
can be realized as the weak constraint equation appearing in the Lagrange multiplier-type mortar method in which λh=C·ε(u)·nε(L2(γ))2 are Lagrange multipliers. The symbol [[w]] denotes the restriction of jump as defined in Eq. (2) for wεΓ. A natural choice for the construction of Lagrange multiplier spaces in nonconforming formulation of the mortar method is to define the Lagrange multiplier basis function locally associated with the discrete nodes. From this point of view, we regard the proposed embedded meshfree method as a degenerated Lagrange multiplier-type mortar method due to the fact that the displacement jump vanishes at the interface-fitted nodes, i.e.,
[[wh]]I=(wI+−wI−)=0∀Iε{tilde over ({tilde over (Z)}2 (43)
Applying the point collocation method to the assembly of discrete constraint equation in mortar method and using Eq. (42) lead to the primal problem presented in Eq. (32).
Despite the vanishing constraint equation in the primal problem, the proposed embedded meshfree method still presents certain boundary quadrature error by different quadrature rules and requires a consistency error estimate to ensure a stable and convergent meshfree discretization. This consistency error estimate resembles the consistency error estimate of meshfree solution in standard Galerkin meshfree method using moving least-squares approximation or reproducing kernel approximation when Dirichlet boundary conditions are imposed point-wisely.
Lemma 3
Assume that uεH2(Ω) be the solution of elasticity interface problem in Eq. (7), there exists a constant cc independent of h and function u such that
|â(u,wh)−{circumflex over (l)}(wh)|≦cch|u|2∥wh∥1,h∀whε{circumflex over (V)}h (45)
Proof.
Let Dn
where u=u+=u−, C=(Cij)2×2 is the material coefficient matrix and n is the outward unit normal vector on Γ as defined in
where we have related the largest nodal support radius h to the maximum length of element edge ei along the interface Γ by defining
∃c1≧1, h=supIεZ
Lemma 3 implies
which indicates that the estimate of consistency error for the embedded meshfree method is of O(h). Combining inequalities (Eqs. (40), (41) and (49)), we obtain the following result for the energy error.
Proposition 1.
Let uεH2(Ω) and uhε{circumflex over (V)}h be respectively the solutions of the weak problem (7) and of the discretized problem (32). Then it holds
∥u−uh∥1,h≦h|u|2 (50)
By Proposition 1, we expect the optimal rate of convergence for embedded meshfree method to be one for u in H2(Ω).
According to one aspect, the present invention is directed towards one or more computer systems capable of carrying out the functionality described herein. An example of a computer system 500 is shown in
Computer system 500 also includes a main memory 508, preferably random access memory (RAM), and may also include a secondary memory 510. The secondary memory 510 may include, for example, one or more hard disk drives 512 and/or one or more removable storage drives 514, representing a floppy disk drive, a magnetic tape drive, an optical disk drive, etc. The removable storage drive 514 reads from and/or writes to a removable storage unit 518 in a well-known manner. Removable storage unit 518, represents a floppy disk, magnetic tape, optical disk, etc. which is read by and written to by removable storage drive 514. As will be appreciated, the removable storage unit 518 includes a computer usable storage medium having stored therein computer software and/or data.
In alternative embodiments, secondary memory 510 may include other similar means for allowing computer programs or other instructions to be loaded into computer system 500. Such means may include, for example, a removable storage unit 522 and an interface 520. Examples of such may include a program cartridge and cartridge interface (such as that found in video game devices), a removable memory chip (such as an Erasable Programmable Read-Only Memory (EPROM), Universal Serial Bus (USB) flash memory, or PROM) and associated socket, and other removable storage units 522 and interfaces 520 which allow software and data to be transferred from the removable storage unit 522 to computer system 500. In general, Computer system 500 is controlled and coordinated by operating system (OS) software, which performs tasks such as process scheduling, memory management, networking and I/O services.
There may also be a communications interface 524 connecting to the bus 502. Communications interface 524 allows software and data to be transferred between computer system 500 and external devices. Examples of communications interface 524 may include a modem, a network interface (such as an Ethernet card), a communications port, a Personal Computer Memory Card International Association (PCMCIA) slot and card, etc. Software and data transferred via communications interface 524 are in the form of signals which may be electronic, electromagnetic, optical, or other signals capable of being received by communications interface 524. The computer 500 communicates with other computing devices over a data network based on a special set of rules (i.e., a protocol). One of the common protocols is TCP/IP (Transmission Control Protocol/Internet Protocol) commonly used in the Internet. In general, the communication interface 524 manages the assembling of a data file into smaller packets that are transmitted over the data network or reassembles received packets into the original data file. In addition, the communication interface 524 handles the address part of each packet so that it gets to the right destination or intercepts packets destined for the computer 500. In this document, the terms “computer program medium”, “computer usable medium”, and “computer readable medium” are used to generally refer to media such as removable storage drive 514, and/or a hard disk installed in hard disk drive 512. These computer program products are means for providing software to computer system 500. The invention is directed to such computer program products.
The computer system 500 may also include an input/output (I/O) interface 530, which provides the computer system 500 to access monitor, keyboard, mouse, printer, scanner, plotter, and alike.
Computer programs (also called computer control logic) are stored as application modules 506 in main memory 508 and/or secondary memory 510. Computer programs may also be received via communications interface 524. Such computer programs, when executed, enable the computer system 500 to perform the features of the present invention as discussed herein. In particular, the computer programs, when executed, enable the processor 504 to perform features of the present invention. Accordingly, such computer programs represent controllers of the computer system 500.
In an embodiment where the invention is implemented using software, the software may be stored in a computer program product and loaded into computer system 500 using removable storage drive 514, hard drive 512, or communications interface 524. The application module 506, when executed by the processor 504, causes the processor 504 to perform the functions of the invention as described herein.
The main memory 508 may be loaded with one or more application modules 506 that can be executed by one or more processors 504 with or without a user input through the I/O interface 530 to achieve desired tasks. In operation, when at least one processor 504 executes one of the application modules 506, the results are computed and stored in the secondary memory 510 (i.e., hard disk drive 512). The status of the computer simulation of embedded bi-material (e.g., meshfree method results) is reported to the user via the I/O interface 530 either in a text or in a graphical representation.
Although the present invention has been described with reference to specific embodiments thereof, these embodiments are merely illustrative, and not restrictive of, the present invention. Various modifications or changes to the specifically disclosed exemplary embodiments will be suggested to persons skilled in the art. For example, whereas only one immersed material has been shown and described, other numbers of immersed materials may be employed, for example, two, three or four immersed materials within a base material. Additionally, the shape of base material has been shown and described as a rectangular and the immersed material as a circular shape, the present invention sets no limit to the shapes (i.e., any enclosed shapes can be used for either material). Furthermore, only two-dimensional examples have been shown and described, the present invention applies to three-dimensional problem. In summary, the scope of the invention should not be restricted to the specific exemplary embodiments disclosed herein, and all modifications that are readily suggested to those of ordinary skill in the art should be included within the spirit and purview of this application and scope of the appended claims
Number | Name | Date | Kind |
---|---|---|---|
6718291 | Shapiro et al. | Apr 2004 | B1 |
7499050 | Wu et al. | Mar 2009 | B2 |
7660480 | Wu et al. | Feb 2010 | B1 |
8479588 | Simkins, Jr. | Jul 2013 | B1 |
Entry |
---|
Sukumar et al, “Modeling Holes and Inclusions by Level Sets in the Extended Finite Element Method”, Comput. Methods Appl. Mech. Engrg. 190, 2001. |
Belytschko et al, “A Coupled Finite Element-Element-Free Galerkin Method”, Computational Mechanics 17, 1995. |
Krongauz et al, “EFG Approximation with Discontinuous Derivatives”, International Journal for Numerical Methods in Engineering, vol. 21, 1998) teaches utilizing meshless methods in modeling bi-materials (Figure 7), Liu et al (“A Fast Boundary Element Method for the Analysis of Fiber-Reinforced Composites Based on a Rigid Inclusion Model”, Journal o. |
Karutz et al, “Automatic Adaptive Generation of a Coupled Finite Element/Element-Free Galerkin Discretization”, Finite Elements in Analysis and Design, 38, pp. 1075-1091, 2002. |
Zhu et al, “Three-Dimensional Numerical Modelling by XFEM of Spring-Layer Imperfect Curved Interfaces with Applications to Linearly Elastic Composite Materials”, International Journal for Numerical Methods in Engineering, Mar. 17, 2011. |
Hegen, D., “Element-Free Galerkin Methods in Combination with Finite Element Approaches”, Comput. Methods Appl. Mech. Engrg. 135, 1996. |
Hu et al, “Locally Enhanced Voronoi Cell Finite Element Model (LE-VCFEM) for Simulating Involving Fracture in Ductile Microstructures Containing Inclusions”, International Journal for Numerical Methods in Engineering, 76, 2008. |
Liu et al, “A Fast Boundary Element Method for the Analysis of Fiber-Reinforced Composites Based on a Rigid Inclusion Model”, Journal of Applied Mechanics, vol. 72, Jan. 2005. |
Number | Date | Country | |
---|---|---|---|
20130085727 A1 | Apr 2013 | US |