diff --git a/notebooks/2015-06-25-elasticity-solver-example.ipynb b/notebooks/2015-06-25-elasticity-solver-example.ipynb index b27ccad..02f5773 100644 --- a/notebooks/2015-06-25-elasticity-solver-example.ipynb +++ b/notebooks/2015-06-25-elasticity-solver-example.ipynb @@ -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" ] }, { diff --git a/src/elements.jl b/src/elements.jl index 63dcb8d..f882ba3 100644 --- a/src/elements.jl +++ b/src/elements.jl @@ -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. diff --git a/verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.comm b/verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.comm new file mode 100644 index 0000000..74ffe1c --- /dev/null +++ b/verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.comm @@ -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() diff --git a/verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.mail b/verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.mail new file mode 100644 index 0000000..9972c1e --- /dev/null +++ b/verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.mail @@ -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 diff --git a/verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.resu b/verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.resu new file mode 100644 index 0000000..4c6650f --- /dev/null +++ b/verification/2015-10-22-plane-stress/cplan_grot_gdep_point_load.resu @@ -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) + + + !------------------------------------------------------------------------------------! + ! ! + ! ! + ! 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 + + FERMETURE DE LA BASE "GLOBALE" EFFECTUEE. + + Arrêt normal dans "FIN". + ARRET NORMAL DANS "FIN" PAR APPEL A "JEFINI". + + MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION : 20.99 Mo + MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION : 27.32 Mo + 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 * + ******************************************************************************** + diff --git a/verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.comm b/verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.comm new file mode 100644 index 0000000..a464a18 --- /dev/null +++ b/verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.comm @@ -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() diff --git a/verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.mail b/verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.mail new file mode 100644 index 0000000..9972c1e --- /dev/null +++ b/verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.mail @@ -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 diff --git a/verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.resu b/verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.resu new file mode 100644 index 0000000..62d1ed0 --- /dev/null +++ b/verification/2015-10-22-plane-stress/cplan_grot_gdep_volume_force.resu @@ -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) + + + !------------------------------------------------------------------------------------! + ! ! + ! ! + ! 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 + + FERMETURE DE LA BASE "GLOBALE" EFFECTUEE. + + Arrêt normal dans "FIN". + ARRET NORMAL DANS "FIN" PAR APPEL A "JEFINI". + + MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION : 20.99 Mo + MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION : 27.32 Mo + 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 * + ******************************************************************************** +