|
#******************************************************************* #** #** v e m b l d e x m 1 1 #** #** the coupling of 2-D linear elastic problem (plain stress) and #** a temperatur diffusion. The mesh is read from an I-DEAS #** universal file. #** #** by L. Grosz Karlsruhe, June 1995 #** #******************************************************************* #** #** The data set of this examples has two parts (search for #** 'cut here'). The first part specifies the problem #** (please copy it to 'vembldexm11.equation') and the second part #** defines the control parameters (please copy it to #** 'vembldexm11.resource'). The FORTRAN code for the solution #** of the problem is generated by entering #** 'vembuild vembldexm11' into your shell. #** #******************************************************************* #>>>>>>> cut here to get vembldexm11.equation <<<<<<<<<<<<<<<<<<<<<<<<< #******************************************************************* #** #** The searched displacement (u1,u2) and temperture distribution #** u3 in the body are the solution of the coupled conservation #** law of the stress and the diffusion equation for the #** temperature. The coupling term is the thermal strain in the #** conservation law of the stress. The body is loaded by a #** surface load p=(p1,p2) and is fixed at two point in both #** directions (=> Dirichlet conditions for component 1 and 2). #** The temperature is prescribed at all surfaces (Dirichlet #** conditions for component 3). #** #******************************************************************* #** #** The geometry in the I-DEAS universal file fins.unv #** describes the following geometry: #** #** surface load #** ~~~~~~~~ #** |------| |------| |------| #** | | | | | | #** | | | | | | #** | \------/ \------/ | #** | /---------\ /--------\ | u3=20 #** | | hole | | hole | | #** | | u3=1000 | body | u3=800 | | #** | \---------/ \--------/ | #** u1=0 >\----------------------------------/< u1=0 #** ^ ^ #** u2=0 u2=0 #** #** The nodes with the Dirichlet conditions for component 1, which #** is the displacement in x-direction, are indexed in restraint #** set 1. The nodes with the Dirichlet conditions for #** component 2, which is the displacement in y-direction, are #** indexed in restraint set 2. The nodes with the Dirichlet #** conditions for component 3, which is the temperature, are #** indexed in restraint set 3. In I-DEAS the displacements in #** x-direction for all these nodes are set to the value we want #** to have for the solution at this location in VECFEM. #** Especially the nodes on the surface of the left hole get the #** value 1000 but on the surface of the right hole the value 800. #** #******************************************************************* #** #** material parameter: #** nu=.3 # poisson's number alpha=0.01 # thermal coefficient of expansion #** #******************************************************************* #** #** the boundary values are defined in the universal file : #** u1=prevalue u2=prevalue u3=prevalue #** #******************************************************************* #** #** the external load works on the surface elements: #** p1=1000 p2=0 #** #******************************************************************* #** #** the strains of the searched displacements: #** eps11=u1x1-alpha*(u3-20) eps22=u2x2-alpha*(u3-20) eps12=(u1x2+u2x1)/2 #** #** the term alpha*u3 considers the thermal expansion of the body. #** #******************************************************************* #** #** the resulting stresses : #** sig11=eps11+nu*eps22 sig22=eps22+nu*eps11 sig12=eps12*(1-nu)/2 #** #******************************************************************* #** #** the conversion equation for the stress: #** line{ p1*v1 + p2*v2} + area{ v1x1*sig11+(v1x2+v2x1)*sig12+v2x2*sig22 + #** #** the diffusion of the temperature #** v3x1*u3x1+v3x2*u3x2} + #** #******************************************************************* >>>>>>>> cut here to get vembldexm11.resource <<<<<<<<<<<<<<<<<<<<<<<<< #******************************************************************* #** #** The problem has a two dimensional domain and three solution #** component: #** NK=3 DIM=2 #** #******************************************************************* #** #** One processor with maximal 20 Mbytes are used. Maximal 4000 #** nodes and 1000 elements are allowed: #** PROCESS_STORAGE=20 PROCESS_MAXNN=4000 PROCESS_MAXNE=1000 #** #******************************************************************* #** #** the is read from the file I-DEAS universal fins.unv: #** MESH_PREP=i-deas MESH_FILEIN=fins.unv #** #** The output format is I-DEAS universal file: #** MESH_POSTP=i-deas #** #******************************************************************* #** #** The problem is a steady problem : #** SOLVER_STEADY=1 #** #******************************************************************* #** #** for this problem it is better to use BICO: #** SOLVER_MS=2 #** #******************************************************************* #** #** The solution component 1,2, which are the displacements, #** are written to file disp.unv with the title 'displacents' and #** the third solution component is written to file temp.unv with #** the title 'temperature'. The error output considers all #** components: #** OUTPUT_INDEX= 110 001 OUTPUT_FILE= disp.unv, temp.unv OUTPUT_TITLE= displacements, temperature OUTPUT_ERRINDEX=111 OUTPUT_ERRFILE=error.unv OUTPUT_ERRELEM=1 # The error indicator is given on the element # centre, so that a mesh adaption can be started. #** #******************************************************************* |