Error with Input Script in Non Homogeneous Elasticity #28896
Unanswered
ppandit95
asked this question in
Q&A General
Replies: 0 comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
Check these boxes if you have followed the posting rules.
Question
Hi Developers
In my attempt to solve steady state mechanical equilibrium equation in a domain with circular precipitate,I ended up with the input script as -
[GlobalParams]
displacements = "ux uy"
[]
[Variables]
[./ux]
family = LAGRANGE
order = FIRST
[../]
[./uy]
family = LAGRANGE
order = FIRST
[../]
[]
[Mesh]
type = GeneratedMesh
dim = 2
nx = 500
ny = 500
elem_type = QUAD4
[]
[AuxVariables]
[./phi]
order = CONSTANT
family = MONOMIAL
[../]
[]
[AuxKernels]
[./phi_aux]
type = FunctionAux
variable = phi
function = circle_function
[../]
[]
[Functions]
[./circle_function]
type = ParsedFunction
expression = "sqrt(pow(x-250,2) + pow(y-250,2)) < 50 ? 1.0 : 0.0"
[../]
[]
[Materials]
[./matrix]
type = ComputeIsotropicElasticityTensor
youngs_modulus = "210000"
poissons_ratio = "0.3"
block = 0
[../]
[./precipitate]
type = ComputeIsotropicElasticityTensor
youngs_modulus = "300000"
poissons_ratio = "0.3"
block = 0
[../]
[./variable_elasticity]
type = ComputeVariableIsotropicElasticityTensor
youngs_modulus = "210000 + (300000 - 210000) * phi"
poissons_ratio = "0.3"
args = "phi"
block = 0
[../]
[]
[Physics/SolidMechanics/QuasiStatic]
[./block1]
strain = SMALL
add_variables = true
incremental = true
generate_output = "stress_xx stress_xy stress_yy"
[../]
[]
[Executioner]
type = Steady
solve_type = PJFNK
petsc_options_iname = "-pc_type"
petsc_options_value = "jacobi"
[]
[Outputs]
[./vtk]
type = VTK
file_base = "test"
[../]
[]
[BCs]
[./left]
type = DirichletBC
variable = "ux"
value = "0"
boundary = "left"
[../]
[./bottom]
type = DirichletBC
variable = "uy"
value = "0"
boundary = "bottom"
[../]
[./right]
type = FunctionDirichletBC
variable = "ux"
function = "disp_x"
boundary = "right"
[../]
[./top]
type = FunctionDirichletBC
variable = "uy"
function = "disp_y"
boundary = "top"
[../]
[]
[Functions]
[./disp_x]
type = ParsedFunction
expression = "0.5"
[../]
[./disp_y]
type = ParsedFunction
expression = "0.5"
[../]
[]
But still I m facing the following error and I m not able to debug it so any help in this regard will be appreciated.
Many Thanks
Pushkar
Beta Was this translation helpful? Give feedback.
All reactions