brazil_2d_nl.ctrl

c++的有限元源码,带说明书和测试数据,是难得的有限元开发资料.

brazil_2d_nl.oofem.out
2D brazilian splitting test
#linearstatic nsteps 1 renumber 1
#
#nonlinearstatic nsteps 100 renumber 1 stiffMode 1 controllmode 0 psi 0.0 steplength 0.0001 initialsteplength 0.0002 rtolv 0.0001 MaxIter 200 reqIterations 30 manrmsteps 10 hpcmode 1 hpc 2 3 1 nonlocalext lstype 4 smtype 7
#
nonlinearstatic nsteps 1 nmsteps 5 renumber 1 nonlocalext lstype 4 smtype 7
#
metastep 1 nsteps 4 stiffMode 2 controllmode 0 psi 0.0 steplength 0.0001 initialsteplength 0.0002 rtolv 0.0001 MaxIter 200 reqIterations 30 manrmsteps 10 hpcmode 1 hpc 2 3 1 keepll nonlocalext lstype 4 smtype 7
metastep 2 nsteps 3 stiffMode 2 controllmode 0 psi 0.0 steplength 0.00002 rtolv 0.0001 MaxIter 200 reqIterations 30 manrmsteps 10 hpcmode 1 hpc 2 3 1 keepll nonlocalext lstype 4 smtype 7
metastep 3 nsteps 4 stiffMode 2 controllmode 0 psi 0.0 steplength 0.00001 rtolv 0.0001 MaxIter 200 reqIterations 30 manrmsteps 10 hpcmode 1 hpc 2 3 1 keepll nonlocalext lstype 4 smtype 7
metastep 4 nsteps 9 stiffMode 2 controllmode 0 psi 0.0 steplength 0.0001 rtolv 0.0001 MaxIter 200 reqIterations 30 manrmsteps 10 hpcmode 1 hpc 2 3 1 keepll nonlocalext lstype 4 smtype 7
metastep 5 nsteps 10 stiffMode 2 controllmode 0 psi 0.0 steplength 0.001 rtolv 0.0001 MaxIter 200 reqIterations 30 manrmsteps 10 hpcmode 1 hpc 2 3 1 keepll nonlocalext lstype 4 smtype 7
#
domain PlaneStrain
OutputManager tstep_all dofman_all element_all
2 2 2 1
SimpleCS 1 thick 1.0 width 1.0
SimpleCS 2 thick 1.0 width 1.0
idmnl1 1 d 1.0 E 37.7e3 n 0.2 talpha 0.0 e0 6.75e-5 ef 3.0e-4 R 5.0 equivstraintype 1 regionMap 2 0 1
IsoLE  2 d 1.0 E 300.e3 n 0.2 talpha 0.0
##idmnl1 1 d 1.0 E 30.e3 n 0.18 talpha 0.0 e0 0.00011666666667 ef 0.0059953333333 R 5.0 equivstraintype 1 regionMap 2 0 1
##IsoLE  2 d 1.0 E 210.e3 n 0.3 talpha 0.0
BoundaryCondition 1 loadTimeFunction 1 prescribedvalue 0.0 
ConstantEdgeLoad 2 loadTimeFunction 1 components 2 0 -1.0 loadType 3 ndofs 2
ConstantFunction 1 f(t) 1.0

#30GPa = 30e9N/m2 = 30e9N/1e6mm2 = 30e3N/mm2
#damage scale 0.05 - 0.95

#%BEGIN_CHECK% tolerance 1.e-3
#LOADLEVEL
#NODE number 7 dof 2 unknown d 
#NODE number 3 dof 1 unknown d 
#%END_CHECK%

vertex 1
nodeprop bc 2 1 1

vertex 2 3
curve 1
nodeprop bc 2 0 1

vertex 5 7
curve 4 7
nodeprop bc 2 1 0

curve 9
boundaryload 2

patch 1
elemprop mat 1 crossSect 1
triatype TrPlaneStrain2d
quadtype Quad1PlaneStrain

patch 2
elemprop mat 2 crossSect 2
triatype TrPlaneStrain2d
quadtype Quad1PlaneStrain

property 0
width 5.0