ANSYS Autodyn User子程序指南

版权申诉

5星 · 超过95%的资源 8 浏览量更新于2024-07-13 收藏 2.47MB PDF 举报

"Autodyn User's Subroutines Guide.pdf" 该文档是ANSYS Autodyn用户子程序指南，适用于版本19.0。Autodyn是ANSYS公司的一款强大的非线性动力学分析软件，常用于模拟复杂的工程问题，如碰撞、爆炸、冲击等。此指南详细介绍了如何利用用户自定义子程序来扩展Autodyn的功能，以满足特定的计算需求。在Autodyn中，用户子程序允许用户直接介入求解过程，实现更高级别的定制化。这些子程序通常由FORTRAN语言编写，可以控制材料模型、边界条件、积分方法等多个方面，使分析更加精确和灵活。例如，你可以通过编写用户子程序来实现自定义的材料响应，处理非线性行为，如塑性变形、蠕变或损伤积累。文档中可能涵盖了以下关键知识点： 1. **用户子程序接口**：介绍了如何将用户编写的FORTRAN代码与Autodyn内部算法集成，包括调用和返回数据的机制，以及如何定义和使用用户子程序的输入输出参数。 2. **材料模型开发**：详细解释了如何定义和实现自定义的材料模型，这可能是通过编写USERMAT子程序来实现，以模拟复杂的材料行为，如多晶材料、黏弹性材料或者非均匀材料。 3. **边界条件和加载**：用户可以通过USERSPC和USERNSUB子程序控制约束和载荷，比如创建特殊的边界条件，模拟复杂接触问题或动态加载情况。 4. **积分方法**：用户可能需要修改或扩展积分算法，这可能涉及USERINT子程序，以适应特定的数值稳定性和精度要求。 5. **后处理和结果输出**：用户可能需要自定义结果报告和图形输出，这可以通过USERPOST子程序来完成，提供个性化的结果分析功能。 6. **调试和测试**：文档会指导用户如何调试和验证用户子程序，确保它们正确地与Autodyn软件交互。 7. **实例和应用**：指南可能包含一些示例，展示了如何使用用户子程序解决实际工程问题，帮助用户更好地理解和应用这些概念。 8. **支持和资源**：文档还会提供关于获取技术支持和查找其他相关资源的信息，比如ANSYS社区、在线帮助和更新等。请注意，由于未提供具体内容，以上描述基于一般Autodyn用户子程序指南的典型结构和内容进行推测。实际文档应包含更详尽的技术细节和具体操作步骤。

Additionally, the module file fsrc\materials\mdusersub_call.f90 has been supplied. The

module in this file should not be modified.

On the Windows operating systems, the following Intel/.NET project files should be used:

…..\AUTODYN\<platform>

\ad_usrsub.sln Intel/.NET Solution

\ad_usrsub.vfproj Intel/.NET Fortran project

3.3. Formal Parameters

Most of the user subroutines have formal parameters. There are comments at the head of the subroutines

defining each parameter, stating whether the parameter is an input parameter, an output parameter,

or both. In the simplest case, the writing of a user subroutine might consist of computing the output

parameter(s) from the input parameter(s) as in the MD_STR_USER_1 example in How to Invoke User

Subroutines (p. 3). If the user subroutine requires the use of other variables that are not formal para-

meters then we must obtain them in some other manner. If a variable is global in nature (for example,

cycle number, material data, time) this information is obtained through the USE statement, referencing

the appropriate Autodyn module. If the desired variable is associated with a Part (for example, cell

pressures, velocities, and so on), these are obtained most readily by using Autodyn supplied functions

to retrieve Part data.

3.4. Autodyn Modules (Non-Parameter Data)

If you need data in addition to that supplied as formal parameters to a user subroutine, you can usually

get this data directly from the Autodyn modules, which follow the naming conventions *.mod. The

Autodyn modules for Autodyn-2D and Autodyn-3D are described in Autodyn User Subroutine Mod-

ules (p. 31). To use a particular module in a user subroutine insert the statement:

USE xxxxxx

where xxxxxx is the name of the appropriate module.

The dummy subroutines have USE statements already defined for the most often required global data.

Also, the comments indicate which variables in these modules are likely to be of use.

Caution

Be extremely careful if you modify variables other than the user variables (VAR01 through

VAR20), and that the modified values are consistent with the use of the variables. Global

variables are used by other routines and assigning bad values to them could cause Autodyn

to terminate unexpectedly.

3.5. Autodyn Variables

Autodyn uses the dynamic array allocation features available through Fortran 90/95. As such, there are

no fixed limits on the size of the model that you can generate. However, your computing time may be

greatly limited by such factors as your machine memory.

There are two basic types of data storage and access used in Autodyn depending on whether you are

using the IJK based or Unstructured solvers. These are described in the sections:

3.5.1. IJK Solvers

3.5.2. Unstructured Solvers/Elements

of ANSYS, Inc. and its subsidiaries and affiliates.12

Writing Your Own User Subroutines

3.5.3. Notes About Autodyn Variables

3.5.1. IJK Solvers

For the IJK based solvers (Lagrange, ALE, Shell, Beam, Euler, Euler-FCT, Multi-Material Euler), variables

in the model are accessed by pointers to dynamically allocated arrays. To assist in the retrieving and

storing of grid data a number of Autodyn functions are provided. Examples are also given in Subroutine

EXEDIT - 2D (p. 98) and Subroutine EXEDIT3 – 3D (p. 100) for Autodyn-2D and 3D respectively.

Grid variable arrays and pointers are defined in module mdgrid (p. 55) for 2D and mdgrid3 (p. 58)

for 3D. The listings for Autodyn-2D and Autodyn-3D give all grid variables and their associated pointer

names. Autodyn Variable Listings (p. 73) provides a definition of these variables. The grid arrays and

pointer names are used to reference all grid variables either:

• by direct reference to arrays

• through use of the supplied Autodyn functions

By way of illustration, if we want to access the x-velocity at a particular (I,J,K) in the currently processed

Part in Autodyn-3D, we would include the mdgrid3 module in our subroutine:

USE mdgrid3

For node (I,J,K) of the current Part, we can obtain the index, IJK, for this node by including the statement:

IJK = IJKSET3(I,J,K)

Note

• For a shell Part, set I=1.

• For a beam Part, set I=1 and J=1.

• For an SPH Part, set I=1 and K=1.

• For 2D, use IJK=IJSET(I,J) instead.

Then we can obtain the required velocity, XVEL, either directly through the array reference or by using

a function:

Direct: XVEL = UXN(IJK)

Function: XVEL = GV3(NUXN,IJK)

Where NUXN is the pointer index for the X-velocity array. GV3 (p. 21) is an Autodyn-3D function which

retrieves the specified variable value from the current Part (the equivalent Autodyn-2D function is simply

GV (p. 21)). The various Autodyn functions such as IJKSET3 (p. 20) and GV3 (p. 21) are described in

detail in Autodyn Utility Functions/Subroutines (p. 20).

Normally, the direct reference approach is recommended. This works best when there is a current Part

(NSUB) (see Timing of Calls to User Subroutines (p. 17)). For references to Parts other than current, use

of the functions is recommended.

In the Euler Godunov processor in 3D, the access to the dynamic memory management is slightly dif-

ferent than the other processors. While, the function GV3 (p. 21) can be called as described before, the

of ANSYS, Inc. and its subsidiaries and affiliates.

Autodyn Variables

direct access method uses Fortran90 pointers. Thus, to reference the x velocity for a 3D Euler Godunov

cell:

Direct: XVEL = MTSUB(IJK)%V(NNUXN)

The variable NNUXN is contained in the module mdgrid3 (p. 58). Any specific Euler Godunov variable

index is defined by NNname, where name is the Autodyn-3D internal variable name (see Autodyn-3D

– Structured (IJK) Solvers (p. 78)).

3.5.2. Unstructured Solvers/Elements

For the unstructured solvers, variables associated with nodes and elements can be accessed and stored

using specific functions developed for the user:

CALL GET_ELEM_VAR(index_elem,index_layer)

CALL GET_NODE_VAR(index_node,index_material)

where

index_elem is the internal index for a given element.

index_node is the internal index for a given node.

index_layer is the layer (integration point) number for each element. For single point integrated

solid elements index_layer can be selected as either 0 or 1. For shell elements index_layer defines the

shell sublayer number.

index_material indicates the material that is being referred to on a given node. For all nodes, except

NBS tetrahedral nodes, index_material can be selected as either 0 or 1. For NBS tetrahedral nodes in-

dex_material defines the material number.

These functions copy data from the underlying data structures to local arrays. These arrays contain the

values of all possible variables. Only the variables used by the element/node identified in the above

calls will actually be set. To access these local arrays you will need to include

USE mdvar_all

RVL(index_var) will contain all real variables for the element/node

IVL(index_var) will contain all integer variables for the element/node

The values of index_var required to access a particular variable can be can be obtained through the

Output, Save, Review variables, Review all unstructured variables option in the interface.

For example, to retrieve the internal energy for an element, you could use the statement

INT_ENERGY = RVL(IVR_INT_ENERGY)

Only variables that are actually used in the analysis can be retrieved using the RVL(...) option. The variables

available depend on the element types and materials used in an analysis, and can be reviewed through

the Output, Save, Review variables, Review used unstructured variables option in the interface.

The values of the variables in the array RVL and IVL can be updated and stored by first setting the up-

dated value in the RVL or IVL array: for example,

RVL(IVR_INT_ENERGY)= NEWVAL

Then use the functions

CALL PUT_ELEM_VAR(index_elem,index_layer)

CALL PUT_NODE_VAR(index_node,index_material)

of ANSYS, Inc. and its subsidiaries and affiliates.14

Writing Your Own User Subroutines

This will place all the new data set for variables in the array RVL and IVL back into permanent main

storage for this element or node respectively.

Note

• The values of user variables VAR01 to VAR20 can also be updated and stored by first setting the

updated value in the RVL array and then use the function PUT_ELEM_VAR to place the new values

back in permanent main storage for the element.

• See Access to Nodal Variables within an Element (p. 104) for an elaborated example on how to

retrieve and store element and nodal variables.

3.5.3. Notes About Autodyn Variables

Do not change variables other than user variables 1 through to 20 because these changes will not be

stored to underlying data structures when PUT_ELEM_VAR (p. 23)/PUT_NODE_VAR (p. 23) are called.

Variables passed into the material modeling subroutines as arguments can be changed. Calls to

GET_ELEM_VAR (p. 22)/GET_NODE_VAR (p. 23) and PUT_ELEM_VAR (p. 23)/PUT_NODE_VAR (p. 23)

should not be used in the material user subroutines for materials used to fill ANP and NBS tetrahedra,

with the exception of the user erosion subroutine in mdero_user_1.f90. An example of using these

calls in mdero_user_1.f90 for NBS tetrahedra is given in Access to Nodal Variables for NBS Tetra-

hedral Elements (p. 105).

If GET_ELEM_VAR (p. 22) is called for an element other than the one currently being processed in the

solver (or for a shell sublayer that differs from the current one being processed), it will not be possible

to tell whether the values in RVL for that element are from the current or the previous computational

cycle.

The above methods for retrieving and updating unstructured element or node data require the global

index_elem or index_node as input. For material modeling user subroutines, the index of the current

element will be available via

USE mdstring

Index_elem = ELEM_NOW

For other cases, there are a number of ways in which one can obtain the internal index for an element:

for example, direct access, looping over Parts, looping over a Group. Examples of various types of ele-

ment/node access are given by example in Unstructured Element Data Access (p. 101).

3.6. Accessing Multiple Material Data – All IJK Solvers

Depending upon the Autodyn processor (solver) being used, it is possible for a single cell (element) to

contain more than one type of material. This is primarily applicable to Euler and Euler Godunov pro-

cessors. In Autodyn-2D, for multiple material cells, the standard cell values of compression, internal

energy, temperature, and alpha are mass-weighted averages of multiple material values. In Autodyn-

3D, plotting options allow you to plot multiple material values for a specific material or for mass-weighted

values as in 2D. Normally, if you wish to access the individual cell values of these variables for each

material in the cell, you must make the following subroutine call:

CALL GETMLT(IJK,0) or GETMLT3(IJK,0)

Where IJK is the index of the cell. After the call is made, the cell variables may be accessed according

to the following table:

of ANSYS, Inc. and its subsidiaries and affiliates.

Accessing Multiple Material Data – All IJK Solvers

ML IndexInternal ArrayMulti-material variable

NCVF*CVF(1,matno)Material volume fraction

NCMS*CMS(1,matno)Material mass

NCMU*CMU(1,matno)Material compression

NCEN*CEN(1,matno)Material internal energy

NCTP*CTP(1,matno)Material temperature

NCAL*CAL(1,matno)Material alpha

NCBF*CBF(1,matno)*Material burn fraction

NCDM*CDM(1,matno)*Material damage

NCPS*CPS(1,matno)*Material plastic strain

*Autodyn-3D only

where

matno is the material number for a given material. These material numbers are assigned in sequence

when defining your problem, starting with 1.

You can determine what the material number is for a given material name by testing against the MA-

TERIALS(matno)%NAME array (see the examples in User Subroutine Examples (p. 91)).

In Autodyn-3D, the GETV3 (p. 22) function provides a functional method for accessing multi-material

variables without use of GETMLT3 (p. 24). This is discussed in Autodyn Utility Functions/Sub-

routines (p. 20).

Note that to save updated values of multi-material variables, a call PUTMLT3(IJK,0) must be made after

setting the data in the internal arrays listed above.

For Lagrange, ALE, Shell, SPH and Beam solvers of Autodyn-3D it is also possible to access and set data

in the material arrays directly using the following procedure: after setting IJK for each cell, use the call

ML => MTSUB(IJK)%V(1:NUMMLV)

This sets-up the pointer ML to look at the material data for the current cell. This data can be accessed

and set directly by addressing the appropriate index of the ML pointer array. For example, to set mater-

ial damage to one and internal energy to zero, we could now use

ML(NCDM) = 1.0

ML(NCEN) = 0.0

This direct method of access is significantly more efficient than using the GETMLT3 (p. 24)/PUTM-

LT3 (p. 24) calls in Autodyn-3D.

3.7. Material Modeling User Subroutines

The main material modeling user subroutines (MD_EOS_USER_1.F90, MD_STR_USER_1.F90,

MD_FAI_USER_1.F90, MD_ERO_USER_1.F90) have been updated and modularized to allow more

flexibility. The table below shows the general description of the usage of the routines, and also the four

basic components of each routine:

of ANSYS, Inc. and its subsidiaries and affiliates.16

Writing Your Own User Subroutines

剩余113页未读，继续阅读

少年小鱼

粉丝: 31
资源: 528

ANSYS Autodyn User子程序指南

Autodyn Composite Modeling Guide.pdf

Autodyn Parallel Processing Guide.pdf

ANSYS_Fluent_Users_Guide.pdf

Ansys_Fluent_Users_Guide.pdf

ANSYS Autodyn Users Manual.pdf

Workbench Users Guide.pdf

DesignXplorer Users Guide.pdf

System Coupling Users Guide.pdf

FE Modeler Users Guide.pdf

ANSYS Forte Users Guide.pdf

最新资源