要用到的用户子程序.doc

1.1.3 DFLUX: User subroutine to define nonuniform distributed flux in a heat transfer analysis.应用于定义非均匀分布热流密度。Product: Abaqus/StandardReferences “Thermal loads,” Section 30.4.4 of the Abaqus Analysis Users Manual “Mass diffusion analysis,” Section 6.9.1 of the Abaqus Analysis Users Manual *DFLUX *DSFLUX “DFLUX,” Section 4.1.1 of the Abaqus Verification ManualOverviewUser subroutine DFLUX: can be used to define a nonuniform distributed flux as a function of position, time, temperature,element number, integration point number, etc. in a heat transfer or mass diffusion analysis; will be called at each flux integration point for each element-based or surface-based (heat transferonly) nonuniform distributed flux definition in the analysis; ignores any amplitude references that may appear with the associated nonuniform distributed fluxdefinition; and uses the nodes as flux integration points for first-order heat transfer, first-order coupled temperaturedisplacement,and mass diffusion elements.User subroutine interfaceSUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,1 JLTYP,TEMP,PRESS,SNAME)CINCLUDE ABA_PARAM.INCCDIMENSION FLUX(2), TIME(2), COORDS(3)CHARACTER*80 SNAMEuser coding to define FLUX(1) and FLUX(2)RETURNEND1.1.31DFLUXVariables to be definedFLUX(1)Magnitude of flux flowing into the model at this point. In heat transfer cases the units are JT1L2 forsurface fluxes and JT1L3 for body flux. In transient heat transfer cases where a non-default amplitudeis used to vary the applied fluxes, the time average flux over the time increment must be defined ratherthan the value at the end of the time increment. In mass diffusion cases the units are PLT1 for surfacefluxes and PT1 for body flux.FLUX(1) will be passed into the routine as the magnitude of the flux specified as part of theelement-based or surface-based flux definition. If the magnitude is not defined, FLUX(1) will bepassed in as zero.This flux is not available for output purposes.FLUX(2)In heat transfer cases: , the rate of change of the flux with respect to the temperature at this point.The units are JT1L2 1 for surface fluxes and JT1L3 1 for body flux.Inmass diffusion cases: , the rate of change of the flux with respect to themass concentrationat this point. The units are LT1 for surface fluxes and T1 for body flux.The convergence rate during the solution of the nonlinear equations in an increment is improvedby defining this value, especially when the flux is a strong function of temperature in heat transferanalysis or concentration in mass diffusion analysis.Variables passed in for informationSOLEstimated value of the solution variable (temperature in a heat transfer analysis or concentration in amass diffusion analysis) at this time at this point.KSTEPStep number.KINCIncrement number.TIME(1)Current value of step time (defined only in transient analysis).TIME(2)Current value of total time (defined only in transient analysis).NOELElement number.1.1.32DFLUXNPTIntegration point number in the element or on the elements surface. The integration scheme dependson whether this is a surface or a body flux.COORDSAn array containing the coordinates of this point. These are the current coordinates if geometricnonlinearity is accounted for during the step (see “Procedures: overview,” Section 6.1.1 of the AbaqusAnalysis Users Manual); otherwise, the array contains the original coordinates of the point.JLTYPIdentifies the flux type for which this call to DFLUX is being made. The flux type may be a body flux,a surface-based flux, or an element-based surface flux. For element-based surface fluxes, this variableidentifies the element face for which this call to DFLUX is being made. This information is useful whenseveral different nonuniform distributed fluxes are being imposed on an element at the same time. SeePart VI, “Elements,” of the Abaqus Analysis Users Manual for element face identification. The keyis as follows:JLTYP Flux type0 Surface-based flux1 BFNU11 S1NU (SNEGNU for heat transfer shells)12 S2NU (SPOSNU for heat transfer shells)13 S3NU14 S4NU15 S5NU16 S6NUTEMPCurrent value of temperature at this integration point (defined only for a mass diffusion analysis).Temperature for a heat transfer analysis is passed in as variable SOL.PRESSCurrent value of the equivalent pressure stress at this integration point (defined only for amass diffusionanalysis).SNAMESurface name for a surface-based flux definition (JLTYP=0). For a body flux or an element-basedsurface flux the surface name is passed in as blank.1.1.5 DLOAD: User subroutine to specify nonuniform distributed loads.用于定义非均匀分布载荷Product: Abaqus/StandardReferences “Distributed loads,” Section 30.4.3 of the Abaqus Analysis Users Manual *DLOAD *DSLOAD “Nonuniform crack-face loading and J -integrals,” Section 1.16.7 of the Abaqus BenchmarksManual “Pure bending of a cylinder: CAXA elements,” Section 1.3.33 of the Abaqus Verification Manual “Cylinder subjected to asymmetric pressure loads: CAXA elements,” Section 1.3.35 of the AbaqusVerification Manual “Patch test for axisymmetric elements,” Section 1.5.4 of the Abaqus Verification Manual “Transient internal pressure loading of a viscoelastic cylinder,” Section 2.2.8 of the AbaqusVerification Manual “DLOAD,” Section 4.1.3 of the Abaqus Verification ManualOverviewUser subroutine DLOAD: can be used to define the variation of the distributed load magnitude as a function of position, time,element number, load integration point number, etc.; will be called at each load integration point for each element-based or surface-based nonuniformdistributed load definition during stress analysis; will be called at each stiffness integration point for computing the effective axial force, ESF1, forpipe elements subjected to nonuniform load types PENU and PINU; cannot be used in mode-based procedures to describe the time variation of the load; and ignores any amplitude references that may appear with the associated step definition or nonuniformdistributed load definition.User subroutine interfaceSUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,1 COORDS,JLTYP,SNAME)CINCLUDE ABA_PARAM.INCC1.1.51DLOADDIMENSION TIME(2), COORDS (3)CHARACTER*80 SNAMEuser coding to define FRETURNENDVariable to be definedFMagnitude of the distributed load. Units are FL2 for surface loads and FL3 for body forces. F will bepassed into the routine as the magnitude of the load specified as part of the element-based or surfacebaseddistributed load definition. If the magnitude is not defined, F will be passed in as zero. For astatic analysis that uses the modified Riks method (“Static stress analysis,” Section 6.2.2 of the AbaqusAnalysis Users Manual) F must be defined as a function of the load proportionality factor, . Thedistributed load magnitude is not available for output purposes.Variables passed in for informationKSTEPStep number.KINCIncrement number.TIME(1)Current value of step time or current value of the load proportionality factor, , in a Riks step.TIME(2)Current value of total time.NOELElement number.NPTLoad integration point number within the element or on the elements surface, depending on theload type. (Stiffness integration point number while computing effective axial force, ESF1, for pipeelements subjected to load types PENU and PINU.)LAYERLayer number (for body forces in layered solids).KSPTSection point number within the current layer.1.1.52DLOADCOORDSAn array containing the coordinates of the load integration point. These are the current coordinates ifgeometric nonlinearity is accounted for during the step (see “Procedures: overview,” Section 6.1.1of the Abaqus Analysis Users Manual); otherwise, the array contains the original coordinates ofthe point. For axisymmetric elements that allow nonaxisymmetric deformation, COORDS(3) is theangular position of the integration point, in degrees.JLTYPIdentifies the load type for which this call to DLOAD is being made. The load type may be a body force,a surface-based load, or an element-based surface load. For element-based surface loads, this variableidentifies the element face for which this call to DLOAD is being made. This information is useful whenseveral different nonuniform distributed loads are being imposed on an element at the same time. SeePart VI, “Elements,” of the Abaqus Analysis Users Manual for element face identification. The keyis as follows:JLTYP Load type0 Surface-based load1 BXNU1 BRNU2 BYNU (except for axisymmetric elements)2 BZNU (for axisymmetric elements only)3 BZNU (for three-dimensional elements and asymmetric-axisymmetricelements)20 PNU21 P1NU22 P2NU23 P3NU24 P4NU25 P5NU26 P6NU27 PINU28 PENU41 PXNU42 PYNU43 PZNU1.1.53DLOADSNAMESurface name for a surface-based load definition (JLTYP=0). For a body force or an element-basedsurface load the surface name is passed in as blank.1.1.47 UTRACLOAD: User subroutine to specify nonuniform traction loads.用于定义非均匀分布切向载荷。Product: Abaqus/StandardReferences “Distributed loads,” Section 30.4.3 of the Abaqus Analysis Users Manual *DLOAD *DSLOAD “Distributed traction and edge loads,” Section 1.4.17 of the Abaqus Verification ManualOverviewUser subroutine UTRACLOAD: can be used to define the variation of the distributed traction loadmagnitude as a function of position,time, element number, load integration point number, etc.; if needed, can be used to define the initial loading direction for the distributed traction load as afunction of position, element number, load integration point number, etc.; will be called at each load integration point for each element-based, edge-based, or surface-basednonuniform distributed traction load definition during stress analysis; cannot be used in mode-based procedures to describe the time variation of the load; and ignores any amplitude references that may appear with the associated step definition or nonuniformdistributed traction load definition.User subroutine interfaceSUBROUTINE UTRACLOAD(ALPHA,T_USER,KSTEP,KINC,TIME,NOEL,NPT,1 COORDS,DIRCOS,JLTYP,SNAME)CINCLUDE ABA_PARAM.INCCDIMENSION T_USER(3), TIME(2), COORDS(3), DIRCOS(3,3)CHARACTER*80 SNAMEuser coding to define ALPHA and T_USERRETURNEND1.1.471UTRACLOADVariables to be definedALPHAMagnitude of the distributed traction load. Units are FL2 for surface loads, FL1 for edge loads, and Ffor edge moments. ALPHA is passed into the routine as the magnitude of the load specified as part ofthe element-based or surface-based distributed load definition. If the magnitude is not defined, ALPHAis passed in as zero. For a static analysis that uses the modified Riks method (“Unstable collapse andpostbuckling analysis,” Section 6.2.4 of the Abaqus Analysis Users Manual) ALPHA must be definedas a function of the load proportionality factor, . The distributed load magnitude is not available foroutput purposes.T_USERLoading direction of the distributed traction load. T_USER is passed into the routine as the loaddirection specified as part of the element-based or surface-based distributed load definition. Thevector T_USER passed out of the subroutine is used as the initial loading direction discussed in“Distributed loads,” Section 30.4.3 of the Abaqus Analysis Users Manual. The direction of T_USERas defined by the subroutine should not change during a step. If it does, convergence difficultiesmight arise. Load directions are needed only for a nonuniform general surface traction, shear surfacetraction, and general edge traction. If a direction is defined for the nonuniform normal edge traction,shear edge traction, transverse edge traction, or edge moment, it will be ignored. See “Distributedloads,” Section 30.4.3 of the Abaqus Analysis Users Manual, for details.Variables passed in for informationKSTEPStep number.KINCIncrement number.TIME(1)Current value of step time or current value of the load proportionality factor, , in a Riks step.TIME(2)Current value of total time.NOELUser-defined element number.NPTLoad integration point number within the element or on the elements surface, depending on the loadtype.1.1.472UTRACLOADCOORDSAn array containing the coordinates of the load integration point. These are the current coordinatesif geometric nonlinearity is accounted for during the step (see “Procedures: overview,” Section 6.1.1of the Abaqus Analysis Users Manual); otherwise, the array contains the original coordinates of thepoint.DIRCOSOrientation of the face or edge in the reference configuration. For three-dimensional facets the first andsecond columns are the normalized local directions in the plane of the surface, and the third columnis the normal to the face. For solid elements the normal points inward, which is the negative of whatis defined in “Conventions,” Section 1.2.2 of the Abaqus Analysis Users Manual; for shell elementsthe normal definition is consistent with the convention. For two-dimensional facets the first columnis the normalized tangent, the second column is the facet normal, and the third column is not used.For three-dimensional shell edges the first column is the tangent to the shell edge (shear direction),the second column is the in-plane normal (normal direction), and the third column is the normal to theplane of the shell (transverse direction).JLTYPIdentifies the load type for which this call to UTRA