要用到的用户子程序.doc

《要用到的用户子程序.doc》由会员分享，可在线阅读，更多相关《要用到的用户子程序.doc（12页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、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 *DS

2、FLUX “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

3、 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 f

4、or 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

5、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 t

6、he 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

7、 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 JT1

8、L2 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

9、 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 concentrati

10、on 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

11、 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

12、 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

13、 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

14、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).Temperat

15、ure 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

16、 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

17、.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 t

18、he 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

19、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 th

20、e 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

21、 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 fo

22、r 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 met

23、hod (“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)Cu

24、rrent 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 co

25、mputing 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 curre

26、nt 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

27、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 w

28、hich 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 Su

29、rface-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

30、 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

31、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 integratio

32、n 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-basednonu

33、niform 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 subr

34、outine 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

35、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

36、 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 p

37、urposes.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 discuss

38、ed 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 surface

39、traction, 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 in

40、formationKSTEPStep 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, de

41、pending 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

42、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

43、 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

44、 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 dir

45、ection).JLTYPIdentifies the load type for which this call to UTRACLOAD is being made. The load type may be anelement-based surface load, an edge-based load, or a surface-based load. This variable identifies theelement face or edge for which this call to UTRACLOAD is being made. This information is u

46、sefulwhen several different nonuniform distributed loads are being imposed on an element at the sametime. See Part VI, “Elements,” of the Abaqus Analysis Users Manual for element face and edgeidentification. The load labels are shown in Table 1.1.471. For surface- or edge-based loading(TRSHRNU, TRVE

47、CNU, EDLDNU, EDNORNU, EDSHRNU, EDTRANU, EDMOMNU), j in theload type identifies the face or edge of the element underlying the surface.Table 1.1.471 JLTYP values for surface traction and edge load labels.Load Label JLTYP Load Label JLTYP Load Label JLTYPTRSHRNU 510+j EDLDNU 540+j EDTRANU 570+jTRSHR1NU 511 EDLD1NU 543 EDTRANU 573TRSHR2NU 512 EDLD2NU 544 EDTRANU 574TRSHR3NU 513 EDLD3NU 545 EDTRANU 575TRSHR4NU 514 EDLD4NU 546 EDTRANU 576TRSHR5NU 515 EDNORNU 550+j EDMOMNU 580+jTRSHR6NU 516 EDNOR1NU 55