Home

Computational Material Failure

1. MATERIAL followed by the material label that is associated with the section e NEW_SECTION Keyword which opens the field with data that defines a section of the database e END_NEW_SECTION Key word here closes the field with data that defines a section of the database Rules This field allows you to communicate with the Section database containing data associated with the element type e g material associated with an element of the mesh if the problem is a two dimensional plane strain etc The following commands may be repeated as necessary such as DELETE NEW_SECTION and END_NEW_SECTION ASSOCIATE option should be used when the material associated with the section is removed and then re created In these circumstances the association of the section breaks at the time when the material is removed it is necessary to re link the material for the purpose of the program and is considered as a new material Example BEGIN_SECTION DELETE section_1 DELETE section_2 ASSOCIATE SECTION section_1 MATERIAL material_1 NEW_SECTION END_NEW_SECTION NEW_SECTION END_NEW_SECTION END_SECTION 15 commatfail user guide 2 10 Creating the sections NEW_SECTION and END_NEW_SECTION are the key words that identify the beginning and end respectively of the data set associated with each of the sections involved in the mesh for example if a body is composed of two materials it would have at least two sectio
2. Chapter 5 Strain injection implementation in commatfail n 0 arctan Ny Establish values of y for each corner node according to their position with respect to the local axes IF z gt 0 gt 6 1 ELSE 6 0 Calculate the characteristic length for the actual integration point nnode zl LO E n DE e n cos 0 a n sin 0 a i 1 Compute vecrtor v to compare with n in the next time step Vu VN Au Vel yVu 1 1 Vu 1 2 Vu 2 1 Vu 2 2 Valle Vu 1 1 Vu 1 2 Vully Vu 2 1 Vu 2 2 j Vu Vu min Vulle Vully 5 4 Subroutine smoth var 5 4 1 Description Performs the double smooth procedure for the computation of the crack path 5 4 2 Procedure 1 Compute elemental mass matrix and assemble the global mass matrix M A NOT NY a2 e Q Qe 2 Begin loop over elements to make the first smooth a Compute elemental mass matrices MW NOT N dQ Qe b Calculate local smoothing variable at element nodes 53 commatfail user guide ale NOT QO dO nte c Assemble to calculate first smooth variable Y M A m OAMI a End loop for first smooth variable 3 Begin loop over elements to compute derivatives and make the second smooth process a Compute normal derivative of the first smooth variable make it the second variable to smooth SE TNO pO nl b Compute elemental mass matrices M
3. e TITLE is the command or keyword used to indicate the title of the problem described by the data file Must be preceded by a string or sentence e SCRLEVEL is the command or keyword used to indicate the level of messages going to be appeared on screen and in the scr file This command must be preceded by a positive integer that indicates the level of messages per screen the higher the number the greater number of messages will appear By default the value is considered 0 e RSNLEVEL is the command or keyword used to indicate the level of messages in the file rsn This command must be preceded by a positive integer that indicates the level of messages per screen the higher the number the greater number of messages will appear By default the value is considered 0 e NDIME is the command or keyword used to indicate the dimension of the problem This command must be preceded by a positive integer that indicates the dimension of the problem 1 to 1D problems there are currently no scheduled element for 1D 2 and 3 for 2D problems and for 3D problems e STAGE_TYPE is the command or keyword to indicate the type of stage e g an incremental integration scheme in time implied a re accommodation or body and should be preceded by a keyword predefined by the program according to the types of stages that has been implemented Currently the program takes into consideration the following IMPLICIT_NR implicit integration scheme using
4. 3 commatfail problemtype for GiD 3 1 Problemtype instalation ee 3 2 Loading the commatfail problemtype 020000500004 3 3 Using the commatfail problemtype e 3 3 1 Problem Data window 0 020002 eee ee eee 3 3 2 Interval Data window 0 020002 eee ee 3 3 3 Conditions window s t a iya aoei Et e aoe PR Sone doe a EE 3 3 4 Materials window 0 02 ee ee ee 4 User example for commatfail 4 1 Problem to solve cor a asna 24 De De Ae ee ee ee ale he SME 4 2 Generating the input fle 0 0 0 0 e ASS Input file a aval we 6 Meh en eos edi 4 4 Processing the problem 0 00 eee ee 4 5 Post processing the problem 2 2 a ill 29 29 29 30 30 31 32 33 commatfail user guide 4 6 Showing the results ee A 61 Curve plots gto ek a ee o Bat wee Bg A ae Hates 436 2 Gildroutputs ico a ha eh a ath oe idk oak mn ee ee A a Strain injection implementation in commatfail 5 1 Subroutine calcimpstg 2 e bdd Description a EE E Se BR a AA ho BS Hale PLOT tna be BO Moh oto te ee hehe Le Sie A BO Ae 5 2 Subroutine cfint_e2dstd_sssdlt_wsda eee a 5 2 L a Description td 5 2 8 Ml a ee a a le ee bet oe AS Hi2s2 PROCECUEC ss 2 s eh ne Ble Sel tie Bee da Soe blo A se ae 53 gt Subroutine Le calci iri aer a el nd cB Met se eae Bee gk fee kk Didi L Description Eee ee lien ale ee aS 32 Proceed Ure an S58 dg ee A AA
5. 33 1 25 0 34 1 5 0 25 41 commatfail user guide 35 36 37 38 39 40 41 42 43 44 45 END_PARAMETERS END_NEW_NODE_SET END_NODE_SET 1 75 1 1 75 0 75 1 5 0 1 75 0 5 1 75 0 25 2 1 1 75 0 2 0 75 2 0 5 2 0 25 2 0 Elements sets definition BEGIN_ELEMENT_SET New element set NEW_ELEMENT_SET NAME Elem_set_1 element set label TYPE E2DSTD element set type existing types gt esetyp FORMULATION SMALL_STRA_SMALL_DISP_LT_WSDA SUBTYPE LINEAR_QUAD BEGIN_PARAMETERS 1 18 12 11 16 Sec_1 2 21 14 12 18 Sec_1 3 2320 14 21 Sec_i 4 28 25 20 23 Sec_2 5 33 29 25 28 Sec_i 6 37 34 29 33 Sec_1 7 41 39 34 37 Sec_1 8 45 44 39 41 Sec_i 9 12 76 11 Sec_1 10 149 712 Sec_1 11 20 15 9 14 Sec_i 12 25 22 15 20 Sec_1 13 29 27 22 25 Sec_1 14 34 32 27 29 Sec_1 15 39 38 32 34 Sec_1 16 44 43 38 39 Sec_1 17 743 6 Sec_1 18 9847 Sec_1 19 15 13 8 9 Sec_1 20 22 19 13 15 Sec_i 21 27 26 19 22 Sec_1 22 32 31 26 27 Sec_1 23 38 36 31 32 Sec_1 24 43 42 36 38 Sec_1 25 42 1 3 Sec_1 26 85 2 4 Sec_i 27 13 10 5 8 Sec_1 28 19 17 10 13 Sec_1 29 26 24 17 19 Sec_1 30 31 30 24 26 Sec_1 31 36 35 30 31 Sec_1 32 42 40 35 36 Sec_1 END_PARAMETERS END_NEW_ELEMENT_SET END_ELEMENT_SET 42 element set subtype existing types gt etype element set formulation existing types Chapter 4 User example for commatfail Boundary conditions BEGIN_BOUNDARY_CONDITIONS BEGIN_DOF_CONDITIONS Constra
6. BEGIN_PARAMETERS YOUNG 10000 young modulus POISSON 0 45 poisson ratio END_PARAMETERS END_NEW_SUBMATERIAL NEW_SUBMATERIAL NAME Submat_p1_2 submaterial label TYPE J2_IMPLEX submaterial type existing types gt sbmat stype BEGIN_PARAMETERS SIGMAU 80 0 elastic limit HARD_TYPE LINEAR hardening type GFVAL 8 0 plastic modulus END_PARAMETERS END_NEW_SUBMATERIAL END_NEW_MATERIAL END_MATERIAL Sections definition 40 Chapter 4 User example for commatfail BEGIN_SECTION New section NEW_SECTION NAME Sec_1 section label TYPE PLANE_STRAIN section type existing types gt sec stype BEGIN_PARAMETERS MATERIAL Mat_1 material label associated END_PARAMETERS END_NEW_SECTION New section NEW_SECTION NAME Sec_2 section label TYPE PLANE_STRAIN section type existing types gt sec stype BEGIN_PARAMETERS MATERIAL Mat_2 material label associated END_PARAMETERS END_NEW_SECTION END_SECTION Nodes sets definition BEGIN_NODE_SET New node set NEW_NODE_SET NAME Nodes_coord node set label BEGIN_PARAMETERS 1 0 1 2 0 25 1 3 0 0 75 4 0 25 0 75 5 0 5 1 6 0 0 5 7 0 25 0 5 8 0 5 0 75 9 0 5 0 5 10 0 75 1 11 0 0 25 12 0 25 0 25 13 0 75 0 75 14 0 5 0 25 15 0 75 0 5 16 0 0 17 1 1 18 0 25 0 19 1 0 75 20 0 75 0 25 21 0 5 0 22 1 0 5 23 0 75 0 24 1 25 1 25 1 0 25 26 1 25 0 75 27 1 25 0 5 28 1 0 29 1 25 0 25 30 1 5 1 31 1 5 0 75 32 1 5 0 5
7. label to be removed from the database of material Must be preceded by the label for the material to be removed e NEW_MATERIAL Keyword which opens the data field that defines a material database see scope of materials e END_NEW_MATERIAL Key word which closes the field with the data that define a material database see scope of materials Rules This field allows you to communicate with the database of materials used in the calculation creating and removing the materials from the database The following commands may be repeated as necessary DELETE NEW_MATERIAL and END_NEW_MATERIAL Example BEGIN_MATERIAL DELETE material_1 DELETE material_2 NEW_MATERIAL END_NEW_MATERIAL NEW_MATERIAL END_NEW_MATERIAL END_MATERIAL 2 6 Creating the material properties NEW_MATERIAL and END_NEW_MATERIAL are respectively the key words that identify the beginning and the end of the data set that defines a material NEW_MATERIAL NAME lt material label gt TYPE lt type of_material gt NEW_SUBMATERIAL lt properties of the material gt END_NEW_SUBMATERIAL END_NEW_MATERIAL e NEW_MATERIAL is the keyword that indicates the beginning of the data that characterize the material e END_NEW_MATERIAL is the keyword that indicates the end of the data that characterize the material 11 commatfail user guide e NAME is the keyword to specify the name of material label that will identify the da
8. o Compute the band width for Strong Discontinuity analysis k parak xI input dat gt gt PARAK 0 01 e Va lt Qsp inj HO 4 ko Koo 1 5 exp agate Vasp ing lt 4 lt oo q dko kas Vq gt Goo or for Weak Discontinuity analysis 51 commatfail user guide pO 18 p Write output elemental results output_vector output_scalar En loop over elements 2 Perform the double smoothing process to determine the crack path smoth_var Y u M7 1 ecn M7 E MO MOV a 1 A aro 5 3 Subroutine le calc 5 3 1 Description Performs the double smooth procedure for the computation of the crack path 5 3 2 Procedure 1 Choose vector normal to discontinuity n IF n is to be chosen from to one more collinear to v Vu Vul a Get the two normals to the discontinuity obtained from previous bifurcation analysis we and select one comparing with v Vu Vu n n3 bifurcation di n v dg N2 V IF d gt da aa m n2 IF d lt da ee m n Rotate the normal vector if needed after comparing with a reference normal introduced by the user input dat gt gt NORMAL REFERENCE NX 0 70710678 NY 0 70710678 select N Nref IF select lt 0 gt f TaN m m ELSE call n previously chosen as the one more collinear with v V u V after double smoothing process smoth_var of the last time step 2 Compute angle from normal components 52
9. remesh b Initialize residual variable and vector I r 0 R 0 c Modify DOF by prescribed DOF in first iteration of the actual time step actualization of displacements calc_odof Un 1 Un 1 du d Compute forces and residual vectors for initial displacement calc_resid_2 47 commatfail user guide Pet A NT bayi dQ ia IN i ar calc_extload eE ate een Ft A BOT a Aun 1 a r Aun 1 cfint_e2dstd_sssdlt_wsda ate R Fet_ pit calc forctoresid e Begin loop over time step iterations k until the error gets lower than certain tolerance or a maximum number of iterations is complete IIr 0 i Calculate stiffness matrix calc_stiff ili K A x KO K KR e stiff e2dstd sssdlt wsda Calculate condensate residual conds_resid_e2dstd_c2 R R A x KaK iar Solve the general system of equations for du and u IF arclength solver is used Jul K u R arclen actualizedof ELSE A gt k u K u R solver actualizedof lul K Tajla r kjuju u act_condsvar_e2dstd_c2 Update displacement and increment uk Uni Uppi OU k AUn 1 Un Un Compute forces and residual vectors for actual iteration calc_resid_2 Pet A NT bya dQ Le N to ar calc_extload eE ate F A on BOT o Aun 1 a r Aun 1 cfint_e2dstd_sssdlt_wsda R F F calc forctoresid 48 Chapter 5 Strain in
10. Description sname character Section user label secnum integer Section number stype character Section type props real Section properties mat tp_mat Material data next tp_sec Pointer to next section B 5 Material variables General material variables are stored in variable mat of type tp_mat Table B 5 Material variables tp_mat Variable Type Description mname character Material user label matnum integer Material number mtype character Material type nsbmat integer Number of sub material model heads lasts tp_sbmat List of submaterial models next tp_mat Pointer to next material 62 Appendix B Variables B 6 Sub Material variables General sub material variables are stored in variable heads and lasts of type tp_sbmat Table B 6 Sub Material variables tp_sbmat Variable Type Description sname character Sub Material user label stype character Sub Material type smodel character Sub Material models props real Sub Material real properties ncurve integer Number of curves headc lastc nsbmat heads lasts next tp_crv List of associated curves integer Number of sub material level tp_sbmat List of sub material models inside sub material tp_sbmat Pointer to next sub material B 7 Elemental real variables for elements of type e2dsta Table B 7 Element real variables for e2dstd ev elmr Variable Type Description elmr 1 19 Characteristic length of the element e
11. NGAUS lt number of integration points gt MGAUS lt number of mass matrix integr points gt BEGIN_PARAMETERS lt node_1 gt lt node_2 gt lt node_n gt lt section gt lt node_i gt lt node_2 gt lt node_n gt lt section gt lt node_i gt lt node_2 gt lt node_n gt lt section gt END_PARAMETERS END_NEW_ELEMENT_SET e NEW_ELEMENT_SET is the keyword that indicates the beginning of the data that characterize the set of elements e END_NEW_ELEMENT_SET is the keyword that indicates the end of the data that characterize the set of elements e NAME is the keyword to specify the name of the set of elements tag that will identify the set in the database set of elements Must be preceded by a word that will gain the category label It is a mandatory parameter and is recommended to be located at the beginning of the field e TYPE is the keyword to specify the type of elements of the set and therefore that data that is needed and should be interpreted Must be preceded by a keyword predefined by the program depending on the set of elements It is a mandatory parameter and is recommended to be located at the beginning of the field with the only restriction to be declared before the BEGIN_PARAMETERS field and other parameters for the type of elements The elements that are currently available are E2DSTD 2D standard isoparametric element 21 commatfail user guide E2DPFEM 2D iso
12. NOT NW dQ nte c Assemble to calculate second smooth variable e m A momy in dn End loop for second smoothing 4 Compute the function vy which is going to take unit or null values for nodes belonging to Qt and 07 respectively It is calculated searching for the sides crossed by the discontinuity One can know which side is crossed looking for a change of sign in the variable u between the nodes that define the side IP sign 1 ne A sign 1 7 gt side crossed then the function y may take values as 0100 0010 1011 eee O 0001 54 Chapter 5 Strain injection implementation in commatfail 5 Calculate the gradient of the function Ve YNOT pls 6 Begin loop over elements to choose vector normal to discontinuity n only IF n is to be chosen from to one more collinear to v V u V pl a Compute the derivative of the smooth variable and the vector to select n Vuy VN y Vel VVa 1 Vu 2 Ve v PV Get the two normals to the discontinuity obtained from previous bifurcation analysis and select one comparing with v V u V al n1 n2 bifurcation dj nj v dg N2 Uv IF d gt da Sd Eo m ma IF d lt da T m ni Rotate the normal vector if needed after comparing with a reference normal introduced by the user input dat gt gt NORMAL_REFERENCE NX 0 70710678 NY 0 70710678 select N Nref IF select lt 0 gt f tet gt m m E
13. a common property which is referenced throughout instead of a group e g a set of nodes that have the same displacement When a subset of nodes is cleared during a stage with the command DELETE_SUBSET then during this stage it can not generate a subset of nodes with the same name as the deleted one It avoids ambiguities in interpreting the subset of nodes The nodes can belong to several subsets of nodes subset or none Example 7 BEGIN_NODE_SET DELETE_SET node_set_1 NEW_NODE_SET END_NEW_NODE_SET NEW_NODE_SUBSET END_NEW_NODE_SUBSET END_NODE_SET 2 12 Creating the node set NEW_NODE_SET and END_NEW_NODE_SET are the key words that identify the beginning and end respec tively of the data set associated with each of the sets of nodes that make up the mesh e g the discretization of several bodies could be represented as a set of nodes for each body NEW_NODE_SET NAME lt node set label gt BEGIN_PARAMETERS lt node_1 gt lt coord_1 gt lt coord_2 gt lt coord_3 gt lt node_1 gt lt coord_1 gt lt coord_2 gt lt coord_3 gt lt node_1 gt lt coord_1 gt lt coord_2 gt lt coord_3 gt END_PARAMETERS END_NEW_NODE_SET e NEW_NODE_SET is the keyword that indicates the beginning of the data that characterize the set of nodes e END_NEW_NODE_SET is the keyword that indicates the end of the data that characterize the set of nodes e NAME is the keyword to specify the name of the set of nodes t
14. eS lt A oe Fe 5 4 Subroutine smoth_var iasa aae a a a Taa a aa A e ra i bE aa e To a 54 1 Description sasar A a A a ds a 54 2 E aE dee Ged EE ee an Be hd ve seh Bee ee BE Interface subroutines A Interne Cale stress lt x ea A ek wh ee ee ae Ae ee E N Au Ta Description w Yay fa geese le nnd a ah ek anton eee a ee ACA a He ash des Aly Die ee pe RAIS E do M eget Mig dot kat cae a Oe ate A 2 Interface calc_stress_sssd eri iaaea E a a a a a a Oa ADV SDSSCrIptiGM anarag s a n So a oie eG ay Ee hed i ah eae E r sels ae Dd ASES o nn oot ape Mela deeds Winky Do Sevan ty e adria O ei e AJB Intertac e2dstdes lt 0 a a A a a eile ALB Description sasis acia des a e Be G TR Cases E aoe ne a Sk ol Se A ee es De AA Jnterface elemset cmo a Sct Ree dh AAA Description ose oa A E GR ei Aa ho AR aS ATADOS gt a eh ay AVE eS AOE es De A ie Ae O BO Sok Ao Interface CalecCelas ia aon a A eed Cede ee Eas Cee a Pa la A A 5 1 Description A A EI Mee bite Se te ah aN OR de il Sk to Ne gk at ats A 6 Interface stiff_e2dstd ooe a kay a a N A a a e e E a A ALG Description 6 are a Seog mala ie bee SEY Ahd vores A e ely A622 i Casese a te hs wick ed ee A ee he en S Ee ee eee os o Variables B 1 Variables of the e2dstd set of elements 2 2 e Bi2 Gauss point variables a o sier aea ade ew oR eh te Pe eee oe AS a s B 3 Element variables gt cit Soke aS be A NO Oe ee Oe ae A BA Section Variables 4 25 ek ceded atte aa a
15. field Rules This field can contain any command Any command on the information of Stage located outside the field will not be interpreted as a figure of the Stage A transaction may not contain field to another field Stage However it can contain other fields such as field Stage Basics Material field Section field the Set of Elements field etc Example NEW_STAGE BEGIN_STAGE_BASICS END_STAGE_BASICS BEGIN_MATERIAL END_MATERIAL END_NEW_STAGE 2 4 Creating the stage basics BEGIN_STAGE_BASICS and END_STAGE_BASICS are the key words that identify the beginning and the end respectively of the general data of a stage as it includes information on stage name size or time integration scheme the solver the tolerance the number of increments or the time of completion of the stage BEGIN_STAGE_BASICS TITLE lt Title of the stage gt SCRLEVEL lt enter number gt RSNLEVEL lt enter number gt STAGE_TYPE lt Type of stage gt SOLVER_TYPE lt Type of solver gt MAX_ITERATION lt Max no of iterations per step gt TOLERANCE lt tolerance of convergence gt NUMBER_STEPS lt No of increments of the stage gt NUMBER_POST lt No of steps for post processing gt PERIOD_POST lt time period for post processing gt END_TIME lt final time of the stage gt DYNAMIC_TYPE lt Type of dynamic integration gt MESH_POST lt YES NO gt END_STAGE_BASICS e BEGIN_STAGE_BASICS is the command
16. field can contain Stage Basics However when a new stage is set NEW_STAGE END_NEW_STAGE it is necessary to define a field as STAGE_BASICS which defines the characteristics of the new stage The command SCRLEVEL is optional and can be located anywhere in the Stage Basics field how ever it is strongly recommended to write it in the beginning or otherwise some messages might be lost during reading The declared value in this field will overwrite the value that was defined in the Problem Basics see Problem Basics The whole number into SCRLEVEL indicates the level of messages you want to get in the screen The lowest integer is 0 and says there is no upper limit The higher the level the greater the number of messages per screen Thus if the level is 0 during the execution only the minimum number of messages will be written The command RSNLEVEL is optional and can be located anywhere in the Stage Basics field how ever it is strongly recommended to write it in the beginning or otherwise some messages might be lost during the reading The declared value in this field will overlap the value that was defined in the Problem Basics see field Problem Basics The whole number into RSNLEVEL indicates the level of messages you want to get into the file rsn The lowest integer is 0 and says there is no upper limit The higher the level the greater the number of messages sent to the file rsn Thus if the level is 0 during the execution only the minimum n
17. file It is the field that encompasses all other fields in the data file that can be defined and interpreted by the program However it only allows defining a single data input field This field may contain general commands such as the import file ImportFile the title of the problem Begin_Problem_Basics or End_Problem_basics or the definition of stages New_Stage or End_New_Stage Any command outside this field is not interpreted by the program Example BEGIN_PROBLEM BEGIN_PROBLEM_BASICS END_PROBLEM_BASICS NEW_STAGE END_NEW_STAGE END_PROBLEM 2 2 Creating the problem basics BEGIN_PROBLEM_BASICS and END_PROBLEM_BASICS are the key words that identify the beginning and the end respectively of the general data of the problem such as the name title of the problem size or the default time integration scheme the solver and the tolerance of the stages BEGIN_PROBLEM_BASICS TITLE lt Problem Title gt SCRLEVEL lt number gt RSNLEVEL lt number gt NDIME lt Problem Dimension gt STAGE_TYPE lt Type of stage gt commatfail user guide SOLVER_TYPE lt Solver Type gt MAX_ITERATION lt Max no of iterations per time step gt TOLERANCE lt Convergence tolerance gt MESH_POST Option lt YES NO gt END_PROBLEM_BASICS e BEGIN_PROBLEM_BASICS is the command or keyword that identifies the beginning of the field e END_PROBLEM is the command or keyword that identifies the end of the field
18. freedom such as a standard displacement of solids e END_DOF_CONDITIONS is the keyword that indicates the end of the data corresponding to the definition of the Dirichlet conditions for specific degrees of freedom e BEGIN_ADD_NODE_SUBSET is the keyword that indicates the beginning of the data to manage the subset of nodes of the corresponding node database e END_ADD_NODE_SUBSET is the keyword that indicates the end of data to manage the subset of nodes of the corresponding node database 23 commatfail user guide e DELETE_SUBSET is the keyword to specify the name of the subset of nodes label to be removed from the database Must be preceded by the label for the subset of nodes to remove This parameter must be contained in the BEGIN_ADD_NODE_SUBSET and END_ADD_NODE_SUBSET e NEW_NODE_SUBSET Keyword which opens the field with the data defining a group of nodes in the database see section Creating the Subset of Nodes This field must be contained between BEGIN_ADD_NODE_SUBSET and END_ADD_NODE_SUBSET e END_NEW_NODE_SUBSET Key word here closes the field with the data defining a group of nodes see section Creating the Subset of Nodes This field must be contained between BEGIN_ADD_NODE_SUBSET and END_ADD_NODE_SUBSET Rules This field allows you to communicate with the database that defines the contour allowing you to create modify or delete the contour of each stage For that reason this field may
19. or keyword that identifies the beginning of the field e END_STAGE_BASICS is the command or keyword that identifies the end of the field e TITLE is the command or keyword used to indicate the title of the stage Must be preceded by a string or sentence e SCRLEVEL is the command or keyword used to indicate the level of messages on screen and in the scr file in the stage This command must be preceded by a positive integer that indicates the level of messages per screen the higher the value the greater number of messages By default the value is considered 0 Chapter 2 Input data file for commatfail RSNLEVEL is the command or keyword used to indicate the level of messages in the file rsn during the stage This command must be preceded by a positive integer that indicates the level of messages per screen the higher the value the greater number of messages By default the value is considered 0 STAGE_TYPE is the command or keyword to indicate the type of stage and should be preceded by a keyword predefined by the program according to the types of Stages implemented Currently the program takes into consideration the following IMPLICIT_NR implicit integration scheme using Newton Raphson scheme IMPLICIT_NR_V2 implicit integration scheme using Newton Raphson scheme similar to the preceding but which converges in one step if the problem linear default option EXPLICIT explicit integration scheme not implemented y
20. quit El 0 0009951 Command Figure 4 5 Mesh of the problem 7 Assign the point constraints to some points Data gt conditions Point Constraints condi tion a x 0 0 and y 0 0 37 commatfail user guide b Select the 5 left points 1 3 6 11 16 8 Assign known displacement to some points Data gt conditions Point Displacement con dition a Specify the properties of the displacement curve gt IPoint Displacements Curve type 2 points v Curve name horiz_disp 2 Displacement direction DISP Me NZ ordinate point 1 0 0 abscissa point 1 0 0 ordinate point 2 abscissa point 2 0 3 0 Assign Entities Draw Unassign Close Figure 4 6 Displacement curve properties b Select the 5 right points 40 42 43 44 45 9 Assign a point to create a force displacement plot Data gt conditions Point Plots condi tion a Specify the properties of the plot e 2 Ordinate variable DISPLACEMENT Ord variable direction X Ne Abscissa variable REACTION Abs variable direction x v g Assign Entities Draw Unassign Close Figure 4 7 Plot properties b Assign Select the point 43 10 Materials Data gt Materials a Create two new materials of type elasto plastic with the properties on Figure 4 8 38 Chapter 4 Use
21. such as DELETE_SET or NEW_NODE_SET and END_NEW_NODE_SET may be repeated as necessary The set of nodes should be interpreted as how to declare the program all nodes involved in the problem This statement should be made for at least one set of nodes that may be unique regardless of the type of problem However this statement can subdivide the set of nodes in a number of sets for the sake of clarity e g a set of nodes for each body that composes the problem It is mandatory that the nodes are declared through a set of nodes but a node can not belong simultaneously in two or more different sets of nodes If a calculation for adding new nodes after remeshing is produced it is possible to change the labels of the nodes added when they come into conflict with the labels of nodes generated during remeshing In this case all references to the nodes in the corresponding stage like the table of 17 commatfail user guide connectivities degrees of freedom of a node etc are considered for the new added nodes whose labels have been corrected This field also allows communication with the database Subset of Nodes that contain data asso ciated with any particular group of nodes identified only by their labels This may be repeated as necessary with the commands DELETE_SUBSET or NEW_NODE_SUBSET and END_NEW_NODE_SUBSET The subset of nodes to be interpreted as ways of defining groups of nodes This is especially useful when these nodes share
22. the category label It is a mandatory parameter and is recommended to be located at the beginning of the field e TYPE is the keyword to specify the type of submodel linear elastic elastoplastic etc and must be preceded by a keyword predefined by the program according to the type of material submodel It is a mandatory parameter and is recommended to be located at the begin ning of the field with the only restriction to be observed that it comes before the command BEGIN_PARAMETERS e BEGIN_PARAMETERS Keyword which opens the field to interpret the constants that define the submodel of material such as elastic model surface creep and so on e END_PARAMETERS Key word which declares the closing of the field of material parameters for the submodel e NEW_CURVE Keyword that paves the way to interpret a curve see Definition field curves e END_NEW_CURVE Key word which declares the closing of the field of material curves Rules This field can contain any command or data on the parameters of submodel and any material information located outside shall not be interpreted It is mandatory to define a label identifying the material which is achieved through the NAME command in this field Failure generates an error in reading inputs that will stop the program It is also mandatory to define the type of submodel through the TYPE command the omission of which also generates a reading error that stops the execution To define the properti
23. to be located at the beginning of the field e TYPE is the keyword to specify the type of the curve if a cloud of points or any other family of curves that is scheduled and should be preceded by a keyword predefined by the program according to the type of curve It is a mandatory parameter and is recommended to be located at the beginning of the field with the only restriction to be observed that it comes before the field BEGIN_PARAMETERS Currently the program takes into consideration the following types BY_POINTS curve defined by points e BEGIN_PARAMETERS Keyword that paves the way to interpret the data that define the curve e END_PARAMETERS Key word that closes the field of parameters that define the curve Rules This field can contain any command or data on the curves and any command or data outside of it will not be interpreted It is mandatory to define a label that identifies and distinguishes the rest of the curve which is achieved through the NAME command in this field Failure generates an error in input reading that stops the program Furthermore the label of the curve is restricted to the set of curves defined by the base curve Two different materials for example can have separate curves with identical labels and there is no conflict because the discrimination between the two is done via the label on the material and therefore are two different base curves It is also mandatory to define the type of curve through
24. type E2DSTD 22 Chapter 2 Input data file for commatfail Example NEW_ELEMENT_SET NAME set_elements_i TYPE E2DSTD FORMULATION Small_Stra_Small_Disp_LA SUBTYPE Linear_Triangle BEGIN_PARAMETERS 1 2 4 Section_1 Section_1 Section_1 Section_2 Section_1 Section_1 SECTION_2 SECTION_2 8 10 Section_2 11 10 Section_2 9 11 Section_2 12 11 Section_2 END_PARAMETERS END_NEW_ELEMENT_SET ODWAAAHWMH ooN Na Ss OADAANAOAAATBPWN ND 2 16 Creating the boundary conditions BEGIN_BOUNDARY_CONDITIONS and END_BOUNDARY_CONDITIONS are the key words that identify the be ginning and end respectively of the data set associated with the boundary conditions BEGIN_BOUNDARY_CONDITIONS BEGIN_DOF_CONDITIONS lt Def degrees of freedom imposed Dirichlet conditions gt END_DOF_CONDITIONS BEGIN_ADD_NODE_SUBSET DELETE_SUBSET lt label of subset to eliminate gt NEW_NODE_SUBSET lt body of nodes set with imposed DOF gt END_NEW_NODE_SUBSET END_ADD_NODE_SUBSET END_BOUNDARY_CONDITIONS BEGIN_BOUNDARY_CONDITIONS is the keyword that indicates the beginning of the data corre sponding to the definition of boundary conditions e END_BOUNDARY_CONDITIONS is the keyword that indicates the end of the data corresponding to the definition of boundary conditions e BEGIN_DOF_CONDITIONS is the keyword that indicates the beginning of the data corresponding to the definition of the Dirichlet conditions for specific degrees of
25. Bee a fee B5 Materialovaniables io ii a A ee A A E B 6 Sub Material variables es B 7 Elemental real variables for elements of type e2dstd o lv 47 47 47 47 49 49 49 52 52 52 53 53 53 57 57 57 57 57 57 57 58 58 58 59 59 59 59 59 59 59 59 59 Introduction The preliminary purpose of this document is to aid commatfail users and developers in getting some basic tools when using the code for themselves It was conceived after the research work of many people and as the final work to stop the implementation of new features in this code It can not be viewed neither as a rigorous manual nor an exhaustive description of each line in the code Readers who are a only user of the code and are not interested in the implementation of new features should probably go to read only chapter 2 and 3 Those chapters will guide the user through the rules to keep in mind when writing or creating an input file For developers we recommend a careful reading of the entire document so they have a better idea of how the code is structured Pre requisites The commatfail users do not require to have additional software in their machines Executable binaries are already compiled for Windows platform and should be used in a command line prompt The commatfail developers should have installed Microsoft Visual Studio 2010 and Intel Visual Fortran Composer 2011 Other IDEs compilers are not recommend
26. COM MAT FAIL Computational Material Failure User Guide CIMNE 2013 Contents 1 Overview of commatfail tl Modules a didatabase 2 4 cca fave eee a A dd Bae oe he 112 Stages bacara A ar dou ds A e We ee he Be a a a BS 13 The malo programi i e acras io a a ee Re eee ee 1 4 Description of files managed by commatfail 2 04 1 5 Fil commands 0 22 40e0h 4 ed Pk Re a eee ae de ee 1 6 Gode comments ss ay bies a OR eo et a Gee o ds OR ara 2 Input data file for commatfail 2 1 Creating input datafile e 2 2 Creating the problem basics e 2 3 Creatine the staee ts 4 225 sor A RA A A ee 2 4 Creating the stage basics ee 2 5 Creating the material data base a 2 6 Creating the material properties e 2 7 Creating the sub material properties e 2 8 Defining CUIVES aoa axe E A a as 2 9 Establishment of the basis of sections e a 2 10 Creating the sections ee 2 11 Establishment of the basis of nodes o a 2 12 Creating the node set ee 2 13 Creating the subset of nodes 2 ee 2 14 Establishment of the Basis of Elements oaa aaa a 2 15 Creating the set of elements 2 16 Creating the boundary conditions 0 0 0 0002 eee eee 2 17 Allocation of conditions on the degrees of freedom o 2 18 Creating the localization data module e
27. EW_PROBLEM e NEW_STAGE e RESTART_STAGE e RESTART_CALC e RESTART_NEWMSH 1 3 The main program The main commatfail procedure is done by mainpg It can be divided into three basic parts l Initialization and input data The initialization of the execution parameters is carried out in this phase at the very beginning of the execution of the program The reading precess of the input data is done at the beginning of each stage from the input data text file or the restart file The incremental finite element procedure This is the main body of the program where all the numerical calculations associated to the finite element method are done This procedure is done one time for each stage of the problem after reading the input data Essentially it consist of a time steps integration process with an interface routine calcstg that chooses the calculation method depending on the integration scheme implicit explicit or implicit explicit This comprises all output operations necessary to print converged finite element solutions into the results file and or dump an image of the database into a restart output file Output routines are called from inside the main loop over load increments Write output results This comprises all output operations necessary to print converged finite element solutions into the results or restart files Those routines write in files with the correct format to be read with the post processor commatf
28. Gauss point coordinates posgpm real Gauss point coord for mass matrix weigp real Gauss point weigths weigpm real Gauss point weigths for mass matrix nshape real Shape functions nshapem real Shape functions for mass matrix nderiv real Shape fun derivatives nderivm real Shape fun derivatives for mass matrix dmass real Mass integration point by element gv tp_e2dstd_gp Gauss point internal variables array ev tp_e2dstd_el Element variables next tp_e2dstd Pointer to next set of elements B 2 Gauss point variables Gauss point internal variables are stored in variable gv of type tp_e2dstd_gp 61 commatfail user guide Table B 2 Gauss point internal variables tp_e2dstd_gp Variable Type Description gausv real Internal las converged variables tmpgv tp_rv Temporal step increment variables eout tp_out Post process data by integration point B 3 Element variables Element variables are stored in variable ev of type tp_e2dstd_el Table B 3 Element variables tp_e2dstd_el Variable Type Description sec tp_sec Section label by element elmi integer Integer variables by element elmr real Real variables by element elmc real Condensate variable by element gout tp_gout Post process data by element B 4 Section variables Section variables are stored in variable sec of type tp_sec Table B 4 Section variables tp_sec Variable Type
29. LASTICLINEAR PLANESTRESS CALL calc_isolin_celas_pstress A 6 Interface stiff_e2dstd A 6 1 Description Call specific subroutines to calculate the stiffness matrix for element type E2DSTD according to the formulation of the problem A 6 2 Cases 59 commatfail user guide Table A 6 Options in interface stiff_e2dstd Case Action SMALLSTRASMALLDISPLA CALL stiff_e2dstd_sssdla SMALLSTRASMALLDISPLT CALL stiff_e2dstd_sssdlt gt SMALLSTRASMALLDISPLTBSTD CALL stiff_e2dstd_sssdlt_bstd gt SMALLSTRASMALLDISPLTBBAR CALL stiff_e2dstd_sssdlt_bstd gt SMALLSTRASMALLDISPLTWSDA 60 CALL stiff_e2dstd_sssdlt_wsda Appendix B Variables This chapter is for the users not used to the program variable saving system It will help to know where is located each variable of the set of elements e2dstd B 1 Variables of the e2dstd set of elements Variables for the hole set of elements of type e2dstd This set of variables is stored in one variable of type tp_e2dstd Table B 1 Element set variables tp_e2dstd Variable Type Description ename character Set name etype character Element type eform character Formulation type qtype character Quadrature type nelem integer Number of elements nnode integer Number of nodes element ngaus integer Number of Gauss points mgaus integer Number of Gauss points for mass matrix lnods real Conectivities posgp real
30. NT NODE REACTION NODE 43 F_X 43 F_X END_PLOT END_PLOTS END_NEW_STAGE END_PROBLEM 43 commatfail user guide End of file 4 4 Processing the problem After having the input file the user should have the following files in the same folder in order to make the computations example1 dat commatfail exe commatfailpost The processing phase is launched with the following command in a command prompt commatfail example1 The beginning and ending of this process should look like shown in Figure 4 9 EX Procesador de comandos de Windows gt EX Procesador de comandos de Windows estrada example gt comfm examplel ce RERERER a 7183966 Cond Variable 6 14786845 100 000000A 3 Error 8 3133034390E Cond Variable Error 6 13396868 fully donet routine ends STAGE horizontal load NORMAL END EXECUTION OF example 1 Ended at 11 38 25 on 22 65 2013 iC Users mestrada example gt a Begin of process b End of process Figure 4 9 Processing the problem 4 5 Post processing the problem The post processing phase is launched with the following command in the same command prompt commatfailpost examplel The beginning and ending of this process should look like shown in Figure 4 10 EMM Procesador de comandos de Windows e Loa EM Procesador de comandos de Windows ke ecco fe fat a 96 90008008 ic Users mestrada example gt comf
31. Newton Raphson scheme IMPLICIT_NR_V2 implicit integration scheme using Newton Raphson scheme similar to the preceding but which converges in one step if the problem is linear default option EXPLICIT explicit integration scheme not scheduled yet BODY_MOVE rearranging or moving solid rigid bodies not scheduled yet The default option IMPLICIT_NR is assigned when this command is not explicitly stated in some of the stages However the obligation of the command depends on the type of SOLVER_TYPE e SOLVER_TYPE is the command or keyword to indicate the type of solver to be used and should be preceded by a keyword predefined by the program according to the type of solver that has already been implemented Currently the program takes into consideration the following solvers GAUSSJORDAN direct solver by Gauss Jordan elimination with full matrix Use of this solver is limited for smaller examples LU direct solver with full storage array Use of this solver is not meant to medium large problems so its use is limited only to debug small examples SUPERLU direct solver with storage of the sparse matrix type This constraint solver has the ability to contain the decomposition of the matrix in memory Chapter 2 Input data file for commatfail e MAX_ITERATION is the command or keyword to specify the value for the maximum number of iterations required for convergence per time step and must be preceded by a
32. OnOFWBN 2 13 Creating the subset of nodes NEW_NODE_SUBSET and END_NEW_NODE_SUBSET are the key words that identify the beginning and end respectively of the data set associated with a group of nodes having some common features and will be used in the program later jointly e g nodes with a same contour NEW_NODE_SUBSET NAME lt node subset label gt BEGIN_PARAMETERS lt node_1 gt lt node_2 gt lt node_3 gt END_PARAMETERS END_NEW_NODE_SUBSET e NEW_NODE_SUBSET is the keyword that indicates the beginning of the data that characterize the subset of nodes e END_NEW_NODE_SUBSET is the keyword that indicates the end of the data that characterize the subset of nodes e NAME is the keyword to specify the name of the subset of nodes tag that will identify you in the cluster nodes database Must be preceded by a word that will gain the category label Although it is an optional parameter is recommended and is situated at the beginning of the field 19 commatfail user guide e BEGIN_PARAMETERS Keyword which opens the field to interpret the nodes labels that make up the group e END_PARAMETERS Keyword closing the field used to interpret the nodes labels that make up the group Rules This field is optional in the data file and can be declared at the base of nodes defined by BEGIN_NODE_SUBSET and END_NODE_SUBSET but alternatively there are other fields in the data file where they can be dec
33. TA e BEGIN_LOCALIZATIONDATA is the keyword to specify the beginning of the field of strain injection data e END_LOCALIZATIONDATA is the keyword to specify the end of the field of strain injection data e DIS_INJ_TYPE is the command to indicate the strain injection type Currently the code has the following options WD_inj 26 Chapter 2 Input data file for commatfail SD_inj Sym_SD_inj No_dis_inj e TAU_LOC is the command to indicate the value of tio for the first injection e TAU_MXD is the command to indicate the value of Tnzq for the elements that has not been injected e Q_INJECTION is the command to specify the value of qin which indicate the moment for injection e Q_INF is the command to specify the lowest value of q e NORMAL_REFERENCE is the command to specify the Cartesian components of the vector of reference that will help to choose the direction of discontinuity n NX refers to component z and NY to the component y of the vector e PARAK is the command to specify the band width in Strong Discontinuities Analysis as a fraction of 1 9 Rules All parameters TAU_MXD TAU_LOC Q_INJECTION Q_INF and PARAK must have values between 0 0 and 1 0 Example cr BEGIN_LOCALIZATIONDATA DIS_INJ_TYPE SD_inj TAU_MXD 0 1 TAU_LOC 0 0 Q_INJECTION 0 70 Q_INF 0 1 NORMAL_REFERENCE NX 0 70710678 NY 0 70710678 PARAK 1 00 END_LOCALIZATIONDATA 27 C
34. XISIMMETRIC Elements 2D axisymmetric problems e BEGIN_PARAMETERS Keyword that paves the way to interpret the data that define the section such as the material to which it is associated and so on e END_PARAMETERS Keyword closing the field used to describe the parameters that define the section Rules This field can contain any command or data on the sections and any command or data located outside shall not be interpreted in relation to the database of the sections It is mandatory to define a label that identifies the section which is achieved through the NAME command in this field Failure generates an error in input reading that will stop the program It is also mandatory to define the type of section through the TYPE command the omission of which also generates a read error that stops the execution The section data is introduced into a single field parameter BEGIN_PARAMETERS and END_PARAMETERS Example NEW_SECTION NAME section_1 TYPE PLANE_STRESS 16 Chapter 2 Input data file for commatfail BEGIN_PARAMETERS MATERIAL AL4030 THICK 1 0E 3 END_PARAMETERS END_NEW_SECTION 2 11 Establishment of the basis of nodes BEGIN_NODE_SET and END_NODE_SET are the key words that identify the beginning and end respectively of the data set node labels and their coordinates associated with the set of nodes BEGIN_NODE_SET DELETE_SET lt node set label gt DELETE_SUBSET lt node subset label gt NEW
35. _NODE_SET lt body that defines the set of nodes gt END_NEW_NODE_SET NEW_NODE_SUBSET lt body that defines the subset of nodes gt END_NEW_NODE_SUBSET END_NODE_SET e BEGIN_NODE_SET is the keyword that indicates the beginning of the data for the node e END_NODE_SET is the keyword that indicates the end of the data for the node e DELETE_SET is the keyword to specify the name of the set of nodes label to be removed from the database Must be preceded by the label for the set of nodes to remove e DELETE_SUBSET is the keyword to specify the name of the subset of nodes label to be removed from the database Must be preceded by the label for the subset of nodes to remove e NEW_NODE_SET Keyword which opens the field of the data that define a set of nodes in the database see section Creating the Node Set e END_NEW_NODE_SET Key word here closes the field with the data that define a set of nodes in the database see section Creating the Node Set e NEW_NODE_SUBSET Keyword which opens the field with the data defining a grouping of nodes in the database see section Creating the Subset of Nodes e END_NEW_NODE_SUBSET Key word here closes the field with the data defining a grouping of nodes in the database see section Creating the Subset of Nodes Rules This field allows you to communicate with the database of Set of nodes that contain data associated with nodes tags and coordinates The commands
36. _TYPE SUPER_LU solver type existing types gt solvtype MAX_ITERATION 50 maximum iterations for convergence TOLERANCE 1e 6 tolerance level for convergence END_PROBLEM_BASICS Stage horizontal load NEW_STAGE 39 commatfail user guide Stage basic information BEGIN_STAGE_BASICS TITLE horizontal load NUMBER_STEPS 600 time steps increments NUMBER_POST 100 postprocess time steps END_TIME 100 stage final time END_STAGE_BASICS Materials definition BEGIN_MATERIAL New material NEW_MATERIAL NAME Mat_1 material label TYPE ELASTO_PLASTIC2 material type existing types gt mat mtype NEW_SUBMATERIAL NAME Submat_el_1 submaterial label TYPE ELASTIC_LINEAR submaterial type existing types gt sbmat stype BEGIN_PARAMETERS YOUNG 10000 young modulus POISSON 0 45 poisson ratio END_PARAMETERS END_NEW_SUBMATERIAL NEW_SUBMATERIAL NAME Submat_p1_1i submaterial label TYPE J2_IMPLEX submaterial type existing types gt sbmat stype BEGIN_PARAMETERS SIGMAU 100 0 elastic limit HARD_TYPE LINEAR hardening type GFVAL 8 0 plastic modulus END_PARAMETERS END_NEW_SUBMATERIAL END_NEW_MATERIAL New material NEW_MATERIAL NAME Mat_2 material label TYPE ELASTO_PLASTIC2 material type existing types gt mat mtype NEW_SUBMATERIAL NAME Submat_el_2 submaterial label TYPE ELASTIC_LINEAR submaterial type existing types gt sbmat stype
37. _sssd_el ELASTOPLASTIC CALL calc_stress_sssd_elpl Continue on next page 57 commatfail user guide Table A 2 continued from previous page Case Action gt ELASTOPLASTIC2 CALL calc_stress_sssd_elp12 ELASTICDAMAGE CALL calc_stress_sssd_eldmg gt RANKINEMODEL CALL calc_stress_sssd_rankine gt RANKINEMODEL2 CALL calc_stress_sssd_rankine2 gt RANKINEMODEL3 CALL calc_stress_sssd_rankine3_v3 gt MATTEST CALL calc_stress_sssd_mattest A 3 Interface e2dstd A 3 1 Description Choose the following interface of calls for specific tasks depending on the type of finite element A 3 2 Cases Table A 3 Options in interface e2dstd Case Action RESTART_RD CALL rest_rd_e2dstd RESTART WR CALL rest_wr_e2dstd gt CALCDOF CALL cdof_e2dstd gt CALC_FINT CALL cfint_e2dstd gt CALC_FEXT CALL cfext_e2dstd gt CALC_FDYN CALL cfdyn_e2dstd gt OUTPUT CALL wroutpost_e2dstd gt CALC_STIFF CALL stiff_e2dstd gt CALC_DYN_STIFF CALL stiff_dyn_e2dstd gt GET_REMESH_DATA CALL rmshedata_e2dstd gt TRANSF_REMESH DATA CALL transfedata_e2dstd CALL stdtrianmesh veer gn CALL getlnods_e2dstd gt GET_BOUND CALL getlnodb_e2dstd gt MESH_OUTPUT CALL wroutmsh_e2dstd SEARCH found TRIM lwcase ename TRIM lwcase eset ename ARRGESET CALL ini_e2dstd DELETE CALL de
38. ag that will identify you in the set of nodes database Must be preceded by a word that will gain the category label Although it is an optional parameter it is recommended and situated at the beginning of the field e BEGIN_PARAMETERS Keyword employee who opened the field to interpret the data that define the set of nodes label and coordinates e END_PARAMETERS Keyword closing the field that is used to interpret the data that define the set of nodes 18 Chapter 2 Input data file for commatfail Rules This field should appear in the data file at least once It can contain any command or data on the set of nodes and any command or data located outside shall not be interpreted as a figure corresponding to the database of sets of nodes The label that identifies the set of nodes is introduced through the NAME command in this field Its omission is not recommended and it assigns a default value corresponding to the word NODE_SET followed by a number A node can not belong to different sets of nodes Data from the set of nodes tags and coordinates for the discretization are introduced into single field parameters BEGIN_PARAMETERS and END_PARAMETERS Example NEW_NODE_SET NAME nodeset 1 BEGIN_PARAMETERS 1 0 000 0 000 5 000 0 000 10 00 0 000 0 000 5 000 5 000 5 000 10 00 5 000 0 000 10 00 5 000 10 00 9 10 00 10 00 10 0 000 15 00 11 5 000 15 00 12 10 00 15 00 END_PARAMETERS END_NEW_NODE_SET ON
39. ailpost Output files from commatfailpost are then in the correct format to be imported in GiD 1 4 Description of files managed by commatfail commatfail manages the opening and closing of files through two routines The routine that handles the opening of files is called openfi while the routine that manages the closing of files is called closefi In turn both depend on a common module called filevar_db in which units and extensions are defined The manipulation of files at the programming level required to make co ordinated development of this module and these two routines e commatfail dat This file actually corresponds to the data file that will be used by the program and it is a single file with normal ASCII characters One of the first tasks of commatfail when launching a new calculation is to build this file from the data files provided by the user In case the data file provided by the user is unique commatfail dat is a copy of it Chapter 1 Overview of commatfail e lt name gt scr A type of report file that replicates all the information that comes to the screen e lt name gt rpt A type of report file that provides basic information such as command line launched if the calculation is new or a restart start date and time of release number of elements for the problem and real time calculation of CPU etc and it is useful for a novice user e lt name gt rsn A type of file focused on advanced or ver
40. b To create the first corner of the rectangle write into the command line 0 0 enter c To create the second corner of the rectangle write into the command line 4 2 enter The screen should look like the one shown in Figure 4 2 35 commatfail user guide Aplicaciones s 4 may20 4 44PM 2M GiD x64 Project UNNAMED comfm Files View Geometry Utilities Data Mesh Calculate Help 831SE4341 8519 bea ISL ae Pe ea p JIZ INGE 195 S AN 3 Hl 8 BAAN gt RES SERRS RBNA PAU AAT y Es Created 1 new surface 4 new lines and 4 new points Can continue Enter first corner point ESC to leave Figure 4 2 Geometry of the problem 3 Load the commatfail problemtype Data gt Problem type gt commatfail 4 Establish the problem data Data gt Problem data Problem Data Problem Data 2 9 xa General problem basics Elements set General problem basics Elements set Title problem example 1 Elements type E2DSTD v Screen messages 3 v Elements fomulation SMALL STRA SMALL DISP LT WSDA Ne Dimensions 2 v Elements subtype LINEAR QUAD v Stage type IMPLICIT NR2 w Solver type SUPER LU v Max iteration 50 Tolerance 1e 6 Accept Close Accept Close a General problem basics b Elements set Figure 4 3 Problem data 5 Establish the interval data Data gt Interval data 36 Chapter 4 User example for commatfa
41. be fixed on 0 001 and 0 001 respectively Double Smoth 0 001 000077778 0 00055556 0 00033333 0 00011111 0 00011111 0 00033333 0 00055556 0 00077778 0 001 y ine fen Figure 4 13 Crack path results Double smoth 45 commatfail user guide The bifurcated and injected elements are shown with View results gt Contour fill gt ebif EL ebif 1 0 88889 0 77778 0 66667 0 55556 0 44444 0 33333 y 0 22222 y 0 11111 x 0 x co co a Post process step 8 b Post process step 60 EL ebif 2 1 7778 1 5556 1 3333 1 1111 0 88889 0 66667 0 44444 0 22222 0 Figure 4 14 Injection procedure results EL ebif 46 Chapter 5 Strain injection implementation in commatfail This chapter explain in a more comprehensible language than the original source Fortran 90 files the implementation of the methods for the Strain Injection Techniques in Numerical Modeling of Propagating Material Failure The complete formal explanation of this implementation can be found at Dias et al 2012 Principal subroutines related with the implementation of the strain injection techniques in commatfail are explained here 5 1 Subroutine calcimpstg 2 5 1 1 Description Calculate stage integration for implicit scheme 5 1 2 Procedure 1 Initialize stage variables calcimpstg_init F Se OR S038 2 Begin loop over time steps n a Remesh if required
42. contain other fields such as DOF or others that are planned in the future This field may be absent but can never be repeated As a general rule the condition defined on a contour of the degrees of freedom of a particular problem will continue in subsequent phases unless the condition is explicitly removed or modified to that degree of freedom Example s BEGIN_BOUNDARY_CONDITIONS BEGIN_DOF_CONDITIONS END_DOF_CONDITIONS END_BOUNDARY_CONDITIONS y 2 17 Allocation of conditions on the degrees of freedom BEGIN_DOF_CONDITIONS and END_DOF_CONDITIONS are the key words that identify the beginning and end of all the data associated with the contour on the degrees of freedom BEGIN_DOF_CONDITIONS DELETE CURVE lt label of curve to eliminate gt NEW_CURVE lt bodyefining new curve gt END_NEW_CURVE DELETE NODE lt node gt DOFi DOF2 DOFn DELETE NODE_SUBSET lt label gt DOFi DOF2 DOFn DELETE RESET_ALL_DOF CREATE NODE lt node gt DOF1 REAL lt value gt WORD lt curve label gt FACTOR_DOF1 lt value gt DOF2 REAL lt value gt WORD lt curve label gt FACTOR_DOF2 lt value gt CREATE NODE_SUBSET lt label of node subset gt DOF1 REAL lt value gt WORD lt curve label gt FACTOR_DOF1 lt value gt DOF2 REAL lt value gt WORD lt curve label gt FACTOR_DOF2 lt value gt END_DOF_CONDITIONS e BEGIN_DOF_CONDITIONS is the keyword that i
43. ect the reduced integration e Q infinite Specify the fraction of q to inject the strong discontinuity e Normal reference Specify the two Cartesian components for the reference normal vector e Parameter K Write the fraction of the element size for strong discontinuity band width 3 3 3 Conditions window Access it with the top menu Data gt Conditions el Es e Point constrains o OE o A 2 x X flag Curve type 2 points x Ordinate variable DISPLACEMENT v X displ 0 0 Curve name name_of curve Ord variable direction X F x Y flag Displacement direction DISPX Y Abscissa variable DISPLACEMENT vV Y displ 0 0 ordinate point 1 0 0 Abs variable direction x v Z flag abscissa point 1 0 0 Z displ 0 0 ordinate point 2 1 0 abscissa point 2 1 0 Entities Draw Unassign Assign Entities Draw Unassign Assign Entities Draw Unassign dose Cose close a Constraints b Displacements c Plots Figure 3 4 Conditions window for commatfail problemtype Opening the condition window the user will be able to assign constraints plots or imposed displacements to points of the mesh Point constraints The Point Constraints condition should be used to assign the boundary conditions to the geometry 32 Chapter 3 commatfail problemtype for GiD Point plots The Point Plots condition can be used to specify any curve the user wants to plot after the analysis Point d
44. ed from the database Must be preceded by the label for the set of elements to delete e NEW_ELEMENT_SET Keyword which opens the field to the data that define a set of elements of the database see section Creating the Set of Elements e END_NEW_ELEMENT_SET Key word here closes the field with the data that define a set of elements of the database see section Creating the Set of Elements 20 Chapter 2 Input data file for commatfail Rules This field allows you to communicate with the database set containing elements data asso ciated with elements or equivalent e g mesh type the order of integration design etc This may be repeated as necessary or DELETE commands NEW_ELEMENT_SET and END_NEW_ELEMENT_SET Example BEGIN ELEMENT_SET DELETE element_set_1 DELETE element_set_2 NEW_ELEMENT_SET END_NEW_ELEMENT_SET NEW_ELEMENT_SET END_NEW_ELEMENT_SET END_ELEMENT_SET 2 15 Creating the set of elements NEW_ELEMENT_SET and END_NEW_ELEMENT_SET are the key words that identify the beginning and end respectively of the data set associated with each of the sets of elements of the mesh e g the discretization of several bodies are represented with a set of elements for each body NEW_ELEMENT_SET NAME lt element set label gt TYPE lt type of element set gt FORMULATION lt problem formulation gt SUBTYPE lt subtype of element set gt QUADRATURE lt quadrature type gt
45. ed for a correct functionality of the code Installation There is no need to install commatfail If the pre requisites are met the solution file commatfail sin should open in Microsoft Visual Studio with no problems Running commatfail and commatfailpost Binaries and input file should be in the same working directory Computations should be done in a command line prompt with commatfail lt input_file_name gt If GiD will be used to see the results post process should be done after computations in a command line prompt with commatfailpost lt input_file_name gt After post process results can be seen opening the file lt input_file_name gt _s1 flavia in GiD pre post processor Clarity of the code Finite element codes are no exception when facing the dilemma of clarity versus efficiency Even though those are not opposite concepts an effort to achieve one tends to lack the other The commatfail user guide programming philosophy adopted in the development of commatfail has been the optimization of procedures This does not mean that the creators of the code have done their best to make it as comprehensible as possible Portability Most commatfail routines are written in a standard Fortran 90 programming language but some small details are only valid for Intel Fortran Compiler Currently the code is already configured in a Microsoft Visual Basic 2010 solution and its usage with other compilers and IDE s could
46. eometry IGES DXF etc or design a new one with all the CAD tools usually available in any CAD system 3 3 Using the commatfail problemtype Once the user have loaded the problemtype he will have the options that commatfail will need to perform the calculations Next the problemtype interface is described 3 3 1 Problem Data window Access it with the top menu Data gt Problem Data Problem Data Problem Data General problem basics Elements set Title problem example Screen messages 3 bi Dimensions 2 T wg General problem basics Elements set Elements type E2DSTD vi Elements fomulation SMALL STRA SMALL DISP LT WSDA Elements subtype LINEAR QUAD y RE v Stage type IMPLICITNR2 v Solver type SUPER LU v Max iteration 50 Tolerance 1e 5 Close a General problem basics b Elements set Figure 3 2 Problem data window for commatfail problemtype This window deals with the definition of the analysis to be performed General problem basics and Elements set cards are available 30 Chapter 3 commatfail problemtype for GiD General problem basics e Title problem Field to specify some name or title to the problem It can contain spaces and it is recommended to use only alphanumeric characters e Screen messages This option allows the user to choose the quantity of messages to be shown on the screen while the calculus is running e Dimensions Use this option to spec
47. es of the submodel it is defined by a single field parameter BEGIN_PARAMETERS and END_PARAMETERS but you can define as many fields as necessary curves to define the submodel material It is allowed that a submodel of material can be composed of different submodels of material so that the base material is dynamic and flexible Example NEW_SUBMATERIAL NAME model_elastic_1 TYPE ELASTIC_LINEAR BEGIN_PARAMETERS YOUNG 2 1E 6 POISSON 0 29 END_PARAMETERS NEW_CURVE END_NEW_CURVE END_NEW_SUBMATERIAL 2 8 Defining curves NEW_CURVE and END_NEW_CURVE are the key words that identify the beginning and the end respectively of the data set that define a curve which may be defined by a cloud of points or any other data set and define the interpretation of which is scheduled 13 commatfail user guide NEW_CURVE NAME lt label of the curve gt TYPE lt type of curve gt BEGIN_PARAMETERS lt body hat defines the curve gt END_PARAMETERS END_NEW_CURVE e NEW_CURVE is the keyword that indicates the beginning of the data that characterizes the curve e END_NEW_CURVE is the keyword that indicates the end of the data that characterizes the curve e NAME is the keyword to specify the name of the curve label that will identify the rest of the curves defined on the basis of corresponding curves Must be preceded by a word that will gain the category label It is a mandatory parameter and is recommended
48. et BODY_MOVE rearranging or moving solid rigid bodies not implemented yet SOLVER_TYPE is the command or keyword to indicate the type of Solver to use during the stage and must be preceded by a keyword predefined by the program according to the type of solver implemented see Problem basics MAX_ITERATION is the command or keyword to specify the value for the maximum number of iterations for convergence for the stage and should be preceded by a whole number TOLERANCE is the command or keyword to specify the value of convergence tolerance for the stage and should be preceded by a real number NUMBER_STEPS is the command or keyword used to indicate the number of increments or time steps for the stage This command must be preceded by a positive integer indicating the number of increments NUMBER_POST is the command or keyword used to indicate the number of post processing carried out on the stage following a long constant frequency This command must be preceded by a positive integer indicating the number of increments PERIOD_POST is the command or keyword used to indicate the frequency of post processing carried out on the stage by defining the time period This command must be preceded by a positive real number indicating the number of increments END_TIME is the command or keyword to specify the value of the final moment of the stage and should be preceded by a real number DYNAMIC_TYPE is the command or keyword
49. gradient vector to select n 63 Appendix B Variables Bibliography Ivo F Dias Javier Oliver and Alfredo E Huespe Strain Injection Techniques in Numerical Modeling of Propagating Material Failure Number 134 in Monograph CIMNE International Center for Numerical Methods in Engineering CIMNE 2012 ISBN 9788494024375 65
50. hapter 3 commatfail problemtype for GiD The commatfail problemtype is an user interface specifically designed for the pre processor GiD Using this interface is possible to prepare the data file for the mechanical solver commatfail The result files obtained from the analysis can be loaded in the post processing window of GID Using this interface it is possible to run a mechanical analysis according to the following steps e Geometry definition it is possible to create and or import the CAD geometry that defines the mechanical system to be analyzed e Materials it is possible to select from a material data base the material properties to be as signed to each volume previously defined It is also possible to modify the material properties according to the user specification e Conditions the user should assign corresponding loading and boundary conditions e Finite element mesh a finite element mesh can be generated for both 2D and 3D analysis 3 1 Problemtype instalation Before trying to do anything the user need to have the pre and post processor GiD installed in the computer It can be freely downloaded from http www gidhome com All the configuration files of the problemtype are in the folder commatfail gid The only thing the user have to do to install the commatfail problemtype is place that folder into the corresponding problemtypes folder of GiD Windows c Program Files lt GiD_directory gt lt GiD_version gt problem
51. ies the grouping of nodes subset followed by the word key that identifies the degree of freedom to be eliminated As was the case above it can remove some or all the contour of the subset e CREATE is the keyword to specify that you will create a contour This option must be accompanied by one of the following keywords NODE Keyword to specify the node on which it is required to define the boundary conditions This keyword must be followed by the number of the node followed by the word key that identifies the degree of freedom to be prescribed after the KEY is specified the size of the prescription that can be defined through a numerical value or through a curve In the case of a numeric value type the keyword followed by the REAL value whereas if it is a curve type the keyword followed by the label WORD of the curve Moreover it can be possible not mandatory to specify a factor that modifies the value To do this use the keyword composed of the root FACTOR and the keyword that specifies the degree of freedom followed by the numeric value of the defined factor Also on this node and on the same line of data file one or more of the contour can be created each of which corresponds to the degree of freedom to be restricted NODE_SUBSET Keyword to specify the subset of nodes on which one needs to define the boundary conditions This keyword should be followed by the label that identifies the grouping of nodes subset followed b
52. ify the geometrical dimension of the problem 1D 2D 3D e Stage type Use this option to specify the type of calculation to be used in the analysis e Solver type Use this option to specify the type of solver that the user want to use in the analysis when solving the finite element matrix equations e Max iteration Specify in this field the maximum number of iterations allowed in each time step e Tolerance Specify in this field the accuracy that define a tolerance precision for the non linear solver Elements set e Elements type Choose the type of elements of the problem The type of element should be consistent with the dimension of the test e Elements formulation Choose the finite element formulation to be executed for the problem e Elements subtype Chose the subtype of elements if any 3 3 2 Interval Data window Access it with the top menu Data gt Interval Data Interval Data Interval Data kex ja p Jex we ja Stage basics Strain localization data Stage basics Strain localization data Title stage horizontal load A Injection type SD INJ hgh Process steps 500 Tau mixed 0 1 E Postprocess steps 100 Tau localized 0 0 End time 100 7 O iniection 0 98 Close Close a Stage basics b Strain localization data Figure 3 3 Interval data window for commatfail problemtype This window deals with the definition of the general data of each stage of load to be performed Outside the
53. il Interval Data a O K x a Interval Data Stage basics Strain localization data 1 HS IX R la Injection type SD INJ X Stage basics Strain localization data OL Tau localized 0 0 Title stage horizontal load Q injection 0 95 Process steps 1500 Q infinite 0 1 Postprocess steps 100 Normal reference X 0 70710678 End timeo Normal reference Y 0 70710678 Accept Close ij Accept Close a Stage basics b Strain localization data Figure 4 4 Interval data 6 Create the mesh Use the GiD tools from the top menu Mesh to create a structured mesh that look like the one on Figure 4 5 Establish a structured mesh Mesh gt Structured gt Surfaces gt Assign number of cells Select the created rectangle surface Enter the number 4 Assign d Select the vertical left line e Enter the number 8 Assign f Select the horizontal top line g Generate the mesh Mesh gt Generate mesh GiD x64 Project examplel comfm Files View Geometry Utilities Data Mesh Calculate Help 5 O8931S684616519o pr Sl s719 i Pe 2 ga z re K 1 2 5 10 17 24 30 35 40 q NA as B 4 A 13 19 PA BA BA 42 da lt i a 7 a 15 2 p7 R2 RA 43 E ae pa 12 id pa _ Pa pa Ba pa Ad 22 Oo IN 1A if B n 3 R Bz BZ 41 45 98 X E y aL la pS Pick LEFTMOUSE to desplace view ESC to quit if present mouse wheel zooms x 2 6168 Pick LEFTMOUSE to desplace view ESC to
54. integration scheme representative of the stress at the discontinuity B VN calc_bdef sssdlt_bstd e e e e e e 1 e e Ae BO Au Elo ve uy no uy a Y Vu Asf calc_stress Calculate stress for gauss point 6 elastic process is forced once a discontinuity is in jected BO VN_ calc_bdef sssdlt_bstd Ac BO Au 0 ver uy cate Y viu Ae calc_stress Calculate stress for regular gauss points 1 2 3 4 standard 2 x 2 gauss points for quadri lateral elements B VN calc_bdef sssdlt_bstd Ac BO Au 0 v uy a Y Vu Ae calc_stress Compute the increment of stress Ae 1 Eo 2 o E a y 1 7 ao 7 00 m 9 fon a 00 aon 2 Compute the new stress i Elemental internal forces vector 50 Chapter 5 Strain injection implementation in commatfail ngp int e eT le e pino S5 Be a avg gp 1 3 Compute increment of different stresses we e 17 In e n Ag 0 Ao e n we 10 k If the symmetric strong discontinuity formulation is used r nO ARO de FE AGO avo else if non symmetric AZ APRO Ow oh OLOT nx e ee Wve AS eave 1 Compute residual rk gl ple L k k rasal lras y rh Ot m Assemblage of internal force vector pni _ A Fin eE n Calculation of characteristic length ne E j gt cos 0 sin 0 a le _calc
55. ints CREATE NODE 1 DISP_X REAL 0 0 DISP_Y REAL 0 0 CREATE NODE 3 DISP_X REAL 0 0 DISP_Y REAL 0 0 CREATE NODE 6 DISP_X REAL 0 0 DISP_Y REAL 0 0 CREATE NODE 11 DISP_X REAL 0 0 DISP_Y REAL 0 0 CREATE NODE 16 DISP_X REAL 0 0 DISP_Y REAL 0 0 New curve by points NEW_CURVE NAME horiz_disp curve label TYPE BYPOINTS curve type existing types gt crv ctype BEGIN_PARAMETERS 0 0 0 0 100 0 0 3 END_PARAMETERS END_NEW_CURVE Displacement by curves CREATE NODE 40 DISP_X WORD horiz_disp FACTOR_DISP_X 1 CREATE NODE 42 DISP_X WORD horiz_disp FACTOR_DISP_X 1 CREATE NODE 43 DISP_X WORD horiz_disp FACTOR_DISP_X 1 CREATE NODE 44 DISP_X WORD horiz_disp FACTOR_DISP_X 1 CREATE NODE 45 DISP_X WORD horiz_disp FACTOR_DISP_X 1 END_DOF_CONDITIONS END_BOUNDARY_CONDITIONS Strain localization data BEGIN_LOCALIZATIONDATA DIS_INJ_TYPE SD_INJ injection type existing types gt dis_inj_type TAU_MXD 0 1 tau for regular zone 0 1 TAU_LOC 0 0 tau for localized zone 0 1 Q_INJECTION 0 95 fraction of q to inject reduced integration 0 1 Q_INF 0 1 fraction of q to inject strong discontinuity 0 1 NORMAL_REFERENCE NX 0 70710678 NY 0 70710678 reference normal vector PARAK 0 01 fraction of element size for strong discontinuity band width 0 1 END_LOCALIZATIONDATA Plots definition BEGIN_PLOTS New plot NEW_PLOT X AXIS Y AXIS DISPLACEME
56. isplacements The Point Displacements condition should be used to impose some known displacement to any point in the mesh The displacement can be defined by a 2 3 4 or 5 points curve 3 3 4 Materials window Access it with the top menu Data gt Materials TSG 910 xe LEMME o 0 x w e Material type ELASTIC Material type ELASTO PLASTIC2 m Section type PLANE STRAIN M Section type PLANE STRAIN ba Young modulus 100 Elastic Poisson ratio 0 30 Submaterial type 1 ELASTIC_LINEAR Young modulus Poisson ratio Submaterial type 2 Sigma u Hardening type 100 0 30 Plastic J2 IMPLEX v 100 LINEAR Gf 8 0 Assign Draw Unassign Exchange Assign Draw Unassign Exchange Close Close a Elastic b Elastop plastic Figure 3 5 Materials window for commatfail problemtype 33 Chapter 4 User example for commatfail 4 1 Problem to solve Here we will sole step by step the problem sketched on the Figure 4 1 Figure 4 1 Sketch of the problem to solve 4 2 Generating the input file The following steps will guide the user through the geometry generation and the assignment of the required properties for this example All these steps should be done in GiD after the installation of commatfail problemtype 1 Create a new file File gt New 2 Create the geometry a Geometry gt Create gt Object gt Rectangle
57. jection implementation in commatfail vi Error calculation F PEN pont T punt R VR R 1 F l EF ss lell 2 A End loop over time step iterations if ep lt tol and lt tol Write plot info wroutplot_e2dstd g Actualize converged values actimstg_2 h Write post process info wroutpost i Write restart file restart_wr End loop over time steps when stage is over 3 Write output results of last step wroutpost 5 2 Subroutine cfint_e2dstd_sssdlt_wsda 5 2 1 Description Computes the vector of internal forces of elements of type E2DSTD for the symmetric and non symmetric strong discontinuity approach The superscript e in equations indicate an elemental quantity 5 2 2 Procedure 1 Begin loop over all elements a Read previously calculated values 19 9 Au Ju Vo n b Assign mixed formulation values depending on the injection stage r 0 0 Ve Qings Qing C Qmaa Qdis input dat gt gt TAU_LOC 0 0 ri 1 0 Ve Q Qinj general constitutive models 7 0 1 Ve MQinj isochoric constitutive models input dat gt gt TAU_MXD 0 1 E 0 0 Ve O Opis O 1 0 Ve Opis c Compute integration variables on each of the 6 gauss points 49 commatfail user guide SS NS WN WN NS det J calc_jacob AV wgp det J gp ngp vo Y avg gp 1 y av de e Calculate stress for gauss point 5 reduced
58. ke a copy from them in an auxiliary data file called commatfail dat The auxiliary data file is unique a single file and this is the file that commatfail reads during im plementation However this is a temporary file to be deleted automatically after the execution of commatfail 1 6 Code comments commatfail allows you to enter any comments deemed necessary in order to help or clarify or interpret the data file However comments will be given by the program The comments in the data file are marked by the characters Or Rules The character should be located in the first column which indicates the program that the entire row is a comment The character can be located anywhere on the line after any command telling the program from that position until the final line is a comment soy Obviously the character can also be located in the first column to perform the same functions as the character so this is a more general command Chapter 2 Input data file for commatfail 2 1 Creating input data file BEGIN_PROBLEM and END_PROBLEM identify respectively the beginning and the end of the actual data file BEGIN_PROBLEM lt problem body data gt END_PROBLEM e BEGIN_PROBLEM is the keyword to specify the beginning of the actual data file e END_PROBLEM is the keyword to indicate the ending of the actual data file Rules The definition of this field is mandatory in every data
59. l_e2dstd GET_ELEM_GEOMETRY CALL get_gtype_e2dstd GET_ELEM_TYPE strng eset etype gt GET_SECTION_TYPE CALL get_section_e2dstd gt GET_NAME_LIST CALL inc_vstg REMESH_CHECK EXIT gt CALC_EMASS CALL emass_e2dstd gt DENSITY_CHECK CALL denschk_e2dstd gt GET_NODE_LIST CALL get_nodelist gt GET_NGBR CALL get_ngbr gt wroutplot CALL wroutplot_e2dstd gt CONDS_RESID CALL conds_resid_e2dstd Continue on next page 58 Appendix A Interface subroutines Table A 3 continued from previous page Case Action gt ACT_CONDSVAR CALL act_condsvar_e2dstd PROD_CONDSVAR CALL prod_condsvar_e2dstd A 4 Interface elemset A 4 1 Description Choose the following interface of calls for specific tasks depending on the type of finite element A 4 2 Cases Table A 4 Options in interface elemset Case Action elmtyp 1 E2DSTD CALL e2dstd elmtyp 2 E3DSTD CALL e3dstd elmtyp 3 E2PFEM CALL e2dpfem elmtyp 4 E3PFEM CALL e3dpfem elmtyp 5 E2DMIX CALL e2dmix elmtyp 6 E2DMXD CALL e2dmxd A 5 Interface calc_celas A 5 1 Description Choosses between two subroutines that calculate the constitutive tensor for linear elastic material according to the type of problem to solve A 5 2 Cases Table A 5 Options in interface calc_stress Case Action ELASTICLINEAR gt PLANESTRAIN CALL calc_isolin celas E
60. lared It can contain any command or data on the subset of nodes and any command or data located outside shall not be interpreted as a figure corresponding to the database of the subset of nodes The label that identifies the subset of nodes is introduced through the NAME command within this field is required The inclusion of a node in a subset does not mean that the node should be declared in the corresponding field to the set of nodes database The data for the subset of nodes labels for the discretization are introduced into single field parameters BEGIN_PARAMETERS and END_PARAMETERS Example NEW_NODE_SUBSET NAME node subset 1 BEGIN_PARAMETERS 1 aor WN 6 END_PARAMETERS END_NEW_NODE_SUBSET 2 14 Establishment of the Basis of Elements BEGIN_ELEMENT_SET and END_ELEMENT_SET are the key words that identify the beginning and end of all the additional data associated with the set of elements such as type of element connectivities material etc BEGIN_ELEMENT_SET DELETE lt label for the set of elements gt NEW_ELEMENT_SET lt body that defines all or part of the set of elements gt END_NEW_ELEMENT_SET END_ELEMENT_SET P e BEGIN_ELEMENT_SET is the keyword that indicates the beginning of the data for the base elements e END_ELEMENT_SET is the keyword that indicates the end of the data for the base elements e DELETE This is the keyword to specify the name of the set of elements tag to be remov
61. lmr 2 Qbif elmr 3 ni x2 Component x of the first vector n from bifurcation analysis elmr 4 nly Component y of the first vector n from bifurcation analysis elmr 5 n x Component x of the second vector n from bifurcation analysis elmr 6 nily Component y of the second vector n from bifurcation analysis elmr 7 k Band width for Strong Discontinuity elmr 8 Vul Norm of the displacements gradient elmr 9 Vu z Vuzy Normalized displacement gradient vector to select n elmr 10 Vu y Vuzy Normalized displacement gradient vector to select n elmr 11 de 4V 6 i Used for integration in the internal force calculations elmr 12 nn x Component x of the vector n selected from n and na elmr 13 nn y Component y of the vector n selected from n and na elmr 14 nma 2 Component x of the vector perpendicular to n selected from n and n2 elmr 15 nmyj y Component y of the vector perpendicular to n selected from n and n elmr 16 inj elmr 17 sda elmr 18 qq elmr 19 gn Internal variable to make the double smoothing process elmr 20 V x Component x of function y gradient elmr 21 Vy y Component y of function y gradient elmr 22 This field is not used for any variable elmr 23 GG elmr 24 Vyul Norm of the internal variable gradient elmr 25 Ve 2 Vyl Component x of normalized internal variable gradient vector to select n elmr 26 Ve y VyI Component y of normalized internal variable
62. lowing stages Moreover if the data of previous stages is already stored in a memory once it reads back the data file the information it had in memory is overwritten This field contains an independent basis curves that is managed with the commands DELETE to remove curves and NEW_CURVE END_NEW_CURVE to create curves and you can define as many times as needed However this base curve is only accessible through the field BEGIN_DOF_CONDITIONS and END_DOF_CONDITIONS Example BEGIN_DOF_CONDITIONS NEW_CURVE NAME disp TYPE BYPOINTS BEGIN_PARAMETERS 0 00 0 00 1 00 0 25 2 00 1 00 END_PARAMETERS END_NEW_CURVE CREATE NODE DISP_X CREATE NODE DISP_X DELETE NODE DISP_Y DELETE NODE 15 DISP_X DISP_Y DELETE NODE 16 DISP_Y CREATE NODE 3 DISP_X WORD disp DISP_Y REAL 1 0 FACTOR_DISP_X 2 0 CREATE NODE 9 DISP_X WORD disp END_DOF_CONDITIONS REAL REAL 0 0 DISP_Y REAL 0 0 0 0 ENN 2 18 Creating the localization data module BEGIN_LOCALIZATIONDATA and END_LOCALIZATIONDATA identify respectively the beginning and the end of the strain injection parameters BEGIN_LOCALIZATIONDATA DIS_INJ_TYPE lt injection type gt TAU_MXD lt tau for regular zone gt TAU_LOC lt tau for localization gt Q_INJECTION lt q for injection gt Q_INF lt q at infinity gt NORMAL_REFERENCE NX lt x component gt NY lt y component gt PARAK lt bandwidth fraction gt END_LOCALIZATIONDA
63. matfailpost lt input_file gt The user could write a simple script to execute those commands more efficiently especially when more than one simulation is to be done sequentially 1 1 Modules and database Once the input file is supplied there are some subroutines in commatfail which reads all the problem data and create a temporary database to check whether all required fields are complete If so an internal database is created to be used later while the calculations are done This database is stored in modules and its main role is to communicate among the calculation subroutines and keep the variables in one place These modules can be found in the files whose name ends with _db 90 Some of those modules are e elmset_db f90 stores the information on element sets Note that the input and output files should be generated and read with some external software capable of pre process and post process a finite element problem This implies that variables are global and not local This brings advantages e g not use many variables as input and output arguments in all subroutines commatfail user guide e nodeset_db f90 stores the information on nodes as the coordinates e mat_db f90 stores the information of materials properties 1 2 Stages A simulation in commatfail can be done by stages which may have different conditions of loads or element properties Suppose that we The type of stages allowed are e N
64. mpost examplel 91 99008008 92 99008088 93 99008088 94 90008088 95 o8008088 96 o8008088 97 99008088 98 o8098088 99 99008088 109 08009008 a Begin of process b End of process Figure 4 10 Processing the problem 4 6 Showing the results The additional files created during the process and post process phases contain the information to show the results in a more user friendly way 4 6 1 Curve plots Curve plot data for node 43 are stored in the files example1 c01 The plot of those data is shown in Figure 4 11 44 Chapter 4 User example for commatfail 30 T T T examplel c01 Force 0 li li ji l 0 0 05 0 1 0 15 0 2 0 25 0 3 Displacement Figure 4 11 Force Displacement curve plot for node 43 4 6 2 GiD outputs The user have to change GiD to the post process window in order to see the results This can be done with the top menu Files gt Postprocess or with the corresponding icon in the menu bar Once the flavia file is loaded GiD can show many results Some of them are the ones in the figures below Deformed shape can be shown with the top menu View results gt Deformation gt TOTAL DISP y Le Gp Gp a Undeformed shape b Deformed shape Figure 4 12 Graphic displacement results The crack path is shown by plotting the double smooth variable View results gt Contour lines gt Double smoth Then the minimum and maximum values should
65. nd loop over elements to choose n 55 Appendix A Interface subroutines This chapter shows some subroutines of commatfail that serve as interfaces between general proce dures and specific calculus These interfaces are a comprehensible way to determine which specific procedure will be used For example the computation of stresses in a small strain and small deformation formulation can be different for each constitutive material model Then the interface calc_stress_sssd is used here to indicate the procedure to follow A 1 Interface calc_stress A 1 1 Description Chooses the next interface whether the problem is based on small strain and small displacement or not As commatfail is still in a development stage it has only the option for small strains and small displacements but other formulations can be implemented A 1 2 Cases Table A 1 Options in interface calc_stress Case Action gt SMALLSTRASMALLDISPLA CALL calc_stress_sssd gt SMALLSTRASMALLDISPLT CALL calc_stress_sssd gt SMALLSTRASMALLDISPLTBSTD CALL calc_stress_sssd gt SMALLSTRASMALLDISPLTBBAR CALL calc_stress_sssd gt SMALLSTRASMALLDISPLTWSDA CALL calc_stress_sssd A 2 Interface calc_stress_sssd A 2 1 Description Call specific subroutines to calculate stress for small strains and small dispacements formulation A 2 2 Cases Table A 2 Options in interface calc_stress_sssd Case Action ELASTIC CALL calc_stress
66. ndicates the beginning of the data corresponding to the definition of the conditions imposed on the degrees of freedom e END_DOF_CONDITIONS is the keyword that indicates the end of the data corresponding to the definition of the conditions imposed on the degrees of freedom 24 Chapter 2 Input data file for commatfail e NEW_CURVE Keyword that paves the way to interpret a curve see Definition field curves e END_NEW_CURVE Key word here closes the field of curves which specify the conditions on the contour of degrees of freedom e DELETE This is the keyword to specify that it will eliminate a curve or a contour of the database This option must be accompanied by one of the following keywords CURVE Keyword to specify removal of a curve This word should be followed by the label of the curve you want to delete RESET_ALL_DOF Keyword to specify removal of all conditions imposed on the degrees of freedom defined in steps stages above as well as curves that are associated NODE Keyword to specify removal of conditions imposed on a node This keyword must be followed by the number of the node followed by the word s key that identifies the degree of freedom to be eliminated You can delete all or some of the contour of the node NODE_SUBSET Keyword to specify removal of conditions imposed on all the nodes that comprise the subset of nodes specified This keyword should be followed by the label that identif
67. need some extra configurations for its correct operation General convention for character fonts The meaning associated to specific font styles is given below Standard typos are not explained here and some exceptions could be highlighted in the text and not here Description Example Scalar values and functions a a A Vectors a n Second order tensors Bio Forth order tensors C S Spaces Uu V Fortran codes plain text file names lt input gt dat NEW_MATERIAL Key or button pressing indications Apply Software menus File gt Save as Data gt Problem data vi Chapter 1 Overview of commatfail commatfail is a multipurpose finite element code written in Fortran 90 As many finite element codes it has only the calculus process of a problem That means that all needed information is taken from an input file before the computations start and all results are written in the output files when finished Figure 1 1 shows this basic way of how commatfail works input_file output_file Figure 1 1 commatfail basic operating process To perform a simulations with commatfail 1 Place the input file and the executable commatfail file in the same directory 2 Navigate to that directory through a command window 3 Execute commatfail lt input_file gt After the calculations are done post process can be done by 1 Place the executable file commatfailpost in the same directory 2 Execute in the command window com
68. ns each one of which would be associated respectively with the elements of the mesh and with the material they are made NEW_SECTION NAME lt section label gt TYPE lt type of section gt BEGIN_PARAMETERS MATERIAL lt associated material label gt lt body that define the properties of section gt END_PARAMETERS END_NEW_SECTION e NEW_SECTION is the keyword that indicates the beginning of the data that characterize the section e END_NEW_SECTION is the keyword that indicates the end of the data that characterize the section e NAME is the keyword to specify the section name label that will identify in the database of the sections Must be preceded by a word that will gain the category label It is a mandatory parameter and is recommended to be located at the beginning of the field e TYPE is the keyword to specify the type of section and therefore that data is needed and should be interpreted Must be preceded by a keyword predefined by the program according to the type of section It is a mandatory parameter and is recommended to be located at the beginning of the field with the only restriction to be observed that it is before the field BEGIN_PARAMETERS Currently the code has the following options MATERIAL_SECTION Elements 3D elements in general and that they only need to associate a material PLANE_STRAIN Elements 2D plane strain problems PLANE_STRESS Elements 2D plane strain problems A
69. p the program Example NEW_MATERIAL NAME material_1 TYPE ELASTOPLASTIC NEW_SUBMATERIAL END_NEW_SUBMATERIAL END_NEW_MATERIAL 2 7 Creating the sub material properties NEW_SUBMATERIAL and END_NEW_SUBMATERIAL are respectively the key words that identify the begin ning and the end of the data set that defines the characteristics of the material or the submodel e g the elastic model the surface creep hardening laws etc NEW _SUBMATERIAL NAME lt label of_sub material gt TYPE lt type of_sub material gt BEGIN_PARAMETERS lt properties of the sub material gt END_PARAMETERS NEW_CURVE lt curves that define the points of the sub material gt END_NEW_CURVE NEW_SUBMATERIAL lt defining nested sub material gt END_NEW_SUBMATERIAL END_NEW_SUBMATERIAL e NEW_SUBMATERIAL is the keyword that indicates the beginning of the data that characterize the sub material If this keyword is declared within the field itself it implies declaring another submodel of material nested in the actual submodel 12 Chapter 2 Input data file for commatfail e END_NEW_SUBMATERIAL is the keyword that indicates the end of the data that characterize the sub material If there is another submodel within the stated submodel then this keyword closes the later field e NAME is the keyword to specify the name of the material submodel label that will identify the sub material Must be preceded by a word that will gain
70. parametric particle element PFEM E3DSTD 3D standard isoparametric element E3DPFEM 3D isoparametric particle element PFEM e FORMULATION This is the keyword that indicates the type of formulation used for the set of elements in question and must be preceded by a keyword predefined by the program according to the formulation Currently the program takes into account the following SMALL_STRA_SMALL_DISP_LA Formulation for an updated Lagrangian with small defor mations and small displacements SMALL_STRA_SMALL_DISP_LT Formulation for a total Lagrangian with small deformations and small displacements LARGE_STRA_LARGE_DISP_LA Formulation for an updated Lagrangian with large deforma tions and large displacements SMALL_STRA_SMALL_DISP_LT_WSDA Formulation for a total Lagrangian with small defor mations and small displacements plus the capacity to generate weak and strong discon tinuities e SUBTYPE the keyword to specify the type of element within the category of elements E2DSTD and must be preceded by a keyword predefined by the program according to the type of item Currently the program takes into consideration the following types LINEAR_TRIANGLE 3 node triangular element LINEAR_QUAD Element 4 square nodes e QUADRATURE is the keyword to specify the quadrature type used for integration and should be preceded by a keyword predefined by the program depending on the quadrature rule Cur
71. r example for commatfail Materials Materials EEE p Material 2 Bi lglo xXh a Material 2 BOK 2 Material type ELASTO_PLASTIC2 Material type ELASTO_PLASTIC2 Section type PLANE STRAIN v Section type PLANE STRAIN v Elastic Elastic Submaterial type 1 ELASTIC_LINEAR Submaterial type 1 ELASTIC_LINEAR Young modulus 10000 Young modulus 10000 Poisson ratio 0 45 Poisson ratio 0 45 Plastic Plastic Submaterial type 2 J2 IMPLEX v Submaterial type 2 J2 IMPLEX ad Sigma u 100 Sigma wfo Hardening type LINEAR Hardening type LINEAR Assign Draw Unassign i Exchange Assign Draw Unassign A Exchange Close Close a Material 1 b Material 2 Figure 4 8 Creation of two new material of elasto plastic type b Assign Material 2 to element number 4 c Assign Material 1 to all other elements 11 Generate input file Files gt Export gt Calculation file 4 3 Input file The input file to run this example was generated with commatfail probemtype and is shown below It can also be found as example1 dat in the folder examples that comes with the download File generated by commatfail gid problemtype Input file for problem example 1 BEGIN_PROBLEM Basic problem information BEGIN_PROBLEM_BASICS TITLE example 1 SCRLEVEL 3 screen messages level RSNLEVEL 3 rsn file messages level NDIME 2 problem dimension STAGE_TYPE IMPLICIT_NR2 stage type existing types gt stgtype SOLVER
72. rently the program takes into consideration the following types GAUSS Gauss Legendre Quadrature e NGAUS is the keyword to specify the number of integration points and should be preceded by a positive integer e MGAUS is the keyword to specify the number of integration points of the mass matrix and should be preceded by a positive integer e BEGIN_PARAMETERS Keyword employee who opened the field to interpret the data that define the set of elements such as grid connectivity e END_PARAMETERS Keyword closing the field used to interpret the data that define the set of elements Rules This field can contain any command or data on the set of elements and any command or data located outside shall not be interpreted in relation to a database of sets of elements It is mandatory to define a label that identifies the set of elements which is achieved through the NAME command in this field Failure generates an error in input reading that will stop the program It is also mandatory to define the type of elements set by the TYPE command the omission of which also generates a read error that stops the execution Depending on the type of element some additional commands to be used to define the set of elements Furthermore data from the set of elements for the discretization are introduced into a single field parameters BEGIN_PARAMETERS and END_PARAMETERS The FORMULATION SMALL_STRA_SMALL_DISP_LT_WSDA is only implemented for the element
73. screen The lowest integer is 0 and there is no upper limit The higher the level the greater the number of messages per screen Thus if the level is 0 during the execution minimum number of messages will be written e RSNLEVEL the command is optional and can be located anywhere in the field Problem Basics The declared value in this field will be valid for the whole problem unless explicitly specified in a Stage see Creation Basics Field Stage The whole number into RSNLEVEL indicates the level of messages you want to get in the file rsn The lowest integer is 0 and there is no upper limit for it The higher the level the greater the number of messages sent to the file rsn Thus if the level is O during the execution minimum number of messages will be written e NDIME this command is required Example BEGIN_PROBLEM BASICS TITLE 5061 piece with a sequence of three movements SCRLEVEL 1 RSNLEVEL 1 NDIME 2 STAGE_TYPE IMPLICIT_NR SOLVER SuperLU MAX_ITERATION 100 TOLERANCE 1d 6 END_PROBLEM_BASICS 2 3 Creating the stage NEW_STAGE and END_NEW_STAGE are the key words that identify the beginning and the end of the field respectively with the data describing the Stage commatfail user guide NEW_STAGE lt Stage body data gt END_NEW_STAGE e NEW_STAGE is the command or keyword used to indicate the beginning of the field e END_NEW_STAGE is the command or keyword used to indicate the end of the
74. tabase and must be preceded by a word that will gain the category label It is a mandatory parameter and is recommended to be located at the beginning of the field e TYPE is the keyword to specify the type of material if material is a standard or defined by the user etc and must be preceded by a keyword predefined by the program according to the type of material It is a mandatory parameter and must be located at the beginning of the field before the definition of submaterials e NEW_SUBMATERIAL Keyword which opens the field of sub material which defines at least a portion of the material characteristics such as the elastic model the surface creep and so on see submaterials field e END_NEW_SUBMATERIAL Keyword which closes the sub material field Rules This field can contain any command or data on the parameters of the material and any information on the material located outside the field shall not be interpreted It is mandatory to define a label identifying the material which is achieved through the NAME command in this field Failure generates an error in reading inputs that will stop the program To define the material properties one can define as many sub fields of material as needed It is mandatory to define the type of material before the submaterials which is achieved through the TYPE command at the beginning of this field Omission of this or statement after submaterial generates an error in reading inputs that will sto
75. the TYPE command the omission of which also generates a reading error that stops the execution The data of the curve are introduced in a single field parameter BEGIN_PARAMETERS and END_PARAMETERS Example NEW_CURVE NAME curve_1 TYPE BY_POINTS BEGIN_PARAMETERS Curve defined by 3 points 0 0 0 0 1 0 20 0 2 0 20 0 END_PARAMETERS END_NEW_CURVE 14 Chapter 2 Input data file for commatfail 2 9 Establishment of the basis of sections BEGIN_SECTION and END_SECTION are the key words that identify the beginning and end of all the additional data associated with the element including the reference material and are not directly related to the geometric definition of the mesh BEGIN_SECTION DELETE lt section label gt ASSOCIATE SECTION lt section label gt MATERIAL lt material label gt NEW_SECTION lt body that defines all or part of the base sections gt END_NEW_SECTION END_SECTION e BEGIN_SECTION is the keyword that indicates the beginning of the data for the sections e END_SECTION is the keyword that indicates the end of the data for the sections e DELETE This is the keyword to specify the section name label to be removed from the database Must be preceded by the label of the section to remove e ASSOCIATE is the keyword to specify the link to the material section This option must be accompanied by keywords SECTION followed by the label of the section to associate
76. types Linux lt GiD_directory gt lt GiD_version gt problemtypes MacOS Applications lt GiD_directory gt lt GiD_version gt problemtypes Restart GiD to complete installation 3 2 Loading the commatfail problemtype To load commatfail problemtype interface choose from the top menu Data gt Problem type gt commatfail When pre process is finished the user can generate the commatfail input file by choosing from the top menu Files gt Export gt Calculation Files and save the dat input file 29 commatfail user guide GiD x64 Project example comfm Files View Geometry Utilities Data Mesh Calculate Help New Ctrl x Ctrl n 2 Q Layero m Open Ctrl o Recent projects Y Save Ctrl s 3 Save as Ctrl x Ctrl s Project UNNAMED comfm Import H lities Data Mesh Calculate ANSYS Help IGES m DXF _Probiemtype 0 anes emcee pr a Recent poss mee Rhinoceros Conditions amples 5 RI Materials e Print to file l oa Interval Data SAP2000 a Page pil cari satinga Amelet EBProblem Data Print Interval Transform quit aria Text data report Internet Retrieve Y ASCII project Local axes Load iai K ON layers Unload ae Debugger A a Ge ataate J a load problemtype Figure 3 1 Loading the problemtype and exporting the dat input file b export file Within this pre processor it is either possible to import a CAD g
77. umber of messages will be written The command STAGE_TYPE is mandatory if it has not been declared in the Problem Basics However compulsory SOLVER_ TYPE command depends on the integration scheme employed Commands such as NUMBER_STEPS NUMBER_POST PERIOD_POST and END_TIME are unique to the stage Commands NUMBER_POST and PERIOD_POST are optional and mutually exclusive In case both are set simultaneously the information defined by the command NUMBER_POST precedes over that defined by the command PERIOD_POST In the event that the stage does not define any of these commands only the final stage will be post processed Example BEGIN_STAGE_BASICS TITLE Stage pressing First movement of punches SCRLEVEL 3 RSNLEVEL 3 SOLVER SuperLU MAX_ITERATION 100 TOLERANCE 1d 3 NUMBER_STEPS 100 NUMBER_POST 3 End_time 1d 2 END_STAGE_BASICS 2 5 Creating the material data base BEGIN_MATERIAL and END_MATERIAL are respectively the key words that identify the beginning and the end of the materials database BEGIN_MATERIAL DELETE lt label of the material gt 10 Chapter 2 Input data file for commatfail NEW_MATERIAL lt body defining new material data gt END_NEW_MATERIAL END_MATERIAL e BEGIN_MATERIAL is the keyword indicating the starting of the material data e END_MATERIAL is the key word to indicate the end of the material data e DELETE This is the keyword to specify the name of material
78. used to indicate the type of dynamic integration to be used for the stage and should be preceded by a keyword Currently the code supports the dynamic integration STATIC Keyword to specify a static calculation default option Requires no additional parameters G_ALPHA Keyword to specify a dynamic calculation using the general Alpha method This method requires the definition of additional parameters MESH_POST is the command or keyword to specify the option type in the post processing stage If the option is YES is post processed all meshes although these results are not partners If the option is NO only post processed mesh containing those results in which case it overwrites the filenames that are not associated any results BEGIN_DYNAMIC_PARAMETERS keyword which opens the field to interpret the parameters of the dynamic integration commatfail user guide e END_DYNAMIC_PARAMETERS keyword closing the field which is used to interpret the parameters of the dynamic integration Rules The field of STAGE_BASICS is declared within a Stage Field if declared in a different field will cause an error and can be located anywhere in the field although it is strongly recommended that it is the first stage of the field or in case some messages may be lost during reading It can contain any command of general variables and any command on this field if written outside this field will not be interpreted as Stage Basics No other
79. whole number e TOLERANCE is the command or keyword to specify the value of tolerance required for conver gence and it must be preceded by a real number e MESH_POST is the command or keyword to specify the option of post process for different stage If the option is YES the results will be post processed for all meshes although these results are not partners If the option is NO only post processed mesh containing those results in which case it overwrites the filenames that are not associated any results The default option is NO Rules This is a mandatory field whose absence will create an error that will stop the program This statement is made in the data input field if stated in a different field will cause an error and can be located anywhere in the field although it is recommended to write at the beginning It can contain any command on the general variables of the problem and any command outside this field will not be interpreted as Problem Basics Problem Basics no one field to another field may contain either one Problem Basics Field Stage e SCRLEVEL the command is optional and can be located anywhere in the field Problem Basics although it is recommended to be in the starting of the field The declared value in this field will be valid for the whole problem unless otherwise explicitly specified in a Stage see Creation Basics Field Stage The whole number into SCRLEVEL indicates the level of messages you want to get in the
80. window some actions over stages are e Data gt Interval gt New creates a new load stage e Data gt Interval gt Current defines which load stage is the actual one e Data gt Interval gt Delete deletes a load stage Stage basics and Strain localization data cards are available 31 commatfail user guide Stage basics e Title stage Write in this field a title to the stage It serve the user to remember what happens in that stage or to distinguish it from the others This field can contain spaces and it is strongly recommended to use only alpha numeric characters e Process steps Specify here the number of steps of the stage Increasing this number will increase both the time of the analysis and the quality of the results e Post process steps Specify here how many of the process steps will be in the post process It must never be higher than the process steps e End time Specify the final time for the stage so the analysis can compute the size of each time step in the actual stage Strain localization data e Injection type Choose here the type of injection for the material failure process e Tau mixed White the value of the stabilization parameter for the regular zone The moment of the first strain injection bifurcation time e Tau localized White the value of the stabilization parameter for the localized zone The second injection strain discontinuity e Q injection Specify the fraction of q to inj
81. y advanced user It includes all infor mation to be introduced to monitor the tasks performed by the code You may have even intended to be used in debug mode e lt name gt rsf An internal file in commatfail to save all the information necessary to relaunch the calculation either by an unwanted interruption of the program or to launch a new phase or stage e lt name gt sun It is a binary file that contains all the results of post processing problem in the form of commatfail which must then be interpreted by the module commatfailpost e lt name gt nlb It is an ASCII file that contains for each stage all the changes to the labels of the nodes that overlap with other nodes automatically added by commatfail or added by the user in previous stages e lt name gt c01 lt name gt c02 etc They are text files that contain all data for plotting the plots defined in the input file e lt name gt _s1 flavia This is a binary file that contains all the results in the correct format to be read by pre post process tool GiD It should be open from GiD Rules Although the operating system allows it commatfail does not allow path or file names contain spaces 1 5 File commands e ImportFile lt file_name gt Rules The command ImportFile can be located anywhere in the data file However it must be the only command line One of the first tasks undertaken by commatfail is to interpret the file or data file provided by the user and ma
82. y the word key that identifies the degree of freedom to be prescribed after the KEY is specified the prescription that can be defined through a numerical value or through a curve In the case of a numeric value type the keyword followed by the REAL value whereas if it is a curve type the keyword followed by the label WORD of the curve In addition you can specify a factor not required amending the prescribed value To do this use the keyword composed of the root FACTOR and the keyword that specifies the degree of freedom followed by the numeric value that defines the factor As was the case earlier on this group of nodes subset and on the same line of data file one or more conditions can be created on the contour subset of nodes each of which corresponds to the degree of freedom to be restricted e Currently the program considers the following keywords to identify the degrees of freedom DISP_X Displacement in the X direction or first coordinate axis DISP_Y Displacement in the Y direction or second coordinate axis 25 commatfail user guide DISP_Z Displacement in the Z direction or third coordinate axis Rules This field allows you to communicate with the database that defines the conditions imposed on the contour of degrees of freedom of the problem allowing you to create modify or delete these conditions In this connection it must be borne in mind that if not erased they remain in memory during the fol

Computational Material Failure

Contents

Download Pdf Manuals

Related Search

Related Contents