Test_NeoHookean

 

 

제너럴한 목적은,

Growth가 없을 때, Convergence 를 얻을 수 있느냐 없느냐이다. 

Growth가 없으면, 얻을 수 있어야 한다. 

 

그니까 이런상황이었음. 둘다 내가 쓴 코드고

Growth Code 가	있다.	없다.
M_e_crit	100 (그래서 저 Growth Code skipped)
Proportional force	Batch / DT,,dt / PROP,,1 / End / 2, 1 0    0 / 10   1 / 20   0 / 30   0 / 40   1 / 200  10	Batch / DT,,dt / PROP,,1 / End /  2, 1 0    0 / 10   1 / 20   0 / 30   0 / 40   1 / 200  10
Residual	E-00,E-04,E-09,E-14	E-00,E-04,E-09,E-15
	mat02안에서, F를 받아서, matinv로 detF를 계산.	kine3df에서 받은 detF로 계산.

그래서 이걸보고

아씨발 저 Growth Code 가 뭔가 잘못됐구나 했다. 

그래서 그럼 빌트인 코드에서는 R이 어떻게 계산되나 봤다.

그래서 해본게 아래꺼다.

빌트인에서의 차이점은 E,nu를 필요로 한다는것이다.

해봤더니 심지어 내가쓴 코드가 더 정확하다는것.

내 코드가 쿼드라틱하다.

by_feapcode.tar.gz

6.53MB

Left: 내가 쓴 코드로 neo Hookean 계산한거. Right: Feap built-in 으로 계산한거.

빌트인

내코드

 

 

자, 이제 그럼.

<timestep = 2.d0> 에서는 그랬는데 그럼. 계속 진행되면서도 괜찮은거임?