Modal Submodeling
This post was inspired by CAEAI's blog post on modal analysis sub-modeling and an XANSYS question. Their instructions were not explicit and everything was done in Ansys Classic. The goal here is to show modal submodeling in Ansys Workbench. Here's a walk through of an example:
Step 1
Solve the coarse model Modal Analysis (System A).
Step 2
Create a finely meshed submodel Harmonic Response Model (System B). This will be the submodel.
After creating the desired mesh, add a named selection at the cut boundaries called cut_face. Add the following command snippet in the Setup branch.
!! Creates cut boundary NODE file
! WB will complain but boundarynodes.node file should be created
/prep7
cmsel, s, cut_face
nwrite, boundarynodes, node
finish
/exit
The above step writes out a file listing the nodes and their coordinates which will be used for interpolation. After solving for the model, Ansys Workbench will flash a lot of errors. Ignore them. The key output is boundarynodes.node which is located at the Solver File Directory (RMB Solution > Open Solver File Directory). Copy that file.
Step 3
Open the Solver File Directory of System A (Modal Analysis) and paste the boundarynodes.node file in that directory. In System A, add the following command snippet into the /POST1 snippet:
!! Creates interpolation file
/post1
set, list
set, 1, 7 ! extracts Mode #7 for submodeling later
cbdof, boundarynodes, node,, Dfile, cbdo
After solving, the Dfile.cbdo will be created in the working directory. This contains the interpolated displacements of the coarse model on the fine model's cut boundary nodes. Copy this file and paste it into the working directory of System B (Harmonic Analysis).
Step 4
Open System B (Harmonic Analysis). Firstly suppress the earlier snippet.
Next, add another snippet (named HarmonicSubmodel3 in this example):
! Cut Boundary Conditions on Fine Model
/input, Dfile, cbdo
The one line script above will apply the cut boundary displacements as D commands.
For analysis settings, it is critical Full Solution method is used as non-zero imposed displacements are not supported in Mode Superposition Method.
The Range Maximum should be set at the natural frequency of the mode that was extracted earlier. Only 1 solution is required. Finally, Output Controls should have both Stress and Strain set to "yes".
The results should look as follows:
Step 5
Modal analysis amplitudes are by-themselves meaningless. Scale the mode shape as required by creating an Identifier in the Equivalent Stress result and use a User Defined Result to scale it. For example an Expression could be 3.1415*MySeqv shown below. The scaling number is used to match experimental results or expected deflections.
Associated Files
Archived Ansys Workbench v18.2: File
Step 1
Solve the coarse model Modal Analysis (System A).
Fig 1: Coarse Model Modal Analysis (System A)
Step 2
Create a finely meshed submodel Harmonic Response Model (System B). This will be the submodel.
Fig 2: Project Schematic
After creating the desired mesh, add a named selection at the cut boundaries called cut_face. Add the following command snippet in the Setup branch.
!! Creates cut boundary NODE file
! WB will complain but boundarynodes.node file should be created
/prep7
cmsel, s, cut_face
nwrite, boundarynodes, node
finish
/exit
The above step writes out a file listing the nodes and their coordinates which will be used for interpolation. After solving for the model, Ansys Workbench will flash a lot of errors. Ignore them. The key output is boundarynodes.node which is located at the Solver File Directory (RMB Solution > Open Solver File Directory). Copy that file.
Step 3
Open the Solver File Directory of System A (Modal Analysis) and paste the boundarynodes.node file in that directory. In System A, add the following command snippet into the /POST1 snippet:
!! Creates interpolation file
/post1
set, list
set, 1, 7 ! extracts Mode #7 for submodeling later
cbdof, boundarynodes, node,, Dfile, cbdo
After solving, the Dfile.cbdo will be created in the working directory. This contains the interpolated displacements of the coarse model on the fine model's cut boundary nodes. Copy this file and paste it into the working directory of System B (Harmonic Analysis).
Step 4
Open System B (Harmonic Analysis). Firstly suppress the earlier snippet.
Fig 3: Only Second Snippet Used
Next, add another snippet (named HarmonicSubmodel3 in this example):
! Cut Boundary Conditions on Fine Model
/input, Dfile, cbdo
The one line script above will apply the cut boundary displacements as D commands.
The Range Maximum should be set at the natural frequency of the mode that was extracted earlier. Only 1 solution is required. Finally, Output Controls should have both Stress and Strain set to "yes".
The results should look as follows:
Fig 4: Submodel Stress
Step 5
Modal analysis amplitudes are by-themselves meaningless. Scale the mode shape as required by creating an Identifier in the Equivalent Stress result and use a User Defined Result to scale it. For example an Expression could be 3.1415*MySeqv shown below. The scaling number is used to match experimental results or expected deflections.
Fig 5: Scaled Submodel Stress
Associated Files
Archived Ansys Workbench v18.2: File
How we can insert the desire input acceleration. Also, the magnitude of stress is too bigger than real. how we can solve this problem?
ReplyDeleteThe amplitudes from modal results are not scaled. To get the acceleration that you want, you can determine the equivalent displacements and scale the modal results accordingly. The stresses should then reflect that.
ReplyDelete