## You are here

# Comments

### gaussian distribution or Goldak heat source model in ANSYS?

*In reply to SIMULATION OF WELDING*

dear friends.. i am learning FEA on ANSYS for Arc Welding. can anyone tell me how to model Gaussian or Goldak heat source model in ANSYS?

### Inquiry about your available PHD positions

*In reply to PhD studentship at URV, Tarragona, Spain*

<p>Dar Prof. Huera,</p>

<p>My research fields exactly coincides with your research areas. I sent an email to you yesterday with my CV attached. I will appreciate you if you take a look at my CV.</p>

<p>Thanks in advance</p>

<p>Morteza Alinia</p>

<p>Email: <a href="mailto:m.alinia@me.iut.ac.ir">m.alinia@me.iut.ac.ir</a></p>

<p> </p>

### hello Mr Jahanshahi

*In reply to Displacement Controlled Analysis*

hello Mr Jahanshahi

thank you so much for your attention.

I will do what you said.

M.Jafarzadeh

### hello peto24ap

*In reply to Geometric nonlinearity*

hello peto24ap

Why do you use line 4 method for computing internal force? I think this is the problem. Use ordinary method for internal force, Just like the method introduced in " introduction to finite elements in engineering" ,by chandrupatla or every other fnite element book, to make your internal force vector.

look for this article "COMPOSITE MATERIALS IN COMPRESSION " by Jens Lycke Wind, In this thechnical report you can find an algorithm for Newton raphson.

I used it before for load base problem but i could not write a code for material nonlinearity which is under displacement BC, and displacement control.

Ask me if there is some thing to ask, our ways are similar. maybe we can help each other in the modeling

M.J

### Hi, There are several arc

*In reply to arc length in delamination*

Hi, There are several arc-length methods that could be used. But depending on the level of non-linearity, they may or may not be able to trace snap-back exactly.

Geometrical nonlinearity simple arc-length method should be sufficient. But in case of severe strain localization, you would need a method like the one suggested by

Crisfield arc-length method. Try this one, incase this one does not work, I may have other suggestions.

Best

### Displacement Controlled Analysis

*In reply to displacement control method *

Dear Mr. Jafarzadeh,

I already have answered your question. The answer lies in equation (2) above. In the fisrt iteration of each load step (which is displacement controlled) you apply an incremental displacement. Of course your stiffness matrix should be comouted in each iteration due to material nonlinerity. In subsequent iterations the second term on the right hand side of equation (2) is set to zero and iterations continue to equilibrate out of balance force. The iterative procedure is exactly similar to force-controlled except for equation (20). If you pay a good attention to what I siad you should be able to proceed without any problem.

Good Luck

### displacement control method

*In reply to Displacement Controlled Analysis*

hello Mr Jahanshahi

thank you for your quick answer

I know how to find reaction forces for the case but my problem is how to use displacement control method for displacement BC case, in which the stiffness matrix in each increment and iteration change. In displacement BC by changing stiffness matrix the total force changes too. I should apply a unique force for each increment to use it in all the iteration.

M.Jafarzadeh

### Displacement Controlled Analysis

*In reply to N-R method*

Consider the following partitioned matrix equation (sorry if I could not show it clearly):

|**K**nn ¦ **K**nu||**u**n| |**R**|

|.......¦.......||....| = |...| (1)

|**K**un ¦ **K**uu||**u**u| |**F**|

In this equation *u**n* is the vector of specified (and thus known) displacements. This is the specified displacement at certain degrees of freedom in displacement controlled analysis. ** R** is the vector of corresponding reactions at these degrees of freedom. The vector

*u**u*is the vector of unknown displacements at other degrees of freedom where no displacement boundary condition is specified.

**is the vector of nodal forces at these latter degrees of freedom. From the second row of this equation we can write:**

*F*[**K**uu][**u**u]=[**F**]-[**K**un][**u**n] (2)

which is the equation that should be solved in displacement controlled analysis. Based on this equation, the vector of forces at unconstrained degrees of freedom should be modified by the second term on the right hand side of equation (2). The rest of the analysis is exactly the same as force controlled analysis. Having calculated *u**u* from preceding equation, the vector of reactions at constrained degrees of freedom can be calculated from the first row of equation (1) as

[**R**]=[**K**nn][**u**n]+[**K**nu][**u**u] (3)

For more information you can refer to page 187 of the following reference:

K.J. Bathe, *Finite Element Procedures*, Prentice-Hall, 1996.

### Helpful tips really nice

### D&DT course - Part II on Fatigue & Damage Tolerance

*In reply to Course on Design Principles, Fatigue & Damage Tolerance for Aeronautical Engineers*

Dear iMechanicians,

Part I of the D&DT courese on Design Principles for Aeronautical Engineers starts just after this weekend. Those, who would be more interested in Fatigue & Damage Tolerance topics, please note this part of the course will be held from November 11 till November 13. The applications, if possible, should be filled in till October 27. If interested, see the details on www.pragtic.com/DDT.php.

Best regards

Jan

### N-R method

*In reply to Nonlinear Algorithms*

Hello Mr Jahanshahi

I wrote a code by N-R method for my problem (a double cantilever beam which has interface elements between two layer for simulating delamination ), every thing was good, I saw the softening of the material, but there was some problem.

For verifying my answers, the Reaction-displacement of the total system should be plotted,for this reason i should impose displacement boundary condition instead of force.

As you know in the material nonlinearity problem the stiffness matrix (D and after that K) will change in every increment and iteration corresponding to displacement (in the geometrical nonlinearity the displacement changes alter the geometry and after that the total stiffness matrix changes). by changing in the stiffness matrix also the force value will change , because of displacement BC formulation in finite element method.

but in the N-R method the at the beginnig a load apply to the system and then in each increment 1 part of it will impose. In my case the force change with K.

now here is my question: Can i use N-R method for displacement BC problem or I must use the displacement control method?

I try displacement control before but I didnt succeed. In this method I have same problems. In displacement control method a first force (not displacement) will introduce, and then a part of force impose in increments. But the base is displacement change. Also In this method the K matrix changes and it is main problem.

I want a method appropriate for displacement BC problems.

Would you please help me in this topic too?

M.Jafarazadeh

### You can post your specific

*In reply to could someone tell what should I do to run Abaqus job in High Performance Computing (HPC) systems?*

You can post your specific question here, https://swym.3ds.com/#community:73/iquestions

### Thanks for your help! A

*In reply to This problemlooks like*

Thanks for your help! A follow up question is that when import a model to Solidworks or AutoCAD, Abaqus can export multiple file extension, such as ACIS (.sat), IGES (.igs)and STEP (.stp). Do you know which is better from 3D reconstruction. By the way, our model is a cylindrical shell.

### Postdoc looking for a job

*In reply to Research Associate/ Postdoctoral Research Associate position in Magnesium Alloy Development*

الأخ المحترم الدكتور بلال منصور

السلام عليكم ورحمة الله وبركاته

أود الاستفسار عن فرصة العمل المذكورة اعلاه إذا لا تزال متوفرة أم لا؟

أنا ابحث عن فرصة عمل في مجال دراسة سبائك الماغنيسيوم وطرق تحسين خواص التآكل بها.

أنا خريج جامعة مانشستر لهذه السنة ولدي درجة دكتوراة في هذا المجال وكذلك درجة ماجستير في علوم المواد والتآكل.

ساقوم بتزويدكم بسيرتي الذاتية في حال كانت هذه الوظيفة مازالت شاغرة.

تقبلوا تحياتي مع جزيل الشكر

والسلام عليكم ورحمة الله وبركاته

د. أحمد الشيخي

### These are a number of research papers and textbooks

*In reply to These are a number of research papers and textbooks*

yes, These are a number of research papers and textbooks by Prof. W. Noll. I suppose they definitely deserve to be here.

### This problemlooks like

*In reply to 3D Printing of an imperfect cylindrical shell*

This problemlooks like surface reconstruction problem. I have heard of plugin in solid works "Scanto3D". See

https://www.youtube.com/watch?v=GrVZR1cldu4

This plugin can be activated if it does not exists.

There are other packages which can generate surface from point cloud:

1. Resurf plugin in Rhinoceros : fits NURBS surface to point clouds.

2. Mesh lab : using "voronoi filtering" algorithm generates facets on the given point cloud.

### Hello Mr Jahanshahi

*In reply to Nonlinear Algorithms*

Hello Mr Jahanshahi

Thank you for your attention.

I will read the Owen book certainly and I will apply the Newton Raphson in my code too.

m.jafarzadeh

### hello Mr Sreenivas

*In reply to Hello
*

hello Mr Sreenivas

Today I saw your comment about the book "Development and Application of the FEM based on MATLAB", but the link doesnt work any more.

I tried to find it by google search,but there wasnt any thing.

Would you please guide me to find this book.

by the way, I'm trying to simulate a material nonlinearity problem, Do you know some thing that can help me in this way?

thanks

m.j

### Nonlinear Algorithms

*In reply to displacement control for material nonlinearity problem*

Hello m.j,

1. You have to be careful about arc length method. If lambda is greater than 1 then the results obtained by arc length method might not be correct. Using full Newton-Raphson method with reasonable step size can be a good substitute if you are not sure of what you are doing.

2. You can consider two tolerances for your algorithm for example 10e-10 and 10e-5. You can enforce the first tolerance in your iterations. If the tolerance obtained after the total iterations are exhausted is greater than the first tolerance but less than the second your results are fairly accurate within the second tolerance and you can move to the next iteration. Otherwise your results are not accurate enough and you have to cease the iterations.

3. You can consult the following book. It is full of algorithms and flow charts which suit your needs:

Computational Methods in Plasticity: Theory and Applications, EA de Souza Neto, D Peric and DRJ Owen, John Wiley & Sons, 2008.

Mohsen

## Recent comments