Glassvision/JuliaFEM.jl
1
0
Fork 0
mirror of https://github.com/JuliaFEM/JuliaFEM.jl.git synced 2026-09-19 09:54:55 +00:00
Code Issues Packages Projects Releases Wiki Activity

volume load

Browse Source
This commit is contained in:
Jukka Aho
2015-10-23 01:11:45 +03:00
parent 86e25d5c4e
commit 2c11b2b2b2
8 changed files with 1006 additions and 126 deletions
Download Patch File Download Diff File Expand all files Collapse all files
notebooks/2015-06-25-elasticity-solver-example.ipynb
+178 -126
View File
@@ -4,45 +4,39 @@
"cell_type": "markdown",
"metadata": {},
"source": [
"# Solving elasticity problems using JuliaFEM\n",
"# Solving elasticity equations using JuliaFEM\n",
"\n",
"Author(s): Jukka Aho\n",
"\n",
"**Abstract**: Elasticity equations design notes."
"**Abstract**: Elasticity equations design notes.\n",
"\n",
"##Weak form\n",
"\n",
"Given function spaces\n",
"\\begin{align}\n",
"\\boldsymbol{\\mathcal{U}} & =\\left\\{ \\boldsymbol{u}\\in H^{1}\\left(\\Omega\\right)|\\boldsymbol{u}\\left(\\boldsymbol{X},t\\right)=\\hat{\\boldsymbol{u}}\\left(\\boldsymbol{X},t\\right)\\text{ on }\\Gamma_{\\mathrm{u}}\\right\\} ,\\\\\n",
"\\boldsymbol{\\mathcal{V}} & =\\left\\{ \\delta\\boldsymbol{u}\\in H^{1}\\left(\\Omega\\right)|\\delta\\boldsymbol{u}\\left(\\boldsymbol{X}\\right)=0\\text{ on }\\Gamma_{\\mathrm{u}}\\right\\} ,\n",
"\\end{align}\n",
"find $\\boldsymbol{u}\\in\\boldsymbol{\\mathcal{U}}$ such that\n",
"\\begin{equation}\n",
"\\delta\\mathcal{W}:=\\int_{\\Omega_{0}}\\rho_{0}\\ddot{\\boldsymbol{u}}\\cdot\\delta\\boldsymbol{u}\\,\\mathrm{d}V_{0}+\\int_{\\Omega_{0}}\\boldsymbol{S}:\\delta\\boldsymbol{E}\\,\\mathrm{d}V_{0}-\\int_{\\Omega_{0}}\\hat{\\boldsymbol{b}}_{0}\\cdot\\delta\\boldsymbol{u}\\,\\mathrm{d}V_{0}-\\int_{\\Gamma_{\\sigma}}\\hat{\\boldsymbol{t}}_{0}\\cdot\\delta\\boldsymbol{u}\\,\\mathrm{d}A_{0} =0 \\qquad\\forall\\delta\\boldsymbol{u}\\in\\boldsymbol{\\mathcal{V}}\n",
"\\end{equation}"
]
},
{
"cell_type": "code",
"execution_count": 1,
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING: Base.String is deprecated, use AbstractString instead.\n",
" likely near /home/jukka/.julia/v0.4/Logging/src/Logging.jl:24\n",
"WARNING: Base.String is deprecated, use AbstractString instead.\n",
" likely near /home/jukka/.julia/v0.4/Logging/src/Logging.jl:24\n",
"WARNING: Base.String is deprecated, use AbstractString instead.\n",
" likely near /home/jukka/.julia/v0.4/Logging/src/Logging.jl:24\n",
"WARNING: Base.String is deprecated, use AbstractString instead.\n",
" likely near /home/jukka/.julia/v0.4/Logging/src/Logging.jl:40\n",
"WARNING: Base.String is deprecated, use AbstractString instead.\n",
" likely near /home/jukka/.julia/v0.4/Logging/src/Logging.jl:40\n",
"WARNING: Base.String is deprecated, use AbstractString instead.\n",
" likely near /home/jukka/.julia/v0.4/Logging/src/Logging.jl:40\n"
]
},
{
"data": {
"text/plain": [
"Logger(root,DEBUG,Base.PipeEndpoint(open, 0 bytes waiting),root)"
]
},
"execution_count": 1,
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
@@ -58,7 +52,7 @@
},
{
"cell_type": "code",
"execution_count": 2,
"execution_count": 3,
"metadata": {
"collapsed": false
},
@@ -69,7 +63,7 @@
"CPS4"
]
},
"execution_count": 2,
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
@@ -107,7 +101,7 @@
},
{
"cell_type": "code",
"execution_count": 40,
"execution_count": 4,
"metadata": {
"collapsed": false
},
@@ -115,79 +109,58 @@
{
"data": {
"text/plain": [
"8-element Array{Float64,1}:\n",
" -189.267 \n",
" 26.1073 \n",
" 5.7554e-13 \n",
" 2.6894e-12 \n",
" -8.73968e-13\n",
" 20.0 \n",
" 189.267 \n",
" -46.1073 "
"has_rhs (generic function with 4 methods)"
]
},
"execution_count": 40,
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"function get_lhs_and_rhs(equation::CPS4, ip, time)\n",
" # boilerplate code start\n",
" element = get_element(equation)\n",
" N = FEM.get_basis(element)\n",
" dN = FEM.diff(N)(ip.xi)\n",
" geometry = element[\"geometry\"](time)\n",
" basis = FEM.get_basis(element)\n",
" dbasis = FEM.diff(basis)(ip.xi)\n",
" grad(u) = dbasis*u*inv(dbasis*geometry)\n",
" # boilerplate code end -- replace with a macro?\n",
"\n",
" # fields at time \n",
" X = element[\"geometry\"](time)\n",
" u = element[\"displacement\"](time)\n",
" young = element[\"young\"](time)\n",
" poisson = element[\"poisson\"](time)\n",
" # interpolate fields in temporal dimension\n",
" young = element[\"youngs modulus\"](time)\n",
" poisson = element[\"poissons ratio\"](time)\n",
" displacement = element[\"displacement\"](time)\n",
"\n",
" # material\n",
" young = interpolate(N, young, ip)\n",
" poisson = interpolate(N, poisson, ip)\n",
" # interpolate material in spatial dimension\n",
" young = interpolate(basis, young, ip)\n",
" poisson = interpolate(basis, poisson, ip)\n",
" mu = young/(2*(1+poisson))\n",
" lambda = young*poisson/((1+poisson)*(1-2*poisson))\n",
" lambda = 2*lambda*mu/(lambda + 2*mu)\n",
" lambda = 2*lambda*mu/(lambda + 2*mu) # <- correction for 2d\n",
"\n",
" gradw = dN*inv(dN*X)\n",
" grad(u) = dN*u*inv(dN*(X+u))\n",
" Grad(u) = dN*u*inv(dN*X)\n",
"\n",
" function calc_R(data::Vector)\n",
" # create new field similar to displacement and fill it with data\n",
" f = similar(u, data)\n",
" # kinematics\n",
" F = I + Grad(f) # deformation gradient\n",
" C = F'*F\n",
" E = 1/2*(C - I) # strain\n",
" function W(data::Vector)\n",
" Wint = 0.0\n",
" # create new field u, similar to field displacement, and fill it with data\n",
" u = similar(displacement, data)\n",
" F = I + grad(u) # deformation gradient\n",
" E = 1/2*(F'*F - I) # strain\n",
" S = 2*mu*E + lambda*trace(E)*I # stress\n",
" P = F*S\n",
" return (P*gradw')[:]\n",
" Wint += 1/2*trace(S*E')\n",
"\n",
" Wext = 0.0\n",
" # any volume load?\n",
" if haskey(element, \"displacement volume load\")\n",
" b = interpolate(element, \"displacement volume load\", ip, time)\n",
" δu = interpolate(basis, u, ip)\n",
" Wext += dot(b, δu)\n",
" end\n",
" return Wint - Wext\n",
" end\n",
"\n",
" function calc_Wint(data::Vector)\n",
" # create new field similar to displacement and fill it with data\n",
" δu = similar(u, data)\n",
" # kinematics\n",
" F = I + Grad(δu) # deformation gradient\n",
" C = F'*F\n",
" E = 1/2*(C - I) # strain\n",
" #E = 1/2*(Grad(δu)' + Grad(δu) + Grad(δu)'*Grad(δu))\n",
" S = 2*mu*E + lambda*trace(E)*I # stress\n",
" δF = Grad(δu)\n",
" #δɛ = 1/2*(δF' + δF)\n",
" #δE = F'*δɛ*F\n",
" #δE = 1/2*(F'*δF + δF'*F)\n",
" #δE = 1/2*(Grad(δu)' + Grad(δu) + Grad(δu)'*Grad(δu))\n",
" δE = 1/2*E # ...???? energy 2x\n",
"\n",
" return trace(S*δE')\n",
" end\n",
"\n",
" R = ForwardDiff.gradient(calc_Wint, u[:])\n",
" Kt = ForwardDiff.hessian(calc_Wint, u[:])\n",
" return Kt, -R\n",
" R = ForwardDiff.gradient(W, displacement[:])\n",
" K = ForwardDiff.hessian(W, displacement[:])\n",
" return K, -R\n",
"end\n",
"\n",
"function JuliaFEM.get_lhs(equation::CPS4, ip, time)\n",
@@ -198,14 +171,12 @@
"end\n",
"\n",
"JuliaFEM.has_lhs(equation::CPS4) = true\n",
"JuliaFEM.has_rhs(equation::CPS4) = true\n",
"\n",
"JuliaFEM.integrate_rhs(equation, Inf)"
"JuliaFEM.has_rhs(equation::CPS4) = true"
]
},
{
"cell_type": "code",
"execution_count": 41,
"execution_count": 5,
"metadata": {
"collapsed": false
},
@@ -214,19 +185,90 @@
"name": "stdout",
"output_type": "stream",
"text": [
"testing primary field with point load versus code aster solution\n",
"increment 1, norm = 5.77653, Wint = 0.000, Wext = 81.525, |Wint-Wext| = 81.52516\n",
"increment 2, norm = 0.99988, Wint = 173.895, Wext = 79.037, |Wint-Wext| = 94.85823\n",
"increment 3, norm = 0.28354, Wint = 87.084, Wext = 82.175, |Wint-Wext| = 4.90953\n",
"increment 4, norm = 0.07071, Wint = 82.548, Wext = 83.100, |Wint-Wext| = 0.55141\n",
"increment 5, norm = 0.00082, Wint = 83.116, Wext = 83.109, |Wint-Wext| = 0.00644\n",
"increment 6, norm = 0.00000, Wint = 83.109, Wext = 83.109, |Wint-Wext| = 0.00000\n",
"increment 7, norm = 0.00000, Wint = 83.109, Wext = 83.109, |Wint-Wext| = 0.00000\n",
"increment 8, norm = 0.00000, Wint = 83.109, Wext = 83.109, |Wint-Wext| = 0.00000\n",
"increment 9, norm = 0.00000, Wint = 83.109, Wext = 83.109, |Wint-Wext| = 0.00000\n",
"increment 10, norm = 0.00000, Wint = 83.109, Wext = 83.109, |Wint-Wext| = 0.00000\n",
"elapsed time: 4.257759666 seconds\n",
"1 fact verified.\n"
]
},
{
"data": {
"text/plain": [
"delayed_handler (generic function with 4 methods)"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"facts(\"testing primary field with point load versus code aster solution\") do\n",
" element = Quad4([1, 2, 3, 4])\n",
" push!(element, FieldSet(\"geometry\", [Field(0.0, Vector[[0.0, 0.0], [10.0, 0.0], [10.0, 1.0], [0.0, 1.0]])]))\n",
" push!(element, FieldSet(\"youngs modulus\", [Field(0.0, 500.0)]))\n",
" push!(element, FieldSet(\"poissons ratio\", [Field(0.0, 0.3)]))\n",
" equation = CPS4(element)\n",
"\n",
" u0 = Field(0.0, Vector[[0.0, 0.0], [0.0, 0.0], [0.0, 0.0], [0.0, 0.0]])\n",
" push!(element[\"displacement\"], u0)\n",
" u = zeros(8)\n",
" du = zeros(8)\n",
" fd = [3, 4, 5, 6]\n",
" f = zeros(8)\n",
" f[6] = -20.0\n",
" tic()\n",
" for i=1:10\n",
" A = JuliaFEM.integrate_lhs(equation, 1.0)\n",
" b = JuliaFEM.integrate_rhs(equation, 1.0)\n",
" du[fd] = A[fd,fd] \\ (b[fd]+f[fd])\n",
" u += du\n",
" new_field = similar(u0, u)\n",
" new_field.time = 1.0\n",
" new_field.increment = i\n",
" push!(element[\"displacement\"], new_field)\n",
" Wint = (-b'*u)[1]\n",
" Wext = (f'*u)[1]\n",
" W = abs(Wint-Wext)\n",
" @printf(\"increment %2d, norm = %8.5f, Wint = %8.3f, Wext = %8.3f, |Wint-Wext| = %8.5f\\n\", i, norm(du), Wint, Wext, W)\n",
" end\n",
" toc()\n",
" # verified using Code Aster.\n",
" @fact interpolate(element, \"displacement\", [1.0, 1.0], Inf)[2] --> roughly(-4.15546385452579E+00)\n",
"end"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"testing primary field with volume load versus code aster solution\n",
"increment "
]
},
{
"data": {
"text/plain": [
"Success :: (line:-1) :: fact was true\n",
" Expression: (interpolate(element,\"displacement\",[1.0,1.0],Inf))[2] --> roughly(-4.15546385452579)\n",
" Expected: -4.15546385452579\n",
" Occurred: -4.155463854525783"
"delayed_handler (generic function with 4 methods)"
]
},
"execution_count": 41,
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
@@ -234,47 +276,57 @@
"name": "stdout",
"output_type": "stream",
"text": [
"1, norm = 5.77653, Wint = 0.000, Wext = 81.525, |Wint-Wext| = 81.52516\n",
"increment 2, norm = 0.99988, Wint = 173.895, Wext = 79.037, |Wint-Wext| = 94.85823\n",
"increment 3, norm = 0.28354, Wint = 87.084, Wext = 82.175, |Wint-Wext| = 4.90953\n",
"increment 4, norm = 0.07071, Wint = 82.548, Wext = 83.100, |Wint-Wext| = 0.55141\n",
"increment 5, norm = 0.00082, Wint = 83.116, Wext = 83.109, |Wint-Wext| = 0.00644\n",
"increment 6, norm = 0.00000, Wint = 83.109, Wext = 83.109, |Wint-Wext| = 0.00000\n",
"83.10927709051566\n"
"1, norm = 14.44128, Wint = -509.167, Wext = 0.000, |Wint-Wext| = 509.16667\n",
"increment 2, norm = 4.01742, Wint = 4335.728, Wext = 0.000, |Wint-Wext| = 4335.72801\n",
"increment 3, norm = 1.54645, Wint = 705.913, Wext = 0.000, |Wint-Wext| = 705.91292\n",
"increment 4, norm = 1.12361, Wint = 60.492, Wext = 0.000, |Wint-Wext| = 60.49208\n",
"increment 5, norm = 0.79119, Wint = -1.486, Wext = 0.000, |Wint-Wext| = 1.48555\n",
"increment 6, norm = 0.12733, Wint = 5.743, Wext = 0.000, |Wint-Wext| = 5.74331\n",
"increment 7, norm = 0.00725, Wint = 0.080, Wext = 0.000, |Wint-Wext| = 0.07997\n",
"increment 8, norm = 0.00001, Wint = 0.000, Wext = 0.000, |Wint-Wext| = 0.00045\n",
"increment 9, norm = 0.00000, Wint = 0.000, Wext = 0.000, |Wint-Wext| = 0.00000\n",
"increment 10, norm = 0.00000, Wint = 0.000, Wext = 0.000, |Wint-Wext| = 0.00000\n",
"elapsed time: 0.068575305 seconds\n",
"1 fact verified.\n"
]
}
],
"source": [
"element = Quad4([1, 2, 3, 4])\n",
"push!(element, FieldSet(\"young\", [Field(0.0, 500.0)]))\n",
"push!(element, FieldSet(\"geometry\", [Field(0.0, Vector[[0.0, 0.0], [10.0, 0.0], [10.0, 1.0], [0.0, 1.0]])]))\n",
"push!(element, FieldSet(\"poisson\", [Field(0.0, 0.3)]))\n",
"equation = CPS4(element)\n",
"facts(\"testing primary field with volume load versus code aster solution\") do\n",
" element = Quad4([1, 2, 3, 4])\n",
" push!(element, FieldSet(\"geometry\", [Field(0.0, Vector[[0.0, 0.0], [10.0, 0.0], [10.0, 1.0], [0.0, 1.0]])]))\n",
" push!(element, FieldSet(\"youngs modulus\", [Field(0.0, 500.0)]))\n",
" push!(element, FieldSet(\"poissons ratio\", [Field(0.0, 0.3)]))\n",
" push!(element, FieldSet(\"displacement volume load\",\n",
" [Field(0.0, Vector[[0.0, -10.0], [0.0, -10.0], [0.0, -10.0], [0.0, -10.0]])]))\n",
" equation = CPS4(element)\n",
"\n",
"u0 = Field(0.0, Vector[[0.0, 0.0], [0.0, 0.0], [0.0, 0.0], [0.0, 0.0]])\n",
"push!(element[\"displacement\"], u0)\n",
"u = zeros(8)\n",
"du = zeros(8)\n",
"fd = [3, 4, 5, 6]\n",
"f = zeros(8)\n",
"f[6] = -20.0\n",
"for i=1:6\n",
" A = JuliaFEM.integrate_lhs(equation, 1.0)\n",
" b = JuliaFEM.integrate_rhs(equation, 1.0)\n",
" du[fd] = A[fd,fd] \\ (b[fd]+f[fd])\n",
" u += du\n",
" new_field = similar(u0, u)\n",
" new_field.time = 1.0\n",
" new_field.increment = i\n",
" push!(element[\"displacement\"], new_field)\n",
" Wint = (-b'*u)[1]\n",
" Wext = (f'*u)[1]\n",
" W = abs(Wint-Wext)\n",
" @printf(\"increment %2d, norm = %8.5f, Wint = %8.3f, Wext = %8.3f, |Wint-Wext| = %8.5f\\n\", i, norm(du), Wint, Wext, W)\n",
"end\n",
"println((f'*u)[1])\n",
"@fact interpolate(element, \"displacement\", [1.0, 1.0], Inf)[2] --> roughly(-4.15546385452579E+00)\n",
"#interpolate(element, \"displacement\", [1.0, 1.0], Inf)[2], (f'*u)[1]"
" u0 = Field(0.0, Vector[[0.0, 0.0], [0.0, 0.0], [0.0, 0.0], [0.0, 0.0]])\n",
" push!(element[\"displacement\"], u0)\n",
" u = zeros(8)\n",
" du = zeros(8)\n",
" fd = [3, 4, 5, 6]\n",
" f = zeros(8)\n",
" f[6] = -20.0*0\n",
" tic()\n",
" for i=1:10\n",
" A = JuliaFEM.integrate_lhs(equation, 1.0)\n",
" b = JuliaFEM.integrate_rhs(equation, 1.0)\n",
" du[fd] = A[fd,fd] \\ (b[fd]+f[fd])\n",
" u += du\n",
" new_field = similar(u0, u)\n",
" new_field.time = 1.0\n",
" new_field.increment = i\n",
" push!(element[\"displacement\"], new_field)\n",
" Wint = (-b'*u)[1]\n",
" Wext = (f'*u)[1]\n",
" W = abs(Wint-Wext)\n",
" @printf(\"increment %2d, norm = %8.5f, Wint = %8.3f, Wext = %8.3f, |Wint-Wext| = %8.5f\\n\", i, norm(du), Wint, Wext, W)\n",
" end\n",
" toc()\n",
" # verified using Code Aster.\n",
" @fact interpolate(element, \"displacement\", [1.0, 1.0], Inf)[2] --> roughly(-8.77303119819776E+00)\n",
"end"
]
},
{
src/elements.jl
+5
View File
@@ -186,6 +186,11 @@ function dinterpolate(element::Element, field_name, ip::IntegrationPoint, time::
dinterpolate(element, field_name, ip.xi, time)
end
""" Check does fieldset exist. """
function Base.haskey(element::Element, what)
haskey(element.fields, symbol(what))
end
"""
Get jacobian of element evaluated at point ξ on element in reference configuration.
verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.comm
+56
View File
@@ -0,0 +1,56 @@
DEBUT()
MAIL = LIRE_MAILLAGE()
MO = AFFE_MODELE(MAILLAGE = MAIL,
AFFE = _F(MAILLE='E1', PHENOMENE='MECANIQUE', MODELISATION='C_PLAN')
)
MAT = DEFI_MATERIAU(ELAS = _F(E = 500.0, NU = 0.3))
CHMAT = AFFE_MATERIAU(
MAILLAGE = MAIL,
AFFE = _F(MAILLE = 'E1', MATER = MAT))
BC = AFFE_CHAR_MECA(
MODELE = MO,
DDL_IMPO = (_F(NOEUD = ('N1','N4'), DX=0, DY=0)))
LO = AFFE_CHAR_MECA(
MODELE = MO,
FORCE_NODALE = _F(NOEUD='N3', FY=-20.0))
LIST = DEFI_LIST_REEL(
DEBUT = 0,
INTERVALLE = _F(JUSQU_A=1.0, NOMBRE=1))
STEP = DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0,0,1,1))
RESU = STAT_NON_LINE(
MODELE=MO,
CHAM_MATER=CHMAT,
EXCIT=(
_F(CHARGE=BC),
_F(CHARGE=LO, FONC_MULT=STEP)),
NEWTON=_F(REAC_INCR=1, MATRICE='TANGENTE', REAC_ITER=1),
COMP_ELAS=_F(DEFORMATION='GROT_GDEP'),
INCREMENT=_F(LIST_INST=LIST),
CONVERGENCE=_F(RESI_GLOB_RELA=1.0E-12))
RESU = CALC_CHAMP(
reuse=RESU,
RESULTAT=RESU,
MODELE=MO,
DEFORMATION=('EPSI_ELNO', 'EPSG_ELNO', 'EPSI_ELGA', 'EPSG_ELGA'),
ENERGIE=('ENEL_ELEM', 'ENEL_NOEU', 'ETOT_ELEM', 'ETOT_NOEU'),
FORCE=('FORC_NODA'),
CONTRAINTE='SIGM_ELNO')
IMPR_RESU(
MODELE = MO,
FORMAT = 'RESULTAT',
RESU = _F(RESULTAT = RESU))
FIN()
verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.mail
+13
View File
@@ -0,0 +1,13 @@
COOR_2D
N1 0.0 0.0
N2 10.0 0.0
N3 10.0 1.0
N4 0.0 1.0
FINSF
QUAD4
E1 N1 N2 N3 N4
FINSF
FIN
verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.resu
+343
View File
@@ -0,0 +1,343 @@
-- CODE_ASTER -- VERSION : EXPLOITATION (stable) --
Version 11.4.0 du 05/06/2013
Copyright EDF R&D 1991 - 2015
Exécution du : Thu Oct 22 22:34:12 2015
Nom de la machine : jukka-desktop
Architecture : 64bit
Type de processeur : x86_64
Système d'exploitation : Linux 3.13.0-66-generic
Langue des messages : en (UTF-8)
!------------------------------------------------------------------------------------!
! <A> <SUPERVIS2_2> !
! !
! Vous utilisez une vieille version de Code_Aster. !
! !
! En mettant à jour votre version, vous bénéficierez des dernières améliorations !
! apportées au code depuis 15 mois. !
! Si vous avez des développements privés, vous risquez d'avoir un travail !
! important de portage si vous ne suivez pas les mises à jour. !
! !
! !
! Ceci est une alarme. Si vous ne comprenez pas le sens de cette !
! alarme, vous pouvez obtenir des résultats inattendus ! !
!------------------------------------------------------------------------------------!
Parallélisme MPI : inactif
Parallélisme OpenMP : actif
Nombre de processus utilisés : 1
Version de la librairie HDF5 : 1.8.8
Version de la librairie MED : 3.0.6
Librairie MUMPS : installée
Version de la librairie SCOTCH : 5.1.10
Mémoire limite pour l'exécution : 4096.00 Mo
consommée par l'initialisation : 197.46 Mo
par les objets du jeu de commandes : 1.53 Mo
reste pour l'allocation dynamique : 3897.01 Mo
Taille limite des fichiers d'échange : 48.00 Go
--------------------------------------------------------------------------------
ASTER 11.04.00 CONCEPT RESU CALCULE LE 22/10/2015 A 22:34:12 DE TYPE EVOL_NOLI
======>
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE DEPL
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
NOEUD DX DY
N1 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE EPSI_ELGA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE EPSI_ELNO
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
N1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE SIEF_ELGA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 SIXX SIYY SIZZ SIXY
1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE SIGM_ELNO
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 SIXX SIYY SIZZ SIXY
N1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE FORC_NODA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
NOEUD DX DY
N1 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE EPSG_ELGA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE EPSG_ELNO
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
N1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE VARI_ELGA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 VARI
1 0.00000000000000E+00
2 0.00000000000000E+00
3 0.00000000000000E+00
4 0.00000000000000E+00
------>
CARTE DE NOM SYMBOLIQUE COMPORTEMENT
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE ENEL_NOEU
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
NOEUD TOTALE
N1 0.00000000000000E+00
N2 0.00000000000000E+00
N3 0.00000000000000E+00
N4 0.00000000000000E+00
------>
CHAMP PAR ELEMENT CONSTANT SUR L'ELEMENT DE NOM SYMBOLIQUE ENEL_ELEM
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 TOTALE
0.00000000000000E+00
------>
CHAMP PAR ELEMENT CONSTANT SUR L'ELEMENT DE NOM SYMBOLIQUE ETOT_ELEM
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 TOTALE
0.00000000000000E+00
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE ETOT_NOEU
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
NOEUD TOTALE
N1 0.00000000000000E+00
N2 0.00000000000000E+00
N3 0.00000000000000E+00
N4 0.00000000000000E+00
======>
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE DEPL
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
NOEUD DX DY
N1 -5.15692621618630E-24 3.92910790962688E-24
N2 -1.13385861994525E+00 -4.00638800902936E+00
N3 -5.36656635814853E-01 -4.15546385452579E+00
N4 8.27180612553028E-25 4.13590306276514E-25
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE EPSI_ELGA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
1 -1.00765499102918E-01 -3.15034329846977E-02 5.66866851804068E-02 -1.38792757642670E-01
2 -1.00765499102918E-01 -1.17572412511733E-01 9.35733906919935E-02 3.36046055064611E-02
3 -6.62860264730919E-02 -1.17572412511733E-01 7.87964738506394E-02 2.93011565301093E-02
4 -6.62860264730919E-02 -3.15034329846977E-02 4.19097683390527E-02 -1.43096206619022E-01
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE EPSI_ELNO
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
N1 -1.13385861994525E-01 2.42861286636753E-17 4.85939408547963E-02 -2.00319400451468E-01
N2 -1.13385861994525E-01 -1.49075845496431E-01 1.12483588924695E-01 9.82815916137292E-02
N3 -5.36656635814853E-02 -1.49075845496431E-01 8.68892181762499E-02 9.08277993389076E-02
N4 -5.36656635814853E-02 2.77555756156289E-17 2.29995701063508E-02 -2.07773192726290E-01
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE SIEF_ELGA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 SIXX SIYY SIZZ SIXY
1 -2.31942215770259E+01 2.75472596932288E+01 -1.95740315902075E-15 -4.42065108646296E+01
2 4.24930377564788E+00 -1.23108782596206E+01 -4.54133574898300E-16 9.80137080537897E+00
3 2.34925742118825E+01 -5.80349743042289E+00 0.00000000000000E+00 6.17854003054832E+00
4 -6.46229422958107E+00 3.75333879264263E+01 0.00000000000000E+00 -5.24880405947485E+01
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE SIGM_ELNO
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 SIXX SIYY SIZZ SIXY
N1 -3.90271033126679E+01 3.80151103339302E+01 -3.42549723273148E-15 -6.03193699855925E+01
N2 6.91434921629353E+00 -2.88157722147076E+01 1.31276792038702E-16 3.02715282869615E+01
N3 4.18367990363141E+01 -1.97500954751241E+01 -3.51755104117185E-17 2.69500981067994E+01
N4 -1.16386827590164E+01 5.75170292855130E+01 9.17859217185445E-16 -7.76168970316193E+01
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE FORC_NODA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
NOEUD DX DY
N1 1.87539318526141E+02 -2.60102307081623E+01
N2 -5.27696481572235E-01 1.66617755509419E-01
N3 6.36360553319175E-01 -1.99738176845921E+01
N4 -1.87647982597888E+02 4.58174306372450E+01
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE EPSG_ELGA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
1 -1.41658197107361E-02 -2.30435218639950E-02 1.59468606748847E-02 -1.38790871348685E-01
2 -1.41658197107361E-02 2.58948887482680E-04 5.96008749568003E-03 3.36116452514522E-02
3 2.09461399718909E-02 2.58948887482680E-04 -9.08789522544582E-03 3.79340506020240E-02
4 2.09461399718909E-02 -2.30435218639950E-02 8.98877953758883E-04 -1.40783029624109E-01
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE EPSG_ELNO
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
N1 -2.70176889312325E-02 -3.15728181299796E-02 2.51102173119481E-02 -2.00319400451468E-01
N2 -2.70176889312325E-02 8.78824515346735E-03 7.81261876189937E-03 9.42872447276373E-02
N3 3.37980091923874E-02 8.78824515346735E-03 -1.82512518625092E-02 1.05777143330804E-01
N4 3.37980091923874E-02 -3.15728181299796E-02 -9.53653312460457E-04 -2.07773192726290E-01
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE VARI_ELGA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 VARI
1 0.00000000000000E+00
2 0.00000000000000E+00
3 0.00000000000000E+00
4 0.00000000000000E+00
------>
CARTE DE NOM SYMBOLIQUE COMPORTEMENT
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE ENEL_NOEU
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
NOEUD TOTALE
N1 8.11427653502166E+00
N2 -1.74924544273778E+00
N3 -2.26851577595122E+00
N4 1.26410781912352E+01
------>
CHAMP PAR ELEMENT CONSTANT SUR L'ELEMENT DE NOM SYMBOLIQUE ENEL_ELEM
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 TOTALE
4.18439837689196E+01
------>
CHAMP PAR ELEMENT CONSTANT SUR L'ELEMENT DE NOM SYMBOLIQUE ETOT_ELEM
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 TOTALE
3.64664959689712E+01
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE ETOT_NOEU
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
NOEUD TOTALE
N1 8.79924750872253E+00
N2 -7.85578326392019E-01
N3 -3.54443135219698E+00
N4 1.01173605574550E+01
<I> <FIN> FERMETURE DE LA BASE "GLOBALE" EFFECTUEE.
<FIN> Arrêt normal dans "FIN".
<I> <FIN> ARRET NORMAL DANS "FIN" PAR APPEL A "JEFINI".
<I> <FIN> MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION : 20.99 Mo
<I> <FIN> MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION : 27.32 Mo
<I> <FIN> MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS LORS DE L'EXECUTION : 226.65 Mo
********************************************************************************
* COMMAND : USER : SYSTEM : USER+SYS : ELAPSED *
********************************************************************************
* init (jdc) : 0.15 : 0.02 : 0.17 : 0.17 *
* . compile : 0.00 : 0.00 : 0.00 : 0.01 *
* . exec_compile : 0.06 : 0.00 : 0.06 : 0.06 *
* . report : 0.00 : 0.00 : 0.00 : 0.00 *
* . build : 0.00 : 0.00 : 0.00 : 0.00 *
* DEBUT : 0.02 : 0.00 : 0.02 : 0.03 *
* LIRE_MAILLAGE : 0.00 : 0.00 : 0.00 : 0.00 *
* AFFE_MODELE : 0.01 : 0.00 : 0.01 : 0.00 *
* DEFI_MATERIAU : 0.00 : 0.00 : 0.00 : 0.01 *
* AFFE_MATERIAU : 0.00 : 0.00 : 0.00 : 0.00 *
* AFFE_CHAR_MECA : 0.01 : 0.00 : 0.01 : 0.00 *
* AFFE_CHAR_MECA : 0.00 : 0.00 : 0.00 : 0.00 *
* DEFI_LIST_REEL : 0.00 : 0.00 : 0.00 : 0.01 *
* DEFI_FONCTION : 0.00 : 0.00 : 0.00 : 0.00 *
* STAT_NON_LINE : 0.07 : 0.00 : 0.07 : 0.07 *
* CALC_CHAMP : 0.03 : 0.00 : 0.03 : 0.03 *
* IMPR_RESU : 0.01 : 0.00 : 0.01 : 0.01 *
* FIN : 0.01 : 0.01 : 0.02 : 0.02 *
* . part Superviseur : 0.18 : 0.02 : 0.20 : 0.23 *
* . part Fortran : 0.14 : 0.01 : 0.15 : 0.13 *
********************************************************************************
* TOTAL_JOB : 0.32 : 0.03 : 0.35 : 0.36 *
********************************************************************************
verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.comm
+55
View File
@@ -0,0 +1,55 @@
DEBUT()
MAIL = LIRE_MAILLAGE()
MO = AFFE_MODELE(MAILLAGE = MAIL,
AFFE = _F(MAILLE='E1', PHENOMENE='MECANIQUE', MODELISATION='C_PLAN'))
MAT = DEFI_MATERIAU(ELAS = _F(E = 500.0, NU = 0.3))
CHMAT = AFFE_MATERIAU(
MAILLAGE = MAIL,
AFFE = _F(MAILLE = 'E1', MATER = MAT))
BC = AFFE_CHAR_MECA(
MODELE = MO,
DDL_IMPO = (_F(NOEUD = ('N1','N4'), DX=0, DY=0)))
LO = AFFE_CHAR_MECA(
MODELE = MO,
FORCE_INTERNE = _F(MAILLE='E1', FY=-10.0))
LIST = DEFI_LIST_REEL(
DEBUT = 0,
INTERVALLE = _F(JUSQU_A=1.0, NOMBRE=1))
STEP = DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0,0,1,1))
RESU = STAT_NON_LINE(
MODELE=MO,
CHAM_MATER=CHMAT,
EXCIT=(
_F(CHARGE=BC),
_F(CHARGE=LO, FONC_MULT=STEP)),
NEWTON=_F(REAC_INCR=1, MATRICE='TANGENTE', REAC_ITER=1),
COMP_ELAS=_F(DEFORMATION='GROT_GDEP'),
INCREMENT=_F(LIST_INST=LIST),
CONVERGENCE=_F(RESI_GLOB_RELA=1.0E-12))
RESU = CALC_CHAMP(
reuse=RESU,
RESULTAT=RESU,
MODELE=MO,
DEFORMATION=('EPSI_ELNO', 'EPSG_ELNO', 'EPSI_ELGA', 'EPSG_ELGA'),
ENERGIE=('ENEL_ELEM', 'ENEL_NOEU', 'ETOT_ELEM', 'ETOT_NOEU'),
FORCE=('FORC_NODA'),
CONTRAINTE='SIGM_ELNO')
IMPR_RESU(
MODELE = MO,
FORMAT = 'RESULTAT',
RESU = _F(RESULTAT = RESU))
FIN()
verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.mail
+13
View File
@@ -0,0 +1,13 @@
COOR_2D
N1 0.0 0.0
N2 10.0 0.0
N3 10.0 1.0
N4 0.0 1.0
FINSF
QUAD4
E1 N1 N2 N3 N4
FINSF
FIN
verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.resu
+343
View File
@@ -0,0 +1,343 @@
-- CODE_ASTER -- VERSION : EXPLOITATION (stable) --
Version 11.4.0 du 05/06/2013
Copyright EDF R&D 1991 - 2015
Exécution du : Thu Oct 22 22:40:02 2015
Nom de la machine : jukka-desktop
Architecture : 64bit
Type de processeur : x86_64
Système d'exploitation : Linux 3.13.0-66-generic
Langue des messages : en (UTF-8)
!------------------------------------------------------------------------------------!
! <A> <SUPERVIS2_2> !
! !
! Vous utilisez une vieille version de Code_Aster. !
! !
! En mettant à jour votre version, vous bénéficierez des dernières améliorations !
! apportées au code depuis 15 mois. !
! Si vous avez des développements privés, vous risquez d'avoir un travail !
! important de portage si vous ne suivez pas les mises à jour. !
! !
! !
! Ceci est une alarme. Si vous ne comprenez pas le sens de cette !
! alarme, vous pouvez obtenir des résultats inattendus ! !
!------------------------------------------------------------------------------------!
Parallélisme MPI : inactif
Parallélisme OpenMP : actif
Nombre de processus utilisés : 1
Version de la librairie HDF5 : 1.8.8
Version de la librairie MED : 3.0.6
Librairie MUMPS : installée
Version de la librairie SCOTCH : 5.1.10
Mémoire limite pour l'exécution : 4096.00 Mo
consommée par l'initialisation : 197.46 Mo
par les objets du jeu de commandes : 1.52 Mo
reste pour l'allocation dynamique : 3897.02 Mo
Taille limite des fichiers d'échange : 48.00 Go
--------------------------------------------------------------------------------
ASTER 11.04.00 CONCEPT RESU CALCULE LE 22/10/2015 A 22:40:02 DE TYPE EVOL_NOLI
======>
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE DEPL
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
NOEUD DX DY
N1 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE EPSI_ELGA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE EPSI_ELNO
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
N1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE SIEF_ELGA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 SIXX SIYY SIZZ SIXY
1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE SIGM_ELNO
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 SIXX SIYY SIZZ SIXY
N1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE FORC_NODA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
NOEUD DX DY
N1 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE EPSG_ELGA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE EPSG_ELNO
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
N1 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N2 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N3 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
N4 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00 0.00000000000000E+00
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE VARI_ELGA
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 VARI
1 0.00000000000000E+00
2 0.00000000000000E+00
3 0.00000000000000E+00
4 0.00000000000000E+00
------>
CARTE DE NOM SYMBOLIQUE COMPORTEMENT
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE ENEL_NOEU
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
NOEUD TOTALE
N1 0.00000000000000E+00
N2 0.00000000000000E+00
N3 0.00000000000000E+00
N4 0.00000000000000E+00
------>
CHAMP PAR ELEMENT CONSTANT SUR L'ELEMENT DE NOM SYMBOLIQUE ENEL_ELEM
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 TOTALE
0.00000000000000E+00
------>
CHAMP PAR ELEMENT CONSTANT SUR L'ELEMENT DE NOM SYMBOLIQUE ETOT_ELEM
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
E1 TOTALE
0.00000000000000E+00
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE ETOT_NOEU
NUMERO D'ORDRE: 0 INST: 0.00000000000000E+00
NOEUD TOTALE
N1 0.00000000000000E+00
N2 0.00000000000000E+00
N3 0.00000000000000E+00
N4 0.00000000000000E+00
======>
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE DEPL
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
NOEUD DX DY
N1 3.76141205561023E-27 9.08767762814606E-28
N2 -4.75507425715860E+00 -8.11866307392406E+00
N3 -3.63268240512585E+00 -8.77303119819776E+00
N4 0.00000000000000E+00 -5.04870979341448E-29
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE EPSI_ELGA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
1 -4.51788495009590E-01 -1.38284255787586E-01 2.52888321770218E-01 -2.94252712954232E-01
2 -4.51788495009590E-01 -5.16083868486120E-01 4.14802441498162E-01 2.97539059994448E-02
3 -3.86987171218855E-01 -5.16083868486120E-01 3.87030445587847E-01 1.08639253645180E-02
4 -3.86987171218855E-01 -1.38284255787586E-01 2.25116325859903E-01 -3.13142693589159E-01
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE EPSI_ELNO
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
N1 -4.75507425715860E-01 0.00000000000000E+00 2.03788896735369E-01 -4.05933153696203E-01
N2 -4.75507425715860E-01 -6.54368124273706E-01 4.84232378566957E-01 1.55262772320174E-01
N3 -3.63268240512585E-01 -6.54368124273706E-01 4.36129870622696E-01 1.22544366106489E-01
N4 -3.63268240512585E-01 5.55111512312578E-17 1.55686388791108E-01 -4.38651559909888E-01
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE SIEF_ELGA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 SIXX SIYY SIZZ SIXY
1 -5.53972120726879E+01 1.46497517450993E+02 0.00000000000000E+00 -4.69527399375953E+01
2 1.77246500715589E+01 -1.47433030015102E+01 1.74132385587134E-18 -4.68980318333484E+00
3 4.93609923779767E+01 9.11867756161782E+00 0.00000000000000E+00 -2.26205092586501E+01
4 -3.98761348479309E+01 2.01496650281590E+02 0.00000000000000E+00 -7.42634623737979E+01
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE SIGM_ELNO
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 SIXX SIYY SIZZ SIXY
N1 -8.56837436156448E+01 1.81213086659204E+02 -8.70661927935669E-19 -5.11689463163885E+01
N2 3.07503680891854E+01 -7.83240430825502E+01 3.24935455127179E-18 1.60859153341716E+01
N3 9.57627890025944E+01 -5.67340439140632E+01 -8.70661927935670E-19 -9.02428651896933E+00
N4 -6.90171179472181E+01 2.96214542630100E+02 2.33293160470887E-19 -1.04419197252192E+02
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE FORC_NODA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
NOEUD DX DY
N1 2.80294095592074E+02 -1.65723804086090E+02
N2 -3.74845717706820E+00 -2.18792970730141E+01
N3 4.57136383646265E+00 -2.69285462982814E+01
N4 -2.81117002251468E+02 2.14531647457386E+02
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE EPSG_ELGA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
1 -8.84617686991390E-03 -1.00593604395780E-01 4.69027633995833E-02 -2.90742122161021E-01
2 -8.84617686991390E-03 8.87918398040899E-03 -1.41459044978941E-05 4.28556092042912E-02
3 6.06866379224402E-02 8.87918398040899E-03 -2.98139236726497E-02 6.23955912195591E-02
4 6.06866379224402E-02 -1.00593604395780E-01 1.71029856314315E-02 -2.99334825338670E-01
------>
CHAMP PAR ELEMENT AUX NOEUDS DE NOM SYMBOLIQUE EPSG_ELNO
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 EPXX EPYY EPZZ EPXY
N1 -3.42969534805539E-02 -1.40663425964584E-01 7.49830197622018E-02 -4.05933153696203E-01
N2 -3.42969534805539E-02 4.89490055492121E-02 -6.27945088656781E-03 1.54039658582593E-01
N3 8.61374145330802E-02 4.89490055492121E-02 -5.78941800352682E-02 2.05719307947658E-01
N4 8.61374145330802E-02 -1.40663425964584E-01 2.33682906135015E-02 -4.38651559909888E-01
------>
CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE VARI_ELGA
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 VARI
1 0.00000000000000E+00
2 0.00000000000000E+00
3 0.00000000000000E+00
4 0.00000000000000E+00
------>
CARTE DE NOM SYMBOLIQUE COMPORTEMENT
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE ENEL_NOEU
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
NOEUD TOTALE
N1 3.49877045496888E+01
N2 -9.74514927024053E+00
N3 -1.96612360938433E+01
N4 9.52268401722956E+01
------>
CHAMP PAR ELEMENT CONSTANT SUR L'ELEMENT DE NOM SYMBOLIQUE ENEL_ELEM
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 TOTALE
2.52020398394751E+02
------>
CHAMP PAR ELEMENT CONSTANT SUR L'ELEMENT DE NOM SYMBOLIQUE ETOT_ELEM
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
E1 TOTALE
5.18770804786897E+01
------>
CHAMP AUX NOEUDS DE NOM SYMBOLIQUE ETOT_NOEU
NUMERO D'ORDRE: 1 INST: 1.00000000000000E+00
NOEUD TOTALE
N1 2.02532742038781E+01
N2 -3.75358980693637E-01
N3 -2.88512611050861E+01
N4 2.97241780733775E+01
<I> <FIN> FERMETURE DE LA BASE "GLOBALE" EFFECTUEE.
<FIN> Arrêt normal dans "FIN".
<I> <FIN> ARRET NORMAL DANS "FIN" PAR APPEL A "JEFINI".
<I> <FIN> MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION : 20.99 Mo
<I> <FIN> MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION : 27.32 Mo
<I> <FIN> MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS LORS DE L'EXECUTION : 226.64 Mo
********************************************************************************
* COMMAND : USER : SYSTEM : USER+SYS : ELAPSED *
********************************************************************************
* init (jdc) : 0.16 : 0.02 : 0.18 : 0.18 *
* . compile : 0.00 : 0.00 : 0.00 : 0.00 *
* . exec_compile : 0.06 : 0.01 : 0.07 : 0.06 *
* . report : 0.00 : 0.01 : 0.01 : 0.01 *
* . build : 0.00 : 0.00 : 0.00 : 0.00 *
* DEBUT : 0.02 : 0.01 : 0.03 : 0.03 *
* LIRE_MAILLAGE : 0.00 : 0.00 : 0.00 : 0.01 *
* AFFE_MODELE : 0.00 : 0.00 : 0.00 : 0.00 *
* DEFI_MATERIAU : 0.01 : 0.00 : 0.01 : 0.00 *
* AFFE_MATERIAU : 0.00 : 0.00 : 0.00 : 0.01 *
* AFFE_CHAR_MECA : 0.00 : 0.00 : 0.00 : 0.00 *
* AFFE_CHAR_MECA : 0.01 : 0.00 : 0.01 : 0.00 *
* DEFI_LIST_REEL : 0.00 : 0.00 : 0.00 : 0.00 *
* DEFI_FONCTION : 0.00 : 0.00 : 0.00 : 0.00 *
* STAT_NON_LINE : 0.07 : 0.01 : 0.08 : 0.09 *
* CALC_CHAMP : 0.03 : 0.00 : 0.03 : 0.03 *
* IMPR_RESU : 0.01 : 0.00 : 0.01 : 0.01 *
* FIN : 0.01 : 0.01 : 0.02 : 0.03 *
* . part Superviseur : 0.18 : 0.03 : 0.21 : 0.23 *
* . part Fortran : 0.14 : 0.02 : 0.16 : 0.16 *
********************************************************************************
* TOTAL_JOB : 0.32 : 0.05 : 0.37 : 0.39 *
********************************************************************************