3D axisymmetry to 2D

[sidharthsarmah]

Hi @GiudGiud ,

When I wished to do this extruded 3D strutted in 2D, I used:

[Problem]
  type = FEProblem
  coord_type = XYZ
[]

Could you suggest, what should be it for an axisymmetry case, I don't want an extruded one like in the figure as both have similar projection. I want the one on the right. Does this work? Thanks.

[Problem]
  type = FEProblem
  coord_type = RZ
[]

[GiudGiud]

you specify the symmetry axis (in the Mesh, in the Problem it is a deprecated parameter)

pick the axis that is horizontal in the mesh you show at the bottom

[sidharthsarmah]

Thanks

[sidharthsarmah]

Hi @GiudGiud , I have one confusion, when I apply a point load in 2D, it resembles a line load in 3D like this:

When I apply a line load in the RZ coord, it seems the load in something like this (right), but I want it to be like the other way(left), a load on top for bending:

[Mesh]
  coord_type = RZ
  rz_coord_axis = X
  [./mesh]
    type = FileMeshGenerator
    file = Cyl.inp
  [../]
  [./bounding_box]
    type = BoundingBoxNodeSetGenerator
    input = 'mesh'
    new_boundary = 'nodes_to_convert_to_sides'
    top_right = '25.3999996  5.71500015 0'
    bottom_left = '20 5 0'
  [../]
  [./create_sideset]
    type = SideSetsFromNodeSetsGenerator
    input = 'bounding_box'
  [../]
[]

[BCs]
  [./u_top_pull]
    type = ADPressure
    variable = disp_y
    boundary = load
    function = 0.1*t
  []

[GiudGiud]

This is all expected, and there's no way around it. It's by definition of these coordinate systems.

[sidharthsarmah]

Hi @GiudGiud , I agree on this but is there any way to just add a load at a particular theta along the radius? That would be interesting.

[GiudGiud]

nope we unfortunatly do not have R-theta-z implemented, only RZ