Manual mesh moving procedure

Research Highlights:

• ABAQUS subroutine UMESHMOTION CANNOT model ablation for composite laminates.

• A new FEA procedure, Manual Mesh Moving (MMM) procedure, is proposed to model ablation of composite laminates in ABAQUS.

• The new procedure is verified with simulations using existing validated method.

• The (laser) ablation of a CFRP composite laminate is captured using the MMM procedure.

Demonstration of MMM Procedure with An Isotropic Material (2D Plane):

• Model Input - Material: 2D isotropic plane, thermal conductivity 6.578E-4 W/m⸱°C, specific heat 1004.3 J/kg⸱°C, density 1.814E-6 kg/m^3

• Model Input - Heat source: uniform heat flux, 250 watt per square meters

• Model Input - Prescribed ablation rate: -308*EXP(-1.5E5*(t-2.5E-6))+100, when t>=10 microseconds

• Modeling methods: using both traditional UMESHMOTION+ALE method and our proposed MMM procedure

• Model output: compare the temperature histories at A, B, C, D, and E points (see Fig. 1 below) and ablation depth (see Fig. 2 below) obtained from the two methods

Instructions for Running the UMESHMOTION+ALE Code and the MMM Code:

• Source code for the UMESHMOTION+ALE method is in the umesh.zip file, in which it contains the ABAQUS input file and the user subroutine. Instructions to run the code is listed below:

  • Step-1: unzip the umesh.zip file, you will see a folder named umesh containing the ABAQUS input file Umesh_NonLinS.inp and the ABAQUS user subroutine file: Umesh_NonLinS.f

  • Step-2: open Windows terminal in the directory that contains the input file and the user subroutine file, execute command: abaqus job=Umesh_NonLinS.inp user=Umesh_NonLinS.f int

• Source code for the Manual Mesh Moving method is in the mmm.zip file, in which it contains the ABAQUS input file, user subroutines (USDFLD and URDFIL), MATLAB and Python scripts. Instructions to run the code is listed below:

  • Step-1: unzip the mmm.zip file, you will get a folder named mmm containing all the necessary files required to run the code. Note that subfolders DatFiles and InpFiles are initially empty. Files in those subfolders will be generated during the simulation.

  • Step-2: find udisp.f, sdv.f, sdvx.f files in the main folder, replace the directory path (C:\Users\Yeqing Wang\Desktop\webpage\mmm) with your own directory path of the mmm folder (lines 15, 16 in udisp.f, lines 54-57 in sdv.f, lines 46-49 in sdvx.f). Likewise, find GetSdvX.m file in the MFiles folder, replace the same directory path with your own directory path of the mmm folder (lines 65-68).Take special cautious when editing the directory path with line continuation symbols.

  • Step-3: start the simulation by executing the Ablation.m file in MATLAB. This should start the simulation.

  • Note: the result files (Odb Files) will be generated in a good number of org-i.odb and org-static-i.odb (where i=1,2,3,...) files. Also note that the org-i.odb files are the temperature and ablation depth results that are of interest. Please delete the org-static-i.odb files after the simulation.

Google Drive Link to Download the Codes:

https://drive.google.com/open?id=1ub2nfdEUEAWuEiypUI2yE731t6jKf8r4

Note: The MMM codes were removed from this site. Please send me an email (yqgwang@gmail.com) to request the code package.

Publication List:

• Y. Wang & C. L. Pasiliao (2018), Modeling Ablation of Laminated Composites: A Novel Manual Mesh Moving Finite Element Analysis Procedure with ABAQUS, Int. J. of Heat and Mass Transfer, Vol. 116, pp. 306-313.

• Y. Wang, O. I. Zhupanska & C. L. Pasiliao, “Verification of a Manual Mesh Moving Finite Element Analysis Procedure for Modeling Ablation in Laminated Composite Materials,” ASME 2017 International Mechanical Engineering Congress and Exposition, Tampa, FL, 2017.