Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

FE method, how ? (software Comsol)

  1. Jun 23, 2014 #1
    Hi all,

    The plastic flow rule for large strains in a continuum medium can be written as :
    [itex]-\dfrac{1}{2} \Big ( \dfrac{d}{dt}(F_p^{-1} ). ^t F_p^{-1} + F_p^{-1} . ^t \dfrac{d}{dt} ( F_p^{-1} ) \Big)= \lambda F^{-1} \dfrac{\partial f}{\partial \tau} F F_p^{-1} . ^t F_p^{-1}[/itex]
    where [itex]F_p[/itex] is the plastic part of the total deformation gradient [itex]F[/itex], [itex]f[/itex] the yield surface and [itex]\tau[/itex] the Kirchhoff stress tensor.
    According to Comsol notes, we denote [itex]M=F_p^{-1}[/itex] and discretize the law by using 'variations with [itex]\Delta[/itex]' :
    [itex]-\dfrac{1}{2} \Big [ 2M. ^t M -M_{old} . ^t M - M . ^t M_{old} \Big ]= \lambda \Delta t F^{-1} \dfrac{\partial f}{\partial \tau} FM. ^t M [/itex]​
    but I don't understand what is stated after. I quote :

    Few things I don't understand or I'm not sure of :

    1/ Gauss points stand for the mesh nodes, right ?
    2/ What does "For each Gauss point the plastic strain variables are computed" mean ? We have matrices that describe the body state at all his points, that is, there is not a matrix for each point of the body, so I don't understand the beginning of their sentence.
    3/ The total deformation gradient [itex]F[/itex] is also a variable, we don't know its value, so how the computing is done for [itex]M[/itex]

    Last edited: Jun 23, 2014
  2. jcsd
  3. Jun 23, 2014 #2
    For the third point, does F depend on the boundary conditions ? For instance, if we consider a round tensile specimen with prescribed displacement at the top and below surface and with free stress on the lateral surface, can we infer F ?
    Sure, we know its impqct on points at the top and below surface, but can we infer all its components ? If the total displacement vector is known at the two outer surfaces, then yes we can know all components of F, but if only the radial displacement is prescribed, some components can not be determined
    Last edited: Jun 23, 2014
  4. Jun 23, 2014 #3
    No, the Gauss points are no the mesh nodes. The Gauss Points, or Gaussian Integration Points, are internal to the elements. More can be found here (last paragraph of the page marked 11) and here.
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook